The Puerto Rican Infant Metagenomic and Epidemiologic Study of Respiratory Outcomes (PRIMERO): Design and Baseline Characteristics for a Birth Cohort Study of Early-life Viral Respiratory Illnesses and Airway Dysfunction in Puerto Rican Children.

Epidemiologic studies demonstrate an association between early-life respiratory illnesses (RIs) and the development of childhood asthma. However, it remains uncertain whether these children are predisposed to both conditions or if early-life RIs induce alterations in airway function, immune responses, or other human biology that contribute to the development of asthma. Puerto Rican children experience a disproportionate burden of early-life RIs and asthma, making them an important population for investigating this complex interplay. PRIMERO, the Puerto Rican Infant Metagenomics and Epidemiologic Study of Respiratory Outcomes , recruited pregnant women and their newborns to investigate how the airways develop in early life among infants exposed to different viral RIs, and will thus provide a critical understanding of childhood asthma development. As the first asthma birth cohort in Puerto Rico, PRIMERO will prospectively follow 2,100 term healthy infants. Collected samples include post-term maternal peripheral blood, infant cord blood, the child's peripheral blood at the year two visit, and the child's nasal airway epithelium, collected using minimally invasive nasal swabs, at birth, during RIs over the first two years of life, and at annual healthy visits until age five. Herein, we describe the study's design, population, recruitment strategy, study visits and procedures, and primary outcomes.

A common polymorphism in the Intelectin-1 gene influences mucus plugging in severe asthma.

By incompletely understood mechanisms, type 2 (T2) inflammation present in the airways of severe asthmatics drives the formation of pathologic mucus which leads to airway mucus plugging. Here we investigate the molecular role and clinical significance of intelectin-1 (ITLN-1) in the development of pathologic airway mucus in asthma. Through analyses of human airway epithelial cells we find that ITLN1 gene expression is highly induced by interleukin-13 (IL-13) in a subset of metaplastic MUC5AC+ mucus secretory cells, and that ITLN-1 protein is a secreted component of IL-13-induced mucus. Additionally, we find ITLN-1 protein binds the C-terminus of the MUC5AC mucin and that its deletion in airway epithelial cells partially reverses IL-13-induced mucostasis. Through analysis of nasal airway epithelial brushings, we find that ITLN1 is highly expressed in T2-high asthmatics, when compared to T2-low children. Furthermore, we demonstrate that both ITLN-1 gene expression and protein levels are significantly reduced by a common genetic variant that is associated with protection from the formation of mucus plugs in T2-high asthma. This work identifies an important biomarker and targetable pathways for the treatment of mucus obstruction in asthma.

Polygenic risk score for ulcerative colitis predicts immune checkpoint inhibitor-mediated colitis.

Immune checkpoint inhibitor-mediated colitis (IMC) is a common adverse event of treatment with immune checkpoint inhibitors (ICI). We hypothesize that genetic susceptibility to Crohn's disease (CD) and ulcerative colitis (UC) predisposes to IMC. In this study, we first develop a polygenic risk scores for CD (PRSCD) and UC (PRSUC) in cancer-free individuals and then test these PRSs on IMC in a cohort of 1316 patients with ICI-treated non-small cell lung cancer and perform a replication in 873 ICI-treated pan-cancer patients. In a meta-analysis, the PRSUC predicts all-grade IMC (ORmeta=1.35 per standard deviation [SD], 95% CI = 1.12-1.64, P = 2×10-03) and severe IMC (ORmeta=1.49 per SD, 95% CI = 1.18-1.88, P = 9×10-04). PRSCD is not associated with IMC. Furthermore, PRSUC predicts severe IMC among patients treated with combination ICIs (ORmeta=2.20 per SD, 95% CI = 1.07-4.53, P = 0.03). Overall, PRSUC can identify patients receiving ICI at risk of developing IMC and may be useful to monitor patients and improve patient outcomes.

Novel insights into the whole-blood DNA methylome of asthma in ethnically diverse children and youth.

