High-coverage whole-genome sequencing of the expanded 1000 Genomes Project cohort including 602 trios.

The 1000 Genomes Project (1kGP) is the largest fully open resource of whole-genome sequencing (WGS) data consented for public distribution without access or use restrictions. The final, phase 3 release of the 1kGP included 2,504 unrelated samples from 26 populations and was based primarily on low-coverage WGS. Here, we present a high-coverage 3,202-sample WGS 1kGP resource, which now includes 602 complete trios, sequenced to a depth of 30X using Illumina. We performed single-nucleotide variant (SNV) and short insertion and deletion (INDEL) discovery and generated a comprehensive set of structural variants (SVs) by integrating multiple analytic methods through a machine learning model. We show gains in sensitivity and precision of variant calls compared to phase 3, especially among rare SNVs as well as INDELs and SVs spanning frequency spectrum. We also generated an improved reference imputation panel, making variants discovered here accessible for association studies.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Mapping the 17q12-21.1 Locus for Variants Associated with Early-Onset Asthma in African Americans.

Rationale: The 17q12-21.1 locus is one of the most highly replicated genetic associations with asthma. Individuals of African descent have lower linkage disequilibrium in this region, which could facilitate identifying causal variants.Objectives: To identify functional variants at 17q12-21.1 associated with early-onset asthma among African American individuals.Methods: We evaluated African American participants from SAPPHIRE (Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity) (n = 1,940), SAGE II (Study of African Americans, Asthma, Genes and Environment) (n = 885), and GCPD-A (Study of the Genetic Causes of Complex Pediatric Disorders-Asthma) (n = 2,805). Associations with asthma onset at ages under 5 years were meta-analyzed across cohorts. The lead signal was reevaluated considering haplotypes informed by genetic ancestry (i.e., African vs. European). Both an expression-quantitative trait locus analysis and a phenome-wide association study were performed on the lead variant.Measurements and Main Results: The meta-analyzed results from SAPPHIRE, SAGE II, and the GCPD-A identified rs11078928 as the top association for early-onset asthma. A haplotype analysis suggested that the asthma association partitioned most closely with the rs11078928 genotype. Genetic ancestry did not appear to influence the effect of this variant. In the expression-quantitative trait locus analysis, rs11078928 was related to alternative splicing of GSDMB (gasdermin-B) transcripts. The phenome-wide association study of rs11078928 suggested that this variant was predominantly associated with asthma and asthma-associated symptoms.Conclusions: A splice-acceptor polymorphism appears to be a causal variant for asthma at the 17q12-21.1 locus. This variant appears to have the same magnitude of effect in individuals of African and European descent.

Whole-genome bisulfite sequencing with improved accuracy and cost.

DNA methylation patterns in the genome both reflect and help to mediate transcriptional regulatory processes. The digital nature of DNA methylation, present or absent on each allele, makes this assay capable of quantifying events in subpopulations of cells, whereas genome-wide chromatin studies lack the same quantitative capacity. Testing DNA methylation throughout the genome is possible using whole-genome bisulfite sequencing (WGBS), but the high costs associated with the assay have made it impractical for studies involving more than limited numbers of samples. We have optimized a new transposase-based library preparation assay for the Illumina HiSeq X platform suitable for limited amounts of DNA and providing a major cost reduction for WGBS. By incorporating methylated cytosines during fragment end repair, we reveal an end-repair artifact affecting 1%-2% of reads that we can remove analytically. We show that the use of a high (G + C) content spike-in performs better than PhiX in terms of bisulfite sequencing quality. As expected, the loci with transposase-accessible chromatin are DNA hypomethylated and enriched in flanking regions by post-translational modifications of histones usually associated with positive effects on gene expression. Using these transposase-accessible loci to represent the cis-regulatory loci in the genome, we compared the representation of these loci between WGBS and other genome-wide DNA methylation assays, showing WGBS to outperform substantially all of the alternatives. We conclude that it is now technologically and financially feasible to perform WGBS in larger numbers of samples with greater accuracy than previously possible.

Whole-Genome Sequencing Identifies Novel Functional Loci Associated with Lung Function in Puerto Rican Youth.

