Genotyping Array Research Articles

A whole-genome reference panel of 14,393 individuals for East Asian populations accelerates discovery of rare functional variants.

Underrepresentation of non-European (EUR) populations hinders growth of global precision medicine. Resources such as imputation reference panels that match the study population are necessary to find low-frequency variants with substantial effects. We created a reference panel consisting of 14,393 whole-genome sequences including more than 11,000 Asian individuals. Genome-wide association studies were conducted using the reference panel and a population-specific genotype array of 72,298 subjects for eight phenotypes. This panel yields improved imputation accuracy of rare and low-frequency variants within East Asian populations compared with the largest reference panel. Thirty-nine previously unidentified associations were found, and more than half of the variants were East Asian specific. We discovered genes with rare protein-altering variants, including LTBP1 for height and GPR75 for body mass index, as well as putative regulatory mechanisms for rare noncoding variants with cell type-specific effects. We suggest that this dataset will add to the potential value of Asian precision medicine.

Genetic and Metabolic Characteristics of Lean Nonalcoholic Fatty Liver Disease in a Korean Health Examinee Cohort.

The pathophysiology of lean nonalcoholic fatty liver disease (NAFLD) is unclear but has been shown to be associated with more diverse pathogenic mechanisms than that of obese NAFLD. We investigated the characteristics of genetic or metabolic lean NAFLD in a health checkup cohort. This retrospective cross-sectional study analyzed single nucleotide polymorphism data for 6,939 health examinees. Lean individuals were categorized according to a body mass index cutoff of 23 kg/m2. Single nucleotide polymorphisms were analyzed using genotyping arrays. The prevalence of lean NAFLD was 21.6% among all participants with NAFLD, and the proportion of lean NAFLD was 18.5% among lean participants. The prevalence of metabolic syndrome and diabetes among lean patients with NAFLD was 12.4% and 10.4%, respectively. Lean NAFLD appeared to be metabolic-associated in approximately 20.1% of patients. The homozygous minor allele (GG) of PNPLA3 (rs738409) and heterozygous minor alleles (CT, TT) of TM6SF2 (rs58542926) were associated with lean NAFLD. However, the prevalence of fatty liver was not associated with the genetic variants MBOAT7 (rs641738), HSD17B13 (rs72613567), MARC1 (rs2642438), or AGXT2 (rs2291702) in lean individuals. Lean NAFLD appeared to be associated with PNPLA3 or TM6SF2 genetic variation in approximately 32.1% of cases. Multivariate risk factor analysis showed that metabolic risk factors, genetic risk variants, and waist circumference were independent risk factors for lean NAFLD. In a considerable number of patients, lean NAFLD did not appear to be associated with known genetic or metabolic risk factors. Further studies are required to investigate additional risk factors and gain a more comprehensive understanding of lean NAFLD.

Single Nucleotide Polymorphisms of MTHFR (rs1051266) and SLC19A1 (rs1801133) Associated to Genomic Ancestry in Cuban Healthy Population

Background: Several single nucleotide polymorphisms on methotrexate pathway have been implicated with hyperhomocysteinemia, susceptibility to autoimmune diseases and the therapy effectiveness of methotrexate Objective: The present study estimates the ethnogeographic prevalence of rs1801133 (c.665C>T) in methylenetetrahydrofolate reductase and rs1051266 (c.80A>G) in solute carrier family Methods: Genomic data was collected from a dense genome-wide genotyping array analysis of a large sample of individuals from all provinces of Cuba, with a final sample of 946 individuals for rs1801133 and 948 individuals for rs1051266. Results: For rs1801133, T allele and TT genotype were more prevalent in Mayabeque, the province with the highest European (p<0.0001) and the lowest African ancestry proportion (p<0.0001). Whereas, T allele and TT genotype frequency were low in Guantánamo (23.7% and 1.8%), the province with the highest African ancestry proportion (p<0.0001) and the lowest European ancestry proportion (p<0.0001). For rs1051266, the higher frequency of G allele was observed in Villa Clara, Las Tunas, Holguín and Granma and this group was associated with AG and GG genotypes (p=0.0045). This seems to be related to high Native American ancestry proportion in Las Tunas (p<0.0001), Holguín (p<0.0001) and Granma (p<0.0001); with the low African ancestry proportion in Villa Clara (p<0.0001) and with a Native American ancestry-enriched pattern observed for these provinces (p=0.0005). Conclusion: These results provide evidence that ancestry contribution impacts in the ethnogeographic prevalence of rs1801133 (c.665C>T) and rs1051266 (c.80A>G) polymorphisms in healthy Cuban individuals.

