Exploring Population Structure with Admixture Models and Principal Component Analysis

Distinguishing recrudescence from reinfection in lymphatic filariasis

BACKGROUND

The Global Program to Eliminate Lymphatic Filariasis (GPELF) is the largest public health program based on mass drug administration (MDA). Despite decades of MDA, ongoing transmission in some countries remains a challenge. To optimise interventions, it is critical to differentiate between recrudescence and new infections. Since adult filariae are inaccessible in humans, deriving a method that relies on the offspring microfilariae (mf) is necessary.

METHODS

We developed a genome amplification and kinship analysis-based approach using Brugia malayi samples from gerbils, and applied it to analyse Wuchereria bancrofti mf from humans in Côte d'Ivoire. We examined the pre-treatment genetic diversity in 269 mf collected from 18 participants, and further analysed 1-year post-treatment samples of 74 mf from 4 participants. Hemizygosity of the male X-chromosome allowed for direct inference of haplotypes, facilitating robust maternal parentage inference. To enrich parasite DNA from samples contaminated with host DNA, a whole-exome capture panel was created for W. bancrofti.

FINDINGS

By reconstructing and temporally tracking sibling relationships across pre- and post-treatment samples, we differentiated between new and established maternal families, suggesting reinfection in one participant and recrudescence in three participants. The estimated number of reproductively active adult females ranged between 3 and 11 in the studied participants. Population structure analysis revealed genetically distinct parasites in Côte d'Ivoire compared to samples from other countries. Exome capture identified protein-coding variants with ∼95% genotype concordance rate.

INTERPRETATION

We have generated resources to facilitate the development of molecular genetic tools that can estimate adult worm burdens and monitor parasite populations, thus providing essential information for the successful implementation of GPELF.

FUNDING