Rationale: Puerto Ricans have the highest childhood asthma prevalence in the United States (23.6%); however, the etiology is uncertain.Objectives: In this study, we sought to uncover the genetic architecture of lung function in Puerto Rican youth with and without asthma who were recruited from the island (n = 836).Methods: We used admixture-mapping and whole-genome sequencing data to discover genomic regions associated with lung function. Functional roles of the prioritized candidate SNPs were examined with chromatin immunoprecipitation sequencing, RNA sequencing, and expression quantitative trait loci data.Measurements and Main Results: We discovered a genomic region at 1q32 that was significantly associated with a 0.12-L decrease in the lung volume of exhaled air (95% confidence interval, -0.17 to -0.07; P = 6.62 × 10-8) with each allele of African ancestry. Within this region, two SNPs were expression quantitative trait loci of TMEM9 in nasal airway epithelial cells and MROH3P in esophagus mucosa. The minor alleles of these SNPs were associated with significantly decreased lung function and decreased TMEM9 gene expression. Another admixture-mapping peak was observed on chromosome 5q35.1, indicating that each Native American ancestry allele was associated with a 0.15-L increase in lung function (95% confidence interval, 0.08-0.21; P = 5.03 × 10-6). The region-based association tests identified four suggestive windows that harbored candidate rare variants associated with lung function.Conclusions: We identified common and rare genetic variants that may play a critical role in lung function among Puerto Rican youth. We independently validated an inflammatory pathway that could potentially be used to develop more targeted treatments and interventions for patients with asthma.

YES1 amplification is a mechanism of acquired resistance to EGFR inhibitors identified by transposon mutagenesis and clinical genomics.

In ∼30% of patients with EGFR-mutant lung adenocarcinomas whose disease progresses on EGFR inhibitors, the basis for acquired resistance remains unclear. We have integrated transposon mutagenesis screening in an EGFR-mutant cell line and clinical genomic sequencing in cases of acquired resistance to identify mechanisms of resistance to EGFR inhibitors. The most prominent candidate genes identified by insertions in or near the genes during the screen were MET, a gene whose amplification is known to mediate resistance to EGFR inhibitors, and the gene encoding the Src family kinase YES1. Cell clones with transposon insertions that activated expression of YES1 exhibited resistance to all three generations of EGFR inhibitors and sensitivity to pharmacologic and siRNA-mediated inhibition of YES1 Analysis of clinical genomic sequencing data from cases of acquired resistance to EGFR inhibitors revealed amplification of YES1 in five cases, four of which lacked any other known mechanisms of resistance. Preinhibitor samples, available for two of the five patients, lacked YES1 amplification. None of 136 postinhibitor samples had detectable amplification of other Src family kinases (SRC and FYN). YES1 amplification was also found in 2 of 17 samples from ALK fusion-positive lung cancer patients who had progressed on ALK TKIs. Taken together, our findings identify acquired amplification of YES1 as a recurrent and targetable mechanism of resistance to EGFR inhibition in EGFR-mutant lung cancers and demonstrate the utility of transposon mutagenesis in discovering clinically relevant mechanisms of drug resistance.

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma.

RATIONALE: Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. OBJECTIVES: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. METHODS: We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. MEASUREMENTS AND MAIN RESULTS: We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and ß-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings. CONCLUSIONS: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

Differential burden of rare protein truncating variants in Alzheimer's disease patients compared to centenarians.

We compared coding region variants of 53 cognitively healthy centenarians and 45 patients with Alzheimer's disease (AD), all of Ashkenazi Jewish (AJ) ancestry. Despite the small sample size, the known AD risk variant APOE4 reached genome-wide significance, indicating the advantage of utilizing 'super-controls'. We restricted our subsequent analysis to rare variants observed at most once in the 1000 Genomes database and having a minor allele frequency below 2% in our AJ sample. We compared the burden of predicted protein altering variants between cases and controls as normalized by the level of rare synonymous variants. We observed an increased burden among AD subjects for predicted loss-of-function (LoFs) variants defined as stop-gain, frame shift, initiation codon (INIT) and splice site mutations (n = 930, OR = 1.3, P = 1.5×E-5). There was no enrichment across all rare protein altering variants defined as missense plus LoFs, in frame indels and stop-loss variants (n = 13 014, OR = 0.97, P = 0.47). Among LoFs, the strongest burden was observed for INIT (OR = 2.16, P = 0.0097) and premature stop variants predicted to cause non-sense-mediated decay in the majority of transcripts (NMD) (OR = 1.98, P = 0.02). Notably, this increased burden of NMD, INIT and splice variants was more pronounced in a set of 1397 innate immune genes (OR = 4.55, P = 0.0043). Further comparison to additional exomes indicates that the difference in LoF burden originated both from the AD and centenarian sample. In summary, we observed an overall increased burden of rare LoFs in AD subjects as compared to centenarians, and this enrichment is more pronounced for innate immune genes.

