An ancestry informative marker panel design for individual ancestry estimation of Hispanic population using whole exome sequencing data

Exome-wide association analysis identifies novel risk loci for alcohol-associated hepatitis

BACKGROUND AND AIMS

Alcohol-associated hepatitis (AH) is a clinically severe, acute disease that afflicts only a fraction of patients with alcohol use disorder. Genomic studies of alcohol-associated cirrhosis (AC) have identified several genes of large effect, but the genetic and environmental factors that lead to AH and AC, and their degree of genetic overlap, remain largely unknown. This study aims to identify genes and genetic variations that contribute to the development of AH.

APPROACH AND RESULTS

Exome-sequencing of patients with AH (N=784) and heavy drinking controls (N=951) identified an exome-wide significant association for AH at patalin-like phospholipase domain containing 3, as previously observed for AC in genome-wide association study, although with a much lower effect size. Single nucleotide polymorphisms (SNPs) of large effect size at inducible T cell costimulatory ligand ( ICOSLG ) (Chr 21) and TOX4/RAB2B (Chr 14) were also exome-wide significant. ICOSLG encodes a co-stimulatory signal for T-cell proliferation and cytokine secretion and induces B-cell proliferation and differentiation. TOX high mobility group box family member 4 ( TOX4 ) was previously implicated in diabetes and immune system function. Other genes previously implicated in AC did not strongly contribute to AH, and the only prominently implicated (but not exome-wide significant) gene overlapping with alcohol use disorder was alcohol dehydrogenase 1B ( ADH1B ). Polygenic signals for AH were observed in both common and rare variant analysis and identified genes with roles associated with inflammation.

CONCLUSIONS

This study has identified 2 new genes of high effect size with a previously unknown contribution to alcohol-associated liver disease and highlights both the overlap in etiology between liver diseases and the unique origins of AH.

Utility of Candidate Genes From an Algorithm Designed to Predict Genetic Risk for Opioid Use Disorder

Importance

Recently, the US Food and Drug Administration gave premarketing approval to an algorithm based on its purported ability to identify individuals at genetic risk for opioid use disorder (OUD). However, the clinical utility of the candidate genetic variants included in the algorithm has not been independently demonstrated.

Objective

To assess the utility of 15 genetic variants from an algorithm intended to predict OUD risk.

Design, Setting, and Participants

This case-control study examined the association of 15 candidate genetic variants with risk of OUD using electronic health record data from December 20, 1992, to September 30, 2022. Electronic health record data, including pharmacy records, were accrued from participants in the Million Veteran Program across the US with opioid exposure (n = 452 664). Cases with OUD were identified using International Classification of Diseases, Ninth Revision, or International Classification of Diseases, Tenth Revision, diagnostic codes, and controls were individuals with no OUD diagnosis.

Exposures

Number of risk alleles present across 15 candidate genetic variants.

Main Outcome and Measures

Performance of 15 genetic variants for identifying OUD risk assessed via logistic regression and machine learning models.

Results

A total of 452 664 individuals with opioid exposure (including 33 669 with OUD) had a mean (SD) age of 61.15 (13.37) years, and 90.46% were male; the sample was ancestrally diverse (with individuals of genetically inferred European, African, and admixed American ancestries). Using Nagelkerke R2, collectively, the 15 candidate genes accounted for 0.40% of variation in OUD risk. In comparison, age and sex alone accounted for 3.27% of the variation. The ensemble machine learning. The ensemble machine learning model using the 15 variants as predictive factors correctly classified 52.83% (95% CI, 52.07%-53.59%) of individuals in an independent testing sample.

Conclusions and Relevance

Results of this study suggest that the candidate genetic variants included in the approved algorithm do not meet reasonable standards of efficacy in identifying OUD risk. Given the algorithm's limited predictive accuracy, its use in clinical care would lead to high rates of both false-positive and false-negative findings. More clinically useful models are needed to identify individuals at risk of developing OUD.

Integrative multi-omics characterization of hepatocellular carcinoma in Hispanic patients

BACKGROUND

