Rare coding variants pinpoint genes that control human hematological traits

Red Blood Cell-Related Phenotype-Genotype Correlations in Chronic and Acute Critical Illnesses (Traumatic Brain Injury Cohort and COVID-19 Cohort)

Changes in red blood cell (RBC)-related parameters and anemia are common in both severe chronic and acute diseases. RBC-related phenotypes have a heritable component. However, it is unclear whether the contribution of genetic variability is pronounced when hematological parameters are affected by physiological stress. In this study, we analyzed RBC-related phenotypes and phenotype-genotype correlations in two exome-sequenced patient cohorts with or at a high risk for a critical illness: chronic TBI patients admitted for rehabilitation and patients with acute COVID-19. In the analysis of exome data, we focused on the cumulative effects of rare high-impact variants (qualifying variants, QVs) in specific gene sets, represented by Notch signaling pathway genes, based on the results of enrichment analysis in anemic TBI patients and three predefined gene sets for phenotypes of interest derived from GO, GWAS, and HPO resources. In both patient cohorts, anemia was associated with the cumulative effects of QVs in the GO (TBI: p = 0.0003, OR = 2.47 (1.54-4.88); COVID-19: p = 0.0004, OR = 2.12 (1.39-3.25)) and Notch pathway-derived (TBI: p = 0.0017, OR = 2.33 (1.35-4.02); COVID-19: p = 0.0012, OR =8.00 (1.79-35.74)) gene sets. In the multiple linear regression analysis, genetic variables contributed to RBC indices in patients with TBI. In COVID-19 patients, QVs in Notch pathway genes influenced RBC, HGB, and HCT levels, whereas genes from other sets influenced MCHC levels. Thus, in this exploratory study, exome data analysis yielded similar and different results in the two patient cohorts, supporting the view that genetic factors may contribute to RBC-related phenotypic performance in both severe chronic and acute health conditions.

Plasminogen missense variants and their involvement in cardiovascular and inflammatory disease

Human plasminogen (PLG), the zymogen of the fibrinolytic protease, plasmin, is a polymorphic protein with two widely distributed codominant alleles, PLG/Asp453 and PLG/Asn453. About 15 other missense or non-synonymous single nucleotide polymorphisms (nsSNPs) of PLG show major, yet different, relative abundances in world populations. Although the existence of these relatively abundant allelic variants is generally acknowledged, they are often overlooked or assumed to be non-pathogenic. In fact, at least half of those major variants are classified as having conflicting pathogenicity, and it is unclear if they contribute to different molecular phenotypes. From those, PLG/K19E and PLG/A601T are examples of two relatively abundant PLG variants that have been associated with PLG deficiencies (PD), but their pathogenic mechanisms are unclear. On the other hand, approximately 50 rare and ultra-rare PLG missense variants have been reported to cause PD as homozygous or compound heterozygous variants, often leading to a debilitating disease known as ligneous conjunctivitis. The true abundance of PD-associated nsSNPs is unknown since they can remain undetected in heterozygous carriers. However, PD variants may also contribute to other diseases. Recently, the ultra-rare autosomal dominant PLG/K311E has been found to be causative of hereditary angioedema (HAE) with normal C1 inhibitor. Two other rare pathogenic PLG missense variants, PLG/R153G and PLG/V709E, appear to affect platelet function and lead to HAE, respectively. Herein, PLG missense variants that are abundant and/or clinically relevant due to association with disease are examined along with their world distribution. Proposed molecular mechanisms are discussed when known or can be reasonably assumed.

A genome-wide association study of neutrophil count in individuals associated to an African continental ancestry group facilitates studies of malaria pathogenesis

BACKGROUND

'Benign ethnic neutropenia' (BEN) is a heritable condition characterized by lower neutrophil counts, predominantly observed in individuals of African ancestry, and the genetic basis of BEN remains a subject of extensive research. In this study, we aimed to dissect the genetic architecture underlying neutrophil count variation through a linear-mixed model genome-wide association study (GWAS) in a population of African ancestry (N = 5976). Malaria caused by P. falciparum imposes a tremendous public health burden on people living in sub-Saharan Africa. Individuals living in malaria endemic regions often have a reduced circulating neutrophil count due to BEN, raising the possibility that reduced neutrophil counts modulate severity of malaria in susceptible populations. As a follow-up, we tested this hypothesis by conducting a Mendelian randomization (MR) analysis of neutrophil counts on severe malaria (MalariaGEN, N = 17,056).

RESULTS

We carried out a GWAS of neutrophil count in individuals associated to an African continental ancestry group within UK Biobank, identifying 73 loci (r2 = 0.1) and 10 index SNPs (GCTA-COJO loci) associated with neutrophil count, including previously unknown rare loci regulating neutrophil count in a non-European population. BOLT-LMM was reliable when conducted in a non-European population, and additional covariates added to the model did not largely alter the results of the top loci or index SNPs. The two-sample bi-directional MR analysis between neutrophil count and severe malaria showed the greatest evidence for an effect between neutrophil count and severe anaemia, although the confidence intervals crossed the null.