BACKGROUND: The epigenetic mechanisms of asthma remain largely understudied in African Americans and Hispanics/Latinos, two populations disproportionately affected by asthma. We aimed to identify markers, regions and processes with differential patterns of DNA methylation (DNAm) in whole blood by asthma status in ethnically diverse children and youth, and to assess their functional consequences. METHODS: DNAm levels were profiled with the Infinium MethylationEPIC or HumanMethylation450 BeadChip arrays among 1226 African Americans or Hispanics/Latinos and assessed for differential methylation per asthma status at the CpG and region (differentially methylated region (DMR)) level. Novel associations were validated in blood and/or nasal epithelium from ethnically diverse children and youth. The functional and biological implications of the markers identified were investigated by combining epigenomics with transcriptomics from study participants. RESULTS: 128 CpGs and 196 DMRs were differentially methylated after multiple testing corrections, including 92.3% and 92.8% novel associations, respectively. 41 CpGs were replicated in other Hispanics/Latinos, prioritising cg17647904 (NCOR2) and cg16412914 (AXIN1) as asthma DNAm markers. Significant DNAm markers were enriched in previous associations for asthma, fractional exhaled nitric oxide, bacterial infections, immune regulation or eosinophilia. Functional annotation highlighted epigenetically regulated gene networks involved in corticosteroid response, host defence and immune regulation. Several implicated genes are targets for approved or experimental drugs, including TNNC1 and NDUFA12. Many differentially methylated loci previously associated with asthma were validated in our study. CONCLUSIONS: We report novel whole-blood DNAm markers for asthma underlying key processes of the disease pathophysiology and confirm the transferability of previous asthma DNAm associations to ethnically diverse populations.

Polygenic risk score for ulcerative colitis predicts immune checkpoint inhibitor-mediated colitis.

Immune checkpoint inhibitors (ICIs) are a remarkable advancement in cancer therapeutics; however, a substantial proportion of patients develop severe immune-related adverse events (irAEs). Understanding and predicting irAEs is a key to advancing precision immuno-oncology. Immune checkpoint inhibitor-mediated colitis (IMC) is a significant complication from ICI and can have life-threatening consequences. Based on clinical presentation, IMC mimics inflammatory bowel disease, however the link is poorly understood. We hypothesized that genetic susceptibility to Crohn's disease (CD) and ulcerative colitis (UC) may predispose to IMC. We developed and validated polygenic risk scores for CD (PRSCD) and UC (PRSUC) in cancer-free individuals and assessed the role of each of these PRSs on IMC in a cohort of 1,316 patients with non-small cell lung cancer who received ICIs. Prevalence of all-grade IMC in our cohort was 4% (55 cases), and for severe IMC, 2.5% (32 cases). The PRSUC predicted the development of all-grade IMC (HR=1.34 per standard deviation [SD], 95% CI=1.02-1.76, P=0.04) and severe IMC (HR=1.62 per SD, 95% CI=1.12-2.35, P=0.01). PRSCD was not associated with IMC or severe IMC. The association between PRSUC and IMC (all-grade and severe) was consistent in an independent pan-cancer cohort of patients treated with ICIs. Furthermore, PRSUC predicted severe IMC among patients treated with combination ICIs (OR = 2.20 per SD, 95% CI = 1.07-4.53, P=0.03). This is the first study to demonstrate the potential clinical utility of a PRS for ulcerative colitis in identifying patients receiving ICI at high risk of developing IMC, where risk reduction and close monitoring strategies could help improve overall patient outcomes.

Development and testing of a polygenic risk score for breast cancer aggressiveness.

