Population structure of indigenous inhabitants of Arabia

Near-complete Middle Eastern genomes refine autozygosity and enhance disease-causing and population-specific variant discovery

Advances in long-read sequencing have enabled routine complete assembly of human genomes, but much remains to be done to represent broader populations and show impact on disease-gene discovery. Here, we report highly accurate, near-complete and phased genomes from six Middle Eastern (ME) family trios (n = 18) with neurodevelopmental conditions, representing ancestries from Sudan, Jordan, Syria, Qatar and Afghanistan. These genomes revealed 42.2 Mb of new sequence (13.8% impacting known genes), 75 new HLA/KIR alleles and strong signals of inbreeding, with ROH covering up to one-third of chromosomes 6 and 12 in one individual. Using assembly-based variant calling, we identified 23 de novo and recessive variants as strong candidates for causing previously unresolved symptoms in the probands. The ME genomes revealed unique variation relative to existing references, showing enhanced mappability and variant calling. These results underscore the value of de novo assembly for disease variant discovery and the need for sampled ME-specific references to better characterize population-relevant variation.

Validity of the Arabic version of AAOS-foot and ankle outcomes questionnaire in patients with traumatic foot and ankle injuries

BACKGROUND

Arabic-speaking patients are underrepresented in orthopedic clinical studies, particularly in foot and ankle trauma research. The lack of validated Arabic language tools hinders their inclusion, creating a need for culturally and linguistically adapted instruments. The American Academy of Orthopedic Surgeons Foot and Ankle Outcomes Questionnaire (AAOS-FAOQ) is a widely used tool but has not been adapted for Arabic-speaking patients.

AIM

To translate, cross-culturally adapt, and validate the AAOS-FAOQ for Arabic-speaking patients with traumatic foot and ankle injuries.

METHODS

The cross-cultural adaptation followed established guidelines, involving forward and backward translations, expert review, and pre-testing. The final Arabic version was administered alongside the Arabic Short-Form 36 (SF-36) to 100 patients for validity testing. Reliability was assessed through test-retest methods with 20 patients completing the questionnaire twice within 48 hours. Pearson correlation coefficients measured convergent and divergent validity with SF-36 subscales, while Cronbach's alpha and intraclass correlation coefficients (ICC) determined internal consistency and reliability.

RESULTS

Out of 100 patients, 92 completed the first set of questionnaires. The Arabic AAOS-FAOQ showed strong correlations with the SF-36 subscales, particularly in physical function and bodily pain (r > 0.6). Test-retest reliability was robust, with ICCs of 0.69 and 0.66 for the Global Foot and Ankle Scale and Shoe Comfort Scale, respectively. Cronbach's alpha for internal consistency ranged from 0.7 to 0.9.

CONCLUSION

The Arabic version of the AAOS-FAOQ demonstrated validity and reliability for use in Arabic-speaking patients with traumatic foot and ankle injuries. This adaptation will enhance the inclusion of this population in orthopedic clinical studies, improving the generalizability of research findings and patient care.

Arab founder variants: Contributions to clinical genomics and precision medicine

BACKGROUND

Founder variants are ancestral variants shared by individuals who are not closely related. The large effect size of some of these variants in the context of Mendelian disorders offers numerous precision medicine opportunities.

METHODS

Using one of the largest datasets on Mendelian disorders in the Middle East, we identified 2,908 medically relevant founder variants derived from 18,360 exomes and genomes and investigated their contribution to the clinical annotation of the human genome.

FINDINGS

Strikingly, ∼34% of Arab founder variants are absent in gnomAD. We found a strong contribution of Arab founder variants to the identification of novel gene-disease links (n = 224) and the support/dispute (n = 81 support, n = 101 dispute) of previously reported candidate gene-disease links. The powerful segregation evidence generated by Arab founder variants allowed many ClinVar and Human Gene Mutation Database variants to be reclassified. Overall, 39.5% of diagnostic reports from our clinical lab are based on founder variants, and 19.41% of tested individuals carry at least one pathogenic founder variant. The presumptive loss-of-function mechanism that typically underlies autosomal recessive diseases means that Arab founder variants also offer unique opportunities in "druggable genome" research. Arab founder variants were also informative of migration patterns in the Middle East consistent with documented historical accounts.

CONCLUSIONS

