Ancestry informative marker sets for determining continental origin and admixture proportions in common populations in America

Dense single nucleotide polymorphism testing revolutionizes scope and degree of certainty for source attribution in forensic investigations

The field of forensic DNA analysis has experienced significant advancements over the years, such as the advent of DNA fingerprinting, the introduction of the polymerase chain reaction for increased sensitivity, the shift to a primary genetic marker system based on short tandem repeats, and implementation of national DNA databases. Now, the forensics field is poised for another revolution with the advent of dense single nucleotide polymorphisms (SNPs) testing. SNP testing holds the potential to significantly enhance source attribution in forensic cases, particularly those involving low-quantity or low-quality samples. When coupled with genetic genealogy and kinship analysis, it can resolve countless active cases as well as cold cases and cases of unidentified human remains, which are hindered by the limitations of existing forensic capabilities that fail to generate viable investigative leads with DNA. The field of forensic genetic genealogy combined with genome-wide sequencing can associate relatives as distant as the seventh-degree and beyond. By leveraging volunteer-populated databases to locate near and distant relatives, genetic genealogy can effectively narrow the candidates linked to crime scene evidence or aid in determining the identity of human remains. With decreasing DNA sequencing costs and improving sensitivity of detection, forensic genetic genealogy is expanding its capabilities to generate investigative leads from a wide range of biological evidence.

scAI-SNP: a method for inferring ancestry from single-cell data

Collaborative efforts, such as the Human Cell Atlas, are rapidly accumulating large amounts of single-cell data. To ensure that single-cell atlases are representative of human genetic diversity, we need to determine the ancestry of the donors from whom single-cell data are generated. Self-reporting of race and ethnicity, although important, can be biased and is not always available for the datasets already collected. Here, we introduce scAI-SNP, a tool to infer ancestry directly from single-cell genomics data. To train scAI-SNP, we identified 4.5 million ancestry-informative single-nucleotide polymorphisms (SNPs) in the 1000 Genomes Project dataset across 3201 individuals from 26 population groups. For a query single-cell data set, scAI-SNP uses these ancestry-informative SNPs to compute the contribution of each of the 26 population groups to the ancestry of the donor from whom the cells were obtained. Using diverse single-cell data sets with matched whole-genome sequencing data, we show that scAI-SNP is robust to the sparsity of single-cell data, can accurately and consistently infer ancestry from samples derived from diverse types of tissues and cancer cells, and can be applied to different modalities of single-cell profiling assays, such as single-cell RNA-seq and single-cell ATAC-seq. Finally, we argue that ensuring that single-cell atlases represent diverse ancestry, ideally alongside race and ethnicity, is ultimately important for improved and equitable health outcomes by accounting for human diversity.

1,25-Dihydroxyvitamin D3 Suppresses Prognostic Survival Biomarkers Associated with Cell Cycle and Actin Organization in a Non-Malignant African American Prostate Cell Line

Vitamin D3 is a steroid hormone that confers anti-tumorigenic properties in prostate cells. Serum vitamin D3 deficiency has been associated with advanced prostate cancer (PCa), particularly affecting African American (AA) men. Therefore, elucidating the pleiotropic effects of vitamin D on signaling pathways, essential to maintaining non-malignancy, may provide additional drug targets to mitigate disparate outcomes for men with PCa, especially AA men. We conducted RNA sequencing on an AA non-malignant prostate cell line, RC-77N/E, comparing untreated cells to those treated with 10 nM of vitamin D3 metabolite, 1α,25(OH)2D3, at 24 h. Differential gene expression analysis revealed 1601 significant genes affected by 1α,25(OH)2D3 treatment. Pathway enrichment analysis predicted 1α,25(OH)2D3- mediated repression of prostate cancer, cell proliferation, actin cytoskeletal, and actin-related signaling pathways (p < 0.05). Prioritizing genes with vitamin D response elements and associating expression levels with overall survival (OS) in The Cancer Genome Atlas Prostate Adenocarcinoma (TCGA PRAD) cohort, we identified ANLN (Anillin) and ECT2 (Epithelial Cell Transforming 2) as potential prognostic PCa biomarkers. Both genes were strongly correlated and significantly downregulated by 1α,25(OH)2D3 treatment, where low expression was statistically associated with better overall survival outcomes in the TCGA PRAD public cohort. Increased ANLN and ECT2 mRNA gene expression was significantly associated with PCa, and Gleason scores using both the TCGA cohort (p < 0.05) and an AA non-malignant/tumor-matched cohort. Our findings suggest 1α,25(OH)2D3 regulation of these biomarkers may be significant for PCa prevention. In addition, 1α,25(OH)2D3 could be used as an adjuvant treatment targeting actin cytoskeleton signaling and actin cytoskeleton-related signaling pathways, particularly among AA men.

Heterogeneous associations between interleukin-6 receptor variants and phenotypes across ancestries and implications for therapy

The Phenome-Wide Association Study (PheWAS) is increasingly used to broadly screen for potential treatment effects, e.g., IL6R variant as a proxy for IL6R antagonists. This approach offers an opportunity to address the limited power in clinical trials to study differential treatment effects across patient subgroups. However, limited methods exist to efficiently test for differences across subgroups in the thousands of multiple comparisons generated as part of a PheWAS. In this study, we developed an approach that maximizes the power to test for heterogeneous genotype-phenotype associations and applied this approach to an IL6R PheWAS among individuals of African (AFR) and European (EUR) ancestries. We identified 29 traits with differences in IL6R variant-phenotype associations, including a lower risk of type 2 diabetes in AFR (OR 0.96) vs EUR (OR 1.0, p-value for heterogeneity = 8.5 × 10-3), and higher white blood cell count (p-value for heterogeneity = 8.5 × 10-131). These data suggest a more salutary effect of IL6R blockade for T2D among individuals of AFR vs EUR ancestry and provide data to inform ongoing clinical trials targeting IL6 for an expanding number of conditions. Moreover, the method to test for heterogeneity of associations can be applied broadly to other large-scale genotype-phenotype screens in diverse populations.

Modifying Effects of Genetic Variations on the Association Between Dietary Isothiocyanate Exposure and Non-muscle Invasive Bladder Cancer Prognosis in the Be-Well Study

SCOPE

Dietary isothiocyanate (ITC) exposure from cruciferous vegetable (CV) intake may improve non-muscle invasive bladder cancer (NMIBC) prognosis. This study aims to investigate whether genetic variations in key ITC-metabolizing/functioning genes modify the associations between dietary ITC exposure and NMIBC prognosis outcomes.

METHODS AND RESULTS

In the Bladder Cancer Epidemiology, Wellness, and Lifestyle Study (Be-Well Study), a prospective cohort of 1472 incident NMIBC patients, dietary ITC exposure is assessed by self-reported CV intake and measured in plasma ITC-albumin adducts. Using Cox proportional hazards regression models, stratified by single nucleotide polymorphisms (SNPs) in nine key ITC-metabolizing/functioning genes, it is calculate hazard ratios (HRs) and 95% confidence intervals (CIs) for recurrence and progression. The rs15561 in N-acetyltransferase 1 (NAT1) is alter the association between CV intake and progression risk. Multiple SNPs in nuclear factor E2-related factor 2 (NRF2) and nuclear factor kappa B (NFκB) are modify the associations between plasma ITC-albumin adduct level and progression risk (pint < 0.05). No significant association is observed with recurrence risk. Overall, >80% study participants are present with at least one protective genotype per gene, showing an average 65% reduction in progression risk with high dietary ITC exposure.

CONCLUSION

Despite that genetic variations in ITC-metabolizing/functioning genes may modify the effect of dietary ITCs on NMIBC prognosis, dietary recommendation of CV consumption may help improve NMIBC survivorship.

Systematic analyses of AISNPs screening and classification algorithms based on genome-wide data for forensic biogeographic ancestry inference

Identifying the biogeographic ancestral origin of biological sample left at a crime scene can provide important evidence for judicial case, as well as clue for narrowing down suspect. Ancestry informative single nucleotide polymorphism (AISNP) has become one of the most important genetic markers in recent years for screening ancestry information loci and analyzing the population genetic background and structure due to their high number and wide distributions in the human genome. In this study, based on data from 26 populations in the 1000 Genomes Project Phase 3, a Random Forest classification model was constructed with one-vs-rest classification strategy for embedded feature selection in order to obtain a panel with a small number of efficient AISNPs. The research aim was to clarify differentiations of population genetic structures among continents and subregions of East Asia. ADMIXTURE results showed that based on the 58 AISNPs selected by the machine learning algorithm, the 26 populations involved in the study could be categorized into six intercontinental ancestry components: North East Asia, South East Asia, Africa, Europe, South Asia, and America. The 24 continental-specific AISNPs and 34 East Asian-specific AISNPs were finally obtained, and used to construct the ancestry prediction model using XGBoost algorithm, resulting in the Matthews correlation coefficients of 0.94 and 0.89, and accuracies of 0.94 and 0.92, respectively. The machine learning models that we constructed using population-specific AISNPs were able to accurately predict the ancestral origins of continental and intra-East Asian populations. To summarize, screening a set of high-perform AISNPs to infer biogeographical ancestral information using embedded feature selection has potential application in creating a layered inference system that accurately differentiates from intercontinental populations to local subpopulations.