CONCLUSION

Our GWAS of neutrophil count revealed unique loci present in individuals of African ancestry. We note that a small sample-size reduced our power to identify variants with low allele frequencies and/or low effect sizes in our GWAS. Our work highlights the need for conducting large-scale biobank studies in Africa and for further exploring the link between neutrophils and severe malaria.

Gene-based association study of rare variants in children of diverse ancestries implicates TNFRSF21 in the development of allergic asthma

BACKGROUND

Most genetic studies of asthma and allergy have focused on common variation in individuals primarily of European ancestry. Studying the role of rare variation in quantitative phenotypes and in asthma phenotypes in populations of diverse ancestries can provide additional, important insights into the development of these traits.

OBJECTIVE

We sought to examine the contribution of rare variants to different asthma- or allergy-associated quantitative traits in children with diverse ancestries and explore their role in asthma phenotypes.

METHODS

We examined whole-genome sequencing data from children participants in longitudinal studies of asthma (n = 1035; parent-identified as 67% Black and 25% Hispanic) to identify rare variants (minor allele frequency < 0.01). We assigned variants to genes and tested for associations using an omnibus variant-set test between each of 24,902 genes and 8 asthma-associated quantitative traits. On combining our results with external data on predicted gene expression in humans and mouse knockout studies, we identified 3 candidate genes. A burden of rare variants in each gene and in a combined 3-gene score was tested for its associations with clinical phenotypes of asthma. Finally, published single-cell gene expression data in lower airway mucosal cells after allergen challenge were used to assess transcriptional responses to allergen.

RESULTS

Rare variants in USF1 were significantly associated with blood neutrophil count (P = 2.18 × 10-7); rare variants in TNFRSF21 with total IgE (P = 6.47 × 10-6) and PIK3R6 with eosinophil count (P = 4.10 × 10-5) reached suggestive significance. These 3 findings were supported by independent data from human and mouse studies. A burden of rare variants in TNFRSF21 and in a 3-gene score was associated with allergy-related phenotypes in cohorts of children with mild and severe asthma. Furthermore, TNFRSF21 was significantly upregulated in bronchial basal epithelial cells from adults with allergic asthma but not in adults with allergies (but not asthma) after allergen challenge.

CONCLUSIONS

We report novel associations between rare variants in genes and allergic and inflammatory phenotypes in children with diverse ancestries, highlighting TNFRSF21 as contributing to the development of allergic asthma.

BIGKnock: fine-mapping gene-based associations via knockoff analysis of biobank-scale data

We propose BIGKnock (BIobank-scale Gene-based association test via Knockoffs), a computationally efficient gene-based testing approach for biobank-scale data, that leverages long-range chromatin interaction data, and performs conditional genome-wide testing via knockoffs. BIGKnock can prioritize causal genes over proxy associations at a locus. We apply BIGKnock to the UK Biobank data with 405,296 participants for multiple binary and quantitative traits, and show that relative to conventional gene-based tests, BIGKnock produces smaller sets of significant genes that contain the causal gene(s) with high probability. We further illustrate its ability to pinpoint potential causal genes at [Formula: see text] of the associated loci.

The link between allergic rhinitis and chronic rhinosinusitis

PURPOSE OF REVIEW

Allergic rhinitis and chronic rhinosinusitis (CRS) are common disorders affecting millions of people worldwide. Although allergic rhinitis and CRS are distinct clinical entities, certain CRS endotypes share similar pathological mechanisms as those seen in patients with allergic rhinitis. This review assesses the literature behind the similarities and differences seen in patients with CRS and allergic rhinitis, and the role atopy might play in the pathophysiology of CRS.

RECENT FINDINGS

In examining the associations between allergic rhinitis and CRS, most studies have focused primarily on CRS with nasal polyps and type 2 inflammation in CRS. Recent studies have demonstrated the similarities and differences in pathologic mechanisms behind allergic rhinitis and CRS, with an emphasis on patient endotypes, genetics, and the nasoepithelial immunologic barrier. Related immunopathology shared by allergic rhinitis and type 2 inflammation in CRS has allowed for therapeutic overlap with biologic treatments.

SUMMARY

Allergic rhinitis and CRS often present as comorbid conditions, and understanding the relationship between allergic rhinitis and CRS is important when considering treatment options. Advances in understanding the genetics and immunology, as well as biologic and immunotherapeutic treatments have improved outcomes in patients with CRS, especially in the setting of atopy.

The contribution of common and rare genetic variants to variation in metabolic traits in 288,137 East Asians