Somatic cell structural variant mutagenesis and neurologic disease

Detection of organismal mosaic states for variant alleles faces technical and analytical challenges, as does the association of such variant alleles with susceptibility to neurologic disease. In this issue of Cell Genomics, Maury etal.1 reanalyze genotyping arrays of a schizophrenia cohort providing evidence for the contribution of somatic structural variant mutagenesis and rare variant alleles.

Mapping of Major Resistance Gene Qse.xn-2BL for Sharp Eyespot in the Wheat-Psathyrostachys huashanica Introgression Line H83.

Sharp eyespot, a soil-borne disease of wheat (Triticum aestivum L.), is one of the most devastating diseases and severely affects grain production. The most efficient and economical method of controlling the disease is the utilization of genetic resistance. In this study, the wheat-Psathyrostachys huashanica introgression line H83 processed the enhanced resistance to Rhizoctonia cerealis isolate R0301 than its wheat parent 7182. A resistance locus in the 600 to 800 Mb interval of chromosome 2BL was screened using 244 segregation population F2 plants of H83×Huixianhong with bulked segregant analysis and wheat axiom 660K genotyping array. Furthermore, by using 12 kompetitive allele-specific PCR markers, a major resistance gene, designated as Qse.xn-2BL, was identified in a secondary segregating population with 138 F3:4 lines and initially mapped to a 765.6 to 775.5 Mb interval on chromosome 2BL. Molecular cytology analysis revealed that H83 probably has an alien introgression at the distal of chromosome 2BL, where it overlapped with the mapping target gene. Above all, H83 showed great potential to improve wheat resistance to sharp eyespot and can be expected to improve resistance in wheat breeding.

Genome-wide Association Study Identifies Novel Risk Loci for Apical Periodontitis

IntroductionApical periodontitis (AP) is a common consequence of root canal infection leading to periapical bone resorption. Microbial and host genetic factors and their interactions have been shown to play a role in AP development and progression. Variations in a few genes have been reported in association with AP; however, the lack of genome-wide studies has hindered progress in understanding the molecular mechanisms involved. Here, we report the first genome-wide association study of AP in a large and well-characterized population. MethodsMale and female adults (n = 932) presenting with deep caries and AP (cases), or deep caries without AP (controls) were included. Genotyping was performed using the Illumina Expanded Multi-Ethnic Genotyping Array (MEGA). Single-variant association testing was performed adjusting for sex and 5 principal components. Subphenotype association testing, analyses of genetically regulated gene expression, polygenic risk score, and phenome-wide association (PheWAS) analyses were also conducted. ResultsEight loci reached near genome-wide significant association with AP (P < 5 × 10−6); gene-focused analyses replicated 3 previously reported associations (P < 8.9 × 10−5). Sex-specific and subphenotype-specific analyses revealed additional significant associations with variants genome-wide. Functionally oriented gene-based analyses revealed 8 genes significantly associated with AP (P < 5 × 10−5), and PheWAS analysis revealed 33 phecodes associated with AP risk score (P < 3.08 × 10−5). ConclusionsThis study identified novel genes/loci contributing to AP and specific contributions to AP risk in men and women. Importantly, we identified additional systemic conditions significantly associated with AP risk. Our findings provide strong evidence for host-mediated effects on AP susceptibility.