Effect of IL28B genotype on early viral kinetics during interferon-free treatment of patients with chronic hepatitis C.

BACKGROUND & AIMS: Although interleukin 28B (interferon, lambda 3) (IL28B) genotype affects the response of patients with chronic hepatitis C to peginterferon and ribavirin, little is known regarding its effect on response to direct-acting antivirals in interferon-free combinations. We analyzed the effects of IL28B genotype on the viral kinetic (VK) response to an interferon-free combination of the nucleoside polymerase inhibitor mericitabine (RG7128) and the hepatitis C virus (HCV) protease inhibitor danoprevir. METHODS: We performed a double-blind, dose-escalation study of patients with chronic HCV genotype 1 infection who were interferon treatment naive or had not responded to previous therapy with peginterferon and ribavirin. Patients were sequentially assigned to 1 of 7 cohorts then randomly assigned to groups that received up to 13 days of treatment with mericitabine (500 or 1000 mg, twice daily) plus danoprevir (100 or 200 mg, every 8 hours, or 600 or 900 mg, twice daily) or placebo. Eighty-three of 87 patients were genotyped for the IL28B single-nucleotide polymorphism rs12979860. VKs were analyzed only in patients who received 13 days of treatment, at optimal doses, using a biphasic model to describe first- and second-phase slopes of viral decay during therapy. RESULTS: At day 14 (the end of interferon-free treatment), the mean reduction in the serum level of HCV RNA was slightly greater in patients with the CC polymorphism (5.01 log(10) IU/mL) than those without (4.59 log(10) IU/mL). Modeling revealed that patients with the CC polymorphism had slightly better early VKs, most apparent in the ß-phase of viral decay. A mixed effect on the α-phase was observed, which was reduced in magnitude but prolonged in patients with CC, who also had better on-treatment response to peginterferon and ribavirin during follow up. CONCLUSIONS: IL28B genotype appears to affect early VKs in patients with chronic hepatitis C receiving interferon-free treatment.

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants.

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-seq data in GTEx v8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased WGS data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Genetic control of mRNA splicing as a potential mechanism for incomplete penetrance of rare coding variants.

Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

Investigation of single nucleotide polymorphisms and biological pathways associated with response to TNFα inhibitors in patients with rheumatoid arthritis.

OBJECTIVE: Recently, two genome-wide association studies identified single nucleotide polymorphisms (SNPs) significantly associated with the treatment response to tumor necrosis factor α (TNFα) inhibitors in patients with rheumatoid arthritis (RA). We aimed to replicate these results and identify SNPs and the possible biological pathways associated with the treatment response to TNFα inhibitors. METHODS: TNFα-naive patients with RA, who had available DNA and initiated TNFα inhibitor therapy between 1999 and 2008, were identified in the DANBIO registry and genotyped using the Illumina HumanHap550K Duo array. The associations between SNPs and changes in the absolute and the relative Disease Activity Score, and European League Against Rheumatism good versus no response after 14 weeks of treatment were tested. SNP data were combined with two independent cohorts in a meta-analysis. A gene-set enrichment analysis (GSEA) was carried out to identify the biological pathways associated with the treatment response. RESULTS: After genotyping and quality control, 486 450 SNPs were analyzed in 196 Danish patients with moderate to severe RA treated with infliximab (n=142), etanercept (n=12), and adalimumab (n=42). None of the previously identified SNPs were confirmed in our dataset or in meta-analyses of available studies. Other potential SNPs were identified, but none achieved genome-wide significance. A GSEA identified the transforming growth factor ß, TNF, mitogen-activated protein kinase, and mammalian target of rapamycin pathways to have a potential influence on the treatment response. CONCLUSION: In a genome-wide association study of 196 genetically homogenous Danish patients with RA and in a meta-analysis, we found no SNPs associated with treatment response to TNFα inhibitors. A GSEA suggested that the transforming growth factor ß, TNF, mitogen-activated protein kinase, and mammalian target of rapamycin pathways may be associated with treatment response.