We highlight the contribution of founder variants from an under-represented population group to precision medicine and inform future prevention programs. Our study also sheds light on the added value of these variants in supplementing other lines of research in tracing population history.

FUNDING

There is no funding for this work.

Patterns of population structure and genetic variation within the Saudi Arabian population

The Arabian Peninsula is considered the initial site of historic human migration out of Africa. The modern-day indigenous Arabians are believed to be the descendants who remained from the ancient split of the migrants into Eurasia. Here, we investigated how the population history and cultural practices such as endogamy have shaped the genetic variation of the Saudi Arabians. We genotyped 3,352 individuals and identified twelve genetic sub-clusters that corresponded to the geographical distribution of different tribal regions, differentiated by distinct components of ancestry based on comparisons to modern and ancient DNA references. These sub-clusters also showed variation across ranges of the genome covered in runs of homozygosity, as well as differences in population size changes over time. Using 25,488,981 variants found in whole genome sequencing data (WGS) from 302 individuals, we found that the Saudi tend to show proportionally more deleterious alleles than neutral alleles when compared to Africans/African Americans from gnomAD (e.g. a 13% increase of deleterious alleles annotated by AlphaMissense between 0.5 - 5% frequency in Saudi, compared to 7% decrease of the benign alleles; P < 0.001). Saudi sub-clusters with greater inbreeding and lower effective population sizes showed greater enrichment of deleterious alleles as well. Additionally, we found that approximately 10% of the variants discovered in our WGS data are not observed in gnomAD; these variants are also enriched with deleterious annotations. To accelerate studying the population-enriched deleterious alleles and their health consequences in this population, we made available the allele frequency estimates of 25,488,981 variants discovered in our samples. Taken together, our results suggest that Saudi's population history impacts its pattern of genetic variation with potential consequences to the population health. It further highlights the need to sequence diverse and unique populations so to provide a foundation on which to interpret medical- and pharmaco- genomic findings from these populations.

A reference quality, fully annotated diploid genome from a Saudi individual

We have used multiple sequencing approaches to sequence the genome of a volunteer from Saudi Arabia. We use the resulting data to generate a de novo assembly of the genome, and use different computational approaches to refine the assembly. As a consequence, we provide a contiguous assembly of the complete genome of an individual from Saudi Arabia for all chromosomes except chromosome Y, and label this assembly KSA001. We transferred genome annotations from reference genomes to fully annotate KSA001, and we make all primary sequencing data, the assembly, and the genome annotations freely available in public databases using the FAIR data principles. KSA001 is the first telomere-to-telomere-assembled genome from a Saudi individual that is freely available for any purpose.

The impact of consanguinity on the design of iPSC banks

The effect of consanguinity on identifying universal induced pluripotent stem cell (iPSC) donors, i.e., homozygous for the major human leukocyte antigen (HLA) loci, is unknown. The discovery sample size was calculated in a consanguineous population using a method (1qF) based on the inbreeding coefficient. The result was orders of magnitude smaller compared to the standard method.

Founder mutations and rare disease in the Arab world

Founder mutations are disease-causing variants that occur frequently in geographically or culturally isolated groups whose shared ancestor(s) carried the pathogenic variant. While some disease alleles may vanish from the genetic pool due to natural selection, variants with weaker effects may survive for a long time, thereby enhancing the prevalence of some rare diseases. These are predominantly autosomal recessive diseases but can also be autosomal dominant traits with late-onset or mild phenotypes. Cultural practices, such as endogamy and consanguinity, in these isolated groups lead to higher prevalence of such rare diseases compared to the rest of the population and worldwide. In this Perspective, we define population isolates and the underlying genetic mechanisms for accumulating founder mutations. We also discuss the current and potential scientific, clinical and public-health implications of studying founder mutations in population isolates around the world, with a particular focus on the Arab population.

Imputation accuracy across global human populations