Aggressive breast cancers portend a poor prognosis, but current polygenic risk scores (PRSs) for breast cancer do not reliably predict aggressive cancers. Aggressiveness can be effectively recapitulated using tumor gene expression profiling. Thus, we sought to develop a PRS for the risk of recurrence score weighted on proliferation (ROR-P), an established prognostic signature. Using 2363 breast cancers with tumor gene expression data and single nucleotide polymorphism (SNP) genotypes, we examined the associations between ROR-P and known breast cancer susceptibility SNPs using linear regression models. We constructed PRSs based on varying p-value thresholds and selected the optimal PRS based on model r2 in 5-fold cross-validation. We then used Cox proportional hazards regression to test the ROR-P PRS's association with breast cancer-specific survival in two independent cohorts totaling 10,196 breast cancers and 785 events. In meta-analysis of these cohorts, higher ROR-P PRS was associated with worse survival, HR per SD = 1.13 (95% CI 1.06-1.21, p = 4.0 × 10-4). The ROR-P PRS had a similar magnitude of effect on survival as a comparator PRS for estrogen receptor (ER)-negative versus positive cancer risk (PRSER-/ER+). Furthermore, its effect was minimally attenuated when adjusted for PRSER-/ER+, suggesting that the ROR-P PRS provides additional prognostic information beyond ER status. In summary, we used integrated analysis of germline SNP and tumor gene expression data to construct a PRS associated with aggressive tumor biology and worse survival. These findings could potentially enhance risk stratification for breast cancer screening and prevention.

Gene expression in African Americans, Puerto Ricans and Mexican Americans reveals ancestry-specific patterns of genetic architecture.

We explored ancestry-related differences in the genetic architecture of whole-blood gene expression using whole-genome and RNA sequencing data from 2,733 African Americans, Puerto Ricans and Mexican Americans. We found that heritability of gene expression significantly increased with greater proportions of African genetic ancestry and decreased with higher proportions of Indigenous American ancestry, reflecting the relationship between heterozygosity and genetic variance. Among heritable protein-coding genes, the prevalence of ancestry-specific expression quantitative trait loci (anc-eQTLs) was 30% in African ancestry and 8% for Indigenous American ancestry segments. Most anc-eQTLs (89%) were driven by population differences in allele frequency. Transcriptome-wide association analyses of multi-ancestry summary statistics for 28 traits identified 79% more gene-trait associations using transcriptome prediction models trained in our admixed population than models trained using data from the Genotype-Tissue Expression project. Our study highlights the importance of measuring gene expression across large and ancestrally diverse populations for enabling new discoveries and reducing disparities.

Profiling the Somatic Mutational Landscape of Breast Tumors from Hispanic/Latina Women Reveals Conserved and Unique Characteristics.

Somatic mutational profiling is increasingly being used to identify potential targets for breast cancer. However, limited tumor-sequencing data from Hispanic/Latinas (H/L) are available to guide treatment. To address this gap, we performed whole-exome sequencing (WES) and RNA sequencing on 146 tumors and WES of matched germline DNA from 140 H/L women in California. Tumor intrinsic subtype, somatic mutations, copy-number alterations, and expression profiles of the tumors were characterized and compared with data from tumors of non-Hispanic White (White) women in The Cancer Genome Atlas (TCGA). Eight genes were significantly mutated in the H/L tumors including PIK3CA, TP53, GATA3, MAP3K1, CDH1, CBFB, PTEN, and RUNX1; the prevalence of mutations in these genes was similar to that observed in White women in TCGA. Four previously reported Catalogue of Somatic Mutations in Cancer (COSMIC) mutation signatures (1, 2, 3, 13) were found in the H/L dataset, along with signature 16 that has not been previously reported in other breast cancer datasets. Recurrent amplifications were observed in breast cancer drivers including MYC, FGFR1, CCND1, and ERBB2, as well as a recurrent amplification in 17q11.2 associated with high KIAA0100 gene expression that has been implicated in breast cancer aggressiveness. In conclusion, this study identified a higher prevalence of COSMIC signature 16 and a recurrent copy-number amplification affecting expression of KIAA0100 in breast tumors from H/L compared with White women. These results highlight the necessity of studying underrepresented populations. SIGNIFICANCE: Comprehensive characterization of genomic and transcriptomic alterations in breast tumors from Hispanic/Latina patients reveals distinct genetic alterations and signatures, demonstrating the importance of inclusive studies to ensure equitable care for patients. See related commentary by Schmit et al., p. 2443.

Whole exome sequencing and replication for breast cancer among Hispanic/Latino women identifies FANCM as a susceptibility gene for estrogen-receptor-negative breast cancer.