Population-enriched innate immune variants may identify candidate gene targets at the intersection of cancer and cardio-metabolic disease

Both cancer and cardio-metabolic disease disparities exist among specific populations in the US. For example, African Americans experience the highest rates of breast and prostate cancer mortality and the highest incidence of obesity. Native and Hispanic Americans experience the highest rates of liver cancer mortality. At the same time, Pacific Islanders have the highest death rate attributed to type 2 diabetes (T2D), and Asian Americans experience the highest incidence of non-alcoholic fatty liver disease (NAFLD) and cancers induced by infectious agents. Notably, the pathologic progression of both cancer and cardio-metabolic diseases involves innate immunity and mechanisms of inflammation. Innate immunity in individuals is established through genetic inheritance and external stimuli to respond to environmental threats and stresses such as pathogen exposure. Further, individual genomes contain characteristic genetic markers associated with one or more geographic ancestries (ethnic groups), including protective innate immune genetic programming optimized for survival in their corresponding ancestral environment(s). This perspective explores evidence related to our working hypothesis that genetic variations in innate immune genes, particularly those that are commonly found but unevenly distributed between populations, are associated with disparities between populations in both cancer and cardio-metabolic diseases. Identifying conventional and unconventional innate immune genes that fit this profile may provide critical insights into the underlying mechanisms that connect these two families of complex diseases and offer novel targets for precision-based treatment of cancer and/or cardio-metabolic disease.

Targeted enrichment of whole-genome SNPs from highly burned skeletal remains

Genetic assessment of highly incinerated and/or degraded human skeletal material is a persistent challenge in forensic DNA analysis, including identifying victims of mass disasters. Few studies have investigated the impact of thermal degradation on whole-genome single-nucleotide polymorphism (SNP) quality and quantity using next-generation sequencing (NGS). We present whole-genome SNP data obtained from the bones and teeth of 27 fire victims using two DNA extraction techniques. Extracts were converted to double-stranded DNA libraries then enriched for whole-genome SNPs using unpublished biotinylated RNA baits and sequenced on an Illumina NextSeq 550 platform. Raw reads were processed using the EAGER (Efficient Ancient Genome Reconstruction) pipeline, and the SNPs filtered and called using FreeBayes and GATK (v. 3.8). Mixed-effects modeling of the data suggest that SNP variability and preservation is predominantly determined by skeletal element and burn category, and not by extraction type. Whole-genome SNP data suggest that selecting long bones, hand and foot bones, and teeth subjected to temperatures <350°C are the most likely sources for higher genomic DNA yields. Furthermore, we observed an inverse correlation between the number of captured SNPs and the extent to which samples were burned, as well as a significant decrease in the total number of SNPs measured for samples subjected to temperatures >350°C. Our data complement previous analyses of burned human remains that compare extraction methods for downstream forensic applications and support the idea of adopting a modified Dabney extraction technique when traditional forensic methods fail to produce DNA yields sufficient for genetic identification.

Association of ESR1 germline variants with TP53 somatic variants in breast tumors in a genome-wide study

Background

In breast tumors, somatic mutation frequencies in TP53 and PIK3CA vary by tumor subtype and ancestry. HER2 positive and triple negative breast cancers (TNBC) have a higher frequency of TP53 somatic mutations than other subtypes. PIK3CA mutations are more frequently observed in hormone receptor positive tumors. Emerging data suggest tumor mutation status is associated with germline variants and genetic ancestry. We aimed to identify germline variants that are associated with somatic TP53 or PIK3CA mutation status in breast tumors.

Methods

A genome-wide association study was conducted using breast cancer mutation status of TP53 and PIK3CA and functional mutation categories including TP53 gain of function (GOF) and loss of function mutations and PIK3CA activating/hotspot mutations. The discovery analysis consisted of 2850 European ancestry women from three datasets. Germline variants showing evidence of association with somatic mutations were selected for validation analyses based on predicted function, allele frequency, and proximity to known cancer genes or risk loci. Candidate variants were assessed for association with mutation status in a multi-ancestry validation study, a Malaysian study, and a study of African American/Black women with TNBC.

Results

The discovery Germline x Mutation (GxM) association study found five variants associated with one or more TP53 phenotypes with P values <1×10-6, 33 variants associated with one or more TP53 phenotypes with P values <1×10-5, and 44 variants associated with one or more PIK3CA phenotypes with P values <1×10-5. In the multi-ancestry and Malaysian validation studies, germline ESR1 locus variant, rs9383938, was associated with the presence of TP53 mutations overall (P values 6.8×10-5 and 9.8×10-8, respectively) and TP53 GOF mutations (P value 8.4×10-6). Multiple variants showed suggestive evidence of association with PIK3CA mutation status in the validation studies, but none were significant after correction for multiple comparisons.

Conclusions

We found evidence that germline variants were associated with TP53 and PIK3CA mutation status in breast cancers. Variants near the estrogen receptor alpha gene, ESR1, were significantly associated with overall TP53 mutations and GOF mutations. Larger multi-ancestry studies are needed to confirm these findings and determine if these variants contribute to ancestry-specific differences in mutation frequency.

Lack of racial and ethnic diversity in lung cancer cell lines contributes to lung cancer health disparities

Lung cancer is the leading cause of cancer death in the United States and worldwide, and a major source of cancer health disparities. Lung cancer cell lines provide key in vitro models for molecular studies of lung cancer development and progression, and for pre-clinical drug testing. To ensure health equity, it is imperative that cell lines representing different lung cancer histological types, carrying different cancer driver genes, and representing different genders, races, and ethnicities should be available. This is particularly relevant for cell lines from Black men, who experience the highest lung cancer mortality in the United States. Here, we undertook a review of the available lung cancer cell lines and their racial and ethnic origin. We noted a marked imbalance in the availability of cell lines from different races and ethnicities. Cell lines from Black patients were strongly underrepresented, and we identified no cell lines from Hispanic/Latin(x) (H/L), American Indian/American Native (AI/AN), or Native Hawaiian or other Pacific Islander (NHOPI) patients. The majority of cell lines were derived from White and Asian patients. Also missing are cell lines representing the cells-of-origin of the major lung cancer histological types, which can be used to model lung cancer development and to study the effects of environmental exposures on lung tissues. To our knowledge, the few available immortalized alveolar epithelial cell lines are all derived from White subjects, and the race and ethnicity of a handful of cell lines derived from bronchial epithelial cells are unknown. The lack of an appropriately diverse collection of lung cancer cell lines and lung cancer cell-of-origin lines severely limits racially and ethnically inclusive lung cancer research. It impedes the ability to develop inclusive models, screen comprehensively for effective compounds, pre-clinically test new drugs, and optimize precision medicine. It thereby hinders the development of therapies that can increase the survival of minority and underserved patients. The noted lack of cell lines from underrepresented groups should constitute a call to action to establish additional cell lines and ensure adequate representation of all population groups in this critical pre-clinical research resource.

Forensic biogeographical ancestry inference: recent insights and current trends

BACKGROUND

As a powerful complement to the paradigmatic DNA profiling strategy, biogeographical ancestry inference (BGAI) plays a significant part in human forensic investigation especially when a database hit or eyewitness testimony are not available. It indicates one's biogeographical profile based on known population-specific genetic variations, and thus is crucial for guiding authority investigations to find unknown individuals. Forensic biogeographical ancestry testing exploits much of the recent advances in the understanding of human genomic variation and improving of molecular biology.

OBJECTIVE

In this review, recent development of prospective ancestry informative markers (AIMs) and the statistical approaches of inferring biogeographic ancestry from AIMs are elucidated and discussed.

METHODS

We highlight the research progress of three potential AIMs (i.e., single nucleotide polymorphisms, microhaplotypes, and Y or mtDNA uniparental markers) and discuss the prospects and challenges of two methods that are commonly used in BGAI.

CONCLUSION

While BGAI for forensic purposes has been thriving in recent years, important challenges, such as ethics and responsibilities, data completeness, and ununified standards for evaluation, remain for the use of biogeographical ancestry information in human forensic investigations. To address these issues and fully realize the value of BGAI in forensic investigation, efforts should be made not only by labs/institutions around the world independently, but also by inter-lab/institution collaborations.

Neuron-binding antibody responses are associated with Black ethnicity in multiple sclerosis during natalizumab treatment