Accounting for population structure in genetic studies of cystic fibrosis.

CFTR F508del (c.1521_1523delCTT, p.Phe508delPhe) is the most common pathogenic allele underlying cystic fibrosis (CF), and its frequency varies in a geographic cline across Europe. We hypothesized that genetic variation associated with this cline is overrepresented in a large cohort (N > 5,000) of persons with CF who underwent whole-genome sequencing and that this pattern could result in spurious associations between variants correlated with both the F508del genotype and CF-related outcomes. Using principal-component (PC) analyses, we showed that variation in the CFTR region disproportionately contributes to a PC explaining a relatively high proportion of genetic variance. Variation near CFTR was correlated with population structure among persons with CF, and this correlation was driven by a subset of the sample inferred to have European ancestry. We performed genome-wide association studies comparing persons with CF with one versus two copies of the F508del allele; this allowed us to identify genetic variation associated with the F508del allele and to determine that standard PC-adjustment strategies eliminated the significant association signals. Our results suggest that PC adjustment can adequately prevent spurious associations between genetic variants and CF-related traits and are therefore effective tools to control for population structure even when population structure is confounded with disease severity and a common pathogenic variant.

Genomic Epidemiology and Serology Associated with a SARS-CoV-2 R.1 Variant Outbreak in New Jersey.

Examining the neutralizing capacity of monoclonal antibodies (MAbs) used to treat COVID-19, as well as antibodies recovered from unvaccinated, previously vaccinated, and infected individuals, against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) remains critical to study. Here, we report on a SARS-CoV-2 nosocomial outbreak caused by the SARS-CoV-2 R.1 variant harboring the E484K mutation in a 281-bed psychiatric facility in New Jersey among unvaccinated inpatients and health care professionals (HCPs). A total of 81 inpatients and HCPs tested positive for SARS-Cov-2 by reverse transcription (RT)-PCR from 29 October 9 to 30 November 2020. The R.1 variant exhibits partial or complete resistance to two MAbs in clinical use, as well as 2 receptor binding domain MAbs and 4 N-terminal domain (NTD) MAbs. NTD MAbs against pseudovirus harboring single characteristic R.1 mutations highlight the role of S255F in loss of activity. Additionally, we note dampened neutralization capacity by plasma from individuals with previous SARS-CoV-2 infection or sera from vaccinated individuals. The relative resistance of the R.1 variant is likely lower than that of B.1.351 and closer to that of P.1 and B.1.526. The R.1 lineage has been reported in 47 states in the United States and 40 countries. Although high proportions exhibited symptoms (26% and 61% among patients and HCPs, respectively) and relative antibody resistance, we detected only 10 R.1 variants from over 2,900 samples (~0.34%) collected from January to October 2021. Among 3 vaccinated individuals previously infected with R.1, we observed robust neutralizing antibody responses against SARS-CoV-2 wild type and VOCs. IMPORTANCE The neutralizing capacities of monoclonal antibodies used to treat COVID-19 and of those recovered from previously infected and vaccinated individuals against SARS-CoV-2 variants of concern (VOCs) remain important questions. We report on a nosocomial outbreak caused by a SARS-CoV-2 R.1 variant harboring an E484K mutation among 81 unvaccinated inpatients and health care professionals. We note high attack rates with symptoms in nearly 50% of infected individuals, in sharp contrast to an unrelated institutional outbreak caused by the R.1 variant among a vaccinated population. We found little evidence of significant community spillover. This variant exhibits partial or complete resistance to two monoclonal antibodies in clinical use and dampened the neutralization capacity of convalescent-phase plasma from individuals with previous infection or sera from vaccinated individuals. Among three vaccinated individuals previously infected with R.1, we observed robust neutralizing antibody responses against SARS-CoV-2 wild type and VOCs. These findings underscore the importance of vaccination for prevention of symptomatic COVID-19 disease.

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.