Metabolic traits are heritable phenotypes widely-used in assessing the risk of various diseases. We conduct a genome-wide association analysis (GWAS) of nine metabolic traits (including glycemic, lipid, liver enzyme levels) in 125,872 Korean subjects genotyped with the Korea Biobank Array. Following meta-analysis with GWAS from Biobank Japan identify 144 novel signals (MAF ≥ 1%), of which 57.0% are replicated in UK Biobank. Additionally, we discover 66 rare (MAF < 1%) variants, 94.4% of them co-incident to common loci, adding to allelic series. Although rare variants have limited contribution to overall trait variance, these lead, in carriers, substantial loss of predictive accuracy from polygenic predictions of disease risk from common variant alone. We capture groups with up to 16-fold variation in type 2 diabetes (T2D) prevalence by integration of genetic risk scores of fasting plasma glucose and T2D and the I349F rare protective variant. This study highlights the need to consider the joint contribution of both common and rare variants on inherited risk of metabolic traits and related diseases.

Dual Immune Regulatory Roles of Interleukin-33 in Pathological Conditions

Interleukin-33 (IL-33), a member of the IL-1 cytokine family and a multifunctional cytokine, plays critical roles in maintaining host homeostasis and in pathological conditions, such as allergy, infectious diseases, and cancer, by acting on multiple types of immune cells and promoting type 1 and 2 immune responses. IL-33 is rapidly released by immune and non-immune cells upon stimulation by stress, acting as an “alarmin” by binding to its receptor, suppression of tumorigenicity 2 (ST2), to trigger downstream signaling pathways and activate inflammatory and immune responses. It has been recognized that IL-33 displays dual-functioning immune regulatory effects in many diseases and has both pro- and anti-tumorigenic effects, likely depending on its primary target cells, IL-33/sST2 expression levels, cellular context, and the cytokine microenvironment. Herein, we summarize our current understanding of the biological functions of IL-33 and its roles in the pathogenesis of various conditions, including inflammatory and autoimmune diseases, infections, cancers, and cases of organ transplantation. We emphasize the nature of context-dependent dual immune regulatory functions of IL-33 in many cells and diseases and review systemic studies to understand the distinct roles of IL-33 in different cells, which is essential to the development of more effective diagnoses and therapeutic approaches for IL-33-related diseases.

Genetic Influence on Frequencies of Myeloid-Derived Cell Subpopulations in Mouse

Differences in frequencies of blood cell subpopulations were reported to influence the course of infections, atopic and autoimmune diseases, and cancer. We have discovered a unique mouse strain B10.O20 containing extremely high frequency of myeloid-derived cells (MDC) in spleen. B10.O20 carries 3.6% of genes of the strain O20 on the C57BL/10 genetic background. It contains much higher frequency of CD11b⁺Gr1⁺ cells in spleen than both its parents. B10.O20 carries O20-derived segments on chromosomes 1, 15, 17, and 18. Their linkage with frequencies of blood cell subpopulations in spleen was tested in F₂ hybrids between B10.O20 and C57BL/10. We found 3 novel loci controlling MDC frequencies: Mydc1, 2, and 3 on chromosomes 1, 15, and 17, respectively, and a locus controlling relative spleen weight (Rsw1) that co-localizes with Mydc3 and also influences proportion of white and red pulp in spleen. Mydc1 controls numbers of CD11b⁺Gr1⁺ cells. Interaction of Mydc2 and Mydc3 regulates frequency of CD11b⁺Gr1⁺ cells and neutrophils (Gr1⁺Siglec-F^- cells from CD11b⁺ cells). Interestingly, Mydc3/Rsw1 is orthologous with human segment 6q21 that was shown previously to determine counts of white blood cells. Bioinformatics analysis of genomic sequence of the chromosomal segments bearing these loci revealed polymorphisms between O20 and C57BL/10 that change RNA stability and genes’ functions, and we examined expression of relevant genes. This identified potential candidate genes Smap1, Vps52, Tnxb, and Rab44. Definition of genetic control of MDC can help to personalize therapy of diseases influenced by these cells.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

From musk to body odor: Decoding olfaction through genetic variation

The olfactory system combines input from multiple receptor types to represent odor information, but there are few explicit examples relating olfactory receptor (OR) activity patterns to odor perception. To uncover these relationships, we performed genome-wide scans on odor-perception phenotypes for ten odors in 1000 Han Chinese and validated results for six of these odors in an ethnically diverse population (n = 364). In both populations, consistent with previous studies, we replicated three previously reported associations (β-ionone/OR5A1, androstenone/OR7D4, cis-3-hexen-1-ol/OR2J3 LD-band), but not for odors containing aldehydes, suggesting that olfactory phenotype/genotype studies are robust across populations. Two novel associations between an OR and odor perception contribute to our understanding of olfactory coding. First, we found a SNP in OR51B2 that associated with trans-3-methyl-2-hexenoic acid, a key component of human underarm odor. Second, we found two linked SNPs associated with the musk Galaxolide in a novel musk receptor, OR4D6, which is also the first human OR shown to drive specific anosmia to a musk compound. We noticed that SNPs detected for odor intensity were enriched with amino acid substitutions, implying functional changes of odor receptors. Furthermore, we also found that the derived alleles of the SNPs tend to be associated with reduced odor intensity, supporting the hypothesis that the primate olfactory gene repertoire has degenerated over time. This study provides information about coding for human body odor, and gives us insight into broader mechanisms of olfactory coding, such as how differential OR activation can converge on a similar percept.