Genotype imputation is now fundamental for genome-wide association studies but lacks fairness due to the underrepresentation of references from non-European ancestries. The state-of-the-art imputation reference panel released by the Trans-Omics for Precision Medicine (TOPMed) initiative improved the imputation of admixed African-ancestry and Hispanic/Latino samples, but imputation for populations primarily residing outside of North America may still fall short in performance due to persisting underrepresentation. To illustrate this point, we imputed the genotypes of over 43,000 individuals across 123 populations around the world and identified numerous populations where imputation accuracy paled in comparison to that of European-ancestry populations. For instance, the mean imputation r-squared (Rsq) for variants with minor allele frequencies between 1% and 5% in Saudi Arabians (n = 1,061), Vietnamese (n = 1,264), Thai (n = 2,435), and Papua New Guineans (n = 776) were 0.79, 0.78, 0.76, and 0.62, respectively, compared to 0.90-0.93 for comparable European populations matched in sample size and SNP array content. Outside of Africa and Latin America, Rsq appeared to decrease as genetic distances to European-ancestry reference increased, as predicted. Using sequencing data as ground truth, we also showed that Rsq may over-estimate imputation accuracy for non-European populations more than European populations, suggesting further disparity in accuracy between populations. Using 1,496 sequenced individuals from Taiwan Biobank as a second reference panel to TOPMed, we also assessed a strategy to improve imputation for non-European populations with meta-imputation, but this design did not improve accuracy across frequency spectra. Taken together, our analyses suggest that we must ultimately strive to increase diversity and size to promote equity within genetics research.

Pragmatic Approach to Applying Polygenic Risk Scores to Diverse Populations

Polygenic risk scores (PRS) estimate genetic susceptibility of an individual to disease and have the potential of providing utility in multiple clinical contexts. However, their performance, computation, and reporting in diverse populations remain challenging. Here, we present a pragmatic approach to optimize a PRS for a population of interest that leverages publicly available data and methods and consists of seven steps that are easily implemented without the requirement of expertise in complex genetics: step 1, selecting source genome-wide association studies (GWAS) and imputation; step 2, selecting methods to compute polygenic score; step 3, adjusting scores using principal components of genetic ancestry; step 4, selecting the best performing score; step 5, defining percentiles of a population distribution; step 6, validating performance of the optimized polygenic score; and step 7, implementing the optimized polygenic score in clinical practice. © 2023 Wiley Periodicals LLC.

Imputation Accuracy Across Global Human Populations

Genotype imputation is now fundamental for genome-wide association studies but lacks fairness due to the underrepresentation of populations with non-European ancestries. The state-of-the-art imputation reference panel released by the Trans-Omics for Precision Medicine (TOPMed) initiative contains a substantial number of admixed African-ancestry and Hispanic/Latino samples to impute these populations with nearly the same accuracy as European-ancestry cohorts. However, imputation for populations primarily residing outside of North America may still fall short in performance due to persisting underrepresentation. To illustrate this point, we curated genome-wide array data from 23 publications published between 2008 to 2021. In total, we imputed over 43k individuals across 123 populations around the world. We identified a number of populations where imputation accuracy paled in comparison to that of European-ancestry populations. For instance, the mean imputation r-squared (Rsq) for 1-5% alleles in Saudi Arabians (N=1061), Vietnamese (N=1264), Thai (N=2435), and Papua New Guineans (N=776) were 0.79, 0.78, 0.76, and 0.62, respectively. In contrast, the mean Rsq ranged from 0.90 to 0.93 for comparable European populations matched in sample size and SNP content. Outside of Africa and Latin America, Rsq appeared to decrease as genetic distances to European reference increased, as predicted. Further analysis using sequencing data as ground truth suggested that imputation software may over-estimate imputation accuracy for non-European populations than European populations, suggesting further disparity between populations. Using 1496 whole genome sequenced individuals from Taiwan Biobank as a reference, we also assessed a strategy to improve imputation for non-European populations with meta-imputation, which can combine results from TOPMed with smaller population-specific reference panels. We found that meta-imputation in this design did not improve Rsq genome-wide. Taken together, our analysis suggests that with the current size of alternative reference panels, meta-imputation alone cannot improve imputation efficacy for underrepresented cohorts and we must ultimately strive to increase diversity and size to promote equity within genetics research.

Clinical utility of polygenic scores for cardiometabolic disease in Arabs

Arabs account for 5% of the world population and have a high burden of cardiometabolic disease, yet clinical utility of polygenic risk prediction in Arabs remains understudied. Among 5399 Arab patients, we optimize polygenic scores for 10 cardiometabolic traits, achieving a performance that is better than published scores and on par with performance in European-ancestry individuals. Odds ratio per standard deviation (OR per SD) for a type 2 diabetes score was 1.83 (95% CI 1.74-1.92), and each SD of body mass index (BMI) score was associated with 1.18 kg/m2 difference in BMI. Polygenic scores associated with disease independent of conventional risk factors, and also associated with disease severity-OR per SD for coronary artery disease (CAD) was 1.78 (95% CI 1.66-1.90) for three-vessel CAD and 1.41 (95% CI 1.29-1.53) for one-vessel CAD. We propose a pragmatic framework leveraging public data as one way to advance equitable clinical implementation of polygenic scores in non-European populations.

