Report on the effects of fragment size, indexing, and read length on HLA sequencing on the Illumina MiSeq

A national proficiency scheme for human leucocyte antigen typing by next-generation sequencing

BACKGROUND

Accurate human leucocyte antigen (HLA) typing is essential for solid organ transplantation (SOT) and hematopoietic stem cell transplantation (HSCT). Next-generation sequencing (NGS) is increasingly adopted by clinical laboratories performing HLA typing. To ensure reliable NGS-based HLA typing throughout China, a HLA-NGS proficiency testing (PT) scheme was set up.

METHODS

A panel of 10 DNA reference materials for 11 HLA loci were distributed to 24 participants. They were asked to submit NGS typing results for at least six loci including HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQB1, and HLA-DPB1. Information regarding individual laboratory's sample preparation method, NGS platform, and genotyping software was also collected.

RESULTS

Thirteen participants correctly identified HLA alleles in all samples. The total concordance rate for all HLA alleles typed was 99.2% (3390/3420). The overall discordance was 0.2% (3/1440) for Class I, and 1.4% (27/1980) for Class II loci. Further investigation of the typing mistakes pointed to poor sequencing coverage and wrong allele calls. The errors identified reflect the deficient aspect of NGS in HLA typing, which may have implications for quality improvement for relevant laboratories.

CONCLUSIONS

The results of this PT underscore the necessity of regular participation in external quality assessment for clinical laboratories that perform HLA-NGS typing.

Recent Advances of Human Leukocyte Antigen (HLA) Typing Technology Based on High-Throughput Sequencing

The major histocompatibility complex (MHC) in humans is a genetic region consisting of cell surface proteins located on the short arm of chromosome 6. This is also known as the human leukocyte antigen (HLA) region. The HLA region consists of genes that exhibit complex genetic polymorphisms, and are extensively involved in immune responses. Each individual has a unique set of HLAs. Donor-recipient HLA allele matching is an important factor for organ transplantation. Therefore, an established rapid and accurate HLA typing technology is instrumental to preventing graft-verses-host disease (GVHD) in organ recipients. As of recent, high-throughput sequencing has allowed for an increase read length and higher accuracy and throughput, thus achieving complete and high-resolution full-length typing. With more advanced nanotechnology used in high-throughput sequencing, HLA typing is more widely used in third-generation single-molecule sequencing. This review article summarizes some of the most widely used sequencing typing platforms and evaluates the latest developments in HLA typing kits and their clinical applications.

Challenges in the application of NGS in the clinical laboratory

Next-generation sequencing (NGS), also known as massively parallel sequencing, has revolutionized genomic research. The current advances in NGS technology make it possible to provide high resolution, high throughput HLA typing in clinical laboratories. The focus of this review is on the recent development and implementation of NGS in clinical laboratories. Here, we examine the critical role of NGS technologies in clinical immunology for HLA genotyping. Two major NGS platforms (Illumina and Ion Torrent) are characterized including NGS library preparation, data analysis, and validation. Challenges of NGS implementation in the clinical laboratory are also discussed, including sequencing error rate, bioinformatics, result interpretation, analytic sensitivity, as well as large data storage. This review aims to promote the broader applications of NGS technology in clinical laboratories and advocate for the novel applications of NGS to drive future research.

Polishing the Oxford Nanopore long-read assemblies of bacterial pathogens with Illumina short reads to improve genomic analyses