eSCAN: scan regulatory regions for aggregate association testing using whole-genome sequencing data

Multiple statistical methods for aggregate association testing have been developed for whole-genome sequencing (WGS) data. Many aggregate variants in a given genomic window and ignore existing knowledge to define test regions, resulting in many identified regions not clearly linked to genes, and thus, limiting biological understanding. Functional information from new technologies (such as Hi-C and its derivatives), which can help link enhancers to their effector genes, can be leveraged to predefine variant sets for aggregate testing in WGS data. Here, we propose the eSCAN (scan the enhancers) method for genome-wide assessment of enhancer regions in sequencing studies, combining the advantages of dynamic window selection in SCANG (SCAN the Genome), a previously developed method, with the advantages of incorporating putative regulatory regions from annotation. eSCAN, by searching in putative enhancers, increases statistical power and aids mechanistic interpretation, as demonstrated by extensive simulation studies. We also apply eSCAN for blood cell traits using NHLBI Trans-Omics for Precision Medicine WGS data. Results from real data analysis show that eSCAN is able to capture more significant signals, and these signals are of shorter length (indicating higher resolution fine-mapping capability) and drive association of larger regions detected by other methods.

Multi-ethnic genome-wide association analyses of white blood cell and platelet traits in the Population Architecture using Genomics and Epidemiology (PAGE) study

BACKGROUND

Circulating white blood cell and platelet traits are clinically linked to various disease outcomes and differ across individuals and ancestry groups. Genetic factors play an important role in determining these traits and many loci have been identified. However, most of these findings were identified in populations of European ancestry (EA), with African Americans (AA), Hispanics/Latinos (HL), and other races/ethnicities being severely underrepresented.

RESULTS

We performed ancestry-combined and ancestry-specific genome-wide association studies (GWAS) for white blood cell and platelet traits in the ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) Study, including 16,201 AA, 21,347 HL, and 27,236 EA participants. We identified six novel findings at suggestive significance (P < 5E-8), which need confirmation, and independent signals at six previously established regions at genome-wide significance (P < 2E-9). We confirmed multiple previously reported genome-wide significant variants in the single variant association analysis and multiple genes using PrediXcan. Evaluation of loci reported from a Euro-centric GWAS indicated attenuation of effect estimates in AA and HL compared to EA populations.

CONCLUSIONS

Our results highlighted the potential to identify ancestry-specific and ancestry-agnostic variants in participants with diverse backgrounds and advocate for continued efforts in improving inclusion of racially/ethnically diverse populations in genetic association studies for complex traits.

Whole-genome sequencing association analysis of quantitative red blood cell phenotypes: The NHLBI TOPMed program

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

An Asian-specific MPL genetic variant alters JAK-STAT signaling and influences platelet count in the population

Genomic discovery efforts for hematological traits have been successfully conducted through genome-wide association study on samples of predominantly European ancestry. We sought to conduct unbiased genetic discovery for coding variants that influence hematological traits in a Han Chinese population. A total of 5257 Han Chinese subjects from Beijing, China were included in the discovery cohort and analyzed by an Illumina ExomeChip array. Replication analyses were conducted in 3827 independent Chinese subjects. We analyzed 12 hematological traits and identified 22 exome-wide significant single-nucleotide polymorphisms (SNP)-trait associations with 15 independent SNPs. Our study provides replication for two associations previously reported but not replicated. Further, one association was identified and replicated in the current study, of a coding variant in the myeloproliferative leukemia (MPL) gene, c.793C > T, p.Leu265Phe (L265F) with increased platelet count (β = 20.6 109 cells/l, Pmeta-analysis = 2.6 × 10-13). This variant is observed at ~2% population frequency in East Asians, whereas it has not been reported in gnomAD European or African populations. Functional analysis demonstrated that expression of MPL L265F in Ba/F3 cells resulted in enhanced phosphorylation of Stat3 and ERK1/2 as compared with the reference MPL allele, supporting altered activation of the JAK-STAT signal transduction pathway as the mechanism underlying the novel association between MPL L265F and platelet count.

Heterozygous germ line CSF3R variants as risk alleles for development of hematologic malignancies