The incidence and mortality rates of hepatocellular carcinoma among Hispanic individuals in the United States are much higher than in non-Hispanic White people. We conducted multi-omics analyses to elucidate molecular alterations in hepatocellular carcinoma among Hispanic patients.

METHODS

Paired tumor and adjacent nontumor samples were collected from 31 Hispanic hepatocellular carcinomas in South Texas for genomic, transcriptomic, proteomic, and metabolomic profiling. Serum lipids were profiled in 40 Hispanic and non-Hispanic patients with or without clinically diagnosed hepatocellular carcinoma.

RESULTS

Exome sequencing revealed high mutation frequencies of AXIN2 and CTNNB1 in South Texas Hispanic hepatocellular carcinoma patients, suggesting a predominant activation of the Wnt/β-catenin pathway. TERT promoter mutations were also statistically significantly more frequent in the Hispanic cohort (Fisher exact test, P < .05). Cell cycles and liver function were positively and negatively enriched, respectively, with gene set enrichment analysis. Gene sets representing specific liver metabolic pathways were associated with dysregulation of corresponding metabolites. Negative enrichment of liver adipogenesis and lipid metabolism corroborated with a significant reduction in most lipids in serum samples of hepatocellular carcinoma patients (paired t test, P < .0001). Two hepatocellular carcinoma subtypes from our Hispanic cohort were identified and validated with the Cancer Genome Atlas liver cancer cohort. Patients with better overall survival showed higher activity of immune and angiogenesis signatures and lower activity of liver function-related gene signatures. They also had higher levels of immune checkpoint and immune exhaustion markers.

CONCLUSIONS

Our study revealed specific molecular features of Hispanic hepatocellular carcinoma and potential biomarkers for therapeutic management. It provides a unique resource for studying Hispanic hepatocellular carcinoma.

Avanzando Caminos (Leading Pathways): design and procedures of the Hispanic/Latino Cancer Survivorship Study

Avanzando Caminos (Leading Pathways): The Hispanic/Latino Cancer Survivorship Cohort Study aims to examine the influence of sociocultural, medical, stress-related, psychosocial, lifestyle, behavioral, and biological factors on symptom burden, health-related quality of life, and clinical outcomes among Hispanics/Latinos who have been previously treated for cancer. Avanzando Caminos is a prospective, cohort-based study of 3000 Hispanics/Latinos who completed primary cancer treatment within the past 5 years that is representative of the general Hispanic/Latino population in the United States. Participants will complete self-report measures at baseline (time [T] 1), 6 months (T2), 1 year (T3), 2 years (T4), 3 years (T5), 4 years (T6), and 5 years (T7). Blood samples drawn for assessment of leukocyte gene expression, cardiometabolic markers, and genetic admixture will be collected at baseline (T1), 1 year (T3), 3 years (T5), and 5 years (T7). Medical and cancer characteristics and clinical outcomes will be extracted from the electronic medical record and/or state cancer registry at each time point. Data analysis will include general latent variable modeling and latent growth modeling. Avanzando Caminos will fill critical gaps in knowledge in order to guide future secondary and tertiary prevention efforts to mitigate cancer disparities and optimize health-related quality of life among Hispanic/Latino cancer survivors.

Candidate Genes from an FDA-Approved Algorithm Fail to Predict Opioid Use Disorder Risk in Over 450,000 Veterans

Importance

Recently, the Food and Drug Administration gave pre-marketing approval to algorithm based on its purported ability to identify genetic risk for opioid use disorder. However, the clinical utility of the candidate genes comprising the algorithm has not been independently demonstrated.

Objective

To assess the utility of 15 variants in candidate genes from an algorithm intended to predict opioid use disorder risk.

Design

This case-control study examined the association of 15 candidate genetic variants with risk of opioid use disorder using available electronic health record data from December 20, 1992 to September 30, 2022.

Setting

Electronic health record data, including pharmacy records, from Million Veteran Program participants across the United States.

Participants