This work was financially supported by the Bill and Melinda Gates Foundation (https://www.gatesfoundation.org) under grant OPP1201530 (Co-PIs PUF & Gary J. Weil). B. malayi parasite material was generated with support of the Foundation for Barnes Jewish Hospital (PUF). In addition, the development of computational methods was supported by the National Institutes of Health under grants AI144161 (MM) and AI146353 (MM). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Comparing the performances of SSR and SNP markers for population analysis in Theobroma cacao L., as alternative approach to validate a new ddRADseq protocol for cacao genotyping

Proper cacao (Theobroma cacao L.) plant genotyping is mandatory for the conservation and use of the species genetic resources. A set of 15 international standard SSR markers was assumed as universal cacao genotyping system. Recently, different SNPs and SNP genotyping techniques have been exploited in cacao. However, a consensus on which to use has not been reached yet, driving the search for new approaches. To validate a new ddRADseq protocol for cacao genotyping, we compared the performances for population analysis of a dataset with 7,880 SNPs obtained from ddRADseq and the genotypic data from the aforementioned SSR set, using 158 cacao plants from productive farms and gene bank. Four genetic groups were identified with STRUCTURE and ADMIXTURE softwares using SSR and SNP data, respectively. Similarities of cacao ancestries among these groups allowed the identification of analogous pairs of groups of individuals, referred to as: G1SSR/G1SNP, G2SSR/G2SNP, G3SSR/G3SNP, G4SSR/G4SNP, whether SSRs or SNPs were used. Both marker systems identified Amelonado and Criollo as the most abundant cacao ancestries among all samples. Genetic distance matrices from both data types were significantly similar to each other according to Mantel test (p < 0.0001). PCoA and UPGMA clustering mostly confirmed the identified genetic groups. AMOVA and FST pairwise comparison revealed a moderate to very large genetic differentiation among identified groups from SSR and SNP data. Genetic diversity parameters from SSR (Hobs = 0.616, Hexp = 0.524 and PIC = 0.544) were higher than that from SNP data (0.288, 0.264, 0.230). In both cases, genetic groups carrying the highest Amelonado proportion (G1SSR and G1SNP) had the lowest genetic diversity parameters among the identified groups. The high congruence among population analysis results using both systems validated the ddRADseq protocol employed for cacao SNP genotyping. These results could provide new ways for developing a universal SNP-based genotyping system very much needed for cacao genetic studies.

Genetic diversity, population structure and kinship relationships highlight the environmental influence on Uganda's indigenous goat populations

Knowledge about genetic diversity and population structure among goat populations is essential for understanding environmental adaptation and fostering efficient utilization, development, and conservation of goat breeds. Uganda's indigenous goats exist in three phenotypic groups: Mubende, Kigezi, and Small East African. However, a limited understanding of their genetic attributes and population structure hinders the development and sustainable utilization of the goats. Using the Goat Illumina 60k chip International Goat Genome Consortium V2, the whole-genome data for 1,021 indigenous goats sourced from 10 agroecological zones in Uganda were analyzed for genetic diversity and population structure. A total of 49,337 (82.6%) single-nucleotide polymorphism markers were aligned to the ARS-1 goat genome and used to assess the genetic diversity, population structure, and kinship relationships of Uganda's indigenous goats. Moderate genetic diversity was observed. The observed and expected heterozygosities were 0.378 and 0.383, the average genetic distance was 0.390, and the average minor allele frequency was 0.30. The average inbreeding coefficient (Fis) was 0.014, and the average fixation index (Fst) was 0.016. Principal component analysis, admixture analysis, and discriminant analysis of principal components grouped the 1,021 goat genotypes into three genetically distinct populations that did not conform to the known phenotypic populations but varied across environmental conditions. Population 1, comprising Mubende (90%) and Kigezi (8.1%) goats, is located in southwest and central Uganda, a warm and humid environment. Population 2, which is 59% Mubende and 49% Small East African goats, is located along the Nile Delta in northwestern Uganda and around the Albertine region, a hot and humid savannah grassland. Population 3, comprising 78.4% Small East African and 21.1% Mubende goats, is found in northeastern to eastern Uganda, a hot and dry Commiphora woodlands. Genetic diversity and population structure information from this study will be a basis for future development, conservation, and sustainable utilization of Uganda's goat genetic resources.

Genome-wide and pan-genomic analysis reveals rich variants of NBS-LRR genes in a newly developed wild rice line from Oryza alta Swallen

Introduction

Oryza alta Swallen is an allotetraploid perennial wild rice and contains CCDD genome, which may harbor favorable genes for the enrichment of genetic resource.

Methods

A new wild rice line, Huaye 5, was developed from Oryza alta Swallen in our lab. Whole genome re-sequencing and pan-genomic analysis were employed to analyze its genomic variations and novel genes.

Results and Discussion

More than ten million genomic variations were detected when compared with Asian cultivar. Among the variational genes, 724, 197 and 710 genes coded protein kinase, synthetase and transcription factor, respectively. A total of 353, 131 and 135 variational genes were associated with morphological trait, physiological trait, resistance or tolerance, respectively. A total of 62 were NBS-LRR genes were detected, in which 11 NBS-LRR genes expressed in sheath and mature stem, and 26 expressed in young and mature roots expressed. The pan-genome sequences of wild rice species with CCDD genome were constructed by integrating 8 Oryza alta (OA), 2 Oryza grandiglumis (OG) and 18 Oryza latifolia (OL) accessions. A total of 28 non-reference NBS-LRR genes were revealed, and 7 of which were mainly expressed in mature roots. This research demonstrated rich DNA variation in the Oryza alta Swallen that may provide a new germplasm for rice resistance breeding.

Multi-locus genome-wide association studies reveal the dynamic genetic architecture of flowering time in chrysanthemum

KEY MESSAGE

The dynamic genetic architecture of flowering time in chrysanthemum was elucidated by GWAS. Thirty-six known genes and 14 candidate genes were identified around the stable QTNs and QEIs, among which ERF-1 was highlighted. Flowering time (FT) adaptation is one of the major breeding goals in chrysanthemum, a multipurpose ornamental plant. In order to reveal the dynamic genetic architecture of FT in chrysanthemum, phenotype investigation of ten FT-related traits was conducted on 169 entries in 2 environments. The broad-sense heritability of five non-conditional FT traits, i.e., budding (FBD), visible coloring (VC), early opening (EO), full-bloom (OF) and decay period (DP), ranged from 56.93 to 84.26%, which were higher than that of the five derived conditional FT traits (38.51-75.13%). The phenotypic variation coefficients of OF_EO and DP_OF were relatively large ranging from 30.59 to 36.17%. Based on 375,865 SNPs, the compressed variance component mixed linear model 3VmrMLM was applied for a multi-locus genome-wide association study (GWAS). As a result, 313 quantitative trait nucleotides (QTNs) were identified for the non-conditional FT traits in single-environment analysis, while 119 QTNs and 67 QTN-by-environment interactions (QEIs) were identified in multi-environment analysis. As for the conditional traits, 343 QTNs were detected in single-environment analysis, and 119 QTNs and 83 QEIs were identified in multi- environment analysis. Among the genes around stable QTNs and QEIs, 36 were orthologs of known FT genes in Arabidopsis and other plants; 14 candidates were mined by combining the transcriptomics data and functional annotation, including ERF-1, ACA10, and FOP1. Furthermore, the haplotype analysis of ERF-1 revealed six elite accessions with extreme FBD. Our findings contribute to the understanding of dynamic genetic architecture of FT and provide valuable resources for future chrysanthemum molecular breeding programs.

Using ddRADseq to assess the genetic diversity of in-farm and gene bank cacao resources in the Baracoa region, eastern Cuba, for use and conservation purposes

The Baracoa region, eastern Cuba, hosts around 80 % of the country cacao (Theobroma cacao L.) plantations. Cacao plants in farms are diverse in origin and propagation, with grafted and hybrid plants being the more common ones. Less frequent are plants from cuttings, TSH progeny, and traditional Cuban cacao. A national cacao gene bank is also present in Baracoa, with 282 accessions either prospected in Cuba or introduced from other countries. A breeding program associated with the gene bank started in the 1990s based on agro-morphological descriptors. The genetic diversity of cacao resources in Baracoa has been poorly described, except for traditional Cuban cacao, affecting the proper development of the breeding program and the cacao planting policies in the region. To assess the population structure and genetic diversity of cacao resources in Baracoa region, we genotyped plants from both cacao gene bank (CG) and cacao farms (CF) applying a new ddRADseq protocol for cacao. After data processing, two SNPs datasets containing 11,425 and 6,481 high-quality SNPs were generated with 238 CG and 135 CF plants, respectively. SNPs were unevenly distributed along the 10 cacao chromosomes and laid mainly in noncoding regions of the genome. Population structure analysis with these SNP datasets identified seven and four genetic groups in CG and CF samples, respectively. Clustering using UPGMA and principal component analysis mostly agree with population structure results. Amelonado was the predominant cacao ancestry, accounting for 49.22 % (CG) and 57.73 % (CF) of the total. Criollo, Contamana, Iquitos, and Nanay ancestries were detected in both CG and CF samples, while Nacional and Marañon backgrounds were only identified in CG. Genetic differentiation among CG (FST ranging from 0.071 to 0.407) was higher than among CF genetic groups (FST: 0.093-0.282). Genetic diversity parameters showed similar values for CG and CF samples. The CG and CF genetic groups with the lowest genetic diversity parameters had the highest proportion of Amelonado ancestry. These results should contribute to reinforcing the ongoing breeding program and updating the planting policies on cacao farms, with an impact on the social and economic life of the region.

Genome-wide diversity, population structure and signatures of inbreeding in the African buffalo in Mozambique

The African buffalo, Syncerus caffer, is a key species in African ecosystems. Like other large herbivores, it plays a fundamental role in its habitat acting as an ecosystem engineer. Over the last few centuries, African buffalo populations have declined because of range contraction and demographic decline caused by direct or indirect human activities. In Mozambique, historically home to large buffalo herds, the combined effect of colonialism and subsequent civil wars has created a critical situation that urgently needs to be addressed. In this study, we focused on the analysis of genetic diversity of Syncerus caffer caffer populations from six areas of Mozambique. Using genome-wide SNPs obtained from ddRAD sequencing, we examined the population structure across the country, estimated gene flow between areas under conservation management, including national reserves, and assessed the inbreeding coefficients. Our results indicate that all studied populations of Syncerus caffer caffer are genetically depauperate, with a high level of inbreeding. Moreover, buffaloes in Mozambique present a significant population differentiation between southern and central areas. We found an unexpected genotype in the Gorongosa National Park, where buffaloes experienced a dramatic population size reduction, that shares a common ancestry with southern populations of Catuane and Namaacha. This could suggest the past occurrence of a connection between southern and central Mozambique and that the observed population structuring could reflect recent events of anthropogenic origin. All the populations analysed showed high levels of homozygosity, likely due to extensive inbreeding over the last few decades, which could have increased the frequency of recessive deleterious alleles. Improving the resilience of Syncerus caffer caffer in Mozambique is essential for preserving the ecosystem integrity. The most viable approach appears to be facilitating translocations and re-establishing connectivity between isolated herds. However, our results also highlight the importance of assessing intraspecific genetic diversity when considering interventions aimed at enhancing population viability such as selecting suitable source populations.

Patterns of cytonuclear discordance and divergence between subspecies of the scarlet macaw (Ara macao) in Central America

The scarlet macaw, Ara macao, is a neotropical parrot that contains two described subspecies with broadly discrete geographical distributions. One subspecies, A. m. macao, is found from South America north into southwestern Costa Rica, while the second subspecies, A. m. cyanoptera, is found from eastern Costa Rica north into central Mexico. Our previous research using mitochondrial data to examine phylogeographical divergence across the collective range of these two subspecies concluded that they represent distinct evolutionary entities, with minimal contemporary hybridization between them. Here we further examine phylogenetic relationships and patterns of genetic variation between these two subspecies using a dataset of genetic markers derived from their nuclear genomes. Our analyses show clear nuclear divergence between A. m. macao and A. m. cyanoptera in Central America. Collectively however, samples from this region appear genetically more similar to one another than they do to the examined South American (Brazilian) A. m. macao sample. This observation contradicts our previous assessments based on mitochondrial DNA analyses that A. m. macao in Central and South America represent a single phylogeographical group that is evolutionarily distinct from Central American A. m. cyanoptera. Nonetheless, in agreement with our previous findings, ongoing genetic exchange between the two subspecies appears limited. Rather, our analyses indicate that incomplete lineage sorting is the best supported explanation for cytonuclear discordance within these parrots. High-altitude regions in Central America may act as a reproductive barrier, limiting contemporary hybridization between A. m. macao and A. m. cyanoptera. The phylogeographic complexities of scarlet macaw taxa in this region highlight the need for additional evolutionary examinations of these populations.

Genome-wide analysis highlights genetic admixture in exotic germplasm resources of Eucalyptus and unexpected ancestral genomic composition of interspecific hybrids

Eucalyptus is an economically important genus comprising more than 890 species in different subgenera and sections. Approximately twenty species of subgenus Symphyomyrtus account for 95% of the world's planted eucalypts. Discrimination of closely related eucalypt taxa is challenging, consistent with their recent phylogenetic divergence and occasional hybridization in nature. Admixture, misclassification or mislabeling of Eucalyptus germplasm resources maintained as exotics have been suggested, although no reports are available. Moreover, hybrids with increased productivity and traits complementarity are planted worldwide, but little is known about their actual genomic ancestry. In this study we examined a set of 440 trees of 16 different Eucalyptus species and 44 interspecific hybrids of multi-species origin conserved in germplasm banks in Brazil. We used genome-wide SNP data to evaluate the agreement between the alleged phylogenetic classification of species and provenances as registered in their historical records, and their observed genetic clustering derived from SNP data. Genetic structure analyses correctly assigned each of the 16 species to a different cluster although the PCA positioning of E. longirostrata was inconsistent with its current taxonomy. Admixture was present for closely related species' materials derived from local germplasm banks, indicating unintended hybridization following germplasm introduction. Provenances could be discriminated for some species, indicating that SNP-based discrimination was directly proportional to geographical distance, consistent with an isolation-by-distance model. SNP-based genomic ancestry analysis showed that the majority of the hybrids displayed realized genomic composition deviating from the expected ones based on their pedigree records, consistent with admixture in their parents and pervasive genome-wide directional selection toward the fast-growing E. grandis genome. SNP data in support of tree breeding provide precise germplasm identity verification, and allow breeders to objectively recognize the actual ancestral origin of superior hybrids to more realistically guide the program toward the development of the desired genetic combinations.

Genome-wide study of linkage disequilibrium, population structure, and inbreeding in Iranian indigenous sheep breeds

Knowledge of linkage disequilibrium (LD), genetic structure and genetic diversity are some key parameters to study the breeding history of indigenous small ruminants. In this study, the OvineSNP50 Bead Chip array was used to estimate and compare LD, genetic diversity, effective population size (Ne) and genomic inbreeding in 186 individuals, from three Iranian indigenous sheep breeds consisting of Baluchi (n = 96), Lori-Bakhtiari (n = 47) and Zel (n = 47). The results of principal component analysis (PCA) revealed that all animals were allocated to the groups that they sampled and the admixture analysis revealed that the structure within the populations is best explained when separated into three groups (K  =  3). The average r2 values estimated between adjacent single nucleotide polymorphisms (SNPs) at distances up to 10Kb, were 0.388±0.324, 0.353±0.311, and 0.333±0.309 for Baluchi, Lori-Bakhtiari and Zel, respectively. Estimation of genetic diversity and effective population size (Ne) showed that the Zel breed had the highest heterozygosity and Ne, whereas the lowest value was found in Baluchi breed. Estimation of genomic inbreeding using FROH (based on the long stretches of consecutive homozygous genotypes) showed the highest inbreeding coefficient in Baluchi and the lowest in Zel breed that could be due to higher pressure of artificial selection on Baluchi breed. The results of genomic inbreeding and Ne showed an increase in sharing haplotypes in Baluchi, leading to the enlargement of LD and the consequences of linkage disequilibrium and haplotype blocks confirmed this point. Also, the persistence of the LD phase between Zel and Lori-Bakhtiari was highest indicating that these two breeds would be combined in a multi-breed training population in genomic selection studies.

Population genomics of American mink using genotype data

Understanding the genetic structure of the target population is critically important to develop an efficient genomic selection program in domestic animals. In this study, 2,973 American mink of six color types from two farms (Canadian Centre for Fur Animal Research (CCFAR), Truro, NS and Millbank Fur Farm (MFF), Rockwood, ON) were genotyped with the Affymetrix Mink 70K panel to compute their linkage disequilibrium (LD) patterns, effective population size (Ne), genetic diversity, genetic distances, and population differentiation and structure. The LD pattern represented by average r ², decreased to <0.2 when the inter-marker interval reached larger than 350 kb and 650 kb for CCFAR and MFF, respectively, and suggested at least 7,700 and 4,200 single nucleotide polymorphisms (SNPs) be used to obtain adequate accuracy for genomic selection programs in CCFAR and MFF respectively. The Ne for five generations ago was estimated to be 76 and 91 respectively. Our results from genetic distance and diversity analyses showed that American mink of the various color types had a close genetic relationship and low genetic diversity, with most of the genetic variation occurring within rather than between color types. Three ancestral genetic groups was considered the most appropriate number to delineate the genetic structure of these populations. Black (in both CCFAR and MFF) and pastel color types had their own ancestral clusters, while demi, mahogany, and stardust color types were admixed with the three ancestral genetic groups. This study provided essential information to utilize the first medium-density SNP panel for American mink in their genomic studies.

Large Region of Homozygous (ROH) Identified in Indian Patients with Autosomal Recessive Limb-Girdle Muscular Dystrophy with p.Thr182Pro Variant in SGCB Gene

The sarcoglycanopathies are autosomal recessive limb-girdle muscular dystrophies (LGMDs) caused by the mutations in genes encoding the α, β, γ, and δ proteins which stabilizes the sarcolemma of muscle cells. The clinical phenotype is characterized by progressive proximal muscle weakness with childhood onset. Muscle biopsy findings are diagnostic in confirming dystrophic changes and deficiency of one or more sarcoglycan proteins. In this study, we summarized 1,046 LGMD patients for which a precise diagnosis was identified using targeted sequencing. The most frequent phenotypes identified in the patients are LGMDR1 (19.7%), LGMDR4 (19.0%), LGMDR2 (17.5%), and MMD1 (14.5%). Among the reported genes, each of CAPN3, SGCB, and DYSF variants was reported in more than 10% of our study cohort. The most common variant SGCB p.Thr182Pro was identified in 146 (12.5%) of the LGMD patients, and in 97.9% of these patients, the variant was found to be homozygous. To understand the genetic structure of the patients carrying SGCB p.Thr182Pro, we genotyped 68 LGMD patients using a whole genome microarray. Analysis of the array data identified a large ~1 Mb region of homozygosity (ROH) (chr4:51817441-528499552) suggestive of a shared genomic region overlapping the recurrent missense variant and shared across all 68 patients. Haplotype analysis identified 133 marker haplotypes that were present in ~85.3% of the probands as a double allele and absent in all random controls. We also identified 5 markers (rs1910739, rs6852236, rs13122418, rs13353646, and rs6554360) which were present in a significantly higher proportion in the patients compared to random control set ( $n=128$ ) and the population database. Of note, admixture analysis was suggestive of greater proportion of West Eurasian/European ancestry as compared to random controls. Haplotype analysis and frequency in the population database indicate a probable event of founder effect. Further systematic study is needed to identify the communities and regions where the SGCB p.Thr182Pro variant is observed in higher proportions. After identifying these communities and//or region, a screening program is needed to identify carriers and provide them counselling.

Multi-locus genome-wide association study of chickpea reference set identifies genetic determinants of Pratylenchus thornei resistance

Pratylenchus thornei is an economically important species of root-lesion nematode adversely affecting chickpea (Cicer arietinum) yields globally. Integration of resistant crops in farming systems is recognised as the most effective and sustainable management strategy for plant-parasitic nematodes. However, breeding for P. thornei resistance in chickpea is limited by the lack of genetic diversity. We deployed a genome-wide association approach to identify genomic regions and candidate genes associated with P. thornei resistance in 285 genetically diverse chickpea accessions. Chickpea accessions were phenotyped for P. thornei resistance in replicated glasshouse experiments performed for two years (2018 and 2020). Whole genome sequencing data comprising 492,849 SNPs were used to implement six multi-locus GWAS models. Fourteen chickpea genotypes were found to be resistant to P. thornei. Of the six multi-locus GWAS methods deployed, FASTmrMLM was found to be the best performing model. In all, 24 significant quantitative trait nucleotides (QTNs) were identified, of which 13 QTNs were associated with lower nematode population density and 11 QTNs with higher nematode population density. These QTNs were distributed across all of the chickpea chromosomes, except chromosome 8. We identified, receptor-linked kinases (RLKs) on chromosomes 1, 4 and 6, GDSL-like Lipase/Acylhydrolase on chromosome 3, Aspartic proteinase-like and Thaumatin-like protein on chromosome 4, AT-hook DNA-binding and HSPRO2 on chromosome 6 as candidate genes for P. thornei resistance in the chickpea reference set. New sources of P. thornei resistant genotypes were identified that can be harnessed into breeding programs and putative candidate P. thornei resistant genes were identified that can be explored further to develop molecular markers and accelerate the incorporation of improved P. thornei resistance into elite chickpea cultivars.

Dual Domestication, Diversity, and Differential Introgression in Old World Cotton Diploids

Domestication in the cotton genus is remarkable in that it has occurred independently four different times at two different ploidy levels. Relatively little is known about genome evolution and domestication in the cultivated diploid species Gossypium herbaceum and Gossypium arboreum, due to the absence of wild representatives for the latter species, their ancient domestication, and their joint history of human-mediated dispersal and interspecific gene flow. Using in-depth resequencing of a broad sampling from both species, we provide support for their independent domestication, as opposed to a progenitor-derivative relationship, showing that diversity (mean π = 6 × 10-3) within species is similar, and that divergence between species is modest (FST = 0.413). Individual accessions were homozygous for ancestral single-nucleotide polymorphisms at over half of variable sites, while fixed, derived sites were at modest frequencies. Notably, two chromosomes with a paucity of fixed, derived sites (i.e., chromosomes 7 and 10) were also strongly implicated as having experienced high levels of introgression. Collectively, these data demonstrate variable permeability to introgression among chromosomes, which we propose is due to divergent selection under domestication and/or the phenomenon of F2 breakdown in interspecific crosses. Our analyses provide insight into the evolutionary forces that shape diversity and divergence in the diploid cultivated species and establish a foundation for understanding the contribution of introgression and/or strong parallel selection to the extensive morphological similarities shared between species.

Demographic and Selection Histories of Populations Across the Sahel/Savannah Belt

The Sahel/Savannah belt harbors diverse populations with different demographic histories and different subsistence patterns. However, populations from this large African region are notably under-represented in genomic research. To investigate the population structure and adaptation history of populations from the Sahel/Savannah space, we generated dense genome-wide genotype data of 327 individuals-comprising 14 ethnolinguistic groups, including 10 previously unsampled populations. Our results highlight fine-scale population structure and complex patterns of admixture, particularly in Fulani groups and Arabic-speaking populations. Among all studied Sahelian populations, only the Rashaayda Arabic-speaking population from eastern Sudan shows a lack of gene flow from African groups, which is consistent with the short history of this population in the African continent. They are recent migrants from Saudi Arabia with evidence of strong genetic isolation during the last few generations and a strong demographic bottleneck. This population also presents a strong selection signal in a genomic region around the CNR1 gene associated with substance dependence and chronic stress. In Western Sahelian populations, signatures of selection were detected in several other genetic regions, including pathways associated with lactase persistence, immune response, and malaria resistance. Taken together, these findings refine our current knowledge of genetic diversity, population structure, migration, admixture and adaptation of human populations in the Sahel/Savannah belt and contribute to our understanding of human history and health.

The recombination landscapes of spiny lizards (genus Sceloporus)

Despite playing a critical role in evolutionary processes and outcomes, relatively little is known about rates of recombination in the vast majority of species, including squamate reptiles—the second largest order of extant vertebrates, many species of which serve as important model organisms in evolutionary and ecological studies. This paucity of data has resulted in limited resolution on questions related to the causes and consequences of rate variation between species and populations, the determinants of within-genome rate variation, as well as the general tempo of recombination rate evolution on this branch of the tree of life. In order to address these questions, it is thus necessary to begin broadening our phylogenetic sampling. We here provide the first fine-scale recombination maps for two species of spiny lizards, Sceloporus jarrovii and Sceloporus megalepidurus, which diverged at least 12 Mya. As might be expected from similarities in karyotype, population-scaled recombination landscapes are largely conserved on the broad-scale. At the same time, considerable variation exists at the fine-scale, highlighting the importance of incorporating species-specific recombination maps in future population genomic studies.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of cultivated-type tea plant in Guizhou Plateau

BACKGROUND

Tea plants originated in southwestern China. Guizhou Plateau is an original center of tea plants, and is rich in germplasm resources. However, the genetic diversity, population structure and distribution characteristics of cultivated-type tea plants in the region are unknown. In this study, we explored the genetic diversity and geographical distribution of cultivated-type tea accessions in Guizhou Plateau.

RESULTS

We used 112,072 high-quality genotyping-by-sequencing to analyze the genetic diversity, principal components, phylogeny, population structure, and linkage disequilibrium, and develop a core collection of 253 cultivated-type tea plant accessions from Guizhou Plateau. The results showed Genetic diversity of the cultivated-type tea accessions of the Pearl River Basin was significantly higher than that of the cultivated-type tea accessions of the Yangtze River Basin. Three inferred pure groups (CG-1, CG-2 and CG-3) and one inferred admixture group (CG-4), were identified by a population structure analysis, and verified by principal component and phylogenetic analyses. The highest genetic distance and differentiation coefficients were determined for CG-2 vs CG-3. The lower genetic distance and differentiation coefficients were determined for CG-4 vs CG-2 and CG-4 vs CG-3, respectively. We developed a core set and a primary set. The primary and core sets contained 77.0 and 33.6% of all individuals in the initial set, respectively. The primary set may serve as the primary population in genome-wide association studies, while the core collection may serve as the core population in multiple treatment setting studies.

CONCLUSIONS

The present study demonstrated the genetic diversity and geographical distribution characteristics of cultivated-type tea plants in Guizhou Plateau. Significant differences in genetic diversity and evolutionary direction were detected between the ancient landraces of the Pearl River Basin and the those of the Yangtze River Basin. Major rivers and ancient hubs were largely responsible for the genetic exchange between the Pearl River Basin and the Yangtze River Basin ancient landraces as well as the formation of the ancient hubs evolutionary group. Genetic diversity, population structure and core collection elucidated by this study will facilitate further genetic studies, germplasm protection, and breeding of tea plants.

First Glimpse into the Genomic Characterization of People from the Imperial Roman Community of Casal Bertone (Rome, First–Third Centuries AD)

This paper aims to provide a first glimpse into the genomic characterization of individuals buried in Casal Bertone (Rome, first–third centuries AD) to gain preliminary insight into the genetic makeup of people who lived near a tannery workshop, fullonica. Therefore, we explored the genetic characteristics of individuals who were putatively recruited as fuller workers outside the Roman population. Moreover, we identified the microbial communities associated with humans to detect microbes associated with the unhealthy environment supposed for such a workshop. We examined five individuals from Casal Bertone for ancient DNA analysis through whole-genome sequencing via a shotgun approach. We conducted multiple investigations to unveil the genetic components featured in the samples studied and their associated microbial communities. We generated reliable whole-genome data for three samples surviving the quality controls. The individuals were descendants of people from North African and the Near East, two of the main foci for tannery and dyeing activity in the past. Our evaluation of the microbes associated with the skeletal samples showed microbes growing in soils with waste products used in the tannery process, indicating that people lived, died, and were buried around places where they worked. In that perspective, the results represent the first genomic characterization of fullers from the past. This analysis broadens our knowledge about the presence of multiple ancestries in Imperial Rome, marking a starting point for future data integration as part of interdisciplinary research on human mobility and the bio-cultural characteristics of people employed in dedicated workshops.

Multiple evolutionary origins of glyphosate resistance in Lolium multiflorum

The multitude of herbicide resistance patterns that have evolved in different weed species is a remarkable example of the rapid adaptation to anthropogenic‐driven disturbance. Recently, resistance to glyphosate was identified in multiple populations of Lolium multiflorum in Oregon. We used phenotypic approaches, as well as population genomic and gene expression analyses, to determine whether known mechanisms were responsible for glyphosate resistance and whether resistance phenotypes evolved independently in different populations, and to identify potential loci contributing to resistance. We found no evidence of genetic alterations or expression changes at known target and non‐target‐site resistance mechanisms of glyphosate. Population genomic analyses indicated that resistant populations tended to have largely distinct ancestry from one another, suggesting that glyphosate resistance did not spread among populations by gene flow. Rather, resistance appears to have evolved independently on different genetic backgrounds. We also detected potential loci associated with the resistance phenotype, some of which encode proteins with potential effects on herbicide metabolism. Our results suggest that Oregon populations of L. multiflorum evolved resistance to glyphosate due to a novel mechanism. Future studies that characterize the gene or genes involved in resistance will be necessary to confirm this conclusion.

Detection of breeding signatures in wheat using a linkage disequilibrium-corrected mapping approach

Marker assisted breeding, facilitated by reference genome assemblies, can help to produce cultivars adapted to changing environmental conditions. However, anomalous linkage disequilibrium (LD), where single markers show high LD with markers on other chromosomes but low LD with adjacent markers, is a serious impediment for genetic studies. We used a LD-correction approach to overcome these drawbacks, correcting the physical position of markers derived from 15 and 135 K arrays in a diversity panel of bread wheat representing 50 years of breeding history. We detected putative mismapping of 11.7% markers and improved the physical alignment of 5.4% markers. Population analysis indicated reduced genetic diversity over time as a result of breeding efforts. By analysis of outlier loci and allele frequency change over time we traced back the 2NS/2AS translocation of Aegilops ventricosa to one cultivar, "Cardos" (registered in 1998) which was the first among the panel to contain this translocation. A "selective sweep" for this important translocation region on chromosome 2AS was found, putatively linked to plant response to biotic stress factors. Our approach helps in overcoming the drawbacks of incorrectly anchored markers on the wheat reference assembly and facilitates detection of selective sweeps for important agronomic traits.