Genetic Ancestry Estimates within Dutch Family Units and Across Genotyping Arrays: Insights from Empirical Analysis Using Two Estimation Methods.

Accurate inference of genetic ancestry is crucial for population-based association studies, accounting for population heterogeneity and structure. This study analyzes genome-wide SNP data from the Netherlands Twin Register to compare genetic ancestry estimates. The focus is on the comparison of ancestry estimates between family members and individuals genotyped on multiple arrays (Affymetrix 6.0, Affymetrix Axiom, and Illumina GSA). Two conventional methods, principal component analysis and ADMIXTURE, were implemented to estimate ancestry, each serving its specific purpose, rather than for direct comparison. The results reveal that as the degree of genetic relatedness decreases, the Euclidean distances of genetic ancestry estimates between family members significantly increase (empirical p < 0.001), regardless of the estimation method and genotyping array. Ancestry estimates among individuals genotyped on multiple arrays also show statistically significant differences (empirical p < 0.001). Additionally, this study investigates the relationship between the ancestry estimates of non-identical twin offspring with ancestrally diverse parents and those with ancestrally similar parents. The results indicate a statistically significant weak correlation between the variation in ancestry estimates among offspring and differences in ancestry estimates among parents (Spearman's rho: 0.07, p = 0.005). This study highlights the utility of current methods in inferring genetic ancestry, emphasizing the importance of reference population composition in determining ancestry estimates.

How we look: European wild mouflon and feral domestic sheep hybrids.

Hybridization between wild animals and feral individuals from closely related domestic species can occur when they share the same habitat. Such events are generally regarded as a threat to the genetic integrity and survival of established wild populations. The aim of this study is to confirm a hybridization between European mouflon and feral sheep on the Adriatic island Dugi Otok, where a group of individuals phenotypically indicative of a possible introgression was photographed. Using an ovine medium-density SNP array genotypes, we provided genetic evidence of recent gene admixture between the two species. The ADMIXTURE analysis suggests that two hybrid individuals had 82.2% and 94.1% of their ancestry attributable to feral sheep, respectively, with the rest of the mouflon in origin. NEWHYBRIDS analysis estimated the combined probabilities for the putative hybrid individuals belonging to one of the hybrid classes at 87.0% and 90.1%, respectively. The main phenotypic differences in the hybrids were observed in the shape of the horns, a stronger horn circumference than in the mouflon and weaker than for the sheep, pronounced horn rings, a shorter tail, and intermediate color. We therefore recommend constant monitoring of the purity of the mouflon population to maintain its genetic integrity. In addition, all feral sheep should be removed from the wild, and released sheep raised in semi-wild conditions should be better controlled.

Evaluating the potential of (epi)genotype-by-low pass nanopore sequencing in dairy cattle: a study on direct genomic value and methylation analysis

BackgroundGenotype-by-sequencing has been proposed as an alternative to SNP genotyping arrays in genomic selection to obtain a high density of markers along the genome. It requires a low sequencing depth to be cost effective, which may increase the error at the genotype assigment. Third generation nanopore sequencing technology offers low cost sequencing and the possibility to detect genome methylation, which provides added value to genotype-by-sequencing. The aim of this study was to evaluate the performance of genotype-by-low pass nanopore sequencing for estimating the direct genomic value in dairy cattle, and the possibility to obtain methylation marks simultaneously.ResultsLatest nanopore chemistry (LSK14 and Q20) achieved a modal base calling accuracy of 99.55%, whereas previous kit (LSK109) achieved slightly lower accuracy (99.1%). The direct genomic value accuracy from genotype-by-low pass sequencing ranged between 0.79 and 0.99, depending on the trait (milk, fat or protein yield), with a sequencing depth as low as 2 × and using the latest chemistry (LSK114). Lower sequencing depth led to biased estimates, yet with high rank correlations. The LSK109 and Q20 achieved lower accuracies (0.57–0.93). More than one million high reliable methylated sites were obtained, even at low sequencing depth, located mainly in distal intergenic (87%) and promoter (5%) regions.ConclusionsThis study showed that the latest nanopore technology in useful in a LowPass sequencing framework to estimate direct genomic values with high reliability. It may provide advantages in populations with no available SNP chip, or when a large density of markers with a wide range of allele frequencies is needed. In addition, low pass sequencing provided nucleotide methylation status of > 1 million nucleotides at ≥ 10 × , which is an added value for epigenetic studies.

Filling the gap: Genetic risk assessment in hypercholesterolemia using LDL-C and LPA genetic scores.

Routine genetic testing in hypercholesterolemia patients reveals a causative monogenic variant in less than 50% of affected individuals. Incomplete genetic characterization is partly due to polygenic factors influencing low-density-lipoprotein-cholesterol (LDL-C). Additionally, functional variants in the LPA gene affect lipoprotein(a)-associated cholesterol concentrations but are difficult to determine due to the complex structure of the LPA gene. In this study we examined whether complementing standard sequencing with the analysis of genetic scores associated with LDL-C and Lp(a) concentrations improves the diagnostic output in hypercholesterolemia patients. 1.020 individuals including 252 clinically diagnosed hypercholesterolemia patients from the FH Register Austria were analyzed by massive-parallel-sequencing of candidate genes combined with array genotyping, identifying nine novel variants in LDLR. For each individual, validated genetic scores associated with elevated LDL-C and Lp(a) were calculated based on imputed genotypes. Integrating these scores especially the score for Lp(a) increased the proportion of individuals with a clearly defined disease etiology to 68.8% compared to 46.6% in standard genetic testing. The study highlights the major role of Lp(a) in disease etiology in clinically diagnosed hypercholesterolemia patients, of which parts are misclassified. Screening for monogenic causes of hypercholesterolemia and genetic scores for LDL-C and Lp(a) permits more precise diagnosis, allowing individualized treatment.

Schizophrenia-associated somatic copy-number variants from 12,834 cases reveal recurrent NRXN1 and ABCB11 disruptions

While germline copy-number variants (CNVs) contribute to schizophrenia (SCZ) risk, the contribution of somatic CNVs (sCNVs)-present in some but not all cells-remains unknown. We identified sCNVs using blood-derived genotype arrays from 12,834 SCZ cases and 11,648 controls, filtering sCNVs at loci recurrently mutated in clonal blood disorders. Likely early-developmental sCNVs were more common in cases (0.91%) than controls (0.51%, p= 2.68e-4), with recurrent somatic deletions of exons 1-5 of the NRXN1 gene in five SCZ cases. Hi-C maps revealed ectopic, allele-specific loops forming between a potential cryptic promoter and non-coding cis-regulatory elements upon 5' deletions in NRXN1. We also observed recurrent intragenic deletions of ABCB11, encoding a transporter implicated in anti-psychotic response, in five treatment-resistant SCZ cases and showed that ABCB11 is specifically enriched in neurons forming mesocortical and mesolimbic dopaminergic projections. Our results indicate potential roles of sCNVs in SCZ risk.

Bacterial Community Modifies Host Genetics Effect on Early Childhood Caries

By age 5, approximately one-fifth of children have early childhood caries (ECC). Both the oral microbiome and host genetics are thought to influence susceptibility. Whether the oral microbiome modifies genetic susceptibility to ECC has not been tested. We test whether the salivary bacteriome modifies the association of a polygenic score (PGS, a score derived from genomic data that summarizes genetic susceptibility to disease) for primary tooth decay on ECC in the Center for Oral Health Research in Appalachia 2 longitudinal birth cohort. Children were genotyped using the Illumina Multi-Ethnic Genotyping Array and underwent annual dental examinations. We constructed a PGS for primary tooth decay using weights from an independent, genome-wide association meta-analysis. Using Poisson regression, we tested for associations between the PGS (high versus low) and ECC incidence, adjusting for demographic characteristics (n = 783). An incidence-density sampled subset of the cohort (n = 138) had salivary bacteriome data at 24 mo of age. We tested for effect modification of the PGS on ECC case status by salivary bacterial community state type (CST). By 60 mo, 20.69% of children had ECC. High PGS was not associated with an increased rate of ECC (incidence rate ratio, 1.09; 95% confidence interval [CI], 0.83–1.42). However, having a cariogenic salivary bacterial CST at 24 mo was associated with ECC (odds ratio [OR], 7.48; 95% CI, 3.06–18.26), which was robust to PGS adjustment. An interaction existed between the salivary bacterial CST and the PGS on the multiplicative scale (P = 0.04). The PGS was associated with ECC (OR, 4.83; 95% CI, 1.29–18.17) only among individuals with a noncariogenic salivary bacterial CST (n = 70). Genetic causes of caries may be harder to detect when not accounting for cariogenic oral microbiomes. As certain salivary bacterial CSTs increased ECC risk across genetic risk strata, preventing colonization of cariogenic microbiomes would be universally beneficial.

Application of a non-invasive tool and identification of genetic markers in swine to enhance ham quality: a preliminary study.

A genome-wide association study (GWAS) was performed as a preliminary step to identify regions potentially related to ham quality traits. In this research, genomic information was obtained from 238 commercial hybrid pigs utilising the GeneSeek® Genomic Profiler genome-wide porcine genotyping array. Carcasses were tested for hot weight, the thickness of backfat and loin, and lean meat percentage. The corresponding fresh hams were assayed for weight and ultimate pH; the activities of Cathepsin B and Ferrochelatase of Semimembranosus muscle were determined through fluorimetric methods. The lean meat percentage of fresh ham (LMPH), salt absorbed after first (SALT1) and overall salting stages (SALT) were estimated online by the Ham Inspector™ apparatus. Hams were processed in compliance with the procedures established for Protected Designation of Origin Parma ham, and ham weight losses were measured at the main processing stages. Hot carcass weights showed a significant negative correlation with their lean meat percentage and LMPH, while LMPH was correlated positively with carcass lean meat, SALT1, SALT, and weight losses. The GWAS detected genome-wide association for 12 single nucleotide polymorphisms with Ferrochelatase activity. The results obtained in this preliminary study were achieved by combining innovative and non-destructive technologies for screening hams under processing, measures of enzymatic muscle properties relevant to dry-cured ham quality, and genomic information obtained through a GWAS. Additional studies carried out in a larger number of pigs have been planned to investigate the effect of gene variants of Ferrochelatase activity in the dry-cured ham's quality with main reference to colour development and to confirm the GWAS results obtained in this study.

Evaluation of the selection of key individuals for genotype imputation in Chinese yellow-feathered chicken

Genotype imputation is a powerful technique employed by next-generation sequencing (NGS) and genotyping arrays, which can significantly enhance the cost-effectiveness and efficiency of genomic selection. The accuracy of imputation is largely determined by the choice of reference panel, with previous studies generally demonstrating that a closely related population as a reference panel leads to greater accuracy than a more distantly related population. Various strategies have been proposed for selecting desirable individuals via targeted resequencing, but their efficiencies need further improvement. In this study, we present a practical broiler selection methodology for a local Chinese chicken line that integrates established methods based on pedigree, genomics, and random sampling, and leverages genotype and pedigree information from the yellow-plumage dwarf chicken line. The efficacy of these selection strategies was assessed by evaluating their ability to accurately impute masked genotypes from data obtained using a 600K chip. Our findings reveal that the pedigree-based method yields superior accuracy in genotype imputation, whereas the haplotype-based method exhibits greater stability. Nonetheless, the impact of these targeted methods for selecting key individuals is slightly different when initiating a new sequencing project in a production context. Overall, this study highlights the advantages of using the pedigree-based approach as the preferred method for optimizing genotype imputation in broiler chickens.

Detection of genome-wide copy number variation in Murrah buffaloes

Riverine Buffaloes, especially the Murrah breed because of their adaptability to harsh climatic conditions, is farmed in many countries to convert low-quality feed into valuable dairy products and meat. Here, we investigated the copy number variations (CNVs) in 296 Murrah buffalo using the Axiom® Buffalo Genotyping Array 90K (Affymetrix, Santa Clara, CA, USA). The CNVs were detected on the autosomes, using the Copy Number Analysis Module (CNAM) using the univariate analysis. 7937 CNVs were detected in 279 Buffaloes, the average length of the CNVs was 119,048.87 bp that ranged between 7800 and 4,561,030 bp. These CNVs were accounting for 10.33% of the buffalo genome, which was comparable to cattle, sheep, and goat CNV analyses. Further, CNVs were merged and 1541 CNVRs were detected using the Bedtools-mergeBed command. 485 genes were annotated within 196 CNVRs that were identified in at least 10 animals of Murrah population. Out of these, 40 CNVRs contained 59 different genes that were associated with 69 different traits. Overall, the study identified a significant number of CNVs and CNVRs in the Murrah breed of buffalo, with a wide range of lengths and frequencies across the autosomes. The identified CNVRs contained genes associated with important traits related to production and reproduction, making them potentially important targets for future breeding and genetic improvement efforts.

RHAMM regulates MMTV-PyMT-induced lung metastasis by connecting STING-dependent DNA damage sensing to interferon/STAT1 pro-apoptosis signaling

BackgroundRHAMM is a multifunctional protein that is upregulated in breast tumors, and the presence of strongly RHAMM+ve cancer cell subsets associates with elevated risk of peripheral metastasis. Experimentally, RHAMM impacts cell cycle progression and cell migration. However, the RHAMM functions that contribute to breast cancer metastasis are poorly understood.MethodsWe interrogated the metastatic functions of RHAMM using a loss-of-function approach by crossing the MMTV-PyMT mouse model of breast cancer susceptibility with Rhamm−/− mice. In vitro analyses of known RHAMM functions were performed using primary tumor cell cultures and MMTV-PyMT cell lines. Somatic mutations were identified using a mouse genotyping array. RNA-seq was performed to identify transcriptome changes resulting from Rhamm-loss, and SiRNA and CRISPR/Cas9 gene editing was used to establish cause and effect of survival mechanisms in vitro.ResultsRhamm-loss does not alter initiation or growth of MMTV-PyMT-induced primary tumors but unexpectedly increases lung metastasis. Increased metastatic propensity with Rhamm-loss is not associated with obvious alterations in proliferation, epithelial plasticity, migration, invasion or genomic stability. SNV analyses identify positive selection of Rhamm−/− primary tumor clones that are enriched in lung metastases. Rhamm−/− tumor clones are characterized by an increased ability to survive with ROS-mediated DNA damage, which associates with blunted expression of interferon pathway and target genes, particularly those implicated in DNA damage-resistance. Mechanistic analyses show that ablating RHAMM expression in breast tumor cells by siRNA knockdown or CRISPR-Cas9 gene editing blunts interferon signaling activation by STING agonists and reduces STING agonist-induced apoptosis. The metastasis-specific effect of RHAMM expression-loss is linked to microenvironmental factors unique to tumor-bearing lung tissue, notably high ROS and TGFB levels. These factors promote STING-induced apoptosis of RHAMM+ve tumor cells to a significantly greater extent than RHAMM−ve comparators. As predicted by these results, colony size of Wildtype lung metastases is inversely related to RHAMM expression.ConclusionRHAMM expression-loss blunts STING-IFN signaling, which offers growth advantages under specific microenvironmental conditions of lung tissue. These results provide mechanistic insight into factors controlling clonal survival/expansion of metastatic colonies and has translational potential for RHAMM expression as a marker of sensitivity to interferon therapy.

MicroRNA-Related Polymorphism and Their Association with Fibromyalgia.

MicroRNAs are tissue-specific expressed short RNAs that serve post-transcriptional gene regulation. A specific microRNA can bind to mRNAs of different genes and thereby suppress their protein production. In the context of the complex phenotype of fibromyalgia, we used the Axiom miRNA Target Site Genotyping Array to search genome-wide for DNA variations in microRNA genes, their regulatory regions, and in the 3'UTR of protein-coding genes. To identify disease-relevant DNA polymorphisms, a cohort of 176 female fibromyalgia patients was studied in comparison to a cohort of 162 healthy women. The association between 48,329 markers and fibromyalgia was investigated using logistic regression adjusted for population stratification. Results show that 29 markers had p-values < 1 × 10-3, and the strongest association was observed for rs758459 (p-value of 0.0001), located in the Neurogenin 1 gene which is targeted by hsa-miR-130a-3p. Furthermore, variant rs2295963 is predicted to affect binding of hsa-miR-1-3p. Both microRNAs were previously reported to be differentially expressed in fibromyalgia patients. Despite its limited statistical power, this study reports two microRNA-related polymorphisms which may play a functional role in the pathogenesis of fibromyalgia. For a better understanding of the disease pattern, further functional analyses on the biological significance of microRNAs and microRNA-related polymorphisms are required.

Genetic diversity and selection signatures in a gene bank panel of maize inbred lines from Southeast Europe compared with two West European panels

Southeast Europe (SEE) is a very important maize-growing region, comparable to the Corn belt region of the United States, with similar dent germplasm (dent by dent hybrids). Historically, this region has undergone several genetic material swaps, following the trends in the US, with one of the most significant swaps related to US aid programs after WWII. The imported accessions used to make double-cross hybrids were also mixed with previously adapted germplasm originating from several more distant OPVs, supporting the transition to single cross-breeding. Many of these materials were deposited at the Maize Gene Bank of the Maize Research Institute Zemun Polje (MRIZP) between the 1960s and 1980s. A part of this Gene Bank (572 inbreds) was genotyped with Affymetrix Axiom Maize Genotyping Array with 616,201 polymorphic variants. Data were merged with two other genotyping datasets with mostly European flint (TUM dataset) and dent (DROPS dataset) germplasm. The final pan-European dataset consisted of 974 inbreds and 460,243 markers. Admixture analysis showed seven ancestral populations representing European flint, B73/B14, Lancaster, B37, Wf9/Oh07, A374, and Iodent pools. Subpanel of inbreds with SEE origin showed a lack of Iodent germplasm, marking its historical context. Several signatures of selection were identified at chromosomes 1, 3, 6, 7, 8, 9, and 10. The regions under selection were mined for protein-coding genes and were used for gene ontology (GO) analysis, showing a highly significant overrepresentation of genes involved in response to stress. Our results suggest the accumulation of favorable allelic diversity, especially in the context of changing climate in the genetic resources of SEE.

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

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. 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. 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. 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.

South Asian medical cohorts reveal strong founder effects and high rates of homozygosity

The benefits of large-scale genetic studies for healthcare of the populations studied are well documented, but these genetic studies have traditionally ignored people from some parts of the world, such as South Asia. Here we describe whole genome sequence (WGS) data from 4806 individuals recruited from the healthcare delivery systems of Pakistan, India and Bangladesh, combined with WGS from 927 individuals from isolated South Asian populations. We characterize population structure in South Asia and describe a genotyping array (SARGAM) and imputation reference panel that are optimized for South Asian genomes. We find evidence for high rates of reproductive isolation, endogamy and consanguinity that vary across the subcontinent and that lead to levels of rare homozygotes that reach 100 times that seen in outbred populations. Founder effects increase the power to associate functional variants with disease processes and make South Asia a uniquely powerful place for population-scale genetic studies.