A Genomic Portrait of Haplotype Diversity and Signatures of Selection in Indigenous Southern African Populations

We report a study of genome-wide, dense SNP (∼900K) and copy number polymorphism data of indigenous southern Africans. We demonstrate the genetic contribution to southern and eastern African populations, which involved admixture between indigenous San, Niger-Congo-speaking and populations of Eurasian ancestry. This finding illustrates the need to account for stratification in genome-wide association studies, and that admixture mapping would likely be a successful approach in these populations. We developed a strategy to detect the signature of selection prior to and following putative admixture events. Several genomic regions show an unusual excess of Niger-Kordofanian, and unusual deficiency of both San and Eurasian ancestry, which were considered the footprints of selection after population admixture. Several SNPs with strong allele frequency differences were observed predominantly between the admixed indigenous southern African populations, and their ancestral Eurasian populations. Interestingly, many candidate genes, which were identified within the genomic regions showing signals for selection, were associated with southern African-specific high-risk, mostly communicable diseases, such as malaria, influenza, tuberculosis, and human immunodeficiency virus/AIDs. This observation suggests a potentially important role that these genes might have played in adapting to the environment. Additionally, our analyses of haplotype structure, linkage disequilibrium, recombination, copy number variation and genome-wide admixture highlight, and support the unique position of San relative to both African and non-African populations. This study contributes to a better understanding of population ancestry and selection in south-eastern African populations; and the data and results obtained will support research into the genetic contributions to infectious as well as non-communicable diseases in the region.     

Admixture mapping of quantitative traits in Populus hybrid zones: power and limitations

Uncovering the genetic architecture of species differences is of central importance for understanding the origin and maintenance of biological diversity. Admixture mapping can be used to identify the number and effect sizes of genes that contribute to the divergence of ecologically important traits, even in taxa that are not amenable to laboratory crosses because of their long generation time or other limitations. Here, we apply admixture mapping to naturally occurring hybrids between two ecologically divergent Populus species. We map quantitative trait loci for eight leaf morphological traits using 77 mapped microsatellite markers from all 19 chromosomes of Populus. We apply multivariate linear regression analysis allowing the modeling of additive and non-additive gene action and identify several candidate genomic regions associated with leaf morphology using an information-theoretic approach. We perform simulation studies to assess the power and limitations of admixture mapping of quantitative traits in natural hybrid populations for a variety of genetic architectures and modes of gene action. Our results indicate that (1) admixture mapping has considerable power to identify the genetic architecture of species differences if sample sizes and marker densities are sufficiently high, (2) modeling of non-additive gene action can help to elucidate the discrepancy between genotype and phenotype sometimes seen in interspecific hybrids, and (3) the genetic architecture of leaf morphological traits in the studied Populus species involves complementary and overdominant gene action, providing the basis for rapid adaptation of these ecologically important forest trees.     

Genetic dissection of top three leaf traits in rice using progenies from a japonica × indica cross

The size of the top three leaves of rice plants is strongly associated with yield; thus, it is important to consider quantitative traits representing leaf size(e.g.,length and width) when breeding novel rice varieties.It is challenging to measure such traits on a large scale in the field, and little is known about the genetic factors that determine the size of the top three leaves. In the present study, a population of recombinant inbred lines(RILs) and reciprocal single chromosomal segment substitution lines(SSSLs) derived from the progeny of a japonica Asominori ? indica IR24 cross were grown under four diverse environmental conditions. Six morphological traits associated with leaf size were measured,namely length and width of the flag, second and third leaves. In the RIL population, 49 QTLs were identified that clustered in 30 genomic region. Twenty-three of these QTLs were confirmed in the SSSL population. A comparison with previously reported genes/QTLs revealed eight novel genomic regions that contained uncharacterized ORFs associated with leaf size. The QTLs identified in this study can be used for markerassisted breeding and for fine mapping of novel genetic elements controlling leaf size in rice.     