Participants were opioid-exposed individuals enrolled in the Million Veteran Program (n = 452,664). Opioid use disorder cases were identified using International Classification of Disease diagnostic codes, and controls were individuals with no opioid use disorder diagnosis.

Exposures

Number of risk alleles present across 15 candidate genetic variants.

Main Outcome and Measures

Predictive performance of 15 genetic variants for opioid use disorder risk assessed via logistic regression and machine learning models.

Results

Opioid exposed individuals (n=33,669 cases) were on average 61.15 (SD = 13.37) years old, 90.46% male, and had varied genetic similarity to global reference panels. Collectively, the 15 candidate genetic variants accounted for 0.4% of variation in opioid use disorder risk. The accuracy of the ensemble machine learning model using the 15 genes as predictors was 52.8% (95% CI = 52.1 - 53.6%) in an independent testing sample.

Conclusions and Relevance

Candidate genes that comprise the approved algorithm do not meet reasonable standards of efficacy in predicting opioid use disorder risk. Given the algorithm's limited predictive accuracy, its use in clinical care would lead to high rates of false positive and negative findings. More clinically useful models are needed to identify individuals at risk of developing opioid use disorder.

Integrative multi-omics characterization of hepatocellular carcinoma in Hispanic patients

Background

The incidence and mortality rates of hepatocellular carcinoma (HCC) among Hispanics in the United States are much higher than those of non-Hispanic whites. We conducted comprehensive multi-omics analyses to understand molecular alterations in HCC among Hispanic patients.

Methods

Paired tumor and adjacent non-tumor samples were collected from 31 Hispanic HCC in South Texas (STX-Hispanic) for genomic, transcriptomic, proteomic, and metabolomic profiling. Additionally, serum lipids were profiled in 40 Hispanic and non-Hispanic patients with or without clinically diagnosed HCC.

Results

Exome sequencing revealed high mutation frequencies of AXIN2 and CTNNB1 in STX Hispanic HCCs, suggesting a predominant activation of the Wnt/β-catenin pathway. The TERT promoter mutation frequency was also remarkably high in the Hispanic cohort. Cell cycles and liver functions were identified as positively- and negatively-enriched, respectively, with gene set enrichment analysis. Gene sets representing specific liver metabolic pathways were associated with dysregulation of corresponding metabolites. Negative enrichment of liver adipogenesis and lipid metabolism corroborated with a significant reduction in most lipids in the serum samples of HCC patients. Two HCC subtypes from our Hispanic cohort were identified and validated with the TCGA liver cancer cohort. The subtype with better overall survival showed higher activity of immune and angiogenesis signatures, and lower activity of liver function-related gene signatures. It also had higher levels of immune checkpoint and immune exhaustion markers.

Conclusions

Our study revealed some specific molecular features of Hispanic HCC and potential biomarkers for therapeutic management of HCC and provides a unique resource for studying Hispanic HCC.

Genomic profiling of subcutaneous patient-derived xenografts reveals immune constraints on tumor evolution in childhood solid cancer

Subcutaneous patient-derived xenografts (PDXs) are an important tool for childhood cancer research. Here, we describe a resource of 68 early passage PDXs established from 65 pediatric solid tumor patients. Through genomic profiling of paired PDXs and patient tumors (PTs), we observe low mutational similarity in about 30% of the PT/PDX pairs. Clonal analysis in these pairs show an aggressive PT minor subclone seeds the major clone in the PDX. We show evidence that this subclone is more immunogenic and is likely suppressed by immune responses in the PT. These results suggest interplay between intratumoral heterogeneity and antitumor immunity may underlie the genetic disparity between PTs and PDXs. We further show that PDXs generally recapitulate PTs in copy number and transcriptomic profiles. Finally, we report a gene fusion LRPAP1-PDGFRA. In summary, we report a childhood cancer PDX resource and our study highlights the role of immune constraints on tumor evolution.

Comprehensive characterization of patient-derived xenograft models of pediatric leukemia