Multiple sclerosis is an inflammatory degenerative condition of the central nervous system that may result in debilitating disability. Several studies over the past twenty years suggest that multiple sclerosis manifests with a rapid, more disabling disease course among individuals identifying with Black or Latin American ethnicity relative to those of White ethnicity. However, very little is known about immunologic underpinnings that may contribute to this ethnicity-associated discordant clinical severity. Given the importance of B cells to multiple sclerosis pathophysiology, and prior work showing increased antibody levels in the cerebrospinal fluid of Black-identifying, compared to White-identifying multiple sclerosis patients, we conducted a cohort study to determine B cell subset dynamics according to both self-reported ethnicity and genetic ancestry over time. Further, we determined relationships between ethnicity, ancestry, and neuron-binding IgG levels. We found significant associations between Black ethnicity and elevated frequencies of class-switched B cell subsets, including memory B cells; double negative two B cells; and antibody-secreting cells. The frequencies of these subsets positively correlated with West African genetic ancestry. We also observed significant associations between Black ethnicity and increased IgG binding to neurons. Our data suggests significantly heightened T cell-dependent B cell responses exhibiting increased titres of neuron-binding antibodies among individuals with multiple sclerosis identifying with the Black African diaspora. Factors driving this immunobiology may promote the greater demyelination, central nervous system atrophy and disability more often experienced by Black-, and Latin American-identifying individuals with multiple sclerosis.

Exploration of the ancestral inference effectiveness of 126 SNPs and the genetic feature of Inner Mongolian Manchu population

In addition to the validated ancestry-informative single nucleotide polymorphisms (AI-SNPs) in classic panels, there are many new potential AI-SNPs yet to be explored. Moreover, the search for AI-SNPs with high discriminatory power for ancestry inference in inter- and intra-continental populations has become a realistic need. In this study, 126 novel AI-SNPs were selected to distinguish the African, European, Central/South Asian and East Asian populations, and a random forest model was introduced to assess the performance of the AI-SNP set. This set was further used in the genetic analysis of the Manchu people in Inner Mongolia, China, based on 79 reference populations from seven continental regions. Results showed that the 126 AI-SNPs were able to achieve the ancestry inference for African, East Asian, European, and Central/South Asian populations. Genetic analyses indicated that Inner Mongolian Manchus were genetically typical of East Asians and were more closely related to the northern Han Chinese and Japanese than to other Altaic-speaking populations. Overall, this study provided a selection of new promising loci for ancestry inference of major intercontinental populations and intracontinental subgroups, as well as genetic insights and valuable data for dissecting the genetic structure of the Inner Mongolian Manchu group.

Performance of McDonald 2017 multiple sclerosis diagnostic criteria and evaluation of genetic ancestry in patients with a first demyelinating event in Argentina

BACKGROUND

Information on performance of multiple sclerosis (MS) diagnostic criteria is scarce for populations from Latin America, Asia, or the Caribbean.

OBJECTIVE

To assess performance of revised 2017 McDonald criteria as well as evaluate genetic ancestry in a group of MS patients from Argentina experiencing a debut demyelinating event.

METHODS

Demographic and clinical characteristics, cerebrospinal fluid (CSF), and magnetic resonance imaging (MRI) findings and new T2 lesions were recorded at baseline and during relapses. Diagnostic accuracy in predicting conversion to clinically defined MS (CDMS) based on initial imaging applying revised 2017 criteria was evaluated and genetic ancestry-informative markers analyzed.

RESULTS

Of 201 patients experiencing their first demyelinating event (median follow-up 60 months), CDMS was confirmed in 67. We found 2017 diagnostic criteria were more sensitive (84% vs 67%) and less specific (14% vs 33%) than 2010 criteria, especially in a group of patients revised separately, presenting positive oligoclonal bands (88% vs 8%). Genetic testing performed in 128 cases showed 72% of patients were of European ancestry and 27% presented genetic admixture.

CONCLUSION

2017 McDonald criteria showed higher sensitivity and lower specificity compared with 2010 criteria, shortening both time-to-diagnosis and time-to-treatment implementation.

Risk prediction with imperfect survival outcome information from electronic health records

Readily available proxies for the time of disease onset such as the time of the first diagnostic code can lead to substantial risk prediction error if performing analyses based on poor proxies. Due to the lack of detailed documentation and labor intensiveness of manual annotation, it is often only feasible to ascertain for a small subset the current status of the disease by a follow-up time rather than the exact time. In this paper, we aim to develop risk prediction models for the onset time efficiently leveraging both a small number of labels on the current status and a large number of unlabeled observations on imperfect proxies. Under a semiparametric transformation model for onset and a highly flexible measurement error model for proxy onset time, we propose the semisupervised risk prediction method by combining information from proxies and limited labels efficiently. From an initially estimator solely based on the labeled subset, we perform a one-step correction with the full data augmenting against a mean zero rank correlation score derived from the proxies. We establish the consistency and asymptotic normality of the proposed semisupervised estimator and provide a resampling procedure for interval estimation. Simulation studies demonstrate that the proposed estimator performs well in a finite sample. We illustrate the proposed estimator by developing a genetic risk prediction model for obesity using data from Mass General Brigham Healthcare Biobank.

Vitamin D and genetic ancestry are associated with apoptosis rates in benign and malignant prostatic epithelium

PURPOSE

Vitamin D metabolites may be protective against prostate cancer (PCa). We conducted a cross-sectional analysis to evaluate associations between in vivo vitamin D status, genetic ancestry, and degree of apoptosis using prostatic epithelial terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) staining.

EXPERIMENTAL DESIGN

Benign and tumor epithelial punch biopsies of participants with clinically localized PCa underwent indirect TUNEL staining. Serum levels of 25 hydroxyvitamin D [25(OH)D] and 1,25 dihydroxyvitamin D were assessed immediately before radical prostatectomy; levels of prostatic 25(OH)D were obtained from the specimen once the prostate was extracted. Ancestry informative markers were used to estimate the percentage of genetic West African, Native American, and European ancestry.

RESULTS

One hundred twenty-one newly diagnosed men, age 40-79, were enrolled between 2013 and 2018. Serum 25(OH)D correlated positively with both tumor (ρ = 0.17, p = 0.03), and benign (ρ = 0.16, p = 0.04) prostatic epithelial TUNEL staining. Similarly, prostatic 25(OH)D correlated positively with both tumor (ρ = 0.31, p < 0.001) and benign (ρ = 0.20, p = 0.03) epithelial TUNEL staining. Only Native American ancestry was positively correlated with tumor (ρ = 0.22, p = 0.05) and benign (ρ = 0.27, p = 0.02) TUNEL staining. In multivariate regression models, increasing quartiles of prostatic 25(OH)D (β = 0.25, p = 0.04) and Native American ancestry (β = 0.327, p = 0.004) were independently associated with tumor TUNEL staining.

CONCLUSIONS

Physiologic serum and prostatic 25(OH)D levels and Native American ancestry are positively associated with the degree of apoptosis in tumor and benign prostatic epithelium in clinically localized PCa. Vitamin D may have secondary chemoprevention benefits in preventing PCa progression in localized disease.

Verification of insertion-deletion markers (InDels) and microsatellites (STRs) as subsidiary tools for inferring Slavic population ancestry

Genetic markers for the prediction of biogeographical ancestry have proved to be effective tools for law enforcement agencies for many years now. In this study, we attempted to assess the potential of insertion-deletion markers (InDel) and microsatellites (STRs) as subsidiary polymorphisms for inference of Slavic population ancestry. For that purpose, we genotyped Slavic-speaking populations samples from Belarus, the Czech Republic, Poland, Serbia, Ukraine and Russia in 46 InDels and 15 STRs by PCR and capillary electrophoresis and analyzed for between-population differentiation with the use of distance-based methods (FST, principal component analysis and multidimensional scaling). Additionally, we studied a sample from a Polish individual of well-documented genealogy whose biogeographic ancestry had previously been inferred by commercial genomic services using autosomal single nucleotide polymorphisms (SNPs), mitochondrial DNA and Y-SNP markers. For comparative purposes, we used genotype data collected in the “forInDel” browser and allele frequencies from previously published papers. The results obtained for InDels and STRs show that the Slavic populations constitute a genetically homogeneous group, with the exception of the Czechs differing clearly from the other tested populations. The analysis of the known Polish sample in the Snipper application proves the usefulness of the InDel markers on the continental level only. Conversely, microsatellites not only improve prediction, but are also informative if considered as an independent set of ancestry markers. Weryfikacja markerów insercyjno-delecyjnych (InDels) i mikrosatelitarnych (STR) jako narzędzi pomocniczych do wnioskowania o pochodzeniu populacji słowiańskiej Markery genetyczne do przewidywania pochodzenia biogeograficznego od wielu lat okazują się skutecznymi narzędziami dla organów ścigania. W tym badaniu podjęliśmy próbę oceny potencjału markerów insercyjno-delecyjnych (InDel) i mikrosatelitarnych (STR) jako pomocniczych polimorfizmów do wnioskowania o pochodzeniu populacji słowiańskiej. W tym celu genotypowaliśmy próbki populacji słowiańskojęzycznych z Białorusi, Czech, Polski, Serbii, Ukrainy i Rosji w w zakresie 46 markerów InDel oraz 15 loci STR za pomocą PCR i elektroforezy kapilarnej oraz analizowaliśmy pod kątem różnicowania między populacjami za pomocą metod bazujących na dystansach genetycznych (FST, analiza głównych składowych i skalowanie wielowymiarowe). Dodatkowo zbadaliśmy próbkę mężczyzny z populacji polskiej o dobrze udokumentowanej genealogii, którego pochodzenie biogeograficzne zostało wcześniej ustalone przez komercyjne usługi genomiczne przy użyciu autosomalnych polimorfizmów pojedynczych nukleotydów (SNP), mitochondrialnego DNA i markerów Y-SNP. Do celów porównawczych wykorzystaliśmy dane genotypowe zebrane w przeglądarce „forInDel” i częstości alleli z wcześniej opublikowanych artykułów. Uzyskane wyniki dla InDels i STR wskazują, że populacje słowiańskie stanowią grupę genetycznie jednorodną, z wyjątkiem Czechów wyraźnie różniących się od pozostałych badanych populacji. Analiza znanej polskiej próbki w aplikacji Snipper dowodzi przydatności markerów InDel jedynie na poziomie kontynentalnym. Z kolei, mikrosatelity nie tylko poprawiają wyniki predykcji, ale są informatywne jako niezależny zestaw markerów pochodzenia biogeograficznego.