A transcriptomic signature of X chromosome overdosage in Saudi Klinefelter syndrome induced pluripotent stem cells

OBJECTIVE

The transcriptional landscape of Klinefelter syndromeduring early embryogenesis remains elusive. This study aimed to evaluate the impact of X chromosome overdosage in 47,XXY males induced pluripotent stem cells (iPSCs) obtained from patients with different genomic backgrounds and ethnicities.

DESIGN AND METHOD

We derived and characterized 15 iPSC lines from four Saudi 47,XXY KS patients and one Saudi 46,XY male. We performed a comparative transcriptional analysis using the Saudi KS-iPSCs and a cohort of European and North American KS-iPSCs.

RESULTS

We identified a panel of X-linked and autosomal genes commonly dysregulated in Saudi and European/North American KS-iPSCs vs 46,XY controls. Our findings demonstrate that seven PAR1 and nine non-PAR escape genes are consistently dysregulated and mostly display comparable transcriptional levels in both groups. Finally, we focused on genes commonly dysregulated in both iPSC cohorts and identified several gene-ontology categories highly relevant to KS physiopathology, including aberrant cardiac muscle contractility, skeletal muscle defects, abnormal synaptic transmission, and behavioral alterations.

CONCLUSIONS

Our results indicate that a transcriptomic signature of X chromosome overdosage in KS is potentially attributable to a subset of X-linked genes sensitive to sex chromosome dosage and escaping X inactivation, regardless of the geographical area of origin, ethnicity, and genetic makeup.

Portability of Polygenic Risk Scores for Sleep Duration, Insomnia and Chronotype in 33,493 Individuals

Polygenic risk scores (PRSs) estimate genetic liability for diseases and traits. However, the portability of PRSs in sleep traits has remained elusive. We generated PRSs for self-reported insomnia, chronotype and sleep duration using summary data from genome-wide association studies (GWASs) performed in 350,000 to 697,000 European-ancestry individuals. We then projected the scores in two independent Finnish population cohorts (N = 33,493) and tested whether the PRSs were associated with their respective sleep traits. We observed that all the generated PRSs were associated with their corresponding traits (p < 0.05 in all cases). Furthermore, we found that there was a 22.2 min difference in reported sleep between the 5% tails of the PRS for sleep duration (p < 0.001). Our findings indicate that sleep-related PRSs show portability across cohorts. The findings also demonstrate that sleep measures using PRSs for sleep behaviors may provide useful instruments for testing disease and trait associations in cohorts where direct sleep parameters have not yet been measured.

EPIDEMIOLOGICAL DETERMINANTS FOR THE SPREAD OF COVID-19 IN RIYADH PROVINCE OF SAUDI ARABIA

Background and objective

Coronavirus Disease 2019 (COVID-19) has affected millions of individuals all over the world. In addition to the patients' compelling indications, various sociodemographic characteristics were identified to influence infection complications. The purpose of this study was to assess the impact of the aforementioned parameters on the dissemination of COVID-19 among residents of Saudi Arabia's Riyadh region.

Materials and Methods

In the Saudi Arabian province of Riyadh, a cross-sectional retrospective analysis of COVID-19 incidences, recoveries, and case-fatality ratio (CFR) was undertaken. The study was carried out by gathering daily COVID-19 records from the ministry of health's official websites between October 2020 and September 2021. The influencing factors were obtained from the statistical authority. Using the SPSS IBM 25 software, the data was examined. The association between demographic factors as well as the presence of comorbidity on the COVID-19 outcome was determined using Spearman's correlation and regression tests. P<0.05 was considered to indicate the significance of the results.

Results

The data from the study indicated that the highest number of COVID-19 cases were recorded in June 2021, and peak recovery was observed in July 2021. The CFR declined progressively from October 2020 to just over 1, even when the cases peaked. A significant (p<0.05) correlation between diabetes and COVID-19 incidences was observed . The recovery rate had a significant (p<0.05) association with the literacy rate and those aged 14 - 49 years old. Presences of co-morbidities such as Dyslipidemia, hypertension, diabetes, asthma, stroke and heart failure have negatively affected the recovery from COVID-19 in the population. The CFR is significantly (p< 0.05) associated with people over 60, hypertensive patients, and asthma patients. Regression analysis suggested that the risk of complications due to COVID-19 infection is more in males, people above 60 years age and those suffering from co-morbidities .