Native American Admixture in the Quebec Founder Population

For years, studies of founder populations and genetic isolates represented the mainstream of genetic mapping in the effort to target genetic defects causing Mendelian disorders. The genetic homogeneity of such populations as well as relatively homogeneous environmental exposures were also seen as primary advantages in studies of genetic susceptibility loci that underlie complex diseases. European colonization of the St-Lawrence Valley by a small number of settlers, mainly from France, resulted in a founder effect reflected by the appearance of a number of population-specific disease-causing mutations in Quebec. The purported genetic homogeneity of this population was recently challenged by genealogical and genetic analyses. We studied one of the contributing factors to genetic heterogeneity, early Native American admixture that was never investigated in this population before. Consistent admixture estimates, in the order of one per cent, were obtained from genome-wide autosomal data using the ADMIXTURE and HAPMIX software, as well as with the fastIBD software evaluating the degree of the identity-by-descent between Quebec individuals and Native American populations. These genomic results correlated well with the genealogical estimates. Correlations are imperfect most likely because of incomplete records of Native founders’ origin in genealogical data. Although the overall degree of admixture is modest, it contributed to the enrichment of the population diversity and to its demographic stratification. Because admixture greatly varies among regions of Quebec and among individuals, it could have significantly affected the homogeneity of the population, which is of importance in mapping studies, especially when rare genetic susceptibility variants are in play.     

Markers informative for ancestry demonstrate consistent megabase-length linkage disequilibrium in the African American population