Patient-derived xenografts (PDX) remain valuable models for understanding the biology and for developing novel therapeutics. To expand current PDX models of childhood leukemia, we have developed new PDX models from Hispanic patients, a subgroup with a poorer overall outcome. Of 117 primary leukemia samples obtained, successful engraftment and serial passage in mice were achieved in 82 samples (70%). Hispanic patient samples engrafted at a rate (51/73, 70%) that was similar to non-Hispanic patient samples (31/45, 70%). With a new algorithm to remove mouse contamination in multi-omics datasets including methylation data, we found PDX models faithfully reflected somatic mutations, copy-number alterations, RNA expression, gene fusions, whole-genome methylation patterns, and immunophenotypes found in primary tumor (PT) samples in the first 50 reported here. This cohort of characterized PDX childhood leukemias represents a valuable resource in that germline DNA sequencing has allowed the unambiguous determination of somatic mutations in both PT and PDX.

Association of rare, recurrent nonsynonymous variants in the germline of prostate cancer patients of African ancestry

BACKGROUND

Although men of African ancestry (AA) have the highest mortality rate from prostate cancer (PCa), relatively little is known about the germline variants that are associated with PCa risk in AA men. The goal of this study is to systematically evaluate rare, recurrent nonsynonymous variants across the exome for their association with PCa in AA men.

METHODS

Whole exome sequencing (WES) of germline DNA in two AA PCa patient cohorts of Johns Hopkins Hospital (N = 960) and Wayne State University (N = 747) was performed. All nonsynonymous variants present in both case cohorts, with a carrier rate between 0.5% and 1%, were identified. Their carrier rates were compared with rates from 8128 African/African American (AFR) control subjects from The Genome Aggregation Database (gnomAD) using Fisher's exact test. Significant variants, defined as false discovery rate (FDR) adjusted p-value ≤ 0.05, were further evaluated in AA PCa cases (N = 132) and controls (N = 1184) from the UK Biobank (UKB).

RESULTS

Two variants reached a pre-specified statistical significance level. The first was p.R14Q in GPRC5C (found in 0.47% of PCa cases and 0.01% of population controls); odds ratio (OR) for PCa was 37.46 (95% confidence interval CI 4.68-299.72), pexact  = 7.01E-06, FDR-adjusted p-value = 0.05. The second was p.R511Q in IGF1R (found in 0.53% of PCa cases and 0.01% of population controls); OR for PCa was 21.54 (95%CI 4.65-99.76), pexact  = 5.51E-06, FDR-adjusted p-value = 0.05. The mean percentage of African ancestry was similar between variant carriers and noncarriers of each variant, p > 0.05. In the UKB AA men, GPRC5C R14Q was 0.76% and 0.08% in cases and controls, respectively, OR for PCa was 9.00 (95%CI 0.56-145.23), pexact  = 0.19. However, IGF1R R511Q was not found in cases or controls.

CONCLUSIONS

This WES study identified two rare, recurrent nonsynonymous PCa risk-associated variants in AA. Confirmation in additional large populations of AA PCa cases and controls is required.

Frequencies of variants in genes associated with dyslipidemias identified in Costa Rican genomes

Dyslipidemias are risk factors in diseases of significant importance to public health, such as atherosclerosis, a condition that contributes to the development of cardiovascular disease. Unhealthy lifestyles, the pre-existence of diseases, and the accumulation of genetic variants in some loci contribute to the development of dyslipidemia. The genetic causality behind these diseases has been studied primarily on populations with extensive European ancestry. Only some studies have explored this topic in Costa Rica, and none have focused on identifying variants that can alter blood lipid levels and quantifying their frequency. To fill this gap, this study focused on identifying variants in 69 genes involved in lipid metabolism using genomes from two studies in Costa Rica. We contrasted the allelic frequencies with those of groups reported in the 1000 Genomes Project and gnomAD and identified potential variants that could influence the development of dyslipidemias. In total, we detected 2,600 variants in the evaluated regions. However, after various filtering steps, we obtained 18 variants that have the potential to alter the function of 16 genes, nine variants have pharmacogenomic or protective implications, eight have high risk in Variant Effect Predictor, and eight were found in other Latin American genetic studies of lipid alterations and the development of dyslipidemia. Some of these variants have been linked to changes in blood lipid levels in other global studies and databases. In future studies, we propose to confirm at least 40 variants of interest from 23 genes in a larger cohort from Costa Rica and Latin American populations to determine their relevance regarding the genetic burden for dyslipidemia. Additionally, more complex studies should arise that include diverse clinical, environmental, and genetic data from patients and controls and functional validation of the variants.