Conclusions

The findings of the study indicate an association between several of the characteristics studied, such as gender, age, and comorbidity, and the spread of infection, recovery, and mortality. To restrict the spread of COVID-19 and prevent its complications, effective measures are required to control the modifiable risk factors.

Expression of Endogenous Retroviral RNA in Prostate Tumors has Prognostic Value and Shows Differences among Americans of African Versus European/Middle Eastern Ancestry

Simple Summary

Endogenous retroviruses (ERVs) are viral sequences that have been incorporated into the human genome over millions of years via integrations in germ-line cells. In this study, we investigated whether the expression of ERVs was associated with two different aspects of prostate cancer (PCa). First, Black American men have a higher incidence and poorer outcome of PCa compared to White men. We identified differences in ERV expression among prostate tumors between men of primarily African versus primarily European or Middle Eastern ancestry, which may be associated with differences in the mechanism of cancer progression in patients of these distinct ancestries. Second, we determined whether ERV expression might be correlated with the progression of disease, regardless of ancestry. We identified the ERV expression signatures that correlated with biochemical relapse among PCa patients of all ancestries, indicating that ERVs may be useful for identifying cancer patients at greatest risk of progression. The utility of ERV expression for studying cancer progression may extend to other cancers.

Abstract

Endogenous retroviruses (ERVs) are abundant, repetitive elements dispersed across the human genome and are implicated in various diseases. We investigated two potential roles for ERVs in prostate cancer (PCa). First, the PCa of Black Americans (BA) is diagnosed at an earlier median age and at a more advanced stage than the PCa of White Americans (WA). We used publicly available RNA-seq data from tumor-enriched samples of 27 BA and 65 WA PCa patients in order to identify 12 differentially expressed ERVs (p_adj < 0.1) and used a tissue microarray of the PCa cores from an independent set of BA and WA patients to validate the differential protein expression of one of these ERVs, ERV3-1 (p = 2.829 × 10⁻⁷). Second, we used 57 PCa tumors from patients of all ancestries from one hospital as a training set to identify the ERVs associated with time to biochemical relapse. A 29-ERV prognostic panel was then tested and validated on 35 separate PCa tumors from patients obtained in two different hospitals with a dramatic increase in prognostic power relative to clinical parameters alone (p = 7.4 × 10⁻¹¹). In summary, ERV RNA expression differences in the prostate tumors of patients of different ancestries may be associated with dissimilarities in the mechanism of cancer progression. In addition, the correlation of expression of certain ERVs in prostate tumors with the risk of biochemical relapse indicates a possible role for ERV expression in cancer progression.

Thousands of Qatari genomes inform human migration history and improve imputation of Arab haplotypes

Arab populations are largely understudied, notably their genetic structure and history. Here we present an in-depth analysis of 6,218 whole genomes from Qatar, revealing extensive diversity as well as genetic ancestries representing the main founding Arab genealogical lineages of Qahtanite (Peninsular Arabs) and Adnanite (General Arabs and West Eurasian Arabs). We find that Peninsular Arabs are the closest relatives of ancient hunter-gatherers and Neolithic farmers from the Levant, and that founder Arab populations experienced multiple splitting events 12–20 kya, consistent with the aridification of Arabia and farming in the Levant, giving rise to settler and nomadic communities. In terms of recent genetic flow, we show that these ancestries contributed significantly to European, South Asian as well as South American populations, likely as a result of Islamic expansion over the past 1400 years. Notably, we characterize a large cohort of men with the ChrY J1a2b haplogroup (n = 1,491), identifying 29 unique sub-haplogroups. Finally, we leverage genotype novelty to build a reference panel of 12,432 haplotypes, demonstrating improved genotype imputation for both rare and common alleles in Arabs and the wider Middle East.

Arab populations are relatively understudied, especially their genetic architecture and historical relationship with early founders of the ancient Near East. Here, the authors examine 6,218 Qatari whole genomes, revealing insights on migration, population history and genetic structure of populations across the Middle Eastern region.