Introduction: Breast cancer (BC) is one of the most common cancers globally. Genetic testing can facilitate screening and risk-reducing recommendations, and inform use of targeted treatments. However, genes included in testing panels are from studies of European-ancestry participants. We sequenced Hispanic/Latina (H/L) women to identify BC susceptibility genes. Methods: We conducted a pooled BC case-control analysis in H/L women from the San Francisco Bay area, Los Angeles County, and Mexico (4,178 cases and 4,344 controls). Whole exome sequencing was conducted on 1,043 cases and 1,188 controls and a targeted 857-gene panel on the remaining samples. Using ancestry-adjusted SKAT-O analyses, we tested the association of loss of function (LoF) variants with overall, estrogen receptor (ER)-positive, and ER-negative BC risk. We calculated odds ratios (OR) for BC using ancestry-adjusted logistic regression models. We also tested the association of single variants with BC risk. Results: We saw a strong association of LoF variants in FANCM with ER-negative BC (p=4.1×10-7, OR [CI]: 6.7 [2.9-15.6]) and a nominal association with overall BC risk. Among known susceptibility genes, BRCA1 (p=2.3×10-10, OR [CI]: 24.9 [6.1-102.5]), BRCA2 (p=8.4×10-10, OR [CI]: 7.0 [3.5-14.0]), and PALB2 (p=1.8×10-8, OR [CI]: 6.5 [3.2-13.1]) were strongly associated with BC. There were nominally significant associations with CHEK2, RAD51D, and TP53. Conclusion: In H/L women, LoF variants in FANCM were strongly associated with ER-negative breast cancer risk. It previously was proposed as a possible susceptibility gene for ER-negative BC, but is not routinely tested in clinical practice. Our results demonstrate that FANCM should be added to BC gene panels.

Multi-omic approach associates blood methylome with bronchodilator drug response in pediatric asthma.

BACKGROUND: Albuterol is the drug most widely used as asthma treatment among African Americans despite having a lower bronchodilator drug response (BDR) than other populations. Although BDR is affected by gene and environmental factors, the influence of DNA methylation is unknown. OBJECTIVE: This study aimed to identify epigenetic markers in whole blood associated with BDR, study their functional consequences by multi-omic integration, and assess their clinical applicability in admixed populations with a high asthma burden. METHODS: We studied 414 children and young adults (8-21 years old) with asthma in a discovery and replication design. We performed an epigenome-wide association study on 221 African Americans and replicated the results on 193 Latinos. Functional consequences were assessed by integrating epigenomics with genomics, transcriptomics, and environmental exposure data. Machine learning was used to develop a panel of epigenetic markers to classify treatment response. RESULTS: We identified 5 differentially methylated regions and 2 CpGs genome-wide significantly associated with BDR in African Americans located in FGL2 (cg08241295, P = 6.8 × 10-9) and DNASE2 (cg15341340, P = 7.8 × 10-8), which were regulated by genetic variation and/or associated with gene expression of nearby genes (false discovery rate < 0.05). The CpG cg15341340 was replicated in Latinos (P = 3.5 × 10-3). Moreover, a panel of 70 CpGs showed good classification for those with response and nonresponse to albuterol therapy in African American and Latino children (area under the receiver operating characteristic curve for training, 0.99; for validation, 0.70-0.71). The DNA methylation model showed similar discrimination as clinical predictors (P > .05). CONCLUSIONS: We report novel associations of epigenetic markers with BDR in pediatric asthma and demonstrate for the first time the applicability of pharmacoepigenetics in precision medicine of respiratory diseases.

Admixture mapping of severe asthma exacerbations in Hispanic/Latino children and youth.

