Allele-specific quantification of human leukocyte antigen transcript isoforms by nanopore sequencing

Introduction

While tens of thousands of HLA alleles have been identified by DNA sequencing, the contribution of alternative splicing to HLA diversity is not well characterized. In this study, we sought to determine if long-read sequencing could be used to accurately quantify allele-specific HLA transcripts in primary human lymphocytes.

Methods

cDNA libraries were prepared from peripheral blood lymphocytes from 12 donors and sequenced by nanopore long-read sequencing. HLA reads were aligned to donor-specific reference sequences based on the known type of each donor. Allele-specific exon utilization was calculated as the proportion of reads aligning to each allele containing known exons, and transcript isoforms were quantified based on patterns of exon utilization within individual reads.

Results

Splice variants were rare among class I HLA genes (median exon retention rate 99%-100%), except for several HLA-C alleles with exon 5 spliced out of up to 15% of reads. Splice variants were also rare among class II HLA genes (median exon retention rate 98%-100%), except for HLA-DQB1. Consistent with previous work, exon 5 of HLA-DQB1 was spliced out in alleles with a mutated splice acceptor site at rs28688207. Surprisingly, a 28% loss of exon 5 was also observed in HLA-DQB1 alleles with an intact splice acceptor site at rs28688207.

Discussion

We describe a simple bioinformatic workflow to quantify allele-specific expression of HLA transcript isoforms. Further studies are warranted to characterize the repertoire of HLA transcripts expressed in different cell types and tissues across diverse populations.

Implementation of High-Sensitivity Troponin with a Rapid Diagnostic Algorithm Reduces Emergency Department Length of Stay for Discharged Patients

INTRODUCTION

High sensitivity troponin (hs-cTn) and diagnostic algorithms are used to rapidly triage patients with symptoms of acute myocardial infarction in emergency departments (ED). However, few studies have evaluated the impact of simultaneously implementing hs-cTn and a rapid rule-out algorithm on length of stay (LOS).

METHODS

We assessed the impact of transitioning from contemporary cTnI to hs-cTnI in 59,232 ED encounters over three years. hs-cTnI was implemented with an orderable series that included baseline, two-, four-, and six-hour specimens collected at provider discretion and operationalized with an algorithm to calculate the change in hs-cTnI from baseline and provide interpretations of "insignificant", "significant," or "equivocal." Patient demographics, results, chief complaint, disposition, and ED LOS were captured from the electronic medical record.

RESULTS

cTnI was ordered for 31,875 encounters prior to hs-cTnI implementation and 27,357 after. The proportion of cTnI results above the 99th percentile upper reference limit decreased from 35.0% to 27.0% for men and increased from 27.8% to 34.8% for women. Among discharged patients, the median LOS decreased by 0.6 h (0.5-0.7). LOS among discharged patients with a chief complaint of chest pain decreased by 1.0 h (0.8-1.1) and further decreased by 1.2 h (1.0-1.3) if the initial hs-cTnI was below the limit of quantitation. The rate of acute coronary syndrome upon re-presentation within 30 days did not change post-implementation (0.10% versus 0.07%).

CONCLUSION

Implementation of an hs-cTnI assay with a rapid rule-out algorithm decreased ED LOS among discharged patients, particularly among those with a chief complaint of chest pain.

Clinical Decision Support for Laboratory Testing

BACKGROUND

As technology enables new and increasingly complex laboratory tests, test utilization presents a growing challenge for healthcare systems. Clinical decision support (CDS) refers to digital tools that present providers with clinically relevant information and recommendations, which have been shown to improve test utilization. Nevertheless, individual CDS applications often fail, and implementation remains challenging.

CONTENT

We review common classes of CDS tools grounded in examples from the literature as well as our own institutional experience. In addition, we present a practical framework and specific recommendations for effective CDS implementation.

SUMMARY

CDS encompasses a rich set of tools that have the potential to drive significant improvements in laboratory testing, especially with respect to test utilization. Deploying CDS effectively requires thoughtful design and careful maintenance, and structured processes focused on quality improvement and change management play an important role in achieving these goals.

Genome Sequencing as an Alternative to Cytogenetic Analysis in Myeloid Cancers

BACKGROUND

Genomic analysis is essential for risk stratification in patients with acute myeloid leukemia (AML) or myelodysplastic syndromes (MDS). Whole-genome sequencing is a potential replacement for conventional cytogenetic and sequencing approaches, but its accuracy, feasibility, and clinical utility have not been demonstrated.

METHODS

We used a streamlined whole-genome sequencing approach to obtain genomic profiles for 263 patients with myeloid cancers, including 235 patients who had undergone successful cytogenetic analysis. We adapted sample preparation, sequencing, and analysis to detect mutations for risk stratification using existing European Leukemia Network (ELN) guidelines and to minimize turnaround time. We analyzed the performance of whole-genome sequencing by comparing our results with findings from cytogenetic analysis and targeted sequencing.