A Genetic Polymorphism in the WDR72 Gene is Associated With Calcium Nephrolithiasis in the Chinese Han Population

A previous genome-wide association study (GWAS) reported several novel loci for nephrolithiasis in British and Japanese population, some of which were predicted to influence CaSR signaling. In this study, we aimed to evaluate the association of these loci with calcium nephrolithiasis in Chinese Han population. We performed a case-control association analysis involving 691 patients with calcium nephrolithiasis and 1008 control subjects. We were able to genotype a total of 17 single-nucleotide polymorphisms (SNPs), which were previously reported to be significantly associated with nephrolithiasis in GWAS. rs578595 at WDR72 was significantly associated with calcium nephrolithiasis in Chinese Han population (p < 0.001, OR = 0.617). Moreover, rs12654812 at SLC34A1 (p = 0.0427, OR = 1.170), rs12539707 at HIBADH (p = 0.0179, OR = 0.734), rs1037271 at DGKH (p = 0.0096, OR = 0.828) and rs12626330 at CLDN14 (p = 0.0080, OR = 1.213) indicated suggestive associations with calcium nephrolithiasis. Our results elucidated the significance of genetic variation at WDR72, DGKH, CLDN14, SLC34A1, and HIBADH in Chinese patients with nephrolithiasis. Since polymorphisms of WDR72, DGKH, and CLDN14 are predicted to influence in CaSR signaling, our results emphasized the role of abnormal calcium homeostasis in calcium nephrolithiasis.

Disparities in Cancer Genetic Testing and Variants of Uncertain Significance in the Hispanic Population of South Texas

PURPOSE

Racial and ethnic disparities have included a lack of access to both genetic testing and research, resulting in poor understanding of the genomic architecture in under-represented populations. The South Texas population is primarily of Hispanic background and has been largely devoid of genetic services. We extended access to this underserved population and uncovered genetic variants previously not observed, emphasizing the need to continually improve both genomic databases and clarification of variant significance to provide meaningful patient counseling.

METHODS

This study consisted of a retrospective cohort review of patients seen through a cancer genetics education and service program across 24 counties in South Texas. In total, 1,595 individuals were identified as appropriate for cancer genetic counseling and 1,377 completed genetic testing.

RESULTS

Eighty percent of those receiving genetic counseling self-identified as Hispanic, 16% as non-Hispanic White (NHW), 3% as African American, and 1% as other race/ethnicity. Of reported variants, 18.8% were pathogenic and 13.7% were reported as a variant of uncertain significance (VUS). VUS was reported in 17.2% of the Hispanic individuals compared with 9% NHW (P = .005).

CONCLUSION

Individuals of Hispanic ethnicity were significantly more likely to harbor a VUS compared with NHW. The extended reach into our regional communities revealed a gap in the ability to accurately interpret genomic variation with implications for advising patients on screening, prevention, and management strategies. A higher percentage of VUS also emphasizes the challenge of continued follow-up amid existing barriers that led to disparities in access. As understanding of the variants develops, hopefully gaps in knowledge of the genomic landscape will be lessened with increased clarity to provide accurate cancer risk assessment and recommendations for implementing prevention initiatives.

Genetic Ancestry Inference and Its Application for the Genetic Mapping of Human Diseases

Admixed populations arise when two or more ancestral populations interbreed. As a result of this admixture, the genome of admixed populations is defined by tracts of variable size inherited from these parental groups and has particular genetic features that provide valuable information about their demographic history. Diverse methods can be used to derive the ancestry apportionment of admixed individuals, and such inferences can be leveraged for the discovery of genetic loci associated with diseases and traits, therefore having important biomedical implications. In this review article, we summarize the most common methods of global and local genetic ancestry estimation and discuss the use of admixture mapping studies in human diseases.