BACKGROUND: In the USA, genetically admixed populations have the highest asthma prevalence and severe asthma exacerbations rates. This could be explained not only by environmental factors but also by genetic variants that exert ethnic-specific effects. However, no admixture mapping has been performed for severe asthma exacerbations. OBJECTIVE: We sought to identify genetic variants associated with severe asthma exacerbations in Hispanic/Latino subgroups by means of admixture mapping analyses and fine mapping, and to assess their transferability to other populations and potential functional roles. METHODS: We performed an admixture mapping in 1124 Puerto Rican and 625 Mexican American children with asthma. Fine-mapping of the significant peaks was performed via allelic testing of common and rare variants. We performed replication across Hispanic/Latino subgroups, and the transferability to non-Hispanic/Latino populations was assessed in 1001 African Americans, 1250 Singaporeans and 941 Europeans with asthma. The effects of the variants on gene expression and DNA methylation from whole blood were also evaluated in participants with asthma and in silico with data obtained through public databases. RESULTS: Genomewide significant associations of Indigenous American ancestry with severe asthma exacerbations were found at 5q32 in Mexican Americans as well as at 13q13-q13.2 and 3p13 in Puerto Ricans. The single nucleotide polymorphism (SNP) rs1144986 (C5orf46) showed consistent effects for severe asthma exacerbations across Hispanic/Latino subgroups, but it was not validated in non-Hispanics/Latinos. This SNP was associated with DPYSL3 DNA methylation and SCGB3A2 gene expression levels. CONCLUSIONS: Admixture mapping study of asthma exacerbations revealed a novel locus that exhibited Hispanic/Latino-specific effects and regulated DPYSL3 and SCGB3A2.

Analysis of germline-driven ancestry-associated gene expression in cancers.

Differential mRNA expression between ancestry groups can be explained by both genetic and environmental factors. We outline a computational workflow to determine the extent to which germline genetic variation explains cancer-specific molecular differences across ancestry groups. Using multi-omics datasets from The Cancer Genome Atlas (TCGA), we enumerate ancestry-informative markers colocalized with cancer-type-specific expression quantitative trait loci (e-QTLs) at ancestry-associated genes. This approach is generalizable to other settings with paired germline genotyping and mRNA expression data for a multi-ethnic cohort. For complete details on the use and execution of this protocol, please refer to Carrot-Zhang et al. (2020), Robertson et al. (2021), and Sayaman et al. (2021).

Genome-wide association study in minority children with asthma implicates DNAH5 in bronchodilator responsiveness.

Variability in response to short-acting ß2-agonists (e.g., albuterol) among patients with asthma from diverse racial/ethnic groups may contribute to asthma disparities. We sought to identify genetic variants associated with bronchodilator response (BDR) to identify potential mechanisms of drug response and risk factors for worse asthma outcomes. Genome-wide association studies of bronchodilator response (BDR) were performed using TOPMed Whole Genome Sequencing data of the Asthma Translational Genomic Collaboration (ATGC), which corresponded to 1136 Puerto Rican, 656 Mexican and 4337 African American patients with asthma. With the population-specific GWAS results, a trans-ethnic meta-analysis was performed to identify BDR-associated variants shared across the three populations. Replication analysis was carried out in three pediatric asthma cohorts, including CAMP (Childhood Asthma Management Program; n = 560), GACRS (Genetics of Asthma in Costa Rica Study; n = 967) and HPR (Hartford-Puerto Rico; n = 417). A genome-wide significant locus (rs35661809; P = 3.61 × 10-8) in LINC02220, a non-coding RNA gene, was identified in Puerto Ricans. While this region was devoid of protein-coding genes, capture Hi-C data showed a distal interaction with the promoter of the DNAH5 gene in lung tissue. In replication analysis, the GACRS cohort yielded a nominal association (1-tailed P < 0.05). No genetic variant was associated with BDR at the genome-wide significant threshold in Mexicans and African Americans. Our findings help inform genetic underpinnings of BDR for understudied minority patients with asthma, but the limited availability of genetic data for racial/ethnic minority children with asthma remains a paramount challenge.

Association between Ancestry-Specific 6q25 Variants and Breast Cancer Subtypes in Peruvian Women.