RESULTS

Whole-genome sequencing detected all 40 recurrent translocations and 91 copy-number alterations that had been identified by cytogenetic analysis. In addition, we identified new clinically reportable genomic events in 40 of 235 patients (17.0%). Prospective sequencing of samples obtained from 117 consecutive patients was performed in a median of 5 days and provided new genetic information in 29 patients (24.8%), which changed the risk category for 19 patients (16.2%). Standard AML risk groups, as defined by sequencing results instead of cytogenetic analysis, correlated with clinical outcomes. Whole-genome sequencing was also used to stratify patients who had inconclusive results by cytogenetic analysis into risk groups in which clinical outcomes were measurably different.

CONCLUSIONS

In our study, we found that whole-genome sequencing provided rapid and accurate genomic profiling in patients with AML or MDS. Such sequencing also provided a greater diagnostic yield than conventional cytogenetic analysis and more efficient risk stratification on the basis of standard risk categories. (Funded by the Siteman Cancer Research Fund and others.).

Enhancer transcription identifies cis-regulatory elements for photoreceptor cell types

Identification of cell type-specific cis-regulatory elements (CREs) is crucial for understanding development and disease, although identification of functional regulatory elements remains challenging. We hypothesized that context-specific CREs could be identified by context-specific non-coding RNA (ncRNA) profiling, based on the observation that active CREs produce ncRNAs. We applied ncRNA profiling to identify rod and cone photoreceptor CREs from wild-type and mutant mouse retinas, defined by presence or absence, respectively, of the rod-specific transcription factor (TF) NrlNrl-dependent ncRNA expression strongly correlated with epigenetic profiles of rod and cone photoreceptors, identified thousands of candidate rod- and cone-specific CREs, and identified motifs for rod- and cone-specific TFs. Colocalization of NRL and the retinal TF CRX correlated with rod-specific ncRNA expression, whereas CRX alone favored cone-specific ncRNA expression, providing quantitative evidence that heterotypic TF interactions distinguish cell type-specific CRE activity. We validated the activity of novel Nrl-dependent ncRNA-defined CREs in developing cones. This work supports differential ncRNA profiling as a platform for the identification of cell type-specific CREs and the discovery of molecular mechanisms underlying TF-dependent CRE activity.

A massively parallel reporter assay reveals context-dependent activity of homeodomain binding sites in vivo

Cone-rod homeobox (CRX) is a paired-like homeodomain transcription factor (TF) and a master regulator of photoreceptor development in vertebrates. The in vitro DNA binding preferences of CRX have been described in detail, but the degree to which in vitro binding affinity is correlated with in vivo enhancer activity is not known. In addition, paired-class homeodomain TFs can bind DNA cooperatively as both homodimers and heterodimers at inverted TAAT half-sites separated by 2 or 3 nucleotides. This dimeric configuration is thought to mediate target specificity, but whether monomeric and dimeric sites encode distinct levels of activity is not known. Here, we used a massively parallel reporter assay to determine how local sequence context shapes the regulatory activity of CRX binding sites in mouse photoreceptors. We assayed inactivating mutations in more than 1700 TF binding sites and found that dimeric CRX binding sites act as stronger enhancers than monomeric CRX binding sites. Furthermore, the activity of dimeric half-sites is cooperative, dependent on a strict 3-bp spacing, and tuned by the identity of the spacer nucleotides. Saturating single-nucleotide mutagenesis of 195 CRX binding sites showed that, on average, changes in TF binding site affinity are correlated with changes in regulatory activity, but this relationship is obscured when considering mutations across multiple cis-regulatory elements (CREs). Taken together, these results demonstrate that the activity of CRX binding sites is highly dependent on sequence context, providing insight into photoreceptor gene regulation and illustrating functional principles of homeodomain binding sites that may be conserved in other cell types.

Massively parallel cis-regulatory analysis in the mammalian central nervous system

Cis-regulatory elements (CREs, e.g., promoters and enhancers) regulate gene expression, and variants within CREs can modulate disease risk. Next-generation sequencing has enabled the rapid generation of genomic data that predict the locations of CREs, but a bottleneck lies in functionally interpreting these data. To address this issue, massively parallel reporter assays (MPRAs) have emerged, in which barcoded reporter libraries are introduced into cells, and the resulting barcoded transcripts are quantified by next-generation sequencing. Thus far, MPRAs have been largely restricted to assaying short CREs in a limited repertoire of cultured cell types. Here, we present two advances that extend the biological relevance and applicability of MPRAs. First, we adapt exome capture technology to instead capture candidate CREs, thereby tiling across the targeted regions and markedly increasing the length of CREs that can be readily assayed. Second, we package the library into adeno-associated virus (AAV), thereby allowing delivery to target organs in vivo. As a proof of concept, we introduce a capture library of about 46,000 constructs, corresponding to roughly 3500 DNase I hypersensitive (DHS) sites, into the mouse retina by ex vivo plasmid electroporation and into the mouse cerebral cortex by in vivo AAV injection. We demonstrate tissue-specific cis-regulatory activity of DHSs and provide examples of high-resolution truncation mutation analysis for multiplex parsing of CREs. Our approach should enable massively parallel functional analysis of a wide range of CREs in any organ or species that can be infected by AAV, such as nonhuman primates and human stem cell–derived organoids.