Colony-stimulating factor 3 receptor (CSF3R) encodes the receptor for granulocyte colony-stimulating factor (G-CSF), a cytokine vital for granulocyte proliferation and differentiation. Acquired activating heterozygous variants in CSF3R are the main cause of chronic neutrophilic leukemia, a hyperproliferative disorder. In contrast, biallelic germ line hypomorphic variants in CSF3R are a rare cause of severe congenital neutropenia, a hypoproliferative condition. The impact of heterozygous germ line CSF3R variants, however, is unknown. We identified CSF3R as a new germ line hematologic malignancy predisposition gene through analysis of 832 next-generation sequencing tests conducted in 632 patients with hematologic malignancies. Among germ line CSF3R variants, 3 were abnormal in functional testing, indicating their deleterious nature. p.Trp547* was identified in 2 unrelated men with myelodysplastic syndromes diagnosed at 76 and 33 years of age, respectively. p.Trp547* is a loss-of-function nonsense variant in the extracellular domain that results in decreased CSF3R messenger RNA expression and abrogation of CSF3R surface expression and proliferative responses to G-CSF. p.Ala119Thr is a missense variant found in 2 patients with multiple myeloma and acute lymphoblastic leukemia, respectively. This variant is located between the extracellular immunoglobulin-like and cytokine receptor homology domains and results in decreased G-CSF sensitivity. p.Pro784Thr was identified in a 67-year-old man with multiple myeloma. p.Pro784Thr is a missense variant in the cytoplasmic domain that inhibits CSF3R internalization, producing a gain-of-function phenotype and G-CSF hypersensitivity. Our findings identify germ line heterozygous CSF3R variants as risk factors for development of myeloid and lymphoid malignancies.

Integrative Genomic Analysis Reveals Four Protein Biomarkers for Platelet Traits

RATIONALE

Mean platelet volume (MPV) and platelet count (PLT) are platelet measures that have been linked to cardiovascular disease (CVD) and mortality risk. Identifying protein biomarkers for these measures may yield insights into CVD mechanisms.

OBJECTIVE

We aimed to identify causal protein biomarkers for MPV and PLT among 71 CVD-related plasma proteins measured in FHS (Framingham Heart Study) participants.

METHODS AND RESULTS

We conducted integrative analyses of genetic variants associated with PLT/MPV with protein quantitative trait locus variants associated with plasma proteins followed by Mendelian randomization to infer causal relations of proteins for PLT/MPV. We also tested protein-PLT/MPV association in FHS participants. Using induced pluripotent stem cell-derived megakaryocyte clones that produce functional platelets, we conducted RNA-sequencing and analyzed expression differences between low- and high-platelet producing clones. We then performed small interfering RNA gene knockdown experiments targeting genes encoding proteins with putatively causal platelet effects in megakaryocyte clones to examine effects on platelet production. In protein-trait association analyses, ten proteins were associated with MPV and 31 with PLT. Mendelian randomization identified 4 putatively causal proteins for MPV and 4 for PLT. GP-5 (Glycoprotein V), GRN (granulin), and MCAM (melanoma cell adhesion molecule) were associated with PLT, while MPO (myeloperoxidase) showed significant association with MPV in both analyses. RNA-sequencing analysis results were directionally concordant with observed and Mendelian randomization-inferred associations for GP-5, GRN, and MCAM. In siRNA gene knockdown experiments, silencing GP-5, GRN, and MPO decreased PLTs. Genome-wide association study results suggest several of these may be linked to CVD risk.

CONCLUSIONS

We identified 4 proteins that are causally linked to PLTs. These proteins may also have roles in the pathogenesis of CVD via a platelet/blood coagulation-based mechanism.

Transcriptomic signatures across human tissues identify functional rare genetic variation

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits.

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Estimation of non-null SNP effect size distributions enables the detection of enriched genes underlying complex traits

Traditional univariate genome-wide association studies generate false positives and negatives due to difficulties distinguishing associated variants from variants with spurious nonzero effects that do not directly influence the trait. Recent efforts have been directed at identifying genes or signaling pathways enriched for mutations in quantitative traits or case-control studies, but these can be computationally costly and hampered by strict model assumptions. Here, we present gene-ε, a new approach for identifying statistical associations between sets of variants and quantitative traits. Our key insight is that enrichment studies on the gene-level are improved when we reformulate the genome-wide SNP-level null hypothesis to identify spurious small-to-intermediate SNP effects and classify them as non-causal. gene-ε efficiently identifies enriched genes under a variety of simulated genetic architectures, achieving greater than a 90% true positive rate at 1% false positive rate for polygenic traits. Lastly, we apply gene-ε to summary statistics derived from six quantitative traits using European-ancestry individuals in the UK Biobank, and identify enriched genes that are in biologically relevant pathways.

Use of >100,000 NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium whole genome sequences improves imputation quality and detection of rare variant associations in admixed African and Hispanic/Latino populations

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count <10 (including singletons) in the imputation target samples, average information content rescued was >86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

The Cytokines of Asthma