BACKGROUND: Breast cancer incidence in the United States is lower in Hispanic/Latina (H/L) compared with African American/Black or Non-Hispanic White women. An Indigenous American breast cancer-protective germline variant (rs140068132) has been reported near the estrogen receptor 1 gene. This study tests the association of rs140068132 and other polymorphisms in the 6q25 region with subtype-specific breast cancer risk in H/Ls of high Indigenous American ancestry. METHODS: Genotypes were obtained for 5,094 Peruvian women with (1,755) and without (3,337) breast cancer. Associations between genotype and overall and subtype-specific risk for the protective variant were tested using logistic regression models and conditional analyses, including other risk-associated polymorphisms in the region. RESULTS: We replicated the reported association between rs140068132 and breast cancer risk overall [odds ratio (OR), 0.53; 95% confidence interval (CI), 0.47-0.59], as well as the lower odds of developing hormone receptor negative (HR-) versus HR+ disease (OR, 0.77; 95% CI, 0.61-0.97). Models, including HER2, showed further heterogeneity with reduced odds for HR+HER2+ (OR, 0.68; 95% CI, 0.51-0.92), HR-HER2+ (OR, 0.63; 95% CI, 0.44-0.90) and HR-HER2- (OR, 0.77; 95% CI, 0.56-1.05) compared with HR+HER2-. Inclusion of other risk-associated variants did not change these observations. CONCLUSIONS: The rs140068132 polymorphism is associated with decreased risk of breast cancer in Peruvians and is more protective against HR- and HER2+ diseases independently of other breast cancer-associated variants in the 6q25 region. IMPACT: These results could inform functional analyses to understand the mechanism by which rs140068132-G reduces risk of breast cancer development in a subtype-specific manner. They also illustrate the importance of including diverse individuals in genetic studies.

Clotting factor genes are associated with preeclampsia in high-altitude pregnant women in the Peruvian Andes.

Preeclampsia is a multi-organ complication of pregnancy characterized by sudden hypertension and proteinuria that is among the leading causes of preterm delivery and maternal morbidity and mortality worldwide. The heterogeneity of preeclampsia poses a challenge for understanding its etiology and molecular basis. Intriguingly, risk for the condition increases in high-altitude regions such as the Peruvian Andes. To investigate the genetic basis of preeclampsia in a population living at high altitude, we characterized genome-wide variation in a cohort of preeclamptic and healthy Andean families (n = 883) from Puno, Peru, a city located above 3,800 meters of altitude. Our study collected genomic DNA and medical records from case-control trios and duos in local hospital settings. We generated genotype data for 439,314 SNPs, determined global ancestry patterns, and mapped associations between genetic variants and preeclampsia phenotypes. A transmission disequilibrium test (TDT) revealed variants near genes of biological importance for placental and blood vessel function. The top candidate region was found on chromosome 13 of the fetal genome and contains clotting factor genes PROZ, F7, and F10. These findings provide supporting evidence that common genetic variants within coagulation genes play an important role in preeclampsia. A selection scan revealed a potential adaptive signal around the ADAM12 locus on chromosome 10, implicated in pregnancy disorders. Our discovery of an association in a functional pathway relevant to pregnancy physiology in an understudied population of Native American origin demonstrates the increased power of family-based study design and underscores the importance of conducting genetic research in diverse populations.

Nasal airway transcriptome-wide association study of asthma reveals genetically driven mucus pathobiology.

To identify genetic determinants of airway dysfunction, we performed a transcriptome-wide association study for asthma by combining RNA-seq data from the nasal airway epithelium of 681 children, with UK Biobank genetic association data. Our airway analysis identified 95 asthma genes, 58 of which were not identified by transcriptome-wide association analyses using other asthma-relevant tissues. Among these genes were MUC5AC, an airway mucin, and FOXA3, a transcriptional driver of mucus metaplasia. Muco-ciliary epithelial cultures from genotyped donors revealed that the MUC5AC risk variant increases MUC5AC protein secretion and mucus secretory cell frequency. Airway transcriptome-wide association analyses for mucus production and chronic cough also identified MUC5AC. These cis-expression variants were associated with trans effects on expression; the MUC5AC variant was associated with upregulation of non-inflammatory mucus secretory network genes, while the FOXA3 variant was associated with upregulation of type-2 inflammation-induced mucus-metaplasia pathway genes. Our results reveal genetic mechanisms of airway mucus pathobiology.

Epigenome-wide association study of lung function in Latino children and youth with asthma.