Clonal architecture of secondary acute myeloid leukemia defined by single-cell sequencing

Next-generation sequencing has been used to infer the clonality of heterogeneous tumor samples. These analyses yield specific predictions—the population frequency of individual clones, their genetic composition, and their evolutionary relationships—which we set out to test by sequencing individual cells from three subjects diagnosed with secondary acute myeloid leukemia, each of whom had been previously characterized by whole genome sequencing of unfractionated tumor samples. Single-cell mutation profiling strongly supported the clonal architecture implied by the analysis of bulk material. In addition, it resolved the clonal assignment of single nucleotide variants that had been initially ambiguous and identified areas of previously unappreciated complexity. Accordingly, we find that many of the key assumptions underlying the analysis of tumor clonality by deep sequencing of unfractionated material are valid. Furthermore, we illustrate a single-cell sequencing strategy for interrogating the clonal relationships among known variants that is cost-effective, scalable, and adaptable to the analysis of both hematopoietic and solid tumors, or any heterogeneous population of cells.

Human cancers are genetically diverse populations of cells that evolve over the course of their natural history or in response to the selective pressure of therapy. In theory, it is possible to infer how this variation is structured into related populations of cells based on the frequency of individual mutations in bulk samples, but the accuracy of these models has not been evaluated across a large number of variants in individual cells. Here, we report a strategy for analyzing hundreds of variants within a single cell, and we apply this method to assess models of tumor clonality derived from bulk samples in three cases of leukemia. The data largely support the predicted population structure, though they suggest specific refinements. This type of approach not only illustrates the biological complexity of human cancer, but it also has the potential to inform patient management. That is, precise knowledge of which variants are present in which populations of cells may allow physicians to more effectively target combinations of mutations and predict how patients will respond to therapy.

Excess congenital non-synonymous variation in leukemia-associated genes in MLL- infant leukemia: a Children's Oncology Group report

Infant leukemia (IL) is a rare sporadic cancer with a grim prognosis. Although most cases are accompanied by MLL rearrangements and harbor very few somatic mutations, less is known about the genetics of the cases without MLL translocations. We performed the largest exome-sequencing study to date on matched non-cancer DNA from pairs of mothers and IL patients to characterize congenital variation that may contribute to early leukemogenesis. Using the COSMIC database to define acute leukemia-associated candidate genes, we find a significant enrichment of rare, potentially functional congenital variation in IL patients compared with randomly selected genes within the same patients and unaffected pediatric controls. IL acute myeloid leukemia (AML) patients had more overall variation than IL acute lymphocytic leukemia (ALL) patients, but less of that variation was inherited from mothers. Of our candidate genes, we found that MLL3 was a compound heterozygote in every infant who developed AML and 50% of infants who developed ALL. These data suggest a model by which known genetic mechanisms for leukemogenesis could be disrupted without an abundance of somatic mutation or chromosomal rearrangements. This model would be consistent with existing models for the establishment of leukemia clones in utero and the high rate of IL concordance in monozygotic twins.

Genetic variation in MKL2 and decreased downstream PCTAIRE1 expression in extreme, fatal primary human microcephaly

The genetic mechanisms driving normal brain development remain largely unknown. We performed genomic and immunohistochemical characterization of a novel, fatal human phenotype including extreme microcephaly with cerebral growth arrest at 14-18 weeks gestation in three full sisters born to healthy, non-consanguineous parents. Analysis of index cases and parents included familial exome sequencing, karyotyping, and genome-wide single nucleotide polymorphism (SNP) array. From proband, control and unrelated microcephalic fetal cortical tissue, we compared gene expression of RNA and targeted immunohistochemistry. Each daughter was homozygous for a rare, non-synonymous, deleterious variant in the MKL2 gene and heterozygous for a private 185 kb deletion on the paternal allele, upstream and in cis with his MKL2 variant allele, eliminating 24 CArG transcription factor binding sites and MIR4718. MKL1 was underexpressed in probands. Dysfunction of MKL2 and its transcriptional coactivation partner, serum response factor (SRF), was supported by a decrease in gene and protein expression of PCTAIRE1, a downstream target of MKL2:SRF heterodimer transcriptional activation, previously shown to result in severe microcephaly in murine models. While disruption of the MKL2:SRF axis has been associated with severe microcephaly and disordered brain development in multiple model systems, the role of this transcription factor complex has not been previously demonstrated in human brain development.