A custom hybridisation enrichment forensic intelligence panel to infer biogeographic ancestry, hair and eye colour, and Y chromosome lineage

Massively parallel sequencing can provide genetic data for hundreds to thousands of loci in a single assay for various types of forensic testing. However, available commercial kits require an initial PCR amplification of short-to-medium sized targets which limits their application for highly degraded DNA. Development and optimisation of large PCR multiplexes also prevents creation of custom panels that target different suites of markers for identity, biogeographic ancestry, phenotype, and lineage markers (Y-chromosome and mtDNA). Hybridisation enrichment, an alternative approach for target enrichment prior to sequencing, uses biotinylated probes to bind to target DNA and has proven successful on degraded and ancient DNA. We developed a customisable hybridisation capture method, that uses individually mixed baits to allow tailored and targeted enrichment to specific forensic questions of interest. To allow collection of forensic intelligence data, we assembled and tested a custom panel of hybridisation baits to infer biogeographic ancestry, hair and eye colour, and paternal lineage (and sex) on modern male and female samples with a range of self-declared ancestries and hair/eye colour combinations. The panel correctly estimated biogeographic ancestry in 9/12 samples (75%) but detected European admixture in three individuals from regions with admixed demographic history. Hair and eye colour were predicted correctly in 83% and 92% of samples respectively, where intermediate eye colour and blond hair were problematic to predict. Analysis of Y-chromosome SNPs correctly assigned sex and paternal haplogroups, the latter complementing and supporting biogeographic ancestry predictions. Overall, we demonstrate the utility of this hybridisation enrichment approach to forensic intelligence testing using a combined suite of biogeographic ancestry, phenotype, and Y-chromosome SNPs for comprehensive biological profiling.

Comparative evaluation of the MAPlex, Precision ID Ancestry Panel, and VISAGE Basic Tool for biogeographical ancestry inference

