Improving the accuracy and efficiency of identity-by-descent detection in population data

Genome-wide variation in the Angolan Namib Desert reveals unique pre-Bantu ancestry

Ancient DNA studies reveal the genetic structure of Africa before the expansion of Bantu-speaking agriculturalists; however, the impact of now extinct hunter-gatherer and herder societies on the genetic makeup of present-day African groups remains elusive. Here, we uncover the genetic legacy of pre-Bantu populations from the Angolan Namib Desert, where we located small-scale groups associated with enigmatic forager traditions, as well as the last speakers of the Khoe-Kwadi family's Kwadi branch. By applying an ancestry decomposition approach to genome-wide data from these and other African populations, we reconstructed the fine-scale histories of contact emerging from the migration of Khoe-Kwadi-speaking pastoralists and identified a deeply divergent ancestry, which is exclusively shared between groups from the Angolan Namib and adjacent areas of Namibia. The unique genetic heritage of the Namib peoples shows how modern DNA research targeting understudied regions of high ethnolinguistic diversity can complement ancient DNA studies in probing the deep genetic structure of the African continent.

Recombination map tailored to Native Hawaiians improves robustness of genomic scans for positive selection

Recombination events establish the patterns of haplotypic structure in a population and estimates of recombination rates are used in several downstream population and statistical genetic analyses. Using suboptimal maps from distantly related populations may reduce the efficacy of genomic analyses, particularly for underrepresented populations such as the Native Hawaiians. To overcome this challenge, we constructed recombination maps using genome-wide array data from two study samples of Native Hawaiians: one reflecting the current admixed state of Native Hawaiians (NH map), and one based on individuals of enriched Polynesian ancestries (PNS map) with the potential to be used for less admixed Polynesian populations such as the Samoans. We found the recombination landscape to be less correlated with those from other continental populations (e.g. Spearman's rho = 0.79 between PNS and CEU (Utah residents with Northern and Western European ancestry) compared to 0.92 between YRI (Yoruba in Ibadan, Nigeria) and CEU at 50 kb resolution), likely driven by the unique demographic history of the Native Hawaiians. PNS also shared the fewest recombination hotspots with other populations (e.g. 8% of hotspots shared between PNS and CEU compared to 27% of hotspots shared between YRI and CEU). We found that downstream analyses in the Native Hawaiian population, such as local ancestry inference, imputation, and IBD segment and relatedness detections, would achieve similar efficacy when using the NH map compared to an omnibus map. However, for genome scans of adaptive loci using integrated haplotype scores, we found several loci with apparent genome-wide significant signals (|Z-score| > 4) in Native Hawaiians that would not have been significant when analyzed using NH-specific maps. Population-specific recombination maps may therefore improve the robustness of haplotype-based statistics and help us better characterize the evolutionary history that may underlie Native Hawaiian-specific health conditions that persist today.

A computational approach for positive genetic identification and relatedness detection from low-coverage shotgun sequencing data

Several methods exist for detecting genetic relatedness or identity by comparing DNA information. These methods generally require genotype calls, either single-nucleotide polymorphisms or short tandem repeats, at the sites used for comparison. For some DNA samples, like those obtained from bone fragments or single rootless hairs, there is often not enough DNA present to generate genotype calls that are accurate and complete enough for these comparisons. Here, we describe IBDGem, a fast and robust computational procedure for detecting genomic regions of identity-by-descent by comparing low-coverage shotgun sequence data against genotype calls from a known query individual. At less than 1× genome coverage, IBDGem reliably detects segments of relatedness and can make high-confidence identity detections with as little as 0.01× genome coverage.

The genomic history of the indigenous people of the Canary Islands

The indigenous population of the Canary Islands, which colonized the archipelago around the 3rd century CE, provides both a window into the past of North Africa and a unique model to explore the effects of insularity. We generate genome-wide data from 40 individuals from the seven islands, dated between the 3rd-16rd centuries CE. Along with components already present in Moroccan Neolithic populations, the Canarian natives show signatures related to Bronze Age expansions in Eurasia and trans-Saharan migrations. The lack of gene flow between islands and constant or decreasing effective population sizes suggest that populations were isolated. While some island populations maintained relatively high genetic diversity, with the only detected bottleneck coinciding with the colonization time, other islands with fewer natural resources show the effects of insularity and isolation. Finally, consistent genetic differentiation between eastern and western islands points to a more complex colonization process than previously thought.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

Multiple Origins and Genomic Basis of Complex Traits in Sighthounds

Sighthounds, a distinctive group of hounds comprising numerous breeds, have their origins rooted in ancient artificial selection of dogs. In this study, we performed genome sequencing for 123 sighthounds, including one breed from Africa, six breeds from Europe, two breeds from Russia, and four breeds and 12 village dogs from the Middle East. We gathered public genome data of five sighthounds and 98 other dogs as well as 31 gray wolves to pinpoint the origin and genes influencing the morphology of the sighthound genome. Population genomic analysis suggested that sighthounds originated from native dogs independently and were comprehensively admixed among breeds, supporting the multiple origins hypothesis of sighthounds. An additional 67 published ancient wolf genomes were added for gene flow detection. Results showed dramatic admixture of ancient wolves in African sighthounds, even more than with modern wolves. Whole-genome scan analysis identified 17 positively selected genes (PSGs) in the African population, 27 PSGs in the European population, and 54 PSGs in the Middle Eastern population. None of the PSGs overlapped in the three populations. Pooled PSGs of the three populations were significantly enriched in "regulation of release of sequestered calcium ion into cytosol" (gene ontology: 0051279), which is related to blood circulation and heart contraction. In addition, ESR1, JAK2, ADRB1, PRKCE, and CAMK2D were under positive selection in all three selected groups. This suggests that different PSGs in the same pathway contributed to the similar phenotype of sighthounds. We identified an ESR1 mutation (chr1: g.42,177,149 T > C) in the transcription factor (TF) binding site of Stat5a and a JAK2 mutation (chr1: g.93,277,007 T > A) in the TF binding site of Sox5. Functional experiments confirmed that the ESR1 and JAK2 mutation reduced their expression. Our results provide new insights into the domestication history and genomic basis of sighthounds.

Inferring the population history of Kyrgyz in Xinjiang, Northwest China from genome-wide array genotyping

OBJECTIVES

Xinjiang plays a vital role in the trans-Eurasian population migration, language diffusion, and culture and technology exchange. However, the underrepresentation of Xinjiang's genomes has hindered a more comprehensive understanding of Xinjiang's genetic structure and population history.

MATERIALS AND METHODS

We collected and genotyped 70 southern Xinjiang's Kyrgyz (SXJK) individuals and combined the data with modern and ancient Eurasians published. We used allele-frequency methods, including PCA, ADMIXTURE, f-statistics, qpWave/qpAdm, ALDER, Treemix, and haplotype-shared methods including shared-IBD segments, fineSTRUCTURE, and GLOBETROTTER to unveil the fine-scale population structure and reconstruct admixture history.

RESULTS

We identified genetic substructure within the SXJK population with subgroups showing different genetic affinities to West and East Eurasians. All SXJK subgroups were suggested to have close genetic relationships with surrounding Turkic-speaking groups that is, Uyghur, Kyrgyz from north Xinjiang and Tajikistan, and Chinese Kazakh, suggesting a shared ancestry among those populations. Outgroup-f3 and symmetrical f4 statistics showed a high genetic affinity of SXJK to present-day Tungusic, Mongolic-speaking populations and Ancient Northeast Asian (ANA) related groups. Allele sharing and haplotype sharing profiles revealed the east-west admixture pattern of SXJK. The qpAdm-based admixture models showed that SXJK derived ancestry from East Eurasian (ANA and East Asian, 42.7%-83.3%) and West Eurasian (Western Steppe herders and Central Asian, 16.7%-57.3%), the recent east-west admixture event could be traced to 1000 years ago based on ALDER and GLOBETROTTER analysis.

DISCUSSION

The high genetic affinity of SXJK to present-day Tungusic and Mongolic-speaking populations and short-shared IBD segments indicated their shared common ancestry. SXJK harbored a close genetic affinity to ANA-related populations, indicating the Northeast Asian origin of SXJK. The West and East Eurasian admixture models observed in SXJK further provided evidence of the dynamic admixture history in Xinjiang. The east-west admixture pattern and the identified ancestral makeup of SXJK suggested a genetic continuity from some Iron Age Xinjiang populations to present-day SXJK.

Drivers of genomic landscapes of differentiation across a Populus divergence gradient

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

Fast inference of genetic recombination rates in biobank scale data

Although rates of recombination events across the genome (genetic maps) are fundamental to genetic research, the majority of current studies only use one standard map. There is evidence suggesting population differences in genetic maps, and thus estimating population-specific maps, are of interest. Although the recent availability of biobank-scale data offers such opportunities, current methods are not efficient at leveraging very large sample sizes. The most accurate methods are still linkage disequilibrium (LD)-based methods that are only tractable for a few hundred samples. In this work, we propose a fast and memory-efficient method for estimating genetic maps from population genotyping data. Our method, FastRecomb, leverages the efficient positional Burrows-Wheeler transform (PBWT) data structure for counting IBD segment boundaries as potential recombination events. We used PBWT blocks to avoid redundant counting of pairwise matches. Moreover, we used a panel-smoothing technique to reduce the noise from errors and recent mutations. Using simulation, we found that FastRecomb achieves state-of-the-art performance at 10-kb resolution, in terms of correlation coefficients between the estimated map and the ground truth. This is mainly because FastRecomb can effectively take advantage of large panels comprising more than hundreds of thousands of haplotypes. At the same time, other methods lack the efficiency to handle such data. We believe further refinement of FastRecomb would deliver more accurate genetic maps for the genetics community.

Extensive genetic admixture between Tai-Kadai-speaking people and their neighbours in the northeastern region of the Yungui Plateau inferred from genome-wide variations

BACKGROUND

Yungui Plateau in Southwest China is characterized by multi-language and multi-ethnic communities and is one of the regions with the wealthiest ethnolinguistic, cultural and genetic diversity in East Asia. There are numerous Tai-Kadai (TK)-speaking populations, but their detailed evolutionary history and biological adaptations are still unclear.

RESULTS

Here, we genotyped genome-wide SNP data of 77 unrelated TK-speaking Zhuang and Dong individuals from the Yungui Plateau and explored their detailed admixture history and adaptive features using clustering patterns, allele frequency differentiation and sharing haplotype patterns. TK-speaking Zhuang and Dong people in Guizhou are closely related to geographically close TK and Hmong-Mien (HM)-speaking populations. Besides, we identified that Guizhou TK-speaking people have a close genetic relationship with Austronesian (AN)-speaking Atayal and Paiwan people, which is supported by the common origin of the ancient Baiyue tribe. We additionally found subtle genetic differences among the newly studied TK people and previously reported Dais via the fine-scale genetic substructure analysis based on the shared haplotype chunks. Finally, we identified specific selection candidate signatures associated with several essential human immune systems and neurological disorders, which could provide evolutionary evidence for the allele frequency distribution pattern of genetic risk loci.

CONCLUSIONS

Our comprehensive genetic characterization of TK people suggested the strong genetic affinity within TK groups and extensive gene flow with geographically close HM and Han people. We also provided genetic evidence that supported the common origin hypothesis of TK and AN people. The best-fitted admixture models further suggested that ancestral sources from northern millet farmers and southern inland and coastal people contributed to the formation of the gene pool of the Zhuang and Dong people.

Native American genetic ancestry and pigmentation allele contributions to skin color in a Caribbean population

Our interest in the genetic basis of skin color variation between populations led us to seek a Native American population with genetically African admixture but low frequency of European light skin alleles. Analysis of 458 genomes from individuals residing in the Kalinago Territory of the Commonwealth of Dominica showed approximately 55% Native American, 32% African, and 12% European genetic ancestry, the highest Native American genetic ancestry among Caribbean populations to date. Skin pigmentation ranged from 20 to 80 melanin units, averaging 46. Three albino individuals were determined to be homozygous for a causative multi-nucleotide polymorphism OCA2NW273KV contained within a haplotype of African origin; its allele frequency was 0.03 and single allele effect size was -8 melanin units. Derived allele frequencies of SLC24A5A111T and SLC45A2L374F were 0.14 and 0.06, with single allele effect sizes of -6 and -4, respectively. Native American genetic ancestry by itself reduced pigmentation by more than 20 melanin units (range 24-29). The responsible hypopigmenting genetic variants remain to be identified, since none of the published polymorphisms predicted in prior literature to affect skin color in Native Americans caused detectable hypopigmentation in the Kalinago.

The genetic history of the Southern Andes from present-day Mapuche ancestry

The southernmost regions of South America harbor some of the earliest evidence of human presence in the Americas. However, connections with the rest of the continent and the contextualization of present-day indigenous ancestries remain poorly resolved. In this study, we analyze the genetic ancestry of one of the largest indigenous groups in South America: the Mapuche. We generate genome-wide data from 64 participants from three Mapuche populations in Southern Chile: Pehuenche, Lafkenche, and Huilliche. Broadly, we describe three main ancestry blocks with a common origin, which characterize the Southern Cone, the Central Andes, and Amazonia. Within the Southern Cone, ancestors of the Mapuche lineages differentiated from those of the Far South during the Middle Holocene and did not experience further migration waves from the north. We find that the deep genetic split between the Central and Southern Andes is followed by instances of gene flow, which may have accompanied the southward spread of cultural traits from the Central Andes, including crops and loanwords from Quechua into Mapudungun (the language of the Mapuche). Finally, we report close genetic relatedness between the three populations analyzed, with the Huilliche characterized additionally by intense recent exchanges with the Far South. Our findings add new perspectives on the genetic (pre)history of South America, from the first settlement through to the present-day indigenous presence. Follow-up fieldwork took these results back to the indigenous communities to contextualize the genetic narrative alongside indigenous knowledge and perspectives.

Addressing the feasibility of people of African descent finding living African relatives using direct-to-consumer genetic testing

People of African descent use direct-to-consumer genomics services such as 23andMe and AncestryDNA for various family histories and health reasons, including identifying and interacting with the previously unknown living African genetic relatives. In this commentary, I argue that it is reasonable to consider that cousin pairs consisting of an African person and a descendant of an African person enslaved in the Americans during the Transatlantic Slave Trade (i.e., a person of African descent) have genealogical ancestors recent enough to be detected using autosomal DNA testing where the pair has shared ancestors in the range of 20-6 generations ago from the present.

Dispersal inference from population genetic variation using a convolutional neural network

The geographic nature of biological dispersal shapes patterns of genetic variation over landscapes, making it possible to infer properties of dispersal from genetic variation data. Here, we present an inference tool that uses geographically distributed genotype data in combination with a convolutional neural network to estimate a critical population parameter: the mean per-generation dispersal distance. Using extensive simulation, we show that our deep learning approach is competitive with or outperforms state-of-the-art methods, particularly at small sample sizes. In addition, we evaluate varying nuisance parameters during training-including population density, demographic history, habitat size, and sampling area-and show that this strategy is effective for estimating dispersal distance when other model parameters are unknown. Whereas competing methods depend on information about local population density or accurate inference of identity-by-descent tracts, our method uses only single-nucleotide-polymorphism data and the spatial scale of sampling as input. Strikingly, and unlike other methods, our method does not use the geographic coordinates of the genotyped individuals. These features make our method, which we call "disperseNN," a potentially valuable new tool for estimating dispersal distance in nonmodel systems with whole genome data or reduced representation data. We apply disperseNN to 12 different species with publicly available data, yielding reasonable estimates for most species. Importantly, our method estimated consistently larger dispersal distances than mark-recapture calculations in the same species, which may be due to the limited geographic sampling area covered by some mark-recapture studies. Thus genetic tools like ours complement direct methods for improving our understanding of dispersal.

Cycles of satellite and transposon evolution in Arabidopsis centromeres

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation^1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox^3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

Fast integration and accumulation of beneficial breeding alleles through an AB-NAMIC strategy in wheat

Wheat (Triticum aestivum) is among the most important staple crops for safeguarding the food security of the growing world population. To bridge the gap between genebank diversity and breeding programs, we developed an advanced backcross-nested association mapping plus inter-crossed population (AB-NAMIC) by crossing three popular wheat cultivars as recurrent founders to 20 germplasm lines from a mini core collection. Selective backcrossing combined with selection against undesirable traits and extensive crossing within and between sub-populations created new opportunities to detect unknown genes and increase the frequency of beneficial alleles in the AB-NAMIC population. We performed phenotyping of 590 AB-NAMIC lines and a natural panel of 476 cultivars for six consecutive growing seasons and genotyped these 1066 lines with a 660K SNP array. Genome-wide association studies of both panels for plant development and yield traits demonstrated improved power to detect rare alleles and loci with medium genetic effects in AB-NAMIC. Notably, genome-wide association studies in AB-NAMIC detected the candidate gene TaSWEET6-7B (TraesCS7B03G1216700), which has high homology to the rice SWEET6b gene and exerts strong effects on adaptation and yield traits. The commercial release of two derived AB-NAMIC lines attests to its direct applicability in wheat improvement. Valuable information on genome-wide association study mapping, candidate genes, and their haplotypes for breeding traits are available through WheatGAB. Our research provides an excellent framework for fast-tracking exploration and accumulation of beneficial alleles stored in genebanks.

The genomic diversity of Taiwanese Austronesian groups: Implications for the "Into- and Out-of-Taiwan" models

The origin and dispersal of the Austronesian language family, one of the largest and most widespread in the world, have long attracted the attention of linguists, archaeologists, and geneticists. Even though there is a growing consensus that Taiwan is the source of the spread of Austronesian languages, little is known about the migration patterns of the early Austronesians who settled in and left Taiwan, i.e. the "Into-Taiwan" and "out-of-Taiwan" events. In particular, the genetic diversity and structure within Taiwan and how this relates to the into-/out-of-Taiwan events are largely unexplored, primarily because most genomic studies have largely utilized data from just two of the 16 recognized Highland Austronesian groups in Taiwan. In this study, we generated the largest genome-wide data set of Taiwanese Austronesians to date, including six Highland groups and one Lowland group from across the island and two Taiwanese Han groups. We identified fine-scale genomic structure in Taiwan, inferred the ancestry profile of the ancestors of Austronesians, and found that the southern Taiwanese Austronesians show excess genetic affinities with the Austronesians outside of Taiwan. Our findings thus shed new light on the Into- and Out-of-Taiwan dispersals.

The Newfoundland and Labrador mosaic founder population descends from an Irish and British diaspora from 300 years ago

The founder population of Newfoundland and Labrador (NL) is a unique genetic resource, in part due to its geographic and cultural isolation, where historical records describe a migration of European settlers, primarily from Ireland and England, to NL in the 18th and 19th centuries. Whilst its historical isolation, and increased prevalence of certain monogenic disorders are well appreciated, details of the fine-scale genetic structure and ancestry of the population are lacking. Understanding the genetic origins and background of functional, disease causing, genetic variants would aid genetic mapping efforts in the Province. Here, we leverage dense genome-wide SNP data on 1,807 NL individuals to reveal fine-scale genetic structure in NL that is clustered around coastal communities and correlated with Christian denomination. We show that the majority of NL European ancestry can be traced back to the south-east and south-west of Ireland and England, respectively. We date a substantial population size bottleneck approximately 10-15 generations ago in NL, associated with increased haplotype sharing and autozygosity. Our results reveal insights into the population history of NL and demonstrate evidence of a population conducive to further genetic studies and biomarker discovery.

Exome sequencing of affected duos and trios uncovers PRUNE2 as a novel prostate cancer predisposition gene

BACKGROUND

Prostate cancer (PrCa) is one of the most hereditable human cancers, however, only a small fraction of patients has been shown to carry deleterious variants in known cancer predisposition genes.

METHODS

Whole-exome sequencing was performed in multiple affected members of 45 PrCa families to select the best candidate genes behind part of the PrCa missing hereditability. Recurrently mutated genes were prioritised, and further investigated by targeted next-generation sequencing in the whole early-onset and/or familial PrCa series of 462 patients.

RESULTS

PRUNE2 stood out from our analysis when also considering the available data on its association with PrCa development. Ten germline pathogenic/likely pathogenic variants in the PRUNE2 gene were identified in 13 patients. The most frequent variant was found in three unrelated patients and identical-by-descent analysis revealed that the haplotype associated with the variant is shared by all the variant carriers, supporting the existence of a common ancestor.

DISCUSSION

This is the first report of pathogenic/likely pathogenic germline variants in PRUNE2 in PrCa patients, namely in those with early-onset/familial disease. Importantly, PRUNE2 was the most frequently mutated gene in the whole series, with a deleterious germline variant identified in 2.8% of the patients, representing a novel prostate cancer predisposition gene.

Mating allocations in Holstein combining genomic information and linear programming optimization at the herd level

In this study, we explored mating allocation in Holstein using genomic information for 24,333 Holstein females born in Denmark, Finland, and Sweden. We used 2 data sets of bulls: the top 50 genotyped bulls and the top 25 polled genotyped bulls on the Nordic total merit scale. We used linear programming to optimize economic scores within each herd, considering genetic level, genetic relationship, semen cost, the economic impact of genetic defects, polledness, and β-casein. We found that it was possible to reduce genetic relationships and eliminate expression of genetic defects with minimal effect on the genetic level in total merit index. Compared with maximizing only Nordic total merit index, the relative frequency of polled offspring increased from 13.5 to 22.5%, and that of offspring homozygous for β-casein (A2A2) from 66.7 to 75.0% in one generation, without any substantial negative impact on other comparison criteria. Using only semen from polled bulls, which might become necessary if dehorning is banned, considerably reduced the genetic level. We also found that animals carrying the polled allele were less likely to be homozygous for β-casein (A2A2) and more likely to be carriers of the genetic defect HH1. Hence, adding economic value to a monogenic trait in the economic score used for mating allocation sometimes negatively affected another monogenetic trait. We recommend that the comparison criteria used in this study be monitored in a modern genomic mating program.

The multifaceted genomic history of Ashaninka from Amazonian Peru

Despite its crucial location, the western side of Amazonia between the Andes and the source(s) of the Amazon River is still understudied from a genomic and archaeogenomic point of view, albeit possibly harboring essential information to clarify the complex genetic history of local Indigenous groups and their interactions with nearby regions,^{1,2,3,4,5,6,7,8} including central America and the Caribbean.^9,10,11,12 Focusing on this key region, we analyzed the genome-wide profiles of 51 Ashaninka individuals from Amazonian Peru, observing an unexpected extent of genomic variation. We identified at least two Ashaninka subgroups with distinctive genomic makeups, which were differentially shaped by the degree and timing of external admixtures, especially with the Indigenous groups from the Andes and the Pacific coast. On a continental scale, Ashaninka ancestors probably derived from a south-north migration of Indigenous groups moving into the Amazonian rainforest from a southeastern area with contributions from the Southern Cone and the Atlantic coast. These ancestral populations diversified in the variegated geographic regions of interior South America, on the eastern side of the Andes, differentially interacting with surrounding coastal groups. In this complex scenario, we also revealed strict connections between the ancestors of present-day Ashaninka, who belong to the Arawakan language family,¹³ and those Indigenous groups that moved further north into the Caribbean, contributing to the early Ceramic (Saladoid) tradition in the islands.^14,15.

Deciphering the genetic basis of resistance to soybean cyst nematode combining IBD and association mapping

IBD analysis clarified the dynamics of chromosomal recombination during the ZP pedigree breeding process and identified ten genomic regions resistant to SCN race3 combining association mapping. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating pathogens for soybean production worldwide. The cultivar Zhongpin03-5373 (ZP), derived from SCN-resistant progenitor parents, Peking, PI 437654 and Huipizhi Heidou, is an elite line with high resistance to SCN race3. In the current study, a pedigree variation map was generated for ZP and its ten progenitors using 3,025,264 high-quality SNPs identified from an average of 16.2 × re-sequencing for each genome. Through identity by decent (IBD) tracking, we showed the dynamic change of genome and detected important IBD fragments, which revealed the comprehensively artificial selection of important traits during ZP breeding process. A total of 2,353 IBD fragments related to SCN resistance including SCN-resistant genes rhg1, rhg4 and NSF_RAN07 were identified based on the resistant-related genetic paths. Moreover, 23 genomic regions underlying resistance to SCN race3 were identified by genome-wide association study (GWAS) in 481 re-sequenced cultivated soybeans. Ten common loci were found by both IBD tracking and GWAS analysis. Haplotype analysis of 16 potential candidate genes suggested a causative SNP (C/T, - 1065) located in the promoter of Glyma.08G096500 and encoding a predicted TIFY5b-related protein on chr8 was highly correlated with SCN race3 resistance. Our results more thoroughly elucidated the dynamics of genomic fragments during ZP pedigree breeding and the genetic basis of SCN resistance, which will provide useful information for gene cloning and the development of resistant soybean cultivars using a marker-assisted selection approach.

ancIBD - Screening for identity by descent segments in human ancient DNA

Long DNA sequences shared between two individuals, known as Identical by descent (IBD) segments, are a powerful signal for identifying close and distant biological relatives because they only arise when the pair shares a recent common ancestor. Existing methods to call IBD segments between present-day genomes cannot be straightforwardly applied to ancient DNA data (aDNA) due to typically low coverage and high genotyping error rates. We present ancIBD, a method to identify IBD segments for human aDNA data implemented as a Python package. Our approach is based on a Hidden Markov Model, using as input genotype probabilities imputed based on a modern reference panel of genomic variation. Through simulation and downsampling experiments, we demonstrate that ancIBD robustly identifies IBD segments longer than 8 centimorgan for aDNA data with at least either 0.25x average whole-genome sequencing (WGS) coverage depth or at least 1x average depth for in-solution enrichment experiments targeting a widely used aDNA SNP set ('1240k'). This application range allows us to screen a substantial fraction of the aDNA record for IBD segments and we showcase two downstream applications. First, leveraging the fact that biological relatives up to the sixth degree are expected to share multiple long IBD segments, we identify relatives between 10,156 ancient Eurasian individuals and document evidence of long-distance migration, for example by identifying a pair of two approximately fifth-degree relatives who were buried 1410km apart in Central Asia 5000 years ago. Second, by applying ancIBD, we reveal new details regarding the spread of ancestry related to Steppe pastoralists into Europe starting 5000 years ago. We find that the first individuals in Central and Northern Europe carrying high amounts of Steppe-ancestry, associated with the Corded Ware culture, share high rates of long IBD (12-25 cM) with Yamnaya herders of the Pontic-Caspian steppe, signaling a strong bottleneck and a recent biological connection on the order of only few hundred years, providing evidence that the Yamnaya themselves are a main source of Steppe ancestry in Corded Ware people. We also detect elevated sharing of long IBD segments between Corded Ware individuals and people associated with the Globular Amphora culture (GAC) from Poland and Ukraine, who were Copper Age farmers not yet carrying Steppe-like ancestry. These IBD links appear for all Corded Ware groups in our analysis, indicating that individuals related to GAC contexts must have had a major demographic impact early on in the genetic admixtures giving rise to various Corded Ware groups across Europe. These results show that detecting IBD segments in aDNA can generate new insights both on a small scale, relevant to understanding the life stories of people, and on the macroscale, relevant to large-scale cultural-historical events.

The dogs of Chernobyl: Demographic insights into populations inhabiting the nuclear exclusion zone

The 1986 Chernobyl nuclear disaster initiated a series of catastrophic events resulting in long-term and widespread environmental contamination. We characterize the genetic structure of 302 dogs representing three free-roaming dog populations living within the power plant itself, as well as those 15 to 45 kilometers from the disaster site. Genome-wide profiles from Chernobyl, purebred and free-breeding dogs, worldwide reveal that the individuals from the power plant and Chernobyl City are genetically distinct, with the former displaying increased intrapopulation genetic similarity and differentiation. Analysis of shared ancestral genome segments highlights differences in the extent and timing of western breed introgression. Kinship analysis reveals 15 families, with the largest spanning all collection sites within the radioactive exclusion zone, reflecting migration of dogs between the power plant and Chernobyl City. This study presents the first characterization of a domestic species in Chernobyl, establishing their importance for genetic studies into the effects of exposure to long-term, low-dose ionizing radiation.

The Legacy of Infectious Disease Exposure on the Genomic Diversity of Indigenous Southern Mexicans

To characterize host risk factors for infectious disease in Mesoamerican populations, we interrogated 857,481 SNPs assayed using the Affymetrix 6.0 genotyping array for signatures of natural selection in immune response genes. We applied three statistical tests to identify signatures of natural selection: locus-specific branch length (LSBL), the cross-population extended haplotype homozygosity (XP-EHH), and the integrated haplotype score (iHS). Each of the haplotype tests (XP-EHH and iHS) were paired with LSBL and significance was determined at the 1% level. For the paired analyses, we identified 95 statistically significant windows for XP-EHH/LSBL and 63 statistically significant windows for iHS/LSBL. Among our top immune response loci, we found evidence of recent directional selection associated with the major histocompatibility complex (MHC) and the peroxisome proliferator-activated receptor gamma (PPAR-γ) signaling pathway. These findings illustrate that Mesoamerican populations' immunity has been shaped by exposure to infectious disease. As targets of selection, these variants are likely to encode phenotypes that manifest themselves physiologically and therefore may contribute to population-level variation in immune response. Our results shed light on past selective events influencing the host response to modern diseases, both pathogenic infection as well as autoimmune disorders.

Genome-wide association study of 17 serum biochemical indicators in a chicken F2 resource population

BACKGROUND

Serum biochemical indicators are often regarded as direct reflections of animal metabolism and health. The molecular mechanisms underlying serum biochemical indicators metabolism of chicken (Gallus Gallus) have not been elucidated. Herein, we performed a genome-wide association study (GWAS) to identify the variation associated with serum biochemical indicators. The aim of this research was to broaden the understanding of the serum biochemical indicators in chickens.

RESULTS

A GWAS of serum biochemical indicators was carried out on 734 samples from an F2 Gushi× Anka chicken population. All chickens were genotyped by sequencing, 734 chickens and 321,314 variants were obtained after quality control. Based on these variants, a total of 236 single-nucleotide polymorphisms (SNPs) on 9 chicken chromosomes (GGAs) were identified to be significantly (-log10(P) > 5.72) associated with eight of seventeen serum biochemical indicators. Ten novel quantitative trait locis (QTLs) were identified for the 8 serum biochemical indicator traits of the F2 population. Literature mining revealed that the ALPL, BCHE, GGT2/GGT5 genes at loci GGA24, GGA9 and GGA15 might affect the alkaline phosphatase (AKP), cholinesterase (CHE) and γ-glutamyl transpeptidase (GGT) traits, respectively.

CONCLUSION

The findings of the present study may contribute to a better understanding of the molecular mechanisms of chicken serum biochemical indicator regulation and provide a theoretical basis for chicken breeding programs.

On the origin and structure of haplotype blocks

The term "haplotype block" is commonly used in the developing field of haplotype-based inference methods. We argue that the term should be defined based on the structure of the Ancestral Recombination Graph (ARG), which contains complete information on the ancestry of a sample. We use simulated examples to demonstrate key features of the relationship between haplotype blocks and ancestral structure, emphasizing the stochasticity of the processes that generate them. Even the simplest cases of neutrality or of a "hard" selective sweep produce a rich structure, often missed by commonly used statistics. We highlight a number of novel methods for inferring haplotype structure, based on the full ARG, or on a sequence of trees, and illustrate how they can be used to define haplotype blocks using an empirical data set. While the advent of new, computationally efficient methods makes it possible to apply these concepts broadly, they (and additional new methods) could benefit from adding features to explore haplotype blocks, as we define them. Understanding and applying the concept of the haplotype block will be essential to fully exploit long and linked-read sequencing technologies.

Blocks identical by descent in the genomes of the indigenous population of Siberia demonstrate genetic links between populations

The gene pool of the indigenous population of Siberia is a unique system for studying population and evolutionary genetic processes, analyzing genetic diversity, and reconstructing the genetic history of populations. High ethnic diversity is a feature of Siberia, as one of the regions of the peripheral settlement of modern human. The vast expanses of this region and the small number of aboriginal populations contributed to the formation of significant territorial and genetic subdivision. About 40 indigenous peoples are settled on the territory of the Siberian historical and ethnographic province. Within the framework of this work, a large-scale population study of the gene pool of the indigenous peoples of Siberia was carried out for the first time at the level of high-density biochips. This makes it possible to fill in a significant gap in the genogeographic picture of the Eurasian population. For this, DNA fragments were analyzed, which had been inherited without recombination by each pair of individuals from their recent common ancestor, that is, segments (blocks) identical by descent (IBD). The distribution of IBD blocks in the populations of Siberia is in good agreement with the geographical proximity of the populations and their linguistic affiliation. Among the Siberian populations, the Chukchi, Koryaks, and Nivkhs form a separate cluster from the main Siberian group, with the Chukchi and Koryaks being more closely related. Separate subclusters of Evenks and Yakuts, Kets and Chulyms are formed within the Siberian cluster. Analysis of SNPs that fell into more IBD segments of the analyzed populations made it possible to compile a list of 5358 genes. According to the calculation results, biological processes enriched with these genes are associated with the detection of a chemical stimulus involved in the sensory perception of smell. Enriched for the genes found, molecular pathways are associated with the metabolism of linoleic, arachidonic, tyrosic acids and by olfactory transduction. At the same time, an analysis of the literature data showed that some of the selected genes, which were found in a larger number of IBD blocks in several populations at once, can play a role in genetic adaptation to environmental factors.

Relationship of the gene pool of the Khants with the peoples of Western Siberia, Cis-Urals and the Altai-Sayan Region according to the data on the polymorphism of autosomic locus and the Y-chromosome

Khanty are indigenous Siberian people living on the territory of Western Siberia, mainly on the territory of the Khanty-Mansiysk and Yamalo-Nenets Autonomous Okrugs. The present study is aimed at a comprehensive analysis of the structure of the Khanty gene pool and their comparison with other populations of the indigenous population of Southern and Western Siberia. To address the issues of genetic proximity of the Khanty with other indigenous peoples, we performed genotyping of a wide genomic set of autosomal markers using high-density biochips, as well as an expanded set of SNP and STR markers of the Y-chromosome in various ethnic groups: Khakas, Tuvans, Southern Altaians, Siberian Tatars, Chulyms (Turkic language family) and Kets (Yeniseian language family). The structure of the gene pool of the Khanty and other West Siberian and South Siberian populations was studied using a genome-wide panel of autosomal single nucleotide polymorphic markers and Y-chromosome markers. The results of the analysis of autosomal SNPs frequencies by various methods, the similarities in the composition of the Y-chromosome haplogroups and YSTR haplotypes indicate that the Khanty gene pool is quite specific. When analyzing autosomal SNPs, the Ugrian genetic component completely dominates in both samples (up to 99-100 %). The samples of the Khanty showed the maximum match in IBD blocks with each other, with a sample of the Kets, Chulyms, Tuvans, Tomsk Tatars, Khakas, Kachins, and Southern Altaians. The degree of coincidence of IBD blocks between the Khanty, Kets, and Tomsk Tatars is consistent with the results of the distribution of allele frequencies and common genetic components in these populations. According to the composition of the Y-chromosome haplogroups, the two samples of the Khanty differ significantly from each other. A detailed phylogenetic analysis of various Y-chromosome haplogroups made it possible to describe and clarify the differences in the phylogeny and structure of individual ethnospecific sublines, to determine their relationship, traces of population expansion in the Khanty gene pool. Variants of different haplogroups of the Y-chromosome in the Khanty, Khakas and Tuvans go back to their common ancestral lines. The results of a comparative analysis of male samples indicate a close genetic relationship between the Khanty and Nenets, Komi, Udmurts and Kets. The specificity of haplotypes, the discovery of various terminal SNPs confirms that the Khanty did not come into contact with other ethnic groups for a long time, except for the Nenets, which included many Khanty clans.

Genomic divergence and introgression among three Populus species

Genomic divergence with gene flow is very common in both plants and animals. However, divergence and gene flow are two counteracting factors during speciation. Identifying the types of genes that are likely to be introgressed and what genetic factors restrict further effective reproduction of interspecific hybrids is of great interest to biologists. We aimed to address these issues using three related tree species, Populus alba (Pa), P. tremula (Pt), and P. tremuloides (Ps), and the interspecific hybrid of the former two species, P. × canescens (Pc). We collected 105 genomes for these four poplar lineages, including 28 Pa, 38Pt, 21 Ps, and 18 Pc individuals, to reconstruct their evolutionary histories. Our coalescence-based simulations indicated that Pa diverged earliest from Ps and Pt, and asymmetrical gene flow existed between any two lineages, with especially large ancient gene flow occurring between Pa and Pt. The genomic landscape of divergence between pairs of the three species are highly heterogeneous, which may have arisen through both divergent sorting of ancient polymorphisms and ongoing gene flow. We found that extant regions of the genome with introgressed ancestry reduced genetic divergence but elevated recombination rates and accounted for 5.76 % of the total genome. Introgressed genes were functionally associated with stress resistance, including innate immune response, anti-adversity response, and programmed cell death. However, candidate genes underlying postmating barriers of Pc were homozygous and resistant to introgression due to the incompatibility of alleles between loci after hybridization and were associated with endosperm and gamete formation and disease resistance. Our study revealed genomic dynamics during speciation with gene flow and identified regions of the genome that were likely introgressed and adaptive as well as candidate loci responsible for hybrid incompatibility that resulted in the formation of postmating barriers after hybridization.

Structure and origin of Tuvan gene pool according to autosome SNP and Y-chromosome haplogroups

Tuvans are one of the most compactly living peoples of Southern Siberia, settled mainly in the territory of Tuva. The gene pool of the Tuvans is quite isolated, due to endogamy and a very low frequency of interethnic marriages. The structure of the gene pool of the Tuvans and other Siberian populations was studied using a genome-wide panel of autosomal single nucleotide polymorphic markers and Y-chromosome markers. The results of the analysis of the frequencies of autosomal SNPs by various methods, the similarities in the composition of the Y-chromosome haplogroups and YSTR haplotypes show that the gene pool of the Tuvans is very heterogeneous in terms of the composition of genetic components. It includes the ancient autochthonous Yeniseian component, which dominates among the Chulym Turks and Kets, the East Siberian component, which prevails among the Yakuts and Evenks, and the Far Eastern component, the frequency of which is maximum among the Nivkhs and Udeges. Analysis of the composition of IBD-blocks on autosomes shows the maximum genetic relationship of the Tuvans with the Southern Altaians, Khakas and Shors, who were formed during the settlement of the Turkic groups of populations on the territory of the Altai-Sayan region. A very diverse composition of the Tuvan gene pool is shown for various sublines of Y-chromosomal haplogroups, most of which show strong ethnic specificity. Phylogenetic analysis of individual Y-chromosome haplogroups demonstrates the maximum proximity of the gene pool of the Tuvans with the Altaians, Khakas and Shors. Differences in frequencies of Y-chromosome haplogroups between the Todzhans and Tuvans and a change in the frequencies of haplogroups from south to north associated with the East Asian component were found. The majority of the most frequent Y-chromosome haplogroups in the Tuvans demonstrate the founder effect, the formation age of which is fully consistent with the data on their ethnogenesis.

The social lives of isolates (and small language families): the case of the Northwest Amazon

The Americas are home to patches of extraordinary linguistic (genealogical) diversity. These high-diversity areas are particularly unexpected given the recent population of the Americas. In this paper, we zoom in on one such area, the Northwest Amazon, and address the question of how the diversity in this area has persisted to the present. We contrast two hypotheses that claim opposite mechanisms for the maintenance of diversity: the isolation hypothesis suggests that isolation facilitates the preservation of diversity, while the integration hypothesis proposes that conscious identity preservation in combination with contact drives diversity maintenance. We test predictions for both hypotheses across four disciplines: biogeography, cultural anthropology, population genetics and linguistics. Our results show signs of both isolation and integration, but they mainly suggest considerable diversity in how groups of speakers have interacted with their surroundings.

Variability and genetic merits of white Guinea yam landraces in Nigeria

INTRODUCTION

Landraces represent a significant gene pool of African cultivated white Guinea yam diversity. They could, therefore, serve as a potential donor of important traits such as resilience to stresses as well as food quality attributes that may be useful in modern yam breeding. This study assessed the pattern of genetic variability, quantitative trait loci (QTLs), alleles, and genetic merits of landraces, which could be exploited in breeding for more sustainable yam production in Africa.

METHODS

A total of 86 white Guinea yam landraces representing the popular landraces in Nigeria alongside 16 elite clones were used for this study. The yam landraces were genotyped using 4,819 DArTseq SNP markers and profiled using key productivity and food quality traits.

RESULTS AND DISCUSSION

Genetic population structure through admixture and hierarchical clustering methods revealed the presence of three major genetic groups. Genome-wide association scan identified thirteen SNP markers associated with five key traits, suggesting that landraces constitute a source of valuable genes for productivity and food quality traits. Further dissection of their genetic merits in yam breeding using the Genomic Prediction of Cross Performance (GPCP) allowed identifying several landraces with high crossing merit for multiple traits. Thirteen landraces were identified as potential genitors to develop segregating progenies to improve multiple traits simultaneously for desired gains in yam breeding. Results of this study provide valuable insights into the patterns and the merits of local genetic diversity which can be utilized for identifying desirable genes and alleles of interest in yam breeding for Africa.

Chromosome-scale genomes reveal genomic consequences of inbreeding in the South China tiger: A comparative study with the Amur tiger

The South China tiger (Panthera tigris amoyensis, SCT) is the most critically endangered subspecies of tiger due to functional extinction in the wild. Inbreeding depression is observed among the captive population descended from six wild ancestors, resulting in high juvenile mortality and low reproduction. We assembled and characterized the first SCT genome and an improved Amur tiger (P. t. altaica, AT) genome named AmyTig1.0 and PanTig2.0. The two genomes are the most continuous and comprehensive among any tiger genomes yet reported at the chromosomal level. By using the two genomes and resequencing data of 15 SCT and 13 AT individuals, we investigated the genomic signature of inbreeding depression of the SCT. The results indicated that the effective population size of SCT experienced three phases of decline, ~5.0-1.0 thousand years ago, 100 years ago, and since captive breeding in 1963. We found 43 long runs of homozygosity fragments that were shared by all individuals in the SCT population and covered a total length of 20.63% in the SCT genome. We also detected a large proportion of identical-by-descent segments across the genome in the SCT population, especially on ChrB4. Deleterious nonsynonymous single nucleotide polymorphic sites and loss-of-function mutations were found across genomes with extensive potential influences, despite a proportion of these loads having been purged by inbreeding depression. Our research provides an invaluable resource for the formulation of genetic management policies for the South China tiger such as developing genome-based breeding and genetic rescue strategy.

Genome-wide allele and haplotype-sharing patterns suggested one unique Hmong-Mein-related lineage and biological adaptation history in Southwest China

BACKGROUND

Fine-scale genetic structure of ethnolinguistically diverse Chinese populations can fill the gap in the missing diversity and evolutionary landscape of East Asians, particularly for anthropologically informed Chinese minorities. Hmong-Mien (HM) people were one of the most significant indigenous populations in South China and Southeast Asia, which were suggested to be the descendants of the ancient Yangtze rice farmers based on linguistic and archeological evidence. However, their deep population history and biological adaptative features remained to be fully characterized.

OBJECTIVES

To explore the evolutionary and adaptive characteristics of the Miao people, we genotyped genome-wide SNP data in Guizhou HM-speaking populations and merged it with modern and ancient reference populations via a comprehensive population genetic analysis and evolutionary admixture modeling.

RESULTS

The overall genetic admixture landscape of Guizhou Miao showed genetic differentiation between them and other linguistically diverse Guizhou populations. Admixture models further confirmed that Miao people derived their primary ancestry from geographically close Guangxi Gaohuahua people. The estimated identity by descent and effective population size confirmed a plausible population bottleneck, contributing to their unique genetic diversity and population structure patterns. We finally identified several natural selection candidate genes associated with several biological pathways.

CONCLUSIONS

Guizhou Miao possessed a specific genetic structure and harbored a close genetic relationship with geographically close southern Chinese indigenous populations and Guangxi historical people. Miao people derived their major ancestry from geographically close Guangxi Gaohuahua people and experienced a plausible population bottleneck which contributed to the unique pattern of their genetic diversity and structure. Future ancient DNA from Shijiahe and Qujialing will provide new insights into the origin of the Miao people.

Genomic insight into the influence of selection, crossbreeding, and geography on population structure in poultry

BACKGROUND

In poultry, the population structure of local breeds is usually complex mainly due to unrecorded breeding. Local chicken breeds offer an interesting proxy to understand the complexity of population structure in the context of human-mediated development of diverse morphologies and varieties. We studied 37 traditional Dutch chicken breeds to investigate population structure and the corresponding genomic impact using whole-genome sequence data.

RESULTS

Looking at the genetic differences between breeds, the Dutch chicken breeds demonstrated a complex and admixed subdivided structure. The dissection of this complexity highlighted the influence of selection adhering to management purposes, as well as the role of geographic distance within subdivided breed clusters. Identification of signatures of genetic differentiation revealed genomic regions that are associated with diversifying phenotypic selection between breeds, including dwarf size (bantam) and feather color. In addition, with a case study of a recently developed bantam breed developed by crossbreeding, we provide a genomic perspective on the effect of crossbreeding.

CONCLUSIONS

This study demonstrates the complex population structure of local traditional Dutch chicken, and provides insight into the genomic basis and the factors involved in the formation of this complexity.

Distinct traces of mixed ancestry in western commercial pig genomes following gene flow from Chinese indigenous breeds

Studying gene flow between different livestock breeds will benefit the discovery of genes related to production traits and provide insight into human historical breeding. Chinese pigs have played an indispensable role in the breeding of Western commercial pigs. However, the differences in the timing and volume of the contribution of pigs from different Chinese regions to Western pigs are not yet apparent. In this paper, we combine the whole-genome sequencing data of 592 pigs from different studies and illustrate patterns of gene flow from Chinese pigs into Western commercial pigs. We describe introgression patterns from four distinct Chinese indigenous groups into five Western commercial groups. There were considerable differences in the number and length of the putative introgressed segments from Chinese pig groups that contributed to Western commercial pig breeds. The contribution of pigs from different Chinese geographical locations to a given western commercial breed varied more than that from a specific Chinese pig group to different Western commercial breeds, implying admixture within Europe after introgression. Within different Western commercial lines from the same breed, the introgression patterns from a given Chinese pig group seemed highly conserved, suggesting that introgression of Chinese pigs into Western commercial pig breeds mainly occurred at an early stage of breed formation. Finally, based on analyses of introgression signals, allele frequencies, and selection footprints, we identified a ∼2.65 Mb Chinese-derived haplotype under selection in Duroc pigs (CHR14: 95.68-98.33 Mb). Functional and phenotypic studies demonstrate that this PRKG1 haplotype is related to backfat and loin depth in Duroc pigs. Overall, we demonstrate that the introgression history of domestic pigs is complex and that Western commercial pigs contain distinct traces of mixed ancestry, likely derived from various Chinese pig breeds.

FastRecomb: Fast inference of genetic recombination rates in biobank scale data

While rates of recombination events across the genome (genetic maps) are fundamental to genetic research, the majority of current studies only use one standard map. There is evidence suggesting population differences in genetic maps, and thus estimating population-specific maps are of interest. While the recent availability of biobank-scale data offers such opportunities, current methods are not efficient at leveraging very large sample sizes. The most accurate methods are still linkage-disequilibrium (LD)-based methods that are only tractable for a few hundred samples. In this work, we propose a fast and memory-efficient method for estimating genetic maps from population genotyping data. Our method, FastRecomb, leverages the efficient positional Burrows-Wheeler transform (PBWT) data structure for counting IBD segment boundaries as potential recombination events. We used PBWT blocks to avoid redundant counting of pairwise matches. Moreover, we used a panel smoothing technique to reduce the noise from errors and recent mutations. Using simulation, we found that FastRecomb achieves state-of-the-art performance at 10k resolution, in terms of correlation coefficients between the estimated map and the ground truth. This is mainly due to the fact that FastRecomb can effectively take advantage of large panels comprising more than hundreds of thousands of haplotypes. At the same time, other methods lack the efficiency to handle such data. We believe further refinement of FastRecomb would deliver more accurate genetic maps for the genetics community.

Genomic vulnerability to climate change in Quercus acutissima, a dominant tree species in East Asian deciduous forests

Understanding the evolutionary processes that shape the landscape of genetic variation and influence the response of species to future climate change is critical for biodiversity conservation. Here, we sampled 27 populations across the distribution range of a dominant forest tree, Quercus acutissima, in East Asia, and applied genome-wide analyses to track the evolutionary history and predict the fate of populations under future climate. We found two genetic groups (East and West) in Q. acutissima that diverged during Pliocene. We also found a heterogeneous landscape of genomic variation in this species, which may have been shaped by population demography and linked selections. Using genotype-environment association analyses, we identified climate-associated SNPs in a diverse set of genes and functional categories, indicating a model of polygenic adaptation in Q. acutissima. We further estimated three genetic offset metrics to quantify genomic vulnerability of this species to climate change due to the complex interplay between local adaptation and migration. We found that marginal populations are under higher risk of local extinction because of future climate change, and may not be able to track suitable habitats to maintain the gene-environment relationships observed under the current climate. We also detected higher reverse genetic offsets in northern China, indicating that genetic variation currently present in the whole range of Q. acutissima may not adapt to future climate conditions in this area. Overall, this study illustrates how evolutionary processes have shaped the landscape of genomic variation, and provides a comprehensive genome-wide view of climate maladaptation in Q. acutissima.

The power of geohistorical boundaries for modeling the genetic background of human populations: The case of the rural catalan Pyrenees

The genetic variation of the European population at a macro-geographic scale follows genetic gradients which reflect main migration events. However, less is known about factors affecting mating patterns at a micro-geographic scale. In this study we have analyzed 726,718 autosomal single nucleotide variants in 435 individuals from the catalan Pyrenees covering around 200 km of a vast and abrupt region in the north of the Iberian Peninsula, for which we have information about the geographic origin of all grand-parents and parents. At a macro-geographic scale, our analyses recapitulate the genetic gradient observed in Spain. However, we also identified the presence of micro-population substructure among the sampled individuals. Such micro-population substructure does not correlate with geographic barriers such as the expected by the orography of the considered region, but by the bishoprics present in the covered geographic area. These results support that, on top of main human migrations, long ongoing socio-cultural factors have also shaped the genetic diversity observed at rural populations.

Identity by descent mapping of HCV spontaneous clearance in populations of diverse ancestry

Background

Acute infection with hepatitis C virus (HCV) affects millions of individuals worldwide. Host genetics plays a role in spontaneous clearance of the acute infection which occurs in approximately 30% of the individuals. Common variants in GPR158, genes in the interferon lambda (IFNL) cluster, and the MHC region have been associated with HCV clearance in populations of diverse ancestry. Fine mapping of those regions has identified some key variants and amino acids as potential causal variants but the role of rare variants in those regions and in the genome, in general, has not been explored. We aimed to detect haplotypes containing rare variants related to HCV clearance using identity-by-descent (IBD) haplotype sharing between unrelated cases/case pairs and case/controls pairs in 3,608 individuals with European and African ancestry.

Results

We detected 1,711,832 and 5,678,043 and individual pairs of IBD segments in the European and African ancestry individuals, respectively. As expected, individuals of African descent had more, and shorter segments compared to Europeans. We did not detect any significant IBD signals in the known associated gene regions.

Conclusions

IBD is based on sharing of haplotypes and is most powerful in populations with a shared founder or recent common ancestor. For the complex trait of HCV clearance, we used two outbred, global populations that limited our power to detect IBD associations. Overall, in this population-based sample we failed to detect rare variations associated with HCV clearance in individuals of European and African ancestry.

Clustering of Juvenile Canavan disease in an Indian community due to population bottleneck and isolation: genomic signatures of a founder event

Mild/juvenile Canavan disease (M/JCD) is less frequently reported in the literature and little is known about its pathogenetic mechanisms. We report a comprehensive investigation into the pathogenetic mechanism of a novel NM_000049.4(ASPA):c.526G>A variant in two families. The families belong to Telugu Devanga Chettiar community (TDC) from southern India. TDC has a complex history of migration from their historical origin centuries ago with high endogamy. TDC probably has the highest clustering M/JCD recorded historically (around 24 cases). The pathogenic variant was shown to cause non-classical splicing defect resulting in two different transcripts. The splicing aberration, a loss of function mechanism coupled with a milder missense effect can explain the milder phenotype compared to the infantile-onset CD. The high clustering of an extremely rare form of neurodegenerative disorder with reduced fitness, led us to speculate the possibility of a founder event. Genotyping array of TDC and multiple distinct populations of Indian origin for several population genetic parameters was performed. It yielded robust signatures of a founder event in TDC, such as a high fixation index, increased runs of homozygosity and identity-by-descent in the absence of consanguinity; a large haplotype with high linkage disequilibrium among markers comprising the pathogenic variant; a robust population structure; mutation dating, estimating the age of the potential founder of TDC at around 375 years; possibly a high carrier rate in TDC. This study has not only focused its attention on natural history and pathogenetics but also paves way for carrier screening programs in TDC and future therapeutic studies.

Heritability estimates and predictive ability for pig meat quality traits using identity-by-state and identity-by-descent relationships in an F2 population

Genomic relationships can be computed with dense genome-wide genotypes through different methods, either based on identity-by-state (IBS) or identity-by-descent (IBD). The latter has been shown to increase the accuracy of both estimated relationships and predicted breeding values. However, it is not clear whether an IBD approach would achieve greater heritability ( h 2 ) and predictive ability ( r ̂ y , y ̂ ) than its IBS counterpart for data with low-depth pedigrees. Here, we compare both approaches in terms of the estimated of h 2 and r ̂ y , y ̂ , using data on meat quality and carcass traits recorded in experimental crossbred pigs, with a pedigree constrained to only three generations. Three animal models were fitted which differed on the relationship matrix: an IBS model ( G IBS ), an IBD (defined within the known pedigree) model ( G IBD ), and a pedigree model ( A 22 ). In 9 of 20 traits, the range of increase for the estimates of σ u 2 and h 2 was 1.2-2.9 times greater with G IBS and G IBD models than with A 22 . Whereas for all traits, both parameters were similar between genomic models. The r ̂ y , y ̂ of the genomic models was higher compared to A 22 . A scarce increment in r ̂ y , y ̂ was found with G IBS when compared to G IBD , most likely due to the former recovering sizeable relationships among founder F₀ animals.

Out of the Qinghai-Tibet plateau: Genomic biogeography of the alpine monospecific genus Megadenia (Biscutelleae, Brassicaceae)

Numerous high-elevation alpine plants of the Qinghai-Tibet Plateau (QTP) also have disjunct distribution in adjacent low-altitude mountains. The out-of-QTP versus into-the-QTP hypothesis of alpine plants provide strong evidence for the highly disputed assumption of the massive ice sheet developed in the central plateau during the Last Glacial Maximum (LGM). In this study, we sequenced the genomes of most known populations of Megadenia, a monospecific alpine genus of Brassicaceae distributed primarily in the QTP, though rarely found in adjacent low-elevation mountains of north China and Russia (NC-R). All sequenced samples clustered into four geographic genetic groups: one pair was in the QTP and another was in NC-R. The latter pair is nested within the former, and these findings support the out-of-QTP hypothesis. Dating the four genetic groups and niche distribution suggested that Megadenia migrated out of the QTP to adjacent regions during the LGM. The NC-R group showed a decrease in the effective population sizes. In addition, the genes with high genetic divergences in the QTP group were mainly involved in habitat adaptations during low-altitude colonization. These findings reject the hypothesis of development massive ice sheets, and support glacial survival of alpine plants within, as well as further migration out of, the QTP.

The genomic analysis of current-day North African populations reveals the existence of trans-Saharan migrations with different origins and dates

The Sahara Desert has acted as a barrier to human gene-flow between the northern and central parts of Africa since its aridification. Nonetheless, some contacts between both sides of the desert have occurred throughout history, mainly driven by commercial activity. Part of this was the infamous trans-Saharan slave trade, which forcedly brought peoples from south of the Sahara to North Africa from Roman times until the nineteenth century. Although historical records exist, the genetic aspects of these trans-Saharan migrations have not been deeply studied. In the present study, we assess the genetic influence of trans-Saharan migrations in current-day North Africa and characterize its amount, geographical origin, and dates. We confirm the heterogeneous and generally low-frequency presence of genomic segments of sub-Saharan origin in present-day North Africans acquired in recent historical times, and we show evidence of at least two admixture events: one dated around the thirteenth-fourteenth centuries CE between North Africans and a Western-sub-Saharan-like source similar to current-day Senegambian populations, and another one dated around the seventeenth century CE involving Tunisians and an Eastern-sub-Saharan-like source related to current-day south-Sudan and Kenyan populations. Time and location coincide with the peak of trans-Saharan slave-trade activity between Western African empires and North African powers, and are also concordant with the possibility of continuous recent south-to-north gene-flow. These findings confirm the trans-Saharan human genetic contacts, providing new and precise evidence about its possible dates and geographical origins, which are pivotal to understanding the genomic composition of an underrepresented region such as North Africa.

Genome resequencing reveals independent domestication and breeding improvement of naked oat

As an important cereal crop, common oat, has attracted more and more attention due to its healthy nutritional components and bioactive compounds. Here, high-depth resequencing of 115 oat accessions and closely related hexaploid species worldwide was performed. Based on genetic diversity and linkage disequilibrium analysis, it was found that hulled oat (Avena sativa) experienced a more severe bottleneck than naked oat (Avena sativa var. nuda). Combined with the divergence time of ∼51,200 years ago, the previous speculation that naked oat was a variant of hulled oat was rejected. It was found that the common segments that hulled oat introgressed to naked oat cultivars contained 444 genes, mainly enriched in photosynthetic efficiency-related pathways. Selective sweeps during environmental adaptation and breeding improvement were identified in the naked oat genome. Candidate genes associated with smut resistance and the days to maturity phenotype were also identified. Our study provides genomic resources and new insights into naked oat domestication and breeding.

Open-source benchmarking of IBD segment detection methods for biobank-scale cohorts

In the recent biobank era of genetics, the problem of identical-by-descent (IBD) segment detection received renewed interest, as IBD segments in large cohorts offer unprecedented opportunities in the study of population and genealogical history, as well as genetic association of long haplotypes. While a new generation of efficient methods for IBD segment detection becomes available, direct comparison of these methods is difficult: existing benchmarks were often evaluated in different datasets, with some not openly accessible; methods benchmarked were run under suboptimal parameters; and benchmark performance metrics were not defined consistently. Here, we developed a comprehensive and completely open-source evaluation of the power, accuracy, and resource consumption of these IBD segment detection methods using realistic population genetic simulations with various settings. Our results pave the road for fair evaluation of IBD segment detection methods and provide an practical guide for users.

The loss of biodiversity in Madagascar is contemporaneous with major demographic events

Only 400 km off the coast of East Africa, the island of Madagascar is one of the last large land masses to have been colonized by humans. While many questions surround the human occupation of Madagascar, recent studies raise the question of human impact on endemic biodiversity and landscape transformation. Previous genetic and linguistic analyses have shown that the Malagasy population has emerged from an admixture that happened during the last millennium, between Bantu-speaking African populations and Austronesian-speaking Asian populations. By studying the sharing of chromosome segments between individuals (IBD determination), local ancestry information, and simulated genetic data, we inferred that the Malagasy ancestral Asian population was isolated for more than 1,000 years with an effective size of just a few hundred individuals. This isolation ended around 1,000 years before present (BP) by admixture with a small African population. Around the admixture time, there was a rapid demographic expansion due to intrinsic population growth of the newly admixed population, which coincides with extensive changes in Madagascar's landscape and the extinction of all endemic large-bodied vertebrates. Therefore, our approach can provide new insights into past human demography and associated impacts on ecosystems.

Genomic insight into genetic changes and shaping of major inbred rice cultivars in China

The annual planting area of major inbred rice (Oryza sativa) cultivars reach more than half of the total annual planting area of inbred rice cultivars in China. However, how the major inbred rice cultivars changed during decades of genetic improvement and why they can be prevalently cultivated in China remains unclear. Here, we investigated the underlying genetic changes of major inbred cultivars and the contributions of landraces and introduced cultivars during the improvement by resequencing a collection of 439 rice accessions including major inbred cultivars, landraces, and introduced cultivars. The results showed that landraces were the main genetic contribution sources of major inbred Xian (Indica) cultivars, while introduced cultivars were that of major inbred Geng (Japonica) cultivars. Selection scans and haplotype frequency analysis shed light on the reflections of some well-known genes in rice improvement, and breeders had different preferences for the Xian's and Geng's breeding. Six candidate regions associated with agronomic traits were identified by genome-wide association mapping, five of which were under positive selection in rice improvement. Our study provides a comprehensive insight into the development of major inbred rice cultivars and lays the foundation for genomics-based breeding in rice.

Identity-by-descent analysis of a large Tourette's syndrome pedigree from Costa Rica implicates genes involved in neuronal development and signal transduction

Tourette Syndrome (TS) is a heritable, early-onset neuropsychiatric disorder that typically begins in early childhood. Identifying rare genetic variants that make a significant contribution to risk in affected families may provide important insights into the molecular aetiology of this complex and heterogeneous syndrome. Here we present a whole-genome sequencing (WGS) analysis from the 11-generation pedigree (>500 individuals) of a densely affected Costa Rican family which shares ancestry from six founder pairs. By conducting an identity-by-descent (IBD) analysis using WGS data from 19 individuals from the extended pedigree we have identified putative risk haplotypes that were not seen in controls, and can be linked with four of the six founder pairs. Rare coding and non-coding variants present on the haplotypes and only seen in haplotype carriers show an enrichment in pathways such as regulation of locomotion and signal transduction, suggesting common mechanisms by which the haplotype-specific variants may be contributing to TS-risk in this pedigree. In particular we have identified a rare deleterious missense variation in RAPGEF1 on a chromosome 9 haplotype and two ultra-rare deleterious intronic variants in ERBB4 and IKZF2 on the same chromosome 2 haplotype. All three genes play a role in neurodevelopment. This study, using WGS data in a pedigree-based approach, shows the importance of investigating both coding and non-coding variants to identify genes that may contribute to disease risk. Together, the genes and variants identified on the IBD haplotypes represent biologically relevant targets for investigation in other pedigree and population-based TS data.