Admixture mapping is a potentially powerful tool for mapping complex genetic diseases. For application of this method, admixed individuals must have genomes composed of large segments derived intact from each founding population. Such segments are thought to be present in African Americans (AA) and should be demonstrable by examination of linkage disequilibrium (LD). Previous studies using a variety of polymorphic markers have variably reported long-range LD or rapid decay of LD. To further define the extent and characteristics of LD caused by admixture in the AA population, the current study utilized a set of 52 diallelic markers that were selected for large standard variances between putative representatives of the founder populations. LD was examined in over 250 marker-pairs, including linked markers from four different chromosomal regions and an equal number of matched unlinked comparisons. In the representative founder populations, strong LD was not observed for markers separated by more than 10 kb. In contrast, results indicated significant LD ( P<0.001, D'>0.3) in AA over large genomic segments exceeding 10 centiMorgans (cM) and 15 megabases (Mb). Only marginally significant LD was present between unlinked markers in this population, suggesting that choosing appropriate levels of significance for admixture mapping can minimize false positive results. The ability to detect LD for extended chromosomal segments in AA decayed not only as a function of the distance between markers, but also as a function of the standard variance of the markers. This examination of several genomic segments provides strong evidence that appropriate selection of informative markers is a crucial prerequisite for the application of admixture mapping to the AA population.     

Admixture mapping of male nuptial colour and body shape in a recently formed hybrid population of threespine stickleback

Despite recent progress, we still know relatively little about the genetic architecture that underlies adaptation to divergent environments. Determining whether the genetic architecture of phenotypic adaptation follows any predictable patterns requires data from a wide variety of species. However, in many organisms, genetic studies are hindered by the inability to perform genetic crosses in the laboratory or by long generation times. Admixture mapping is an approach that circumvents these issues by taking advantage of hybridization that occurs between populations or species in the wild. Here, we demonstrate the utility of admixture mapping in a naturally occurring hybrid population of threespine sticklebacks (Gasterosteus aculeatus) from Enos Lake, British Columbia. Until recently, this lake contained two species of sticklebacks adapted to divergent habitats within the lake. This benthic–limnetic species pair diverged in a number of phenotypes, including male nuptial coloration and body shape, which were previously shown to contribute to reproductive isolation between them. However, recent ecological disturbance has contributed to extensive hybridization between the species, and there is now a single, admixed population within Enos Lake. We collected over 500 males from Enos Lake and found that most had intermediate nuptial colour and body shape. By genotyping males with nuptial colour at the two extremes of the phenotypic distribution, we identified seven genomic regions on three chromosomes associated with divergence in male nuptial colour. These genomic regions are also associated with variation in body shape, suggesting that tight linkage and/or pleiotropy facilitated adaptation to divergent environments in benthic–limnetic species pairs.     

cis-Regulatory changes in Kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans

Dramatic pigmentation changes have evolved within most vertebrate groups, including fish and humans. Here we use genetic crosses in sticklebacks to investigate the parallel origin of pigmentation changes in natural populations. High-resolution mapping and expression experiments show that light gills and light ventrums map to a divergent regulatory allele of the Kit ligand (Kitlg) gene. The divergent allele reduces expression in gill and skin tissue and is shared by multiple derived freshwater populations with reduced pigmentation. In humans, Europeans and East Asians also share derived alleles at the KITLG locus. Strong signatures of selection map to regulatory regions surrounding the gene, and admixture mapping shows that the KITLG genomic region has a significant effect on human skin color. These experiments suggest that regulatory changes in Kitlg contribute to natural variation in vertebrate pigmentation, and that similar genetic mechanisms may underlie rapid evolutionary change in fish and humans.     

Genetic assessment of safflower (Carthamus tinctorius L.) collection with microsatellite markers acquired via pyrosequencing method

A genetic evaluation of safflower germplasm collections derived from different geographical regions and countries will provide useful information for sustainable conservation and the utilization of genetic diversity. However, the molecular marker information is limited for evaluation of genetic diversity of safflower germplasm. In this study, we acquired 509 putative genomic SSR markers for sufficient genome coverage using next‐generation sequencing methods and characterized thirty polymorphic SSRs in safflower collection composed of 100 diverse accessions. The average allele number and expected heterozygosity were 2.8 and 0.386, respectively. Analysis of population structure and phylogeny based on thirty SSR profiles revealed genetic admixture between geographical regions contrary to genetic clustering. However, the accessions from Korea were genetically conserved in distinctive groups in contrast to other safflower gene pool. In conclusion, these new genomic SSRs will facilitate valuable studies to clarify genetic relationships as well as conduct population structure analyses, genetic map construction and association analysis for safflower.     

Geographic patterns of genome admixture in Latin American Mestizos

The large and diverse population of Latin America is potentially a powerful resource for elucidating the genetic basis of complex traits through admixture mapping. However, no genome-wide characterization of admixture across Latin America has yet been attempted. Here, we report an analysis of admixture in thirteen Mestizo populations (i.e. in regions of mainly European and Native settlement) from seven countries in Latin America based on data for 678 autosomal and 29 X-chromosome microsatellites. We found extensive variation in Native American and European ancestry (and generally low levels of African ancestry) among populations and individuals, and evidence that admixture across Latin America has often involved predominantly European men and both Native and African women. An admixture analysis allowing for Native American population subdivision revealed a differentiation of the Native American ancestry amongst Mestizos. This observation is consistent with the genetic structure of pre-Columbian populations and with admixture having involved Natives from the area where the Mestizo examined are located. Our findings agree with available information on the demographic history of Latin America and have a number of implications for the design of association studies in population from the region.     

Genomic consequences of a recent three‐way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human‐mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make‐up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double‐digested restriction site‐associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple‐way admixture with population genomic data.     

Behavioral genomics and the study of speciation at a porous species boundary

Porous species boundaries are characterized by differential gene flow, where some regions of the genome experience divergent evolution while others experience the homogenizing effects of gene flow. If species can arise or remain distinct despite gene flow between them, speciation can only be understood on a gene by gene level. To understand the genetics of speciation, we therefore must identify the targets of selection that cause divergent evolution and identify the genetic architecture underlying such “speciation phenotypes”. This will enable characterization of genomic regions that are “free to flow” between species, and those that diverge in the face of gene flow. We discuss this problem in the genus Laupala, a morphologically cryptic, flightless group of crickets that has radiated in Hawaii. Because songs are used in courtship and always distinguish close relatives of Laupala as well as species in sympatry, we argue that songs in Laupala are speciation phenotypes. Here, we present our approaches to identify the underlying genomic regions and song genes that differentiate closely related species. We discuss what is known about the genetic basis of this species difference derived from classic quantitative genetics and quantitative trait locus mapping experiments. We also present a model of the molecular expression of cricket song to assist in our goal to identify the genes involved in song variation. As most species are sympatric and exchange genes with congeners, we discuss the importance of understanding the genetic and genomic architecture of song as a speciation phenotype that must be characterized to identify differential patterns of gene flow at porous species boundaries.     

Regional admixture mapping and structured association testing: conceptual unification and an extensible general linear model

Individual genetic admixture estimates, determined both across the genome and at specific genomic regions, have been proposed for use in identifying specific genomic regions harboring loci influencing phenotypes in regional admixture mapping (RAM). Estimates of individual ancestry can be used in structured association tests (SAT) to reduce confounding induced by various forms of population substructure. Although presented as two distinct approaches, we provide a conceptual framework in which both RAM and SAT are special cases of a more general linear model. We clarify which variables are sufficient to condition upon in order to prevent spurious associations and also provide a simple closed form “semiparametric” method of evaluating the reliability of individual admixture estimates. An estimate of the reliability of individual admixture estimates is required to make an inherent errors-in-variables problem tractable. Casting RAM and SAT methods as a general linear model offers enormous flexibility enabling application to a rich set of phenotypes, populations, covariates, and situations, including interaction terms and multilocus models. This approach should allow far wider use of RAM and SAT, often using standard software, in addressing admixture as either a confounder of association studies or a tool for finding loci influencing complex phenotypes in species as diverse as plants, humans, and nonhuman animals.     

Admixture in Mexico City: implications for admixture mapping of Type 2 diabetes genetic risk factors

Admixture mapping is a recently developed method for identifying genetic risk factors involved in complex traits or diseases showing prevalence differences between major continental groups. Type 2 diabetes (T2D) is at least twice as prevalent in Native American populations as in populations of European ancestry, so admixture mapping is well suited to study the genetic basis of this complex disease. We have characterized the admixture proportions in a sample of 286 unrelated T2D patients and 275 controls from Mexico City and we discuss the implications of the results for admixture mapping studies. Admixture proportions were estimated using 69 autosomal ancestry-informative markers (AIMs). Maternal and paternal contributions were estimated from geographically informative mtDNA and Y-specific polymorphisms. The average proportions of Native American, European and, West African admixture were estimated as 65, 30, and 5%, respectively. The contributions of Native American ancestors to maternal and paternal lineages were estimated as 90 and 40%, respectively. In a logistic model with higher educational status as dependent variable, the odds ratio for higher educational status associated with an increase from 0 to 1 in European admixture proportions was 9.4 (95%, credible interval 3.8–22.6). This association of socioeconomic status with individual admixture proportion shows that genetic stratification in this population is paralleled, and possibly maintained, by socioeconomic stratification. The effective number of generations back to unadmixed ancestors was 6.7 (95% CI 5.7–8.0), from which we can estimate that genome-wide admixture mapping will require typing about 1,400 evenly distributed AIMs to localize genes underlying disease risk between populations of European and Native American ancestry. Sample sizes of about 2,000 cases will be required to detect any locus that contributes an ancestry risk ratio of at least 1.5.     

Mapping salinity tolerance during Arabidopsis thaliana germination and seedling growth

To characterize and dissect genetic variation for salinity tolerance, we assessed variation in salinity tolerance during germination and seedling growth for a worldwide sample of Arabidopsis thaliana accessions. By combining QTL mapping, association mapping and expression data, we identified genomic regions involved in salinity response. Among the worldwide sample, we found germination ability within a moderately saline environment (150 mM NaCl) varied considerable, from >90% among the most tolerant lines to complete inability to germinate among the most susceptible. Our results also demonstrated wide variation in salinity tolerance within A. thaliana RIL populations and identified multiple genomic regions that contribute to this variation. These regions contain known candidate genes, but at least four of the regions contain loci not yet associated with salinity tolerance response phenotypes. Our observations suggest A. thaliana natural variation may be an underutilized resource for investigating salinity stress response.     

Breed-specific ancestry studies and genome-wide association analysis highlight an association between the MYH9 gene and heat tolerance in Alaskan sprint racing sled dogs

Alaskan sled dogs are a genetically distinct population shaped by generations of selective interbreeding with purebred dogs to create a group of high-performance athletes. As a result of selective breeding strategies, sled dogs present a unique opportunity to employ admixture-mapping techniques to investigate how breed composition and trait selection impact genomic structure. We used admixture mapping to investigate genetic ancestry across the genomes of two classes of sled dogs, sprint and long-distance racers, and combined that with genome-wide association studies (GWAS) to identify regions that correlate with performance-enhancing traits. The sled dog genome is enhanced by differential contributions from four non-admixed breeds (Alaskan Malamute, Siberian Husky, German Shorthaired Pointer, and Borzoi). A principal components analysis (PCA) of 115,000 genome-wide SNPs clearly resolved the sprint and distance populations as distinct genetic groups, with longer blocks of linkage disequilibrium (LD) observed in the distance versus sprint dogs (7.5-10 and 2.5-3.75?kb, respectively). Furthermore, we identified eight regions with the genomic signal from either a selective sweep or an association analysis, corroborated by an excess of ancestry when comparing sprint and distance dogs. A comparison of elite and poor-performing sled dogs identified a single region significantly associated with heat tolerance. Within the region we identified seven SNPs within the myosin heavy chain 9 gene (MYH9) that were significantly associated with heat tolerance in sprint dogs, two of which correspond to conserved promoter and enhancer regions in the human ortholog.     

Genome-wide association analyses of common wheat (Triticum aestivum L.) germplasm identifies multiple loci for aluminium resistance

Aluminium (Al3+) toxicity restricts productivity and profitability of wheat (Triticum aestivum L.) crops grown on acid soils worldwide. Continued gains will be obtained by identifying superior alleles and novel Al3+ resistance loci that can be incorporated into breeding programs. We used association mapping to identify genomic regions associated with Al3+ resistance using 1055 accessions of common wheat from different geographic regions of the world and 178 polymorphic diversity arrays technology (DArT) markers. Bayesian analyses based on genetic distance matrices classified these accessions into 12 subgroups. Genome-wide association analyses detected markers that were significantly associated with Al3+ resistance on chromosomes 1A, 1B, 2A, 2B, 2D, 3A, 3B, 4A, 4B, 4D, 5B, 6A, 6B, 7A, and 7B. Some of these genomic regions correspond to previously identified loci for Al3+ resistance, whereas others appear to be novel. Among the markers targeting TaALMT1 (the major Al3+-resistance gene located on chromosome 4D), those that detected alleles in the promoter explained most of the phenotypic variance for Al3+ resistance, which is consistent with this region controlling the level of TaALMT1 expression. These results demonstrate that genome-wide association mapping cannot only confirm known Al3+-resistance loci, such as those on chromosomes 4D and 4B, but they also highlight the utility of this technique in identifying novel resistance loci.     

Effect of genomic deficiencies on sexual size dimorphism through modification of developmental time in Drosophila melanogaster

Sexual size dimorphism (SSD), a difference in body size between sexes, is common in many taxa. In insects, females are larger than males in >70% of all taxa in most orders. The fruit fly, Drosophila melanogaster is one prominent model organism to investigate SSD since its clear and representative female-biased SSD and its growth regulation are well studied. Elucidating the number and nature of genetic elements that can potentially influence SSD would be helpful in understanding the evolutionary potential of SSD. Here, we investigated the SSD pattern caused by artificially introduced genetic variation in D. melanogaster, and examined whether variation in SSD was mediated by the sex-specific modification of developmental time. To map the genomic regions that had effects on sexual wing size and/or developmental time differences (SDtD), we reanalyzed previously published genome-wide deficiency mapping data to evaluate the effects of 376 isogenic deficiencies covering a total of ~67% of the genomic regions of the second and third chromosomes of D. melanogaster. We found genetic variation in SSD and SDtD generated by genomic deficiencies, and a negative genetic correlation between size and development time. We also found SSD and SDtD allometries that are not qualitatively congruent, which however overall at best only partly help in explaining the patterns found. We identified several genomic deficiencies with the tendency to either exaggerate or suppress SSD, in agreement with quantitative genetic null expectations of many loci with small effects. These novel findings contribute to a better understanding of the evolutionary potential of sexual dimorphism.     

Mapping QTLs associated with drought resistance in sorghum (Sorghum bicolor L. Moench)

Drought is a major abiotic stress factor limiting crop production. Identification of genetic factors involved in plant responses to drought stress will provide a solid foundation to improve drought resistance. Sorghum is well adapted to hot dry environments and regarded as a model for studying drought resistance among the grasses. Significant progress in genome mapping of this crop has also been made. In sorghum, rapid premature leaf death generally occurs when water is limited during the grain filling period. Premature leaf senescence, in turn, leads to charcoal rot, stalk lodging, and significant yield loss. More than 80% of commercial sorghum hybrids in the United States are grown under non-irrigated conditions and although most of them have pre-flowering drought resistance, many do not have any significant post-flowering drought resistance. Stay-green is one form of drought resistance mechanism, which gives sorghum resistance to premature senescence under soil moisture stress during the post-flowering period. Quantitative trait locus (QTL) studies with recombinant inbred lines (RILs) and near-isogenic lines (NILs) identified several genomic regions associated with resistance to pre-flowering and post-flowering drought stress. We have identified four genomic regions associated with the stay-green trait using a RIL population developed from B35 × Tx7000. These four major stay-green QTLs were consistently identified in all field trials and accounted for 53.5% of the phenotypic variance. We review the progress in mapping stay-green QTLs as a component of drought resistance in sorghum. The molecular genetic dissection of the QTLs affecting stay-green will provide further opportunities to elucidate the underlying physiological mechanisms involved in drought resistance in sorghum and other grasses.     

A powerful regression-based method for admixture mapping of isolation across the genome of hybrids

We propose a novel method that uses natural admixture between divergent lineages (hybridization) to investigate the genetic architecture of reproductive isolation and adaptive introgression. Our method employs multinomial regression to estimate genomic clines and to quantify introgression for individual loci relative to the genomic background (clines in genotype frequency along a genomic admixture gradient). Loci with patterns of introgression that deviate significantly from null expectations based on the remainder of the genome are potentially subject to selection and thus of interest to understanding adaptation and the evolution of reproductive isolation. Using simulations, we show that different forms of selection modify these genomic clines in predictable ways and that our method has good power to detect moderate to strong selection for multiple forms of selection. Using individual-based simulations, we demonstrate that our method generally has a low false positive rate, except when genetic drift is particularly pronounced (e.g. low population size, low migration rates from parental populations, and substantial time since initial admixture). Additional individual-based simulations reveal that moderate selection against heterozygotes can be detected as much as 50 c m away from the focal locus directly experiencing selection, but is not detected at unlinked loci. Finally, we apply our analytical method to previously published data sets from a mouse ( Mus musculus and M. domesticus ) and two sunflower ( Helianthus petiolaris and H. annuus ) hybrid zones. This method should be applicable to numerous species that are currently the focus of research in evolution and ecology and should help bring about new insights regarding the processes underlying the origin and maintenance of biological diversity.