Biogeographical ancestry (BGA) inference from ancestry-informative markers (AIMs) has strong potential to support forensic investigations. Over the past two decades, several forensic panels composed of AIMs have been developed to predict ancestry at a continental scale. These panels typically comprise fewer than 200 AIMs and have been designed and tested with a limited set of populations. How well these panels recover patterns of genetic diversity relative to larger sets of markers, and how accurately they infer ancestry of individuals and populations not included in their design remains poorly understood. The lack of comparative studies addressing these aspects makes the selection of appropriate panels for forensic laboratories difficult. In this study, the model-based genetic clustering tool STRUCTURE was used to compare three popular forensic BGA panels: MAPlex, Precision ID Ancestry Panel (PIDAP), and VISAGE Basic Tool (VISAGE BT) relative to a genome-wide reference set of 10k SNPs. The genotypes for all these markers were obtained for a comprehensive set of 3957 individuals from 228 worldwide human populations. Our results indicate that at the broad continental scale (K=6) typically examined in forensic studies, all forensic panels produced similar genetic structure patterns compared to the reference set (G'≈90%) and had high classification performance across all regions (average AUC-PR > 97%). However, at K= 7 and K= 8, the forensic panels displayed some region-specific clustering deviations from the reference set, particularly in Europe and the region of East and South-East Asia, which may be attributed to differences in the design of the respective panels. Overall, the panel with the most consistent performance in all regions was VISAGE BT with an average weighted AUC̅_W score of 96.26% across the three scales of geographical resolution investigated.

Ancestry resolution of South Brazilians by forensic 165 ancestry-informative SNPs panel

Forensic DNA phenotyping (FDP) includes biogeographic ancestry (BGA) inference and externally visible characteristics (EVCs) prediction directly from an evidential DNA sample as alternatives to provide valuable intelligence when conventional DNA profiling fails to achieve identification. In this context, the application of Massively Parallel Sequencing (MPS) methodologies, which enables simultaneous typing of multiple samples and hundreds of forensic markers, has been gradually implemented in forensic genetic casework. The Precision ID Ancestry Panel (Thermo Fisher Scientific, Waltham, USA) is a forensic multiplex assay consisting of 165 autosomal SNPs designed to provide biogeographic ancestry information. In this work, a sample of 250 individuals from Rio Grande do Sul (RS) State, southern Brazil, apportioned into four main population groups (African-, European-, Amerindian-, and Admixed-derived Gauchos), was evaluated with this panel, to assess the feasibility of this approach in a highly heterogeneous population. Forensic descriptive parameters estimated for each population group revealed that this panel has enough polymorphic and informative SNPs to be used as a supplementary instrument in forensic individual identification and kinship testing regardless of ethnicity. No statistically significant deviation from Hardy-Weinberg equilibrium was observed after Bonferroni correction. However, seven loci pairs displayed linkage disequilibrium in pairwise LD testing (p < 3.70 × 10^-6). Interpopulation comparisons by F_ST analysis, MDS plot, and STRUCTURE analysis among the four RS population groups apart and along with 89 reference worldwide populations demonstrated that Admixed- and African-derived Gauchos present the highest levels of admixture and population stratification, whereas European- and Amerindian-derived exhibit a more homogeneous genetic conformation.

Dietary inflammatory index and cardiovascular disease risk in Hispanic women from the Women's Health Initiative

BACKGROUND

To evaluate the association between the dietary inflammatory index (DII^®) and incident cardiovascular disease (CVD) in Hispanic women from the Women's Health Initiative (WHI), and to determine if body mass index (BMI) interacted with the DII scores.

METHODS

Secondary analysis of baseline dietary data and long-term CVD outcomes among 3,469 postmenopausal women who self-identified as Hispanic enrolled in WHI. DII scores were calculated from self-administered food frequency questionnaires. The CVD outcomes included coronary heart disease (CHD) and stroke. Stratified Cox regression models were used to assess the relationship between DII scores and CVD in women with and without obesity. Models were adjusted for age, lifestyle risk factors, known risk factors, and neighborhood socioeconomic status.

RESULTS

The incidence of CHD was 3.4 and 2.8% for stroke after a median follow-up of 12.9 years. None of the DIIs were associated with CVD risk in this sample of Hispanic women. BMI interacted with the DII (p < 0.20) and stratified models showed that the associations between the DII and CVD were only significant in women with overweight (p < 0.05). In this group, higher DII scores were associated with a higher risk of CHD (HR 1.27; 95% CI: 1.08, 1.51) and a higher risk of stroke (HR 1.32; 95% CI: 1.07, 1.64).

CONCLUSION

Among postmenopausal Hispanic women with overweight, greater adherence to pro-inflammatory diets was associated with higher risk of CVD. Additional research is needed to understand how to promote long-term heart-healthy dietary habits to reduce inflammation and prevent CVD in at-risk Hispanic women.

Ancestry-dependent genetic structure of the Xq28 risk haplotype in the Mexican population and its association with childhood-onset systemic lupus erythematosus

Objective

Here we aimed to investigate the association of the Xq28 risk haplotype (H1) with susceptibility to childhood-onset systemic lupus erythematosus (SLE), and to compare its frequency and genetic structure in the Mexican population with those in other continental populations.

Methods

We genotyped 15 single-nucleotide variants (SNVs) that form the H1 haplotype, using TaqMan real-time PCR. The association analysis [case-control and transmission disequilibrium test (TDT)] included 376 cases and 400 adult controls, all of whom were mestizos (MEZ). To identify risk alleles in Mexican Indigenous individuals, SNVs were imputed from whole-exome sequencing data of 1,074 individuals. The allelic frequencies determined in MEZ and Indigenous individuals were compared with those of the continental populations from the 1,000 Genomes database phase 3. Linkage disequilibrium (LD) analysis of risk alleles was performed on all populations. Interleukin-1 receptor associated kinase 1 (IRAK1) and methyl CpG binding protein 2 (MECP2) mRNA levels were determined using real-time PCR.

Results

Case-control analysis revealed genetic association with childhood-onset SLE for all 15 SNVs (OR = 1.49-1.75; p = 0.0095 to 1.81 × 10^-4) and for the Xq28 risk haplotype (OR = 1.97, p = 4 × 10^-6). Comparing with individuals of European ancestry (0.14-0.16), the frequencies of the risk alleles were significantly higher in the MEZ individuals (0.55-0.68) and even higher in Indigenous individuals (0.57-0.83). LD analysis indicated a differential haplotype structure within the Indigenous groups, which was inherited to the MEZ population as a result of genetic admixture. Individuals homozygous for the Xq28 risk haplotype exhibited decreased levels of both MECP2A and B transcripts.

Conclusion

We found that the H1 risk haplotype differs in its conformation in the Mexican population. This difference could be attributed to positive selection within the Indigenous population, with its inheritance now having an autoimmune health impact in both the Mexican Indigenous and MEZ populations.

Challenges in selecting admixture models and marker sets to infer genetic ancestry in a Brazilian admixed population

The inference of genetic ancestry plays an increasingly prominent role in clinical, population, and forensic genetics studies. Several genotyping strategies and analytical methodologies have been developed over the last few decades to assign individuals to specific biogeographic regions. However, despite these efforts, ancestry inference in populations with a recent history of admixture, such as those in Brazil, remains a challenge. In admixed populations, proportion and components of genetic ancestry vary on different levels: (i) between populations; (ii) between individuals of the same population, and (iii) throughout the individual's genome. The present study evaluated 1171 admixed Brazilian samples to compare the genetic ancestry inferred by tri-/tetra-hybrid admixture models and evaluated different marker sets from those with small numbers of ancestry informative markers panels (AIMs), to high-density SNPs (HDSNP) and whole-genome-sequence (WGS) data. Analyses revealed greater variation in the correlation coefficient of ancestry components within and between admixed populations, especially for minority ancestral components. We also observed positive correlation between the number of markers in the AIMs panel and HDSNP/WGS. Furthermore, the greater the number of markers, the more accurate the tri-/tetra-hybrid admixture models.

Genetic admixture and cardiovascular disease risk in postmenopausal Hispanic women

BACKGROUND

Hispanics are a heterogeneous population with differences in the prevalence of cardiovascular disease (CVD) and its related risk factors among ethnic sub-groups. This study evaluated the association of genetic admixture and CVD in self-identified Hispanic women from the Women's Health Initiative (WHI).

METHODS

Data came from the WHI Observational Study and the Clinical Trial Components conducted among postmenopausal women. The CVD outcomes included coronary heart disease (CHD) and stroke. The proportions of European (EUR), sub-Saharan African (AFR), and Amerindian (AMI) admixture were estimated using 92 ancestry-informative markers. Cox regression models were used to assess the relationship between genetic admixture and CVD adjusting for age, lifestyle risk factors, known risk factors, and neighborhood socioeconomic status.

RESULTS

Among 5195 participants EUR ancestry was associated with a lower CHD risk after adjusting for age (HR 0.41, p = 0.02), and in the fully adjusted model (HR 0.40, p = 0.03). AFR ancestry was associated with a higher CHD risk after adjusting for age (HR 2.91, p = 0.03), but it only showed a trend in in the fully adjusted model (HR 2.46, p = 0.10). AMI ancestry was not statistically significantly associated with CHD and none of the genetic admixture proportions were statistically significantly associated with stroke (p > 0.05).

CONCLUSION

EUR ancestry was associated with a lower risk of CHD in Hispanic women. This highlights the need to account for genetic admixture in future CVD studies to consider different heritage groups to understand the role that genetic, neighborhood socioeconomic status, and environmental factors contribute to CVD health disparities in Hispanic women.

Applications of massively parallel sequencing in forensic genetics

Massively parallel sequencing, also referred to as next-generation sequencing, has positively changed DNA analysis, allowing further advances in genetics. Its capability of dealing with low quantity/damaged samples makes it an interesting instrument for forensics. The main advantage of MPS is the possibility of analyzing simultaneously thousands of genetic markers, generating high-resolution data. Its detailed sequence information allowed the discovery of variations in core forensic short tandem repeat loci, as well as the identification of previous unknown polymorphisms. Furthermore, different types of markers can be sequenced in a single run, enabling the emergence of DIP-STRs, SNP-STR haplotypes, and microhaplotypes, which can be very useful in mixture deconvolution cases. In addition, the multiplex analysis of different single nucleotide polymorphisms can provide valuable information about identity, biogeographic ancestry, paternity, or phenotype. DNA methylation patterns, mitochondrial DNA, mRNA, and microRNA profiling can also be analyzed for different purposes, such as age inference, maternal lineage analysis, body-fluid identification, and monozygotic twin discrimination. MPS technology also empowers the study of metagenomics, which analyzes genetic material from a microbial community to obtain information about individual identification, post-mortem interval estimation, geolocation inference, and substrate analysis. This review aims to discuss the main applications of MPS in forensic genetics.

Influence of Genetic West African Ancestry on Metabolomics among Hypertensive Patients

Patients with higher genetic West African ancestry (GWAA) have hypertension (HTN) that is more difficult to treat and have higher rates of cardiovascular diseases (CVD) and differential responses to antihypertensive drugs than those with lower GWAA. The mechanisms underlying these disparities are poorly understood. Using data from 84 ancestry-informative markers in US participants from the Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) and PEAR-2 trials, the GWAA proportion was estimated. Using multivariable linear regression, the baseline levels of 886 metabolites were compared between PEAR participants with GWAA < 45% and those with GWAA ≥ 45% to identify differential metabolites and metabolic clusters. Metabolites with a false discovery rate (FDR) < 0.2 were used to create metabolic clusters, and a cluster analysis was conducted. Differential clusters were then tested for replication in PEAR-2 participants. We identified 353 differential metabolites (FDR < 0.2) between PEAR participants with GWAA < 45% (n = 383) and those with GWAA ≥ 45% (n = 250), which were used to create 24 metabolic clusters. Of those, 13 were significantly different between groups (Bonferroni p < 0.002). Four clusters, plasmalogen and lysoplasmalogen, sphingolipid metabolism and ceramide, cofactors and vitamins, and the urea cycle, were replicated in PEAR-2 (Bonferroni p < 0.0038) and have been previously linked to HTN and CVD. Our findings may give insights into the mechanisms underlying HTN racial disparities.

EASplex: A panel of 308 AISNPs for East Asian ancestry inference using next generation sequencing

Ancestry inference is useful in many scientific fields, such as forensic genetics, medical genetics, and molecular archaeology. Various ancestry inferring methods have been released for major continental populations. However, few reports refer to sub-populations within the East Asian population using hundreds of ancestry informative SNPs (AISNPs). In this study, we developed a 308-AISNP panel (EASplex NGS DNA panel) using multiplex PCR and next generation sequencing (NGS). This panel included 56 SNPs relevant for the continent-level ancestry inference and 252 Japanese-, Korean-, and/or Han Chinese-specific AISNPs to address the ancestry inference of global populations and regional populations among Japanese (JPT), Korean minority (CHK), and Han Chinese (CHH). A total of 87 CHK and 59 CHH samples were used to check the performance of the EASplex NGS DNA panel. By analyzing 146 profiles of samples with JPT and CHH data from Beijing and South China in 1000 genomes project, the following results were obtained: (1) the 146 tested samples were correctly assigned to the East Asian group; (2) the paired population assignment rate was 99.73% for JPT and CHH, 95% for JPT and CHK, and 90.11% for CHK and CHH; and (3) the whole population assignment was 92.14% for the JPT, CHK, and CHH data. Overall, the EASplex NGS DNA panel displayed informativeness for continental ancestry inference and regional ancestry inference among JPT, CHH, and CHK and has the potential for use in forensic and genetic population studies.

Unveiling forensically relevant biogeographic, phenotype and Y-chromosome SNP variation in Pakistani ethnic groups using a customized hybridisation enrichment forensic intelligence panel

Massively parallel sequencing following hybridisation enrichment provides new opportunities to obtain genetic data for various types of forensic testing and has proven successful on modern as well as degraded and ancient DNA. A customisable forensic intelligence panel that targeted 124 SNP markers (67 ancestry informative markers, 23 phenotype markers from the HIrisplex panel, and 35 Y-chromosome SNPs) was used to examine biogeographic ancestry, phenotype and sex and Y-lineage in samples from different ethnic populations of Pakistan including Pothwari, Gilgit, Baloach, Pathan, Kashmiri and Siraiki. Targeted sequencing and computational data analysis pipeline allowed filtering of variants across the targeted loci. Study samples showed an admixture between East Asian and European ancestry. Eye colour was predicted accurately based on the highest p-value giving overall prediction accuracy of 92.8%. Predictions were consistent with reported hair colour for all samples, using the combined highest p-value approach and step-wise model incorporating probability thresholds for light or dark shade. Y-SNPs were successfully recovered only from male samples which indicates the ability of this method to identify biological sex and allow inference of Y-haplogroup. Our results demonstrate practicality of using hybridisation enrichment and MPS to aid in human intelligence gathering and will open many insights into forensic research in South Asia.