Bloom syndrome in a Mexican American family with rhabdomyosarcoma: evidence of a Mexican founder mutation

Bloom syndrome is a rare autosomal recessive disorder with less than 300 cases reported in the literature. Bloom syndrome is characterized by chromosome instability, physical stigmata, growth deficiency, immunodeficiency, and a predisposition to cancer, most commonly leukemias, although solid tumors are reported as well. Bloom syndrome occurs in multiple ethnic groups with a higher incidence in persons of Ashkenazi Jewish origin. Few patients of Hispanic ethnicity have been reported. We report here a Mexican American family with a BLM pathogenic variant, c.2506_2507delAG, previously reported in a single patient from Mexico. In this family of four siblings, three have phenotypic features of Bloom syndrome, and BLM gene mutation was homozygous in these affected individuals. Our proband developed a rhabdomyosarcoma. Analysis of surrounding markers in the germline DNA revealed a common haplotype, suggesting a previously unrecognized founder mutation in the Hispanic population of Mexican origin.

Genetic markers for treatment-related pancreatitis in a cohort of Hispanic children with acute lymphoblastic leukemia

PURPOSE

Treatment-related pancreatitis (TRP) is a serious complication occurring in children with acute lymphoblastic leukemia (ALL). Those affected are at high risk for severe organ toxicity and treatment delays that can impact outcomes. TRP is associated with asparaginase, a standard therapeutic agent in childhood ALL. Native American ancestry, older age, high-risk leukemia, and increased use of asparaginase are linked to pancreatitis risk. However, dedicated genetic studies evaluating pancreatitis in childhood ALL include few Hispanics. Thus, the genetic basis for higher risk of pancreatitis among Hispanic children with ALL remains unknown.

METHODS

Cases of children with ALL treated in from 1994 through 2013 were reviewed and identified 14, all Hispanic, who developed pancreatitis related to asparaginase therapy. Forty-six controls consisting of Hispanic children treated on the same regimens without pancreatitis were selected for comparison. Total DNA isolated from whole blood was used for targeted DNA sequencing of 23 selected genes, including genes associated with pancreatitis without ALL and genes involved in asparagine metabolism.

RESULTS

Non-synonymous polymorphisms and frameshift deletions were detected in 15 genes. Most children with TRP had variants in ABAT, ASNS, and CFTR. Notably, children with TRP harbored many more CFTR variants (71.4%) compared with controls (39.1%). Among these, V470M (rs213950) was most frequent (OR 4.27, p = 0.025).

CONCLUSIONS

This is the first study of genetic factors in treatment-related pancreatitis in Hispanic children with ALL. Identifying correlative variants in ethnically vulnerable populations may improve screening to identify which patients with ALL are at greatest risk for pancreatitis.

The International Conference on Intelligent Biology and Medicine 2019 (ICIBM 2019): conference summary and innovations in genomics

The goal of this editorial is to summarize the 2019 International Conference on Intelligent Biology and Medicine (ICIBM 2019) conference that took place on June 9–11, 2019 in The Ohio State University, Columbus, OH, and to provide an introductory summary of the seven articles presented in this supplement issue. ICIBM 2019 hosted four keynote speakers, four eminent scholar speakers, five tutorials and workshops, twelve concurrent sessions and a poster session, totaling 23 posters, spanning state-of-the-art developments in bioinformatics, genomics, next-generation sequencing (NGS) analysis, scientific databases, cancer and medical genomics, and computational drug discovery. A total of 105 original manuscripts were submitted to ICIBM 2019, and after careful review, seven were selected for this supplement issue. These articles cover methods and applications for functional annotations of miRNA targeting, clonal evolution of bacterial cells, gene co-expression networks that describe a given phenotype, functional binding site analysis of RNA-binding proteins, normalization of genome architecture mapping data, sample predictions based on multiple NGS data types, and prediction of an individual’s genetic admixture given exonic single nucleotide polymorphisms data.