Oxford Nanopore sequencing has been widely used to achieve complete genomes of bacterial pathogens. However, the error rates of Oxford Nanopore long reads are high. Various polishing algorithms using Illumina short reads to correct the errors in Oxford Nanopore long-read assemblies have been developed. The impact of polishing the Oxford Nanopore long-read assemblies of bacterial pathogens with Illumina short reads on improving genomic analyses was evaluated using both simulated and real reads. Ten species (10 strains) were selected for simulated reads, while real reads were tested on 11 species (11 strains). Oxford Nanopore long reads were assembled with Unicycler to produce a draft assembly, followed by three rounds of polishing with Illumina short reads using two polishing tools, Pilon and NextPolish. One round of NextPolish polishing generated genome completeness and accuracy parameters similar to the reference genomes, whereas two or three rounds of Pilon polishing were needed, though contiguity remained unchanged after polishing. The polished assemblies of Escherichia coli O157:H7, Salmonella Typhimurium, and Cronobacter sakazakii with simulated reads did not provide accurate plasmid identifications. One round of NextPolish polishing was needed for accurately identifying plasmids in Staphylococcus aureus and E. coli O26:H11 with real reads, whereas one and two rounds of Pilon polishing were necessary for these two strains, respectively. Polishing failed to provide an accurate antimicrobial resistance (AMR) genotype for S. aureus with real reads. One round of polishing recovered an accurate AMR genotype for Klebsiella pneumoniae with real reads. The reference genome and draft assembly of Citrobacter braakii with real reads differed, which carried blaCMY-83 and fosA6, respectively, while both genes were present after one round of polishing. However, polishing did not improve the assembly of E. coli O26:H11 with real reads to achieve numbers of virulence genes similar to the reference genome. The draft and polished assemblies showed a phylogenetic tree topology comparable with the reference genomes. For multilocus sequence typing and pan-genome analyses, one round of NextPolish polishing was sufficient to obtain accurate results, while two or three rounds of Pilon polishing were needed. Overall, NextPolish outperformed Pilon for polishing the Oxford Nanopore long-read assemblies of bacterial pathogens, though both polishing strategies improved genomic analyses compared to the draft assemblies.

PROBABILISTIC MODELING AND ANALYSIS OF DNA FRAGMENTATION

Deoxyribonucleic Acid (DNA) sequencing has become indispensable to the modern biological and medicine sciences. DNA fragmentation is a preliminary step in a dominant technique called shotgun sequencing that provides a time and cost effective strategy for the DNA sequencing. In this paper, we propose a probabilistic model for the random DNA fragmentation and derive an average number of fragments with the suitable length along with the probability of covering the entire DNA strand through the de novo assembly or the referenced-based mapping assembly. We formulate the coverage problem in terms of the probability of bond breaking between nucleotides and the number of DNA molecules participating in the fragmentation process, and provide insights into the optimal DNA fragmentation. We obtain the lower bound for the minimum number of suitable fragments required to reconstruct the DNA strand with the specified reliability. We evaluate the derived results with our DNA Fragmentation Tool which demonstrate, the validity of these results based on our model. Finally, we update our model with respect to the fragments’ size distribution of real data.

Long-range PCR libraries and next-generation sequencing for pharmacogenetic studies of patients treated with anti-TNF drugs

Biological therapy with anti-tumor necrosis factor-α (anti-TNF-α) monoclonal antibodies significantly increased the effectiveness of autoimmune disease treatment compared with conventional medicines. However, anti-TNF-α drugs are relatively expensive and a response to the therapy is reported in only 60-70% of patients. Moreover, in up to 5% of patients adverse drug reactions occur. The various effects of biological treatment may be a potential consequence of interindividual genetic variability. Only a few studies have been conducted in this field and which refer to single gene loci. Our aim was to design and optimize a methodology for a broader application of pharmacogenetic studies in patients undergoing anti-TNF-α treatment. Based on the current knowledge, we selected 16 candidate genes: TNFRSF1A, TNFRSF1B, ADAM17, CASP9, FCGR3A, LTA, TNF, FAS, IL1B, IL17A, IL6, MMP1, MMP3, S100A8, S100A9, and S100A12, which are potentially involved in the response to anti-TNF-α therapy. As a research model, three DNA samples from Crohn's disease (CD) patients were used. Targeted genomic regions were amplified in 23 long-range (LR) PCR reactions and after enzymatic fragmentation amplicon libraries were prepared and analyzed by next-generation sequencing (NGS). Our results indicated 592 sequence variations located in all fragments with coverage range of 5-1089. We demonstrate a highly sensitive, flexible, rapid, and economical approach to the pharmacogenetic investigation of anti-TNF-α therapy using amplicon libraries and NGS technology.