Comparison of mixed model based approaches for correcting for population substructure with application to extreme phenotype sampling

Background

Mixed models are used to correct for confounding due to population stratification and hidden relatedness in genome-wide association studies. This class of models includes linear mixed models and generalized linear mixed models. Existing mixed model approaches to correct for population substructure have been previously investigated with both continuous and case-control response variables. However, they have not been investigated in the context of extreme phenotype sampling (EPS), where genetic covariates are only collected on samples having extreme response variable values. In this work, we compare the performance of existing binary trait mixed model approaches (GMMAT, LEAP and CARAT) on EPS data. Since linear mixed models are commonly used even with binary traits, we also evaluate the performance of a popular linear mixed model implementation (GEMMA).

Results

We used simulation studies to estimate the type I error rate and power of all approaches assuming a population with substructure. Our simulation results show that for a common candidate variant, both LEAP and GMMAT control the type I error rate while CARAT’s rate remains inflated. We applied all methods to a real dataset from a Québec, Canada, case-control study that is known to have population substructure. We observe similar type I error control with the analysis on the Québec dataset. For rare variants, the false positive rate remains inflated even after correction with mixed model approaches. For methods that control the type I error rate, the estimated power is comparable.

Conclusions

The methods compared in this study differ in their type I error control. Therefore, when data are from an EPS study, care should be taken to ensure that the models underlying the methodology are suitable to the sampling strategy and to the minor allele frequency of the candidate SNPs.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08297-y).

Genetic European Ancestry and Incident Diabetes in Black Individuals: Insights From the SPRINT Trial

BACKGROUND

Black individuals have high incident diabetes risk, despite having paradoxically lower triglyceride and higher HDL (high-density lipoprotein) cholesterol levels. The basis of this is poorly understood. We evaluated the participants of SPRINT (Systolic Blood Pressure Intervention Trial) to assess the association of estimated European genetic ancestry with the risk of incident diabetes in self-identified Black individuals.

METHODS

Self-identified non-Hispanic Black SPRINT participants free of diabetes at baseline were included. Black participants were stratified into tertiles (T1-T3) of European ancestry proportions estimated using 106 biallelic ancestry informative genetic markers. The multivariable-adjusted association of European ancestry proportion with indices of baseline metabolic syndrome (ie, fasting plasma glucose, triglycerides, HDL cholesterol, body mass index, and blood pressure) was assessed. Multivariable-adjusted Cox regression determined the risk of incident diabetes (fasting plasma glucose ≥126 mg/dL or self-reported diabetes treatment) across tertiles of European ancestry proportion.

RESULTS

Among 2466 Black SPRINT participants, a higher European ancestry proportion was independently associated with higher baseline triglyceride and lower HDL cholesterol levels (P<0.001 for both). European ancestry proportion was not associated with baseline fasting plasma glucose, body mass index, and blood pressure (P>0.05). Compared with the first tertile, those in the second (hazard ratio, 0.64 [95% CI, 0.45-0.90]) and third tertiles (hazard ratio, 0.61 [95% CI, 0.44-0.89]) of the European ancestry proportion had a lower risk of incident diabetes. A 5% point higher European ancestry was associated with a 29% lower risk of incident diabetes (hazard ratio, 0.71 [95% CI, 0.55-0.93]). There was no evidence of a differential association between the European ancestry proportion tertiles and incident diabetes between those randomized to intensive versus standard blood pressure treatment.

CONCLUSIONS

The higher risk of incident diabetes in Black individuals may have genetic determinants in addition to adverse social factors. Further research may help understand the interplay between biological and social determinants of cardiometabolic health in Black individuals. Registration: URL: https://www.clinicaltrials.gov; Unique identifier: NCT01206062.

Peptide ancestry informative markers in uterine neoplasms from women of European, African, and Asian ancestry

Characterization of ancestry-linked peptide variants in disease-relevant patient tissues represents a foundational step to connect patient ancestry with disease pathogenesis. Nonsynonymous single-nucleotide polymorphisms encoding missense substitutions within tryptic peptides exhibiting high allele frequencies in European, African, and East Asian populations, termed peptide ancestry informative markers (pAIMs), were prioritized from 1000 genomes. In silico analysis identified that as few as 20 pAIMs can determine ancestry proportions similarly to >260K SNPs (R² = 0.99). Multiplexed proteomic analysis of >100 human endometrial cancer cell lines and uterine leiomyoma tissues combined resulted in the quantitation of 62 pAIMs that correlate with patient race and genotype-confirmed ancestry. Candidates include a D451E substitution in GC vitamin D-binding protein previously associated with altered vitamin D levels in African and European populations. pAIMs will support generalized proteoancestry assessment as well as efforts investigating the impact of ancestry on the human proteome and how this relates to the pathogenesis of uterine neoplasms.

Assessment of the effectiveness of the EUROFORGEN NAME and Precision ID Ancestry panel markers for ancestry investigations

The EUROFORGEN NAME panel is a regional ancestry panel designed to differentiate individuals from the Middle East, North Africa, and Europe. The first version of the panel was developed for the MassARRAY system and included 111 SNPs. Here, a custom AmpliSeq EUROFORGEN NAME panel with 102 of the original 111 loci was used to sequence 1098 individuals from 14 populations from Europe, the Middle East, North Africa, North-East Africa, and South-Central Asia. These samples were also sequenced with a global ancestry panel, the Precision ID Ancestry Panel. The GenoGeographer software was used to assign the AIM profiles to reference populations and calculate the weight of the evidence as likelihood ratios. The combination of the EUROFORGEN NAME and Precision ID Ancestry panels led to fewer ambiguous assignments, especially for individuals from the Middle East and South-Central Asia. The likelihood ratios showed that North African individuals could be separated from European and Middle Eastern individuals using the Precision ID Ancestry Panel. The separation improved with the addition of the EUROFORGEN NAME panel. The analyses also showed that the separation of Middle Eastern populations from European and South-Central Asian populations was challenging even when both panels were applied.

Genotypes of informative loci from 1000 Genomes data allude evolution and mixing of human populations

Principal Component Analysis (PCA) projects high-dimensional genotype data into a few components that discern populations. Ancestry Informative Markers (AIMs) are a small subset of SNPs capable of distinguishing populations. We integrate these two approaches by proposing an algorithm to identify necessary informative loci whose removal from the data deteriorates the PCA structure. Unlike classical AIMs, necessary informative loci densely cover the genome, hence can illuminate the evolution and mixing history of populations. We conduct a comprehensive analysis to the genotype data of the 1000 Genomes Project using necessary informative loci. Projections along the top seven principal components demarcate populations at distinct geographic levels. Millions of necessary informative loci along each PC are identified. Population identities along each PC are approximately determined by weighted sums of minor (or major) alleles over the informative loci. Variations of allele frequencies are aligned with the history and direction of population evolution. The population distribution of projections along the top three PCs is recapitulated by a simple demographic model based on several waves of founder population separation and mixing. Informative loci possess locational concentration in the genome and functional enrichment. Genes at two hot spots encompassing dense PC 7 informative loci exhibit differential expressions among European populations. The mosaic of local ancestry in the genome of a mixed descendant from multiple populations can be inferred from partial PCA projections of informative loci. Finally, informative loci derived from the 1000 Genomes data well predict the projections of an independent genotype data of South Asians. These results demonstrate the utility and relevance of informative loci to investigate human evolution.

Genome-wide sequencing-based identification of methylation quantitative trait loci and their role in schizophrenia risk

DNA methylation (DNAm) is an epigenetic regulator of gene expression and a hallmark of gene-environment interaction. Using whole-genome bisulfite sequencing, we have surveyed DNAm in 344 samples of human postmortem brain tissue from neurotypical subjects and individuals with schizophrenia. We identify genetic influence on local methylation levels throughout the genome, both at CpG sites and CpH sites, with 86% of SNPs and 55% of CpGs being part of methylation quantitative trait loci (meQTLs). These associations can further be clustered into regions that are differentially methylated by a given SNP, highlighting the genes and regions with which these loci are epigenetically associated. These findings can be used to better characterize schizophrenia GWAS-identified variants as epigenetic risk variants. Regions differentially methylated by schizophrenia risk-SNPs explain much of the heritability associated with risk loci, despite covering only a fraction of the genomic space. We provide a comprehensive, single base resolution view of association between genetic variation and genomic methylation, and implicate schizophrenia GWAS-associated variants as influencing the epigenetic plasticity of the brain.

The authors provide a comprehensive, single base resolution view of association between genetic variation and DNA methylation in human brain. They also show that heritability attributed to schizophrenia GWAS-associated variants reflects the epigenetic plasticity of the brain.