Asthma is a chronic inflammatory airway disease associated with type 2 cytokines interleukin-4 (IL-4), IL-5, and IL-13, which promote airway eosinophilia, mucus overproduction, bronchial hyperresponsiveness (BHR), and immunogloubulin E (IgE) synthesis. However, only half of asthma patients exhibit signs of an exacerbated Type 2 response. "Type 2-low" asthma has different immune features: airway neutrophilia, obesity-related systemic inflammation, or in some cases, few signs of immune activation. Here, we review the cytokine networks driving asthma, placing these in cellular context and incorporating insights from cytokine-targeting therapies in the clinic. We discuss established and emerging paradigms in the context of the growing appreciation of disease heterogeneity and argue that the development of new and improved therapeutics will require understanding the diverse mechanisms underlying the spectrum of asthma pathologies.

Evolutionary history of disease-susceptibility loci identified in longitudinal exome-wide association studies

Background

Our longitudinal exome‐wide association studies previously detected various genetic determinants of complex disorders using ~26,000 single‐nucleotide polymorphisms (SNPs) that passed quality control and longitudinal medical examination data (mean follow‐up period, 5 years) in 4884–6022 Japanese subjects. We found that allele frequencies of several identified SNPs were remarkably different among four ethnic groups. Elucidating the evolutionary history of disease‐susceptibility loci may help us uncover the pathogenesis of the related complex disorders.

Methods

In the present study, we conducted evolutionary analyses such as extended haplotype homozygosity, focusing on genomic regions containing disease‐susceptibility loci and based on genotyping data of our previous studies and datasets from the 1000 Genomes Project.

Results

Our evolutionary analyses suggest that derived alleles of rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs closely located at 12q24.1 associated with type 2 diabetes mellitus, obesity, dyslipidemia, and three complex disorders (hypertension, hyperuricemia, and dyslipidemia), respectively, rapidly expanded after the human dispersion from Africa (Out‐of‐Africa). Allele frequencies of GGA3 and six SNPs at 12q24.1 appeared to have remarkably changed in East Asians, whereas the derived alleles of rs34902660 of SLC17A3 and rs7656604 at 4q13.3 might have spread across Japanese and non‐Africans, respectively, although we cannot completely exclude the possibility that allele frequencies of disease‐associated loci may be affected by demographic events.

Conclusion

Our findings indicate that derived allele frequencies of nine disease‐associated SNPs (rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs at 12q24.1) identified in the longitudinal exome‐wide association studies largely increased in non‐Africans after Out‐of‐Africa.

Identification of rare variants of allergic rhinitis based on whole genome sequencing and gene expression profiling: A preliminary investigation in four families

Background

Despite the success of genome-wide association studies for allergic rhinitis (AR), no definitive causal variants have been identified, and a substantial portion of the heritability of the disease is yet to be discovered.

Methods

Four families, each with at least 1 parent and one child suffering from dust mite (DM) AR, were recruited, and whole-genome sequencing was performed on samples from 9 eligible individuals from these families. Conjoint analysis was performed for existing gene expression profiling data in the literature and the whole genome sequencing data obtained for these individuals; for presence of family-specific variants segregating with AR and the pathways involved. Similar analyses were also performed with data obtained for 96 sporadic house dust mite (HDM) AR patients and 96 healthy controls.

Results

Three rare variants in three genes (FLT1_c.603A ​> ​T; VEGFB_c.322A ​> ​C; and ITGA2_c.502+1G ​> ​A), which are involved in Focal Adhesion pathway, were identified in affected, but not unaffected, subjects in two families. VEGFB_c.322A ​> ​C and/or ITGA2_c.502+1G ​> ​A were further detected in all DM AR patients but not in any healthy individuals in 1 family; which was further investigated for members. The 3 identified variants were not found in any of the sporadic DM AR patients or healthy controls.

Conclusion

Despite the relatively small sample size, this study has identified several potentially functional rare variants in AR candidate genes, and it provides a platform for future work in larger numbers of families and sporadic individuals for a better understanding of the genetic basis of AR.

A large-scale exome array analysis of venous thromboembolism

Although recent Genome-Wide Association Studies have identified novel associations for common variants, there has been no comprehensive exome-wide search for low-frequency variants that affect the risk of venous thromboembolism (VTE). We conducted a meta-analysis of 11 studies comprising 8,332 cases and 16,087 controls of European ancestry and 382 cases and 1,476 controls of African American ancestry genotyped with the Illumina HumanExome BeadChip. We used the seqMeta package in R to conduct single variant and gene-based rare variant tests. In the single variant analysis, we limited our analysis to the 64,794 variants with at least 40 minor alleles across studies (minor allele frequency [MAF] ~0.08%). We confirmed associations with previously identified VTE loci, including ABO, F5, F11, and FGA. After adjusting for multiple testing, we observed no novel significant findings in single variant or gene-based analysis. Given our sample size, we had greater than 80% power to detect minimum odds ratios greater than 1.5 and 1.8 for a single variant with MAF of 0.01 and 0.005, respectively. Larger studies and sequence data may be needed to identify novel low-frequency and rare variants associated with VTE risk.