System-wide transcriptome damage and tissue identity loss in COVID-19 patients.

The molecular mechanisms underlying the clinical manifestations of coronavirus disease 2019 (COVID-19), and what distinguishes them from common seasonal influenza virus and other lung injury states such as acute respiratory distress syndrome, remain poorly understood. To address these challenges, we combine transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues to define body-wide transcriptome changes in response to COVID-19. We then match these data with spatial protein and expression profiling across 357 tissue sections from 16 representative patient lung samples and identify tissue-compartment-specific damage wrought by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, evident as a function of varying viral loads during the clinical course of infection and tissue-type-specific expression states. Overall, our findings reveal a systemic disruption of canonical cellular and transcriptional pathways across all tissues, which can inform subsequent studies to combat the mortality of COVID-19 and to better understand the molecular dynamics of lethal SARS-CoV-2 and other respiratory infections.

Genetic variation in UGT1A1 typical of Gilbert syndrome is associated with unconjugated hyperbilirubinemia in patients receiving tocilizumab.

OBJECTIVE: Tocilizumab, a monoclonal antibody to interleukin-6 receptor, was recently approved for the treatment of moderate-to-severe rheumatoid arthritis. Two patients during clinical development met laboratory, but not clinical, criteria for Hy's law with bilirubin elevations suspected as a result of genetic variation in uridine diphosphoglucose glucuronosyltransferase (UGT1A1) typical of Gilbert syndrome. METHODS: Genotyping of the two cases potentially meeting with Hy's law was performed using commercially available procedures. UGT1A1 single nucleotide polymorphism data were extracted from a genome-wide array database for 1187 patients from tocilizumab trials, and associations of UGT1A1 genotypes with bilirubin elevations were analyzed using logistic regression for associations with baseline and change from baseline in bilirubin levels as continuous variables. RESULTS: Bilirubin elevations were not associated with clinical adverse events. Both patients potentially meeting Hy's law carry homozygous UGT1A1*28 alleles and UGT1A1*60 alleles. UGT1A1*28 and three additional single nucleotide polymorphisms showed odds ratios greater than 25 for associations with elevated bilirubin. The presence of rs6742078 accounted for 32% of the total variance in bilirubin (P=2.2×10). CONCLUSION: Bilirubin increases occurring with tocilizumab appear to be related to anti-inflammatory effects extending to the liver. Thus, in the absence of other signs of hepatic dysfunction, bilirubin elevations after treatment with tocilizumab have a high probability of association with UGT1A1 polymorphism, which should alleviate concerns of serious hepatotoxicity. Our results underscore the value of genotyping in the clinical trial setting to avoid misinterpretations that could lead to terminating development of a promising new agent.

Postvaccination SARS-COV-2 among Health Care Workers in New Jersey: A Genomic Epidemiological Study.