INTRODUCTION: DNA methylation studies have associated methylation levels at different CpG sites or genomic regions with lung function. Moreover, genetic ancestry has been associated with lung function in Latinos. However, no epigenome-wide association study (EWAS) of lung function has been performed in this population. Here, we aimed to identify DNA methylation patterns associated with lung function in pediatric asthma among Latinos. RESULTS: We conducted an EWAS in whole blood from 250 Puerto Rican and 148 Mexican American children and young adults with asthma. A total of five CpGs exceeded the genome-wide significance threshold of p = 1.17 × 10-7 in the combined analyses from Puerto Ricans and Mexican Americans: cg06035600 (MAP3K6, p = 6.13 × 10-8) showed significant association with pre-bronchodilator Tiffeneau-Pinelli index, the probes cg00914963 (TBC1D16, p = 1.04 × 10-7), cg16405908 (MRGPRE, p = 2.05 × 10-8), and cg07428101 (MUC2, p = 5.02 × 10-9) were associated with post-bronchodilator forced vital capacity (FVC), and cg20515679 (KCNJ6) with post-bronchodilator Tiffeneau-Pinelli index (p = 1.13 × 10-8). However, these markers did not show significant associations in publicly available data from Europeans (p > 0.05). A methylation quantitative trait loci analysis revealed that methylation levels at these CpG sites were regulated by genetic variation in Latinos and the Biobank-based Integrative Omics Studies (BIOS) consortium. Additionally, two differentially methylated regions in REXOC and AURKC were associated with pre-bronchodilator Tiffeneau-Pinelli index (adjusted p < 0.05) in Puerto Ricans and Mexican Americans. Moreover, we replicated some of the previous differentially methylated signals associated with lung function in non-Latino populations. CONCLUSIONS: We replicated previous associations of epigenetic markers with lung function in whole blood and identified novel population-specific associations shared among Latino subgroups.

Analytic pipelines to assess the relationship between immune response and germline genetics in human tumors.

Germline genetic variants modulate human immune response. We present analytical pipelines for assessing the contribution of hosts' genetic background to the immune landscape of solid tumors using harmonized data from more than 9,000 patients in The Cancer Genome Atlas (TCGA). These include protocols for heritability, genome-wide association studies (GWAS), colocalization, and rare variant analyses. These workflows are developed around the structure of TCGA but can be adapted to explore other repositories or in the context of cancer immunotherapy. For complete details on the use and execution of this protocol, please refer to Sayaman et al. (2021).

Racial/ethnic differences in eligibility for asthma biologics among pediatric populations.

BACKGROUND: Asthma is a heterogeneous disease. Clinical blood parameters differ by race/ethnicity and are used to distinguish asthma subtypes and inform therapies. Differences in subtypes may explain population-specific trends in asthma outcomes. However, these differences in racial/ethnic minority pediatric populations are unclear. OBJECTIVE: We investigated the association of blood parameters and asthma subtypes with asthma outcomes and examined population-specific eligibility for biologic therapies in minority pediatric populations. METHODS: Using data from 2 asthma case-control studies of pediatric minority populations, we performed case-control (N = 3738) and case-only (N = 2743) logistic regressions to quantify the association of blood parameters and asthma subtypes with asthma outcomes. Heterogeneity of these associations was tested using an interaction term between race/ethnicity and each exposure. Differences in therapeutic eligibility were investigated using chi-square tests. RESULTS: Race/ethnicity modified the association between total IgE and asthma exacerbations. Elevated IgE level was associated with worse asthma outcomes in Puerto Ricans. Allergic asthma was associated with worse outcomes in Mexican Americans, whereas eosinophilic asthma was associated with worse outcomes in Puerto Ricans. A lower proportion of Puerto Ricans met dosing criteria for allergic asthma-directed biologic therapy than other groups. A higher proportion of Puerto Ricans qualified for eosinophilic asthma-directed biologic therapy than African Americans. CONCLUSIONS: We found population-specific associations between blood parameters and asthma subtypes with asthma outcomes. Our findings suggest that eligibility for asthma biologic therapies differs across pediatric racial/ethnic populations. These findings call for more studies in diverse populations for equitable treatment of minority patients with asthma.