14q32 and let-7 microRNAs regulate transcriptional networks in fetal and adult human erythroblasts

In humans, fetal erythropoiesis takes place in the liver whereas adult erythropoiesis occurs in the bone marrow. Fetal and adult erythroid cells are not only produced at different sites, but are also distinguished by their respective transcriptional program. In particular, whereas fetal erythroid cells express γ-globin chains to produce fetal hemoglobin (HbF), adult cells express β-globin chains to generate adult hemoglobin. Understanding the transcriptional regulation of the fetal-to-adult hemoglobin switch is clinically important as re-activation of HbF production in adult erythroid cells would represent a promising therapy for the hemoglobin disorders sickle cell disease and β-thalassemia. We used RNA-sequencing to measure global gene and microRNA (miRNA) expression in human erythroblasts derived ex vivo from fetal liver (n = 12 donors) and bone marrow (n = 12 donors) hematopoietic stem/progenitor cells. We identified 7829 transcripts and 402 miRNA that were differentially expressed (false discovery rate <5%). The miRNA expression patterns were replicated in an independent collection of human erythroblasts using a different technology. By combining gene and miRNA expression data, we developed transcriptional networks which show substantial differences between fetal and adult human erythroblasts. Our analyses highlighted the miRNAs at the imprinted 14q32 locus in fetal erythroblasts and the let-7 miRNA family in adult erythroblasts as key regulators of stage-specific erythroid transcriptional programs. Altogether, our results provide a comprehensive resource to prioritize genes that may modify clinical severity in red blood cell (RBC) disorders, or genes that might be implicated in erythropoiesis by genome-wide association studies of RBC traits.

A loss-of-function variant in ALOX15 protects against nasal polyps and chronic rhinosinusitis

Nasal polyps (NP) are lesions on the nasal and paranasal sinus mucosa and are a risk factor for chronic rhinosinusitis (CRS). We performed genome-wide association studies on NP and CRS in Iceland and the UK (using UK Biobank data) with 4,366 NP cases, 5,608 CRS cases, and >700,000 controls. We found 10 markers associated with NP and 2 with CRS. We also tested 210 markers reported to associate with eosinophil count, yielding 17 additional NP associations. Of the 27 NP signals, 7 associate with CRS and 13 with asthma. Most notably, a missense variant in ALOX15 that causes a p.Thr560Met alteration in arachidonate 15-lipoxygenase (15-LO) confers large genome-wide significant protection against NP (P = 8.0 × 10^-27, odds ratio = 0.32; 95% confidence interval = 0.26, 0.39) and CRS (P = 1.1 × 10^-8, odds ratio = 0.64; 95% confidence interval = 0.55, 0.75). p.Thr560Met, carried by around 1 in 20 Europeans, was previously shown to cause near total loss of 15-LO enzymatic activity. Our findings identify 15-LO as a potential target for therapeutic intervention in NP and CRS.

Whole-Genome Sequencing in Severe Chronic Obstructive Pulmonary Disease

Genome-wide association studies have identified common variants associated with chronic obstructive pulmonary disease (COPD). Whole-genome sequencing (WGS) offers comprehensive coverage of the entire genome, as compared with genotyping arrays or exome sequencing. We hypothesized that WGS in subjects with severe COPD and smoking control subjects with normal pulmonary function would allow us to identify novel genetic determinants of COPD. We sequenced 821 patients with severe COPD and 973 control subjects from the COPDGene and Boston Early-Onset COPD studies, including both non-Hispanic white and African American individuals. We performed single-variant and grouped-variant analyses, and in addition, we assessed the overlap of variants between sequencing- and array-based imputation. Our most significantly associated variant was in a known region near HHIP (combined P = 1.6 × 10-9); additional variants approaching genome-wide significance included previously described regions in CHRNA5, TNS1, and SERPINA6/SERPINA1 (the latter in African American individuals). None of our associations were clearly driven by rare variants, and we found minimal evidence of replication of genes identified by previously reported smaller sequencing studies. With WGS, we identified more than 20 million new variants, not seen with imputation, including more than 10,000 of potential importance in previously identified COPD genome-wide association study regions. WGS in severe COPD identifies a large number of potentially important functional variants, with the strongest associations being in known COPD risk loci, including HHIP and SERPINA1. Larger sample sizes will be needed to identify associated variants in novel regions of the genome.

Deciphering the Emerging Complexities of Molecular Mechanisms at GWAS Loci