Emergence of SARS-CoV-2 with high transmission and immune evasion potential, the so-called variants of concern (VOC), is a major concern. We describe the early genomic epidemiology of SARS-CoV-2 recovered from vaccinated health care professionals (HCP). Our postvaccination COVID-19 symptoms-based surveillance program among HCPs in a 17-hospital network identified all vaccinated HCPs who tested positive for COVID-19 after routine screening or after self-reporting. From 1 January 2021 to 30 April 2021, 23,687 HCPs received either mRNA-1273 or BNT162b2 mRNA vaccine. All available postvaccination SARS-CoV-2 samples and a random collection from nonvaccinated patients during the similar time frame were subjected to VOC screening and whole-genome sequencing (WGS). Sixty-two percent (23,697/37,500) of HCPs received at least one vaccine dose, with 60% (22,458) fully vaccinated. We detected 138 (0.58%, 138/23,697) COVID-19 cases, 105 among partially vaccinated and 33 (0.15%, 33/22,458) among fully vaccinated. Five partially vaccinated required hospitalization, four with supplemental oxygen. VOC screening from 16 fully vaccinated HCPs identified 6 (38%) harboring N501Y and 1 (6%) with E484K polymorphisms; percentage of concurrent nonvaccinated samples was 37% (523/1,404) and 20% (284/1,394), respectively. There was an upward trend from January to April for E484K/Q (3% to 26%) and N501Y (1% to 49%). WGS analysis from vaccinated and nonvaccinated individuals indicated highly congruent phylogenies. We did not detect an increased frequency of any receptor-binding domain (RBD)/N-terminal domain (NTD) polymorphism between groups (P > 0.05). Our results support robust protection by vaccination, particularly among recipients of both doses. Despite VOCs accounting for over 40% of SARS-CoV-2 from fully vaccinated individuals, the genomic diversity appears to proportionally represent VOCs among nonvaccinated populations. IMPORTANCE A number of highly effective vaccines have been developed and deployed to combat the COVID-19 pandemic. The emergence and epidemiological dominance of SARS-CoV-2 mutants with high transmission potential and immune evasion properties, the so-called variants of concern (VOC), continue to be a major concern. Whether these VOCs alter the efficacy of the administered vaccines is of great concern and a critical question to study. We describe the initial genomic epidemiology of SARS-CoV-2 recovered from partial/fully vaccinated health care professionals and probe specifically for VOC enrichment. Our findings support the high level of protection provided by full vaccination despite a steep increase in the prevalence of polymorphisms associated with increased transmission potential (N501Y) and immune evasion (E484K) in the nonvaccinated population. Thus, we do not find evidence of VOC enrichment among vaccinated groups. Overall, the genomic diversity of SARS-CoV-2 recovered postvaccination appears to proportionally represent the observed viral diversity within the community.

Emergence of Multiple SARS-CoV-2 Antibody Escape Variants in an Immunocompromised Host Undergoing Convalescent Plasma Treatment.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs), harboring spike protein N-terminal domain (NTD) or receptor-binding domain (RBD) mutations, exhibit reduced in vitro susceptibility to convalescent-phase serum, commercial antibody cocktails, and vaccine neutralization and have been associated with reinfections. The accumulation of these mutations could be the consequence of intrahost viral evolution due to prolonged infection in immunocompromised hosts. In this study, we document the microevolution of SARS-CoV-2 recovered from sequential tracheal aspirates from an immunosuppressed patient on steroids and convalescent plasma therapy and identify the emergence of multiple NTD and RBD mutations. SARS-CoV-2 genomes from the first swab (day 0) and from three tracheal aspirates (days 7, 21, and 27) were compared at the sequence level. We identified a mixed viral population with five different S protein mutations (141 to 144 deletion, 243 to 244 deletion, E484K, Q493K, and Q493R) at the NTD or RBD region from the second tracheal aspirate sample (day 21) and a predominance of the S protein 141 to 144 LGVY deletion and E484K mutant on day 27. The neutralizing antibodies against various S protein lentiviral pseudovirus mutants, as well as the anti-SARS-CoV-2 total Ig and IgG, showed "U" shape dynamics, in support of the endogenous development of neutralizing antibodies. The patient's compromised immune status, the antirejection regiment, convalescent plasma treatment, and the development of neutralizing antibodies may have resulted in unique selective pressures on the intrahost genomic evolution, and this observation supports the hypotheses that VOCs can independently arise and that immunocompromised patients on convalescent plasma therapy are potential breeding grounds for immune escape mutants. IMPORTANCE Over a year of the COVID-19 pandemic, distinct severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages have arisen in multiple geographic areas around the world. SARS-CoV-2 variants of concern (VOCs), i.e., B.1.1.7 (alpha), B.1.351 (beta), P.1 (gamma), and B.1.617.2 (delta), harboring mutations and/or deletions in spike protein N-terminal domain (NTD) or receptor-binding domain (RBD) regions showed evidence of increased transmissibility and disease severity and possible reduced vaccine efficacy. In this study, we report the emergence of five different NTD and RBD mutations in an uncommon SARS-CoV-2 B.1.369 lineage from an immunosuppressed patient undergoing steroid and convalescent plasma therapy. The observation highlighted that VOCs can independently arise in immunocompromised populations undergoing anti-SARS-CoV-2 therapy, and enhanced measures will be required to reduce the transmission.