The Lineage-Defining Transcription Factors SOX2 and NKX2-1 Determine Lung Cancer Cell Fate and Shape the Tumor Immune Microenvironment.

The major types of non-small-cell lung cancer (NSCLC)-squamous cell carcinoma and adenocarcinoma-have distinct immune microenvironments. We developed a genetic model of squamous NSCLC on the basis of overexpression of the transcription factor Sox2, which specifies lung basal cell fate, and loss of the tumor suppressor Lkb1 (SL mice). SL tumors recapitulated gene-expression and immune-infiltrate features of human squamous NSCLC; such features included enrichment of tumor-associated neutrophils (TANs) and decreased expression of NKX2-1, a transcriptional regulator that specifies alveolar cell fate. In Kras-driven adenocarcinomas, mis-expression of Sox2 or loss of Nkx2-1 led to TAN recruitment. TAN recruitment involved SOX2-mediated production of the chemokine CXCL5. Deletion of Nkx2-1 in SL mice (SNL) revealed that NKX2-1 suppresses SOX2-driven squamous tumorigenesis by repressing adeno-to-squamous transdifferentiation. Depletion of TANs in SNL mice reduced squamous tumors, suggesting that TANs foster squamous cell fate. Thus, lineage-defining transcription factors determine the tumor immune microenvironment, which in turn might impact the nature of the tumor.

Accurate long-read sequencing allows assembly of the duplicated RHD and RHCE genes harboring variants relevant to blood transfusion.

Next-generation sequencing (NGS) technologies have transformed medical genetics. However, short-read lengths pose a limitation on identification of structural variants, sequencing repetitive regions, phasing of distant nucleotide changes, and distinguishing highly homologous genomic regions. Long-read sequencing technologies may offer improvements in the characterization of genes that are currently difficult to assess. We used a combination of targeted DNA capture, long-read sequencing, and a customized bioinformatics pipeline to fully assemble the RH region, which harbors variation relevant to red cell donor-recipient mismatch, particularly among patients with sickle cell disease. RHD and RHCE are a pair of duplicated genes located within an ∼175 kb region on human chromosome 1 that have high sequence similarity and frequent structural variations. To achieve the assembly, we utilized palindrome repeats in PacBio SMRT reads to obtain consensus sequences of 2.1 to 2.9 kb average length with over 99% accuracy. We used these long consensus sequences to identify 771 assembly markers and to phase the RHD-RHCE region with high confidence. The dataset enabled direct linkage between coding and intronic variants, phasing of distant SNPs to determine RHD-RHCE haplotypes, and identification of known and novel structural variations along with the breakpoints. A limiting factor in phasing is the frequency of heterozygous assembly markers and therefore was most successful in samples from African Black individuals with increased heterogeneity at the RH locus. Overall, this approach allows RH genotyping and de novo assembly in an unbiased and comprehensive manner that is necessary to expand application of NGS technology to high-resolution RH typing.

Transcriptome-wide profiling of multiple RNA modifications simultaneously at single-base resolution.

The breadth and importance of RNA modifications are growing rapidly as modified ribonucleotides can impact the sequence, structure, function, stability, and fate of RNAs and their interactions with other molecules. Therefore, knowing cellular RNA modifications at single-base resolution could provide important information regarding cell status and fate. A current major limitation is the lack of methods that allow the reproducible profiling of multiple modifications simultaneously, transcriptome-wide and at single-base resolution. Here we developed RBS-Seq, a modification of RNA bisulfite sequencing that enables the sensitive and simultaneous detection of m5C, Ψ, and m1A at single-base resolution transcriptome-wide. With RBS-Seq, m5C and m1A are accurately detected based on known signature base mismatches and are detected here simultaneously along with Ψ sites that show a 1-2 base deletion. Structural analyses revealed the mechanism underlying the deletion signature, which involves Ψ-monobisulfite adduction, heat-induced ribose ring opening, and Mg2+-assisted reorientation, causing base-skipping during cDNA synthesis. Detection of each of these modifications through a unique chemistry allows high-precision mapping of all three modifications within the same RNA molecule, enabling covariation studies. Application of RBS-Seq on HeLa RNA revealed almost all known m5C, m1A, and ψ sites in tRNAs and rRNAs and provided hundreds of new m5C and Ψ sites in noncoding RNAs and mRNAs. However, our results diverge greatly from earlier work, suggesting ∼10-fold fewer m5C sites in noncoding and coding RNAs and the absence of substantial m1A in mRNAs. Taken together, the approaches and refined datasets in this work will greatly enable future epitranscriptome studies.