Ancestry analysis using autosomal SNPs in northern South America, reveals interpretation differences between an AIM panel and an identification panel

Current human populations are studied to elucidate their ancestry composition and to obtain reference values for an array of genetic markers for forensic practice. This study compared the 79 ancestry informative markers (AIMs) panel with the SNPforID 52plex set used in forensic identification, using samples belonging to Continental Caribbean populations from Colombia with a high percentage of locals self-determined as Native American descendants. The results show a bias in the individual estimation made with the identification markers, which disregards the Native American ancestry component and overestimates the African ancestry component. Also, the analysis made with the Bayesian Classification Algorithm shows better likelihoods for individual assignment with AIMs than with SNPforID 52plex.

Diet Modulates the Effects of Genetic Variants on the Vitamin D Metabolic Pathway and Bone Mineral Density in Mexican Postmenopausal Women

BACKGROUND

Macro- and micronutrients, such as proteins, vitamin D, and calcium (Ca), are important dietary factors that can modify bone mineral density (BMD). Genetic factors can interact with diet, affecting an individual's predisposition to osteoporosis.

OBJECTIVES

This study aimed to evaluate the associations between macro- and micronutrient intakes and BMD in Mexican postmenopausal women, and their interactions with genetic polymorphisms involved in the vitamin D metabolic pathway.

METHODS

We analyzed data from 317 postmenopausal women from the Health Workers Cohort Study, a longitudinal cohort studied in Cuernavaca, Mexico. Postmenopausal women participated in 2 data collection waves (2004-2006 and 2010-2011), with a mean time of 6.4 years. Dietary intake was assessed with a semi-quantitative FFQ. BMD (femoral neck, hip, and lumbar spine) was measured by DXA. Hybrid mixed-effects regression models were used to assess the associations of dietary macro- and micronutrients on BMD, after adjusting for confounding factors and for diet and single nucleotide polymorphism interactions.

RESULTS

At baseline, the median age was 57 years (IQR, 50-64). Mean femoral neck, hip, and lumbar spine BMDs decreased over time. We observed statistically significant longitudinal associations for diet (Ca, vitamin D, magnesium, phosphorus, and protein intake) and BMD. Increases of vitamin D, Ca, and protein intakes by 1 SD were associated with mean increases in the femoral neck BMD (0.083 SD, 0.064 SD, and 0.130 SD, respectively). Multiple significant interactions were identified between several loci (CYP2R1, CYP24A1, CYP27B1, VDR, and DHCR7/NADSYN1) and diet for BMDs (femoral neck, hip, and lumbar spine), mainly for protein intake.

CONCLUSIONS

Our data support associations of vitamin D, Ca, protein, phosphorous, and magnesium consumption with BMD in Mexican postmenopausal women and suggest possible gene-diet interactions. These results could facilitate future personalized nutrition recommendations to help prevent low BMD.

Challenges in providing residual risks in carrier testing

The probability an individual is a carrier for a recessive disorder despite a negative carrier test, referred to as residual risk, has been part of carrier screening for over 2 decades. Residual risks are calculated by subtracting the frequency of carriers of pathogenic variants detected by the test from the carrier frequency in a population, estimated from the incidence of the disease. Estimates of the incidence (and therefore carrier frequency) of many recessive disorders differ among different population groups and are inaccurate or unavailable for many genes on large carrier screening panels for most of the world's populations. The pathogenic variants detected by the test and their frequencies also vary across groups and over time as variants are newly discovered or reclassified, which requires today's residual carrier risks to be continually updated. Even when a residual carrier risk is derived using accurate data obtained in a particular group, it may not apply to many individuals in that group because of misattributed ancestry or unsuspected admixture. Missing or inaccurate data, the challenge of determining meaningful ancestry‐specific risks and applying them appropriately, and a lack of evidence they impact management, suggest that patients be counseled that although carrier screening may miss a small fraction of carriers, residual risks with contemporary carrier screening are well below the risk posed by invasive prenatal diagnosis, even if one member of the couple is a carrier, and that efforts to provide precise residual carrier risks are unnecessary.

What's already known about this topic? What does this study add?

There has been no published discussion of the methods and uncertainties involved in the calculation of residual risk that are discussed here

There has been much discussion of using ancestry in genetic testing but this review highlights the serious problems that arise in calculating and assigning ancestry‐specific residual carrier risks at specific disease loci

The review questions what has not been questioned before: Is there clinical utility to providing what are mostly imprecise residual carrier risks

Analysis of Skin Pigmentation and Genetic Ancestry in Three Subpopulations from Pakistan: Punjabi, Pashtun, and Baloch

Skin pigmentation is one of the most prominent and variable phenotypes in humans. We compared the alleles of 163 SNPs and indels from the Human Pigmentation (HuPi) AmpliSeq™ Custom panel, and biogeographic ancestry with the quantitative skin pigmentation levels on the upper arm, lower arm, and forehead of 299 Pakistani individuals from three subpopulations: Baloch, Pashtun, and Punjabi. The biogeographic ancestry of each individual was estimated using the Precision ID Ancestry Panel. All individuals were mainly of mixed South-Central Asian and European ancestry. However, the Baloch individuals also had an average proportion of Sub-Saharan African ancestry of approximately 10%, whereas it was <1% in the Punjabi and Pashtun individuals. The pairwise genetic distances between the Pashtun, Punjabi, and Baloch subpopulations based on the ancestry markers were statistically significantly different. Individuals from the Pashtun subpopulation had statistically significantly lower skin pigmentation than individuals from the Punjabi and Baloch subpopulations (p < 0.05). The proportions of European and Sub-Saharan African ancestry and five SNPs (rs1042602, rs10831496, rs1426654, rs16891982, and rs12913832) were statistically significantly associated with skin pigmentation at either the upper arm, lower arm or forehead in the Pakistani population after correction for multiple testing (p < 10⁻³). A model based on four of these SNPs (rs1426654, rs1042602, rs16891982, and rs12913832) explained 33% of the upper arm skin pigmentation. The four SNPs and the proportions of European and Sub-Saharan African ancestry explained 37% of the upper arm skin pigmentation. Our results indicate that the four likely causative SNPs, rs1426654, rs1042602, rs16891982, and rs12913832 located in SLC24A5, TYR, SLC45A2, and HERC2, respectively, are essential for skin color variation in the admixed Pakistani subpopulations.

Evaluation of the association of heterozygous germline variants in NTHL1 with breast cancer predisposition: an international multi-center study of 47,180 subjects

Bi-allelic loss-of-function (LoF) variants in the base excision repair (BER) gene NTHL1 cause a high-risk hereditary multi-tumor syndrome that includes breast cancer, but the contribution of heterozygous variants to hereditary breast cancer is unknown. An analysis of 4985 women with breast cancer, enriched for familial features, and 4786 cancer-free women revealed significant enrichment for NTHL1 LoF variants. Immunohistochemistry confirmed reduced NTHL1 expression in tumors from heterozygous carriers but the NTHL1 bi-allelic loss characteristic mutational signature (SBS 30) was not present. The analysis was extended to 27,421 breast cancer cases and 19,759 controls from 10 international studies revealing 138 cases and 93 controls with a heterozygous LoF variant (OR 1.06, 95% CI: 0.82-1.39) and 316 cases and 179 controls with a missense variant (OR 1.31, 95% CI: 1.09-1.57). Missense variants selected for deleterious features by a number of in silico bioinformatic prediction tools or located within the endonuclease III functional domain showed a stronger association with breast cancer. Somatic sequencing of breast cancers from carriers indicated that the risk associated with NTHL1 appears to operate through haploinsufficiency, consistent with other described low-penetrance breast cancer genes. Data from this very large international multicenter study suggests that heterozygous pathogenic germline coding variants in NTHL1 may be associated with low- to moderate- increased risk of breast cancer.

Comprehensive cell type decomposition of circulating cell-free DNA with CelFiE

Circulating cell-free DNA (cfDNA) in the bloodstream originates from dying cells and is a promising noninvasive biomarker for cell death. Here, we propose an algorithm, CelFiE, to accurately estimate the relative abundances of cell types and tissues contributing to cfDNA from epigenetic cfDNA sequencing. In contrast to previous work, CelFiE accommodates low coverage data, does not require CpG site curation, and estimates contributions from multiple unknown cell types that are not available in external reference data. In simulations, CelFiE accurately estimates known and unknown cell type proportions from low coverage and noisy cfDNA mixtures, including from cell types composing less than 1% of the total mixture. When used in two clinically-relevant situations, CelFiE correctly estimates a large placenta component in pregnant women, and an elevated skeletal muscle component in amyotrophic lateral sclerosis (ALS) patients, consistent with the occurrence of muscle wasting typical in these patients. Together, these results show how CelFiE could be a useful tool for biomarker discovery and monitoring the progression of degenerative disease.

Tissue damage and turnover lead to the release of DNA in the blood and can be used to monitor changes in tissue state. Here, the authors developed a tool to accurately estimate the proportion of cell types contributing to cell-free DNA in the blood, with an application to pregnant women and ALS patients.