Genome-wide association studies (GWASs) have identified thousands of loci associated with hundreds of complex diseases and traits, and progress is being made toward elucidating the causal variants and genes underlying these associations. Functional characterization of mechanisms at GWAS loci is a multi-faceted challenge. Challenges include linkage disequilibrium and allelic heterogeneity at each locus, the noncoding nature of most loci, and the time and cost needed for experimentally evaluating the potential mechanistic contributions of genes and variants. As GWAS sample sizes increase, more loci are identified, and the complexities of individual loci emerge. Loci can consist of multiple association signals, each of which can reflect the influence of multiple variants, inseparable by association analyses. Each signal within a locus can influence the same or different target genes. Experimental studies of genes and variants can differ on the basis of cell type, cellular environment, or other context-specific variables. In this review, we describe the complexity of mechanisms at GWAS loci-including multiple signals, multiple variants, and/or multiple genes-and the implications these complexities hold for experimental study design and interpretation of GWAS mechanisms.

Learning by counting blood platelets in population studies: survey and perspective a long way after Bizzozero

Platelet count represents a useful tool in clinical practice to discriminate individuals at higher risk of bleeding. Less obvious is the role of platelet count variability within the normal range of distribution in shaping the individual's disease risk profile. Epidemiological studies have shown that platelet count in the adult general population is associated with a number of health outcomes related to hemostasis and thrombosis. However, recent studies are suggesting a possible role of this platelet index also as an independent risk factor. In this review of adult population studies, we will first focus on known genetic and non-genetic determinants of platelet number variability. Next, we will evaluate platelet count as a marker and/or a predictor of disease risk and its interaction with other risk factors. We will then discuss the role of platelet count variability within the normal distribution range as a contribution to disease and mortality risk. The possibility of considering platelet count as a simple, inexpensive indicator of increased risk of disease and death in general populations could open new opportunities to investigate novel platelet pathophysiological roles as well as therapeutic opportunities. Future studies should also consider platelet count, not only platelet function, as a modulator of disease and mortality risk.

Identification of nine novel loci related to hematological traits in a Japanese population

Recent genome-wide association studies have identified various genetic variants associated with hematological traits. Although it is possible that quantitative data of hematological traits are varied among the years examined, conventional genome-wide association studies have been conducted in a cross-sectional manner that measures traits at a single point in time. To address this issue, we have traced blood profiles in 4,884 Japanese individuals who underwent annual health check-ups for several years. In the present study, longitudinal exome-wide association studies were conducted to identify genetic variants related to 13 hematological phenotypes. The generalized estimating equation model showed that a total of 67 single nucleotide polymorphisms (SNPs) were significantly [false discovery rate (FDR) of <0.01] associated with hematological phenotypes. Of the 67 SNPs, nine SNPs were identified as novel hematological markers: rs4686683 of SENP2 for red blood cell count (FDR = 0.008, P = 5.5 × 10⁻⁶); rs3917688 of SELP for mean corpuscular volume (FDR = 0.005, P = 2.4 × 10⁻⁶); rs3133745 of C8orf37-AS1 for white blood cell count (FDR = 0.003, P = 1.3 × 10⁻⁶); rs13121954 at 4q31.2 for basophil count (FDR = 0.007, P = 3.1 × 10⁻⁵); rs7584099 at 2q22.3 (FDR = 2.6 × 10⁻⁵, P = 8.8 × 10⁻⁸), rs1579219 of HCG17 (FDR = 0.003, P = 2.0 × 10⁻⁵), and rs10757049 of DENND4C (FDR = 0.008, P = 5.6 × 10⁻⁵) for eosinophil count; rs12338 of CTSB for neutrophil count (FDR = 0.007, P = 2.9 × 10⁻⁵); and rs395967 of OSMR-AS1 for monocyte count (FDR = 0.008, P = 3.2 × 10⁻⁵).

Dichotomous function of IL-33 in health and disease: From biology to clinical implications

Interleukin (IL)-33 is a cytokine that is released from epithelial and endothelial cells at barrier surfaces upon tissue stress or damage to operate as an alarmin. IL-33 has been primarily implicated in the induction of T helper (Th) 2 type immune responses. Therefore, IL-33 has attracted a lot of interest as a potential therapeutic target in asthma and other allergic diseases. Over the years, it has become clear that IL-33 has a much broader activity and also contributes to Th1 immunity, expanding the possibilities for therapeutic modulation of IL-33 activity to multiple inflammatory diseases. However, more recently IL-33 has also been shown to mediate immunosuppression and tissue repair by activating regulatory T cells (Treg) and promoting M2 macrophage polarization. These pleiotropic activities of IL-33 illustrate the need for a tight molecular regulation of IL-33 activity, and have to be taken into account when IL-33 or its receptor is targeted for therapeutic modulation. Here we review the multiple molecular mechanisms that regulate IL-33 activity and describe how IL-33 can shape innate and adaptive immune responses by promoting Th1, Th2 and Treg function. Finally, we will discuss the possibilities for therapeutic modulation of IL-33 signaling as well as possible safety issues.