Anti-inflammatory microRNA-146a protects mice from diet-induced metabolic disease.

Identifying regulatory mechanisms that influence inflammation in metabolic tissues is critical for developing novel metabolic disease treatments. Here, we investigated the role of microRNA-146a (miR-146a) during diet-induced obesity in mice. miR-146a is reduced in obese and type 2 diabetic patients and our results reveal that miR-146a-/- mice fed a high-fat diet (HFD) have exaggerated weight gain, increased adiposity, hepatosteatosis, and dysregulated blood glucose levels compared to wild-type controls. Pro-inflammatory genes and NF-κB activation increase in miR-146a-/- mice, indicating a role for this miRNA in regulating inflammatory pathways. RNA-sequencing of adipose tissue macrophages demonstrated a role for miR-146a in regulating both inflammation and cellular metabolism, including the mTOR pathway, during obesity. Further, we demonstrate that miR-146a regulates inflammation, cellular respiration and glycolysis in macrophages through a mechanism involving its direct target Traf6. Finally, we found that administration of rapamycin, an inhibitor of mTOR, was able to rescue the obesity phenotype in miR-146a-/- mice. Altogether, our study provides evidence that miR-146a represses inflammation and diet-induced obesity and regulates metabolic processes at the cellular and organismal levels, demonstrating how the combination of diet and miRNA genetics influences obesity and diabetic phenotypes.

AnthOligo: automating the design of oligonucleotides for capture/enrichment technologies.

SUMMARY: A number of methods have been devised to address the need for targeted genomic resequencing. One of these methods, region-specific extraction (RSE) is characterized by the capture of long DNA fragments (15-20 kb) by magnetic beads, after enzymatic extension of oligonucleotides hybridized to selected genomic regions. Facilitating the selection of the most appropriate capture oligos for targeting a region of interest, satisfying the properties of temperature (Tm) and entropy (ΔG), while minimizing the formation of primer-dimers in a pooled experiment, is therefore necessary. Manual design and selection of oligos becomes very challenging, complicated by factors such as length of the target region and number of targeted regions. Here we describe, AnthOligo, a web-based application developed to optimally automate the process of generation of oligo sequences used to target and capture the continuum of large and complex genomic regions. Apart from generating oligos for RSE, this program may have wider applications in the design of customizable internal oligos to be used as baits for gene panel analysis or even probes for large-scale comparative genomic hybridization array processes. AnthOligo was tested by capturing the Major Histocompatibility Complex (MHC) of a random sample.The application provides users with a simple interface to upload an input file in BED format and customize parameters for each task. The task of probe design in AnthOligo commences when a user uploads an input file and concludes with the generation of a result-set containing an optimal set of region-specific oligos. AnthOligo is currently available as a public web application with URL: http://antholigo.chop.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Assessing the utilization of high-resolution 2-field HLA typing in solid organ transplantation.

HLA typing in solid organ transplantation (SOT) is necessary for determining HLA-matching status between donor-recipient pairs and assessing patients' anti-HLA antibody profiles. Histocompatibility has traditionally been evaluated based on serologically defined HLA antigens. The evolution of HLA typing and antibody identification technologies, however, has revealed many limitations with using serologic equivalents for assessing compatibility in SOT. The significant improvements to HLA typing introduced by next-generation sequencing (NGS) require an assessment of the impact of this technology on SOT. We have assessed the role of high-resolution 2-field HLA typing (HR-2F) in SOT by retrospectively evaluating NGS-typed pre- and post-SOT cases. HR-2F typing was highly instructive or necessary in 41% (156/385) of the cases. Several pre- and posttransplant scenarios were identified as being better served by HR-2F typing. Five different categories are presented with specific case examples. The experience of another center (Temple University Hospital) is also included, whereby 21% of the cases required HR-2F typing by Sanger sequencing, as supported by other legacy methods, to properly address posttransplant anti-HLA antibody issues.