Cancer microenvironment and genomics: evolution in process

Cancer heterogeneity is a result of genetic mutations within the cancer cells. Their proliferation is not only driven by autocrine functions but also under the influence of cancer microenvironment, which consists of normal stromal cells such as infiltrating immune cells, cancer-associated fibroblasts, endothelial cells, pericytes, vascular and lymphatic channels. The relationship between cancer cells and cancer microenvironment is a critical one and we are just on the verge to understand it on a molecular level. Cancer microenvironment may serve as a selective force to modulate cancer cells to allow them to evolve into more aggressive clones with ability to invade the lymphatic or vascular channels to spread to regional lymph nodes and distant sites. It is important to understand these steps of cancer evolution within the cancer microenvironment towards invasion so that therapeutic strategies can be developed to control or stop these processes.

Investigation of triple-negative breast cancer risk alleles in an International African-enriched cohort

Large-scale efforts to identify breast cancer (BC) risk alleles have historically taken place among women of European ancestry. Recently, there are new efforts to verify if these alleles increase risk in African American (AA) women as well. We investigated the effect of previously reported AA breast cancer and triple-negative breast cancer (TNBC) risk alleles in our African-enriched International Center for the Study of Breast Cancer Subtypes (ICSBCS) cohort. Using case-control, case-series and race-nested approaches, we report that the Duffy-null allele (rs2814778) is associated with TNBC risk (OR = 3.814, p = 0.001), specifically among AA individuals, after adjusting for self-indicated race and west African ancestry (OR = 3.368, p = 0.007). We have also validated the protective effect of the minor allele of the ANKLE1 missense variant rs2363956 among AA for TNBC (OR = 0.420, p = 0.005). Our results suggest that an ancestry-specific Duffy-null allele and differential prevalence of a polymorphic gene variant of ANKLE1 may play a role in TNBC breast cancer outcomes. These findings present opportunities for therapeutic potential and future studies to address race-specific differences in TNBC risk and disease outcome.

Developmental validations of a self-developed 39 AIM-InDel panel and its forensic efficiency evaluations in the Shaanxi Han population

Most of insertion/deletion polymorphisms are diallelic molecular markers characterized as small amplicon sizes, high inter-population diversities, and low mutation rates, which make them the promising genetic markers in biogeographic ancestor inference field. The developmental validations of a 39 ancestry informative marker-insertion/deletion (AIM-InDel) panel and the genetic polymorphic investigations of this panel were performed in the Shaanxi Han population of China. The developmental validation included the optimizations of PCR-related indicators, repeatability, reproducibility, precision, accuracy, sensitivity, species specificity, stability of the panel, and the abilities in analyzing degraded, casework, and mixture samples, and the present results demonstrated that this 39 AIM-InDel panel was robust, sensitive, and accurate. For the population diversity analyses, the combined discrimination power value of 38 AIM-InDel loci except for rs36038238 locus was 0.999999999931257, indicating that this novel panel was highly polymorphic, biogeographic informative, and could be also used in individual identifications in the Shaanxi Han population.

Ethical considerations for DNA testing as a proxy for nationality

As nations strengthen borders and restrict refugee admissions, national security officials are screening for fraudulent nationality claims. One tool to investigate nationality claims is DNA testing, either for claimed relationships or for ancestral origins. At the same time, the plight of global statelessness leaves millions without documentation of their nationality, and DNA testing might be the only recourse to provide evidence of heritage or relationships. DNA testing has been used sparsely to date to determine ancestral origin as a proxy for nationality but could increase as border controls tighten. Given the historic lessons in eugenics and the potential for misuse of personal genetic information, it is essential to consider the ethical parameters in order to guide the implementation of genetic data for such purposes. Here, we break down examples of the use of DNA testing for nationality, and the risks and benefits of genetic testing for this purpose. Important ethical considerations discussed include (1) empowerment of stateless individuals with evidence for citizenship proceedings; (2) imprecise correlation between genetic heritage and nationality; (3) effective protection of state interests; and (4) practicalities of DNA testing.

Cytokine gene polymorphisms are associated with response to blinatumomab in B-cell acute lymphoblastic leukemia

Blinatumomab is a bispecific T cell-engaging antibody approved for treatment of relapsed/refractory (r/r) ALL, with 40%-50% complete response (CR)/CR with incomplete count recovery (CRi). Cytokine release syndrome (CRS) as a major adverse effect after blinatumomab therapy. Here, we evaluated the possible association between single-nucleotide polymorphisms (SNPs) in cytokine genes, disease response, and CRS in r/r ALL patients who received blinatumomab between 2012 and 2017 at our center (n = 66), using patients' archived DNA samples. With a median duration of 9.5 months (range: 1-37), 37 patients (56.1%) achieved CR/CRi, 54 (81.8%) experienced CRS (G1: n = 35, G2: n = 14, G3: n = 5), and 9 (13.6%) developed neurotoxicity. By multivariable analysis, after adjusting for high disease burden, one SNP on IL2 (rs2069762), odds ratio (OR) = 0.074 (95% CI: NE-0.43, P = .01) and one SNP on IL17A (rs4711998), OR = 0.28 (95% CI: 0.078-0.92, P = .034) were independently associated with CR/CRi. None of the analyzed SNPs were associated with CRS. To our knowledge, this is the first study demonstrating a possible association between treatment response to blinatumomab and SNPs. Our hypothesis-generated data suggest a potential role for IL-17 and IL-2 in blinatumomab response and justify a larger confirmatory study, which may lead to personalized blinatumomab immunotherapy for B-ALL.

An alternate workflow for preparing Precision ID Ancestry and Identity Panel libraries for Illumina sequencing

Single nucleotide polymorphisms (SNPs) are well-established for forensic applications. Although they are not compatible with existing criminal databases, they offer some advantages over short tandem repeat (STR) markers including smaller amplicons, no stutter artifacts, and biogeographic ancestry and phenotype predictions. The Precision ID NGS System, a commercial workflow by Thermo Fisher Scientific, offers a streamlined solution for genotyping forensically relevant SNPs using next-generation sequencing. The Precision ID Ancestry and Identity Panels combined target 289 SNPs, and their sensitivity, reproducibility, and accuracy have been evaluated by the forensic community. The aim of this study was to develop an alternative workflow to genotype these SNP panels using Illumina chemistry. Commercial genomic DNAs (gDNAs) (n, 3) were amplified using three uracil-tolerant polymerase master mixes. Resulting amplicons were prepared into libraries using the KAPA Hyper Prep Kit (KAPA Biosystems) and sequenced via Illumina's MiniSeq. Reads were analyzed using a published analysis pipeline to compile final genotypes with read depth information. Phusion U Multiplex PCR Master Mix (Thermo Fisher Scientific) statistically outperformed the other master mixes tested (P <0.0001), with respect to the number of SNPs genotyped. To ensure a workflow using Phusion U would be compatible across diverse samples, we optimized PCR cycle number using the same commercial gDNAs (n, 3), reference buccal swabs (n, 3), and environmental (household dust) samples (n, 6). Using the developed workflow, 93.9% of all SNPs were successfully genotyped across sample types. Implementation of the developed workflow should be straightforward for forensic laboratories and suitable for processing reference and casework samples.

Differential Associations of SLCO Transporters with Prostate Cancer Aggressiveness between African Americans and European Americans

BACKGROUND

Androgen receptor signaling is crucial to prostate cancer aggressiveness. Members of the solute carrier family of the organic anion transporting peptides (SLCO) are potential regulators of androgen availability in prostate tissue. It remains unknown whether genetic variations in SLCOs contribute to the differences in prostate cancer aggressiveness in African Americans (AA) and European Americans (EA).

METHODS

SNPs in 11 SLCO members were selected, with addition of 139 potentially functional SNPs and 128 ancestry informative markers. A total of 1,045 SNPs were genotyped and analyzed in 993 AAs and 1,057 EAs from the North Carolina-Louisiana Prostate Cancer Project. Expression and cellular localization of SLCOs were examined using qRT-PCR, IHC, and in situ RNA hybridization in independent sets of prostate cancer cases.

RESULTS

Significant associations with prostate cancer characteristics were found for SNPs in SLCO2A1 and SLCO5A1. The associations differed by race (P _interaction < 0.05). SNPs in SLCO2A1 were associated with reduced tumor aggressiveness and low Gleason score in AAs; whereas, SNPs in SLCO5A1 were associated with high clinical stage in EAs. In prostate tissue, SLCO2A1 and SLCO5A1 were the most expressed SLCOs at the mRNA level and were expressed predominantly in prostate endothelial and epithelial cells at the protein level, respectively.

CONCLUSIONS

SLCO2A1 and SLCO5A1 play important but different roles in prostate cancer aggressiveness in AAs versus EAs.

IMPACT

The finding calls for consideration of racial differences in biomarker studies of prostate cancer and for investigations on functions of SLCO2A1 and SLCO5A1 in prostate cancer.