miR-155 promotes FLT3-ITD-induced myeloproliferative disease through inhibition of the interferon response.

FLT3-ITD+ acute myeloid leukemia (AML) accounts for ∼25% of all AML cases and is a subtype that carries a poor prognosis. microRNA-155 (miR-155) is specifically overexpressed in FLT3-ITD+ AML compared with FLT3 wild-type (FLT3-WT) AML and is critical for the growth of FLT3-ITD+ AML cells in vitro. However, miR-155's role in regulating FLT3-ITD-mediated disease in vivo remains unclear. In this study, we used a genetic mouse model to determine whether miR-155 influences the development of FLT3-ITD-induced myeloproliferative disease. Results indicate that miR-155 promotes FLT3-ITD-induced myeloid expansion in the bone marrow, spleen, and peripheral blood. Mechanistically, miR-155 increases proliferation of the hematopoietic stem and progenitor cell compartments by reducing the growth-inhibitory effects of the interferon (IFN) response, and this involves targeting of Cebpb. Consistent with our observations in mice, primary FLT3-ITD+ AML clinical samples have significantly higher miR-155 levels and a lower IFN response compared with FLT3-WT AML samples. Further, inhibition of miR-155 in FLT3-ITD+ AML cell lines using CRISPR/Cas9, or primary FLT3-ITD+ AML samples using locked nucleic acid antisense inhibitors, results in an elevated IFN response and reduces colony formation. Altogether, our data reveal that miR-155 collaborates with FLT3-ITD to promote myeloid cell expansion in vivo and that this involves a multitarget mechanism that includes repression of IFN signaling.

Nucleotide resolution analysis of TMPRSS2 and ERG rearrangements in prostate cancer.

TMPRSS2-ERG rearrangements occur in approximately 50% of prostate cancers and therefore represent one of the most frequently observed structural rearrangements in all cancers. However, little is known about the genomic architecture of such rearrangements. We therefore designed and optimized a pipeline involving target capture of TMPRSS2 and ERG genomic sequences coupled with paired-end next-generation sequencing to resolve genomic rearrangement breakpoints in TMPRSS2 and ERG at nucleotide resolution in a large series of primary prostate cancer specimens (n = 83). This strategy showed > 90% sensitivity and specificity in identifying TMPRSS2-ERG rearrangements, and allowed identification of intra- and inter-chromosomal rearrangements involving TMPRSS2 and ERG with known and novel fusion partners. Our results indicate that rearrangement breakpoints show strong clustering in specific intronic regions of TMPRSS2 and ERG. The observed TMPRSS2-ERG rearrangements often exhibited complex chromosomal architecture associated with several intra- and inter-chromosomal rearrangements. Nucleotide resolution analysis of breakpoint junctions revealed that the majority of TMPRSS2 and ERG rearrangements (~88%) occurred at or near regions of microhomology or involved insertions of one or more base pairs. This architecture implicates non-homologous end joining (NHEJ) and microhomology-mediated end joining (MMEJ) pathways in the generation of such rearrangements. These analyses have provided important insights into the molecular mechanisms involved in generating prostate cancer-specific recurrent rearrangements.

The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection.

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

A combination of HLA-DP α and ß chain polymorphisms paired with a SNP in the DPB1 3' UTR region, denoting expression levels, are associated with atopic dermatitis.

Introduction: Components of the immune response have previously been associated with the pathophysiology of atopic dermatitis (AD), specifically the Human Leukocyte Antigen (HLA) Class II region via genome-wide association studies, however the exact elements have not been identified. Methods: This study examines the genetic variation of HLA Class II genes using next generation sequencing (NGS) and evaluates the resultant amino acids, with particular attention on binding site residues, for associations with AD. The Genetics of AD cohort was used to evaluate HLA Class II allelic variation on 464 subjects with AD and 384 controls. Results: Statistically significant associations with HLA-DP α and ß alleles and specific amino acids were found, some conferring susceptibility to AD and others with a protective effect. Evaluation of polymorphic residues in DP binding pockets revealed the critical role of P1 and P6 (P1: α31M + (ß84G or ß84V) [protection]; α31Q + ß84D [susceptibility] and P6: α11A + ß11G [protection]) and were replicated with a national cohort of children consisting of 424 AD subjects. Independently, AD susceptibility-associated residues were associated with the G polymorphism of SNP rs9277534 in the 3' UTR of the HLA-DPB1 gene, denoting higher expression of these HLA-DP alleles, while protection-associated residues were associated with the A polymorphism, denoting lower expression. Discussion: These findings lay the foundation for evaluating non-self-antigens suspected to be associated with AD as they potentially interact with particular HLA Class II subcomponents, forming a complex involved in the pathophysiology of AD. It is possible that a combination of structural HLA-DP components and levels of expression of these components contribute to AD pathophysiology.

Genomic characterization of HLA class I and class II genes in ethnically diverse sub-Saharan African populations: A report on novel HLA alleles.

HLA allelic variation has been well studied and documented in many parts of the world. However, African populations have been relatively under-represented in studies of HLA variation. We have characterized HLA variation from 489 individuals belonging to 13 ethnically diverse populations from rural communities from the African countries of Botswana, Cameroon, Ethiopia, and Tanzania, known to practice traditional subsistence lifestyles using next generation sequencing (Illumina) and long-reads from Oxford Nanopore Technologies. We identified 342 distinct alleles among the 11 HLA targeted genes: HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQA1, -DQB1, -DPA1, and -DPB1, with 140 of those alleles containing novel sequences that were submitted to the IPD-IMGT/HLA database. Sixteen of the 140 alleles contained novel content within the exonic regions of the genes, while 110 alleles contained novel intronic variants. Four alleles were found to be recombinants of already described HLA alleles and 10 alleles extended the sequence content of already described alleles. All 140 alleles include complete allelic sequence from the 5' UTR to the 3' UTR that are inclusive of all exons and introns. This report characterizes the HLA allelic variation from these individuals and describes the novel allelic variation present within these specific African populations.

Ancestry adjustment improves genome-wide estimates of regional intolerance.

Genomic regions subject to purifying selection are more likely to carry disease-causing mutations than regions not under selection. Cross species conservation is often used to identify such regions but with limited resolution to detect selection on short evolutionary timescales such as that occurring in only one species. In contrast, genetic intolerance looks for depletion of variation relative to expectation within a species, allowing species-specific features to be identified. When estimating the intolerance of noncoding sequence, methods strongly leverage variant frequency distributions. As the expected distributions depend on ancestry, if not properly controlled for, ancestral population source may obfuscate signals of selection. We demonstrate that properly incorporating ancestry in intolerance estimation greatly improved variant classification. We provide a genome-wide intolerance map that is conditional on ancestry and likely to be particularly valuable for variant prioritization.

Cellular immune selection with hepatitis C virus persistence in humans.

Hepatitis C virus (HCV) infection frequently persists despite substantial virus-specific cellular immune responses. To determine if immunologically driven sequence variation occurs with HCV persistence, we coordinately analyzed sequence evolution and CD8+ T cell responses to epitopes covering the entire HCV polyprotein in subjects who were followed prospectively from before infection to beyond the first year. There were no substitutions in T cell epitopes for a year after infection in a subject who cleared viremia. In contrast, in subjects with persistent viremia and detectable T cell responses, we observed substitutions in 69% of T cell epitopes, and every subject had a substitution in at least one epitope. In addition, amino acid substitutions occurred 13-fold more often within than outside T cell epitopes (P < 0.001, range 5-38). T lymphocyte recognition of 8 of 10 mutant peptides was markedly reduced compared with the initial sequence, indicating viral escape. Of 16 nonenvelope substitutions that occurred outside of known T cell epitopes, 8 represented conversion to consensus (P = 0.015). These findings reveal two distinct mechanisms of sequence evolution involved in HCV persistence: viral escape from CD8+ T cell responses and optimization of replicative capacity.

Large-scale candidate gene analysis of spontaneous clearance of hepatitis C virus.

Human genetic variation is a determinant of recovery from acute hepatitis C virus (HCV) infection; however, to date, single-nucleotide polymorphisms (SNPs) in only a limited number of genes have been studied with respect to HCV clearance. We determined whether SNPs in 112 selected immune response genes are important for HCV clearance, by genotyping 1536 SNPs in a cohort of 343 persons with natural HCV clearance and 547 persons with HCV persistence. PLINK (version 1.05) and Haploview (version 4.1) software packages were used to perform association, permutation, and haplotype analyses stratified by African American and European American race. Of the 1536 SNPs tested, 1426 (92.8%) were successfully genotyped. In African Americans, we identified 18 SNPs located in 11 gene regions that were associated with HCV infection outcome (empirical P value, < .01). In European Americans, there were 20 SNPs located in 8 gene regions associated with HCV infection outcome. Four of the gene regions studied (TNFSF18, TANK, HAVCR1, and IL18BP) contained SNPs for which the empirical P value was <.01 in both of the race groups. In this large-scale analysis of 1426 genotyped SNPs in 112 candidate genes, we identified 4 gene regions that are likely candidates for a role in HCV clearance or persistence in both African Americans and European Americans.

Spatial constrains and information content of sub-genomic regions of the human genome.

Complexity metrics and machine learning (ML) models have been utilized to analyze the lengths of segmental genomic entities of DNA sequences (exonic, intronic, intergenic, repeat, unique) with the purpose to ask questions regarding the segmental organization of the human genome within the size distribution of these sequences. For this we developed an integrated methodology that is based upon the reconstructed phase space theorem, the non-extensive statistical theory of Tsallis, ML techniques, and a technical index, integrating the generated information, which we introduce and named complexity factor (COFA). Our analysis revealed that the size distribution of the genomic regions within chromosomes are not random but follow patterns with characteristic features that have been seen through its complexity character, and it is part of the dynamics of the whole genome. Finally, this picture of dynamics in DNA is recognized using ML tools for clustering, classification, and prediction with high accuracy.

Complex Linkage Disequilibrium Effects in HLA-DPB1 Expression and Molecular Mismatch Analyses of Transplantation Outcomes.

BACKGROUND: HLA molecular mismatch (MM) is a risk factor for de novo donor-specific antibody (dnDSA) development in solid organ transplantation. HLA expression differences have also been associated with adverse outcomes in hematopoietic cell transplantation. We sought to study both MM and expression in assessing dnDSA risk. METHODS: One hundred three HLA-DP-mismatched solid organ transplantation pairs were retrospectively analyzed. MM was computed using amino acids (aa), eplets, and, supplementarily, Grantham/Epstein scores. DPB1 alleles were classified as rs9277534-A (low-expression) or rs9277534-G (high-expression) linked. To determine the associations between risk factors and dnDSA, logistic regression, linkage disequilibrium (LD), and population-based analyses were performed. RESULTS: A high-risk AA:GX (recipient:donor) expression combination (X = A or G) demonstrated strong association with HLA-DP dnDSA (P = 0.001). MM was also associated with HLA-DP dnDSA when evaluated by itself (eplet P = 0.007, aa P = 0.003, Grantham P = 0.005, Epstein P = 0.004). When attempting to determine the relative individual effects of the risk factors in multivariable analysis, only AA:GX expression status retained a strong association (relative risk = 18.6, P = 0.007 with eplet; relative risk = 15.8, P = 0.02 with aa), while MM was no longer significant (eplet P = 0.56, aa P = 0.51). Importantly, these risk factors are correlated, due to LD between the expression-tagging single-nucleotide polymorphism and polymorphisms along HLA-DPB1. CONCLUSIONS: The MM and expression risk factors each appear to be strong predictors of HLA-DP dnDSA and to possess clinical utility; however, these two risk factors are closely correlated. These metrics may represent distinct ways of characterizing a common overlapping dnDSA risk profile, but they are not independent. Further, we demonstrate the importance and detailed implications of LD effects in dnDSA risk assessment and possibly transplantation overall.

A Stat1 bound enhancer promotes Nampt expression and function within tumor associated macrophages.

Tumor associated macrophage responses are regulated by distinct metabolic states that affect their function. However, the ability of specific signals in the local tumor microenvironment to program macrophage metabolism remains under investigation. Here, we identify NAMPT, the rate limiting enzyme in NAD salvage synthesis, as a target of STAT1 during cellular activation by interferon gamma, an important driver of macrophage polarization and antitumor responses. We demonstrate that STAT1 occupies a conserved element within the first intron of Nampt, termed Nampt-Regulatory Element-1 (NRE1). Through disruption of NRE1 or pharmacological inhibition, a subset of M1 genes is sensitive to NAMPT activity through its impact on glycolytic processes. scRNAseq is used to profile in vivo responses by NRE1-deficient, tumor-associated leukocytes in melanoma tumors through the creation of a unique mouse strain. Reduced Nampt and inflammatory gene expression are present in specific myeloid and APC populations; moreover, targeted ablation of NRE1 in macrophage lineages results in greater tumor burden. Finally, elevated NAMPT expression correlates with IFNγ responses and melanoma patient survival. This study identifies IFN and STAT1-inducible Nampt as an important factor that shapes the metabolic program and function of tumor associated macrophages.

Interaction between RANTES promoter variant and CCR5Delta32 favors recovery from hepatitis B.

Recovery from acute hepatitis B virus (HBV) infection occurs in 95% of adult-acquired infections. A 32-bp deletion in CCR5 (CCR5Delta32), which encodes for a nonfunctional receptor, increases the likelihood of recovery. Using 181 subjects with persistent HBV infection and 316 who had recovered, we tested the hypothesis that an epistatic interaction between functional polymorphisms in RANTES (a CCR5 ligand) and CCR5 impacts recovery. Specific models designed to assess individual contributions of compound genotypes demonstrated that the only combination associated with recovery from an HBV infection was RANTES -403A with CCR5Delta32 (odds ratio 0.36, p = 0.02). Because the phenotypic consequence of -403A is reported to be higher levels of RANTES, we propose a model in which excess RANTES in combination with low CCR5 favors recovery from an HBV infection, which will require validation through functional testing.

Hepatitis C virus immune escape via exploitation of a hole in the T cell repertoire.

Hepatitis C virus (HCV) infection frequently persists despite eliciting substantial virus-specific immune responses. Thus, HCV infection provides a setting in which to investigate mechanisms of immune escape that allow for viral persistence. Viral amino acid substitutions resulting in decreased MHC binding or impaired Ag processing of T cell epitopes reduce Ag density on the cell surface, permitting evasion of T cell responses in chronic viral infection. Substitutions in viral epitopes that alter TCR contact residues frequently result in escape, but via unclear mechanisms because such substitutions do not reduce surface presentation of peptide-MHC complexes and would be expected to prime T cells with new specificities. We demonstrate that a known in vivo HCV mutation involving a TCR contact residue significantly diminishes T cell recognition and, in contrast to the original sequence, fails to effectively prime naive T cells. This mutant epitope thus escapes de novo immune recognition because there are few highly specific cognate TCR among the primary human T cell repertoire. This example is the first on viral immune escape via exploitation of a "hole" in the T cell repertoire, and may represent an important general mechanism of viral persistence.

Genomic characterization of MICA gene using multiple next generation sequencing platforms: A validation study.

We have developed a protocol regarding the genomic characterization of the MICA gene by next generation sequencing (NGS). The amplicon includes the full length of the gene and is about 13 kb. A total of 156 samples were included in the study. Ninety-seven of these samples were previously characterized at MICA by legacy methods (Sanger or sequence specific oligonucleotide) and were used to evaluate the accuracy, precision, specificity, and sensitivity of the assay. An additional 59 DNA samples of unknown ethnicity volunteers from the United States were only genotyped by NGS. Samples were chosen to contain a diverse set of alleles. Our NGS approach included a first round of sequencing on the Illumina MiSeq platform and a second round of sequencing on the MinION platform by Oxford Nanopore Technology (ONT), on selected samples for the purpose of either characterizing new alleles or setting phase among multiple polymorphisms to resolve ambiguities or generate complete sequence for alleles that were only partially reported in the IMGT/HLA database. Complete consensus sequences were generated for every allele sequenced with ONT, extending from the 5' untranslated region (UTR) to the 3' UTR of the MICA gene. Thirty-two MICA sequences were submitted to the IMGT/HLA database including either new alleles or filling up the gaps (exonic, intronic and/or UTRs) of already reported alleles. Some of the challenges associated with the characterization of these samples are discussed.