Search

HLA genotyping by next-generation sequencing (NGS) has evolved with significant advancements in the last decade. Here we describe full-length HLA genotyping of 11 loci in 612 individuals comprising a dengue vaccine cohort from Cebu province in the Philippines. The multi-locus individual tagging NGS (MIT-NGS) method that we developed initially for genotyping 4–6 loci in one MiSeq run was expanded to 11 loci including HLA-A, B, C, DPA1, DPB1, DQA1, DQB1, DRB1, and DRB3/4/5. This change did not affect the overall coverage or depth of the sequencing reads. HLA alleles with frequencies greater than 10% were A*11:01:01, A*24:02:01, A*24:07:01, A*34:01:01, B*38:02:01, B*15:35, B*35:05:01, C*07:02:01, C*04:01:01, DPA1*02:02:02, DPB1*05:01:01, DPB1*01:01:01, DQA1*01:02:01, DQA1*06:01:01, DQB1*05:02:01, DQB1*03:01:01, DRB1*15:02:01, DRB1*12:02:01, DRB3*03:01:03, DRB4*01:03:01, and DRB5*01:01:01. Improvements in sequencing library preparation provide uniform and even coverage across all exons and introns. This has led to a marked reduction in allele imbalance and dropout. Furthermore, including more loci, such as DRB3/4/5, decreases cross-mapping and incorrect allele assignment at the DRB1 locus. The increased number of loci sequenced for each sample does not reduce the number of samples that can be multiplexed on a single MiSeq run and is therefore more cost-efficient. We believe that such improvements will help HLA genotyping by NGS to gain momentum over other conventional methods by increasing confidence in the calls. [ABSTRACT FROM AUTHOR]

Abstract The human leukocyte antigen (HLA) genes are highly variable and are known to play an important role in disease outcomes, including infectious diseases. Prior knowledge of HLA polymorphisms in a population usually forms the basis for an effective case-control study design. As a prelude to future disease association analyses, we report HLA class I and II diversity in 334 unrelated donors from a Dengue vaccine efficacy trial conducted in Thailand. Long-range PCR amplification of six HLA loci was performed on DNA extracted from saliva samples. HLA-A, -B, -C, -DPB1, -DQB1 and -DRB1 were genotyped using a next-generation sequencing method presented at the 17th International HLA and Immunogenetics Workshop. In total, we identified 201 HLA alleles, including 35 HLA-A, 57 HLA-B, 28 HLA-C, 24 HLA-DPB1, 21 HLA-DQB1 and 36 HLA-DRB1 alleles. Very common HLA alleles with frequencies greater than 10 percent were A∗11:01:01, A∗33:03:01, A∗24:02:01, B∗46:01:01, C∗07:02:01, C∗01:02:01, C∗08:01:01, DPB1∗05:01:01, DPB1∗13:01:01, DPB1∗04:01:01, DPB1∗02:01:02, DQB1∗03:01:01, DQB1∗05:02:01, DQB1∗03:03:02, DRB1∗12:02:01, DRB1∗09:01:02, and DRB1∗15:02:01. A novel HLA allele, B∗15:450, had a non-synonymous substitution and occurred in more than one donor. Population-based full-length NGS HLA typing is more conclusive and provides a sound foundation for exploring disease association in a given population. [ABSTRACT FROM AUTHOR]

Next generation sequencing ( NGS) methods have been established as an efficient approach for HLA typing because unlike traditional Sanger sequencing, they provide unambiguous results at a reasonable cost. We previously developed a multi-locus index method to genotype four HLA loci (A, B, C, and DRB1) on the Illumina MiSeq platform. We have now expanded this method to include two additional loci, HLA-DPB1 and DQB1. Contiguous full-length amplicons from 5′ UTR through 3′ UTR regions were generated using one long-range PCR reaction per locus for each of the six loci from 96 individuals of different ethnicities. The six amplicons from each donor were pooled, enzymatically fragmented and given a donor-specific index. This approach enabled sequencing of 576 loci from 96 individuals in a single MiSeq run. Donor-specific sequence reads were demultiplexed, and allele calls were generated from FASTQ files using commercially available software. Comparison to HLA genotypes generated from Sanger sequence-based typing ( SBT) identified no discordances among any of the alleles analyzed in this study. Importantly, this method was able to resolve 22 DPB1 and 20 DQB1 alleles that were ambiguous with the SBT method. Furthermore, a novel allele in each of these two loci was identified, with the DQB1*05:01:24 allele having a frequency of greater than five percent. This method was subsequently validated against a blinded panel of 22 samples from the 17th International HLA and Immunogenetics Workshop. The flexibility of the method is further highlighted by successful genotyping of eight loci comprising all classical HLA loci for a subset of the samples. We now present a high-throughput, high-resolution, scalable NGS HLA typing method to accurately and efficiently genotype all classical HLA class I and II loci. [ABSTRACT FROM AUTHOR]

Background: Unambiguous human leukocyte antigen (HLA) typing is important in transplant matching and disease association studies. High-resolution HLA typing that is not restricted to the peptide-binding region can decrease HLA allele ambiguities. Cost and technology constraints have hampered high-throughput and efficient high resolution unambiguous HLA typing. We have developed a method for HLA genotyping that preserves the very high-resolution that can be obtained by next-generation sequencing (NGS) but also achieves substantially increased efficiency. Unambiguous HLA-A, B, C and DRB1 genotypes can be determined for 96 individuals in a single run of the Illumina MiSeq. Results: Long-range amplification of full-length HLA genes from four loci was performed in separate polymerase chain reactions (PCR) using primers and PCR conditions that were optimized to reduce co-amplification of other HLA loci. Amplicons from the four HLA loci of each individual were then pooled and subjected to enzymatic library generation. All four loci of an individual were then tagged with one unique index combination. This multi-locus individual tagging (MIT) method combined with NGS enabled the four loci of 96 individuals to be analyzed in a single 500 cycle sequencing paired-end run of the Illumina-MiSeq. The MIT-NGS method generated sequence reads from the four loci were then discriminated using commercially available NGS HLA typing software. Comparison of the MIT-NGS with Sanger sequence-based HLA typing methods showed that all the ambiguities and discordances between the two methods were due to the accuracy of the MIT-NGS method. Conclusions: The MIT-NGS method enabled accurate, robust and cost effective simultaneous analyses of four HLA loci per sample and produced 6 or 8-digit high-resolution unambiguous phased HLA typing data from 96 individuals in a single NGS run. [ABSTRACT FROM AUTHOR]

A gene signature was previously found to be correlated with mosaic adenovirus 26 vaccine protection in simian immunodeficiency virus and simian-human immunodeficiency virus challenge models in non-human primates. In this report, we investigated the presence of this signature as a correlate of reduced risk in human clinical trials and potential mechanisms of protection. The absence of this gene signature in the DNA/rAd5 human vaccine trial, which did not show efficacy, strengthens our hypothesis that this signature is only enriched in studies that demonstrated protection. This gene signature was enriched in the partially effective RV144 human trial that administered the ALVAC/protein vaccine, and we find that the signature associates with both decreased risk of HIV-1 acquisition and increased vaccine efficacy (VE). Total RNA-seq in a clinical trial that used the same vaccine regimen as the RV144 HIV vaccine implicated antibody-dependent cellular phagocytosis (ADCP) as a potential mechanism of vaccine protection. CITE-seq profiling of 53 surface markers and transcriptomes of 53,777 single cells from the same trial showed that genes in this signature were primarily expressed in cells belonging to the myeloid lineage, including monocytes, which are major effector cells for ADCP. The consistent association of this transcriptome signature with VE represents a tool both to identify potential mechanisms, as with ADCP here, and to screen novel approaches to accelerate the development of new vaccine candidates. [ABSTRACT FROM AUTHOR]

Highlights from the article: Given that there are known HLA associations with VL control in natural infection in the absence of ART, we hypothesized that differences in virus rebound could be due to specific protective HLA alleles. This individual carried the protective HLA-B*57 allele known to associate with HIV-1 VL control in natural infection in the absence of ART. The full genotype included HLA-A*01:01:01+HLA-A*11:01:01 HLA-B*39:09:01+HLA-B*57:01:01 HLA-C*06:02:01+HLA-C*07:02:01 HLA-DPA1*01:03:01+HLA-DPA1*02:01:01 HLA-DPB1*03:01:01+HLA-DPB1*26:01:02 HLA-DQA1*01:04:01+HLA-DQA1*02:01:01 HLA-DQB1*03:03:02+HLA-DQB1*05:02:01 HLA-DRB1*07:01:01+HLA-DRB1*14:54:01.

Aim RNA sequencing (RNA-Seq) is a powerful next-generation sequencing method for discovering, profiling and quantifying RNA transcripts. RNA-Seq can also identify splice junctions and single nucleotide polymorphisms (SNP) in coding exons. In this study, we tested if RNA-Seq reads provide sufficient coverage and depth for high-resolution HLA genotyping. Methods RNA-Seq data was available from 24 individuals from an acutely infected HIV cohort. We used two methods to extract HLA reads from the RNA-Seq data: 1) FR-HIT recruited RNA-Seq reads to the IMGT coding sequence reference database; and 2) the SNP-tolerant aligner GSNAP aligned the reads to a constructed GMAP database consisting of cDNA from HLA genes. The resulting HLA-specific FASTQ reads that were generated were analyzed by commercial HLA typing software. Results We obtained 100% sequence coverage and a median sequence read depth of 331, which was sufficient to genotype 15 HLA loci from each of the donors. HLA genotypes for 8 classical HLA class I and II loci matched genotypes generated from the Sanger sequence-based typing method, as well as targeted next-generation sequencing of genomic DNA. We identified a novel HLA-DQB1 allele in the Thai population with a frequency greater than 5% that was identical to the 05:01:01 allele in all exons except exon 1, where it differed by one SNP at position 81. The DQB1 * 05:01:01 allele, like several other DQB1 alleles, has exon 5 spliced out. This deletion however, cannot be confirmed from the reads that are obtained by sequencing genomic DNA. RNA-Seq confirmed that the novel allele does not express exon 5, as the RNA-Seq data showed a complete absence of exon 5 reads compared to reads from targeted genomic sequences. Conclusions We have successfully genotyped 15 HLA loci and confirmed a splice variant in a novel allele using RNA-Seq. This method highlights the importance of mining RNA-Seq data for high-resolution HLA genotyping. Additional analysis methods may be able to identify allele-specific expression of HLA alleles from RNA-Seq data. [ABSTRACT FROM AUTHOR]

Aim The Multi-locus Individual Tagging – Next Generation Sequencing (MIT-NGS) method was previously shown to efficiently genotype four HLA loci (A, B, C, and DRB1) from 96 donors while maintaining phase, and resolving ambiguities associated with traditional sequence based HLA typing (Ehrenberg et al., BMC Genomics, 2014). Here we expanded MIT-NGS to six loci to include DPB1 and DQB1, and evaluated HLA genotypes for all loci against those obtained from both the HLA sequence-based typing (SBT) method and the four loci MIT-NGS method. Methods Full-length genes of HLA-A, B, C, DPB1, DQB1 and DRB1 loci from 96 donors of Caucasian, African and Asian ethnicities were amplified using long-range PCR. Individual donor pools, each consisting of equimolar inputs corresponding to the target size of the six HLA loci amplicons were enzymatically fragmented and uniquely indexed. This 6-plex approach enabled sequencing of 576 loci from 96 individuals in a single Illumina MiSeq run. Locus specific sequence reads were demultiplexed, aligned, and genotyped using commercially available software and an in-house pipeline. Results There were no observed discordances between the 6-plex MIT-NGS method and the combined SBT or 4-plex NGS methods for all 1152 alleles analyzed in this study. Furthermore, this method was able to resolve 20 DQBI and 22 DPB1 genotypes that were ambiguous with SBT. These alleles included DPB1 ∗ 104:01,DPB1 ∗ 105:01, DPB1 ∗ 107:01, DPB1 ∗ 131:01, DPB1 ∗ 135:01, DPB1 ∗ 296:01, DPB1 ∗ 414:01, DQB1 ∗ 02:02, DQB1 ∗ 03:19 and DQB1 ∗ 06:01:03. All of these ambiguous changes affected the protein coding sequence. Conclusions We have developed a six-locus high resolution HLA typing method to accurately and efficiently genotype all classical HLA class I and II loci while maintaining phase, and resolving ambiguities associated with conventional sequence based HLA typing. [ABSTRACT FROM AUTHOR]

The RV144 HIV vaccine trial in Thailand elicited antibody responses to the envelope of HIV-1, which correlated significantly with the risk of HIV-1 acquisition. Human leukocyte antigen ( HLA) class II molecules are essential in antigen presentation to CD4 T cells for activation of B cells to produce antibodies. We genotyped the classical HLA-DRB1, DQB1, and DPB1 genes in 450 individuals from the placebo arm of the RV144 study to determine the background allele and haplotype frequencies of these genes in this cohort. High-resolution 4 and 6-digit class II HLA typing data was generated using sequencing-based methods. The observed diversity for the HLA loci was 33 HLA-DRB1, 15 HLA-DQB1, and 26 HLA-DPB1 alleles. Common alleles with frequencies greater than 10% were DRB1*07:01, DRB1*09:01, DRB1*12:02, DRB1*15:02, DQB1*02:01/02, DQB1*03:01, DQB1*03:03, DQB1*05:01, DQB1*05:02, DPB1*04:01:01, DPB1*05:01:01, and DPB1*13:01:01. We identified 28 rare alleles with frequencies of less than 1% in the Thai individuals. Ambiguity for HLA- DPB1*28:01 in exon 2 was resolved to DPB1*296:01 by next-generation sequencing of all exons. Multi-locus haplotypes including HLA class I and II loci were reported in this study. This is the first comprehensive report of allele and haplotype frequencies of all three HLA class II genes from a Thai population. A high-resolution genotyping method such as next-generation sequencing avoids missing rare alleles and resolves ambiguous calls. The HLA class II genotyping data generated in this study will be beneficial not only for future disease association/vaccine efficacy studies related to the RV144 study, but also for similar studies in other diseases in the Thai population, as well as population genetics and transplantation studies. [ABSTRACT FROM AUTHOR]

Fragaria vesca , a diploid woodland strawberry with a small and sequenced genome, is an excellent model for studying fruit development. The strawberry fruit is unique in that the edible flesh is actually enlarged receptacle tissue. The true fruit are the numerous dry achenes dotting the receptacle's surface. Auxin produced from the achene is essential for the receptacle fruit set, a paradigm for studying crosstalk between hormone signaling and development. To investigate the molecular mechanism underlying strawberry fruit set, next-generation sequencing was employed to profile early-stage fruit development with five fruit tissue types and five developmental stages from floral anthesis to enlarged fruits. This two-dimensional data set provides a systems-level view of molecular events with precise spatial and temporal resolution. The data suggest that the endosperm and seed coat may play a more prominent role than the embryo in auxin and gibberellin biosynthesis for fruit set. A model is proposed to illustrate how hormonal signals produced in the endosperm and seed coat coordinate seed, ovary wall, and receptacle fruit development. The comprehensive fruit transcriptome data set provides a wealth of genomic resources for the strawberry and Rosaceae communities as well as unprecedented molecular insight into fruit set and early stage fruit development. [ABSTRACT FROM AUTHOR]

The recombinant, live, attenuated, tetravalent dengue vaccine CYD-TDV has shown efficacy against all four dengue serotypes. In this exploratory study (CYD59, NCT02827162), we evaluated potential associations of host human leukocyte antigen (HLA) alleles with dengue antibody responses, CYD-TDV vaccine efficacy, and virologically-confirmed dengue (VCD) cases. Children 4–11 years old, who previously completed a phase 2b efficacy study of CYD-TDV in a single center in Thailand, were included in the study. Genotyping of HLA class I and II loci was performed by next-generation sequencing from DNA obtained from 335 saliva samples. Dengue neutralizing antibody titers (NAb) were assessed as a correlate of risk and protection. Regression analyses were used to assess associations between HLA alleles and NAb responses, vaccine efficacy, and dengue outcomes. Month 13 NAb log geometric mean titers (GMTs) were associated with decreased risk of VCD. In the vaccine group, HLA-DRB1*11 was significantly associated with higher NAb log GMT levels (beta: 0.76; p = 0.002, q = 0.13). Additionally, in the absence of vaccination, HLA associations were observed between the presence of DPB1*03:01 and increased NAb log GMT levels (beta: 1.24; p = 0.005, q = 0.17), and between DPB1*05:01 and reduced NAb log GMT levels (beta: −1.1; p = 0.001, q = 0.07). This study suggests associations of HLA alleles with NAb titers in the context of dengue outcomes. This study was registered with clinicaltrials.gov: NCT02827162. [ABSTRACT FROM AUTHOR]

Flowers are reproductive organs and precursors to fruits and seeds. While the basic tenets of the ABCE model of flower development are conserved in angiosperms, different flowering plants exhibit different and sometimes unique characteristics. A distinct feature of strawberry (Fragaria spp.) flowers is the development of several hundreds of individual apocarpous (unfused) carpels. These individual carpels are arranged in a spiral pattern on the subtending stem tip, the receptacle. Therefore, the receptacle is an integral part of the strawberry flower and is of significant agronomic importance, being the precursor to strawberry fruit. Taking advantage of next-generation sequencing and laser capture microdissection, we generated different tissue- and stage-transcriptomic profiling of woodland strawberry (Fragaria vesca) flower development. Using pairwise comparisons and weighted gene coexpression network analysis, we identified modules of coexpressed genes and hub genes of tissue-specific networks. Of particular importance is the discovery of a developing receptacle-specific module exhibiting similar molecular features to those of young floral meristems. The strawberry homologs of a number of meristem regulators, including LOST MERISTEM and WUSCHEL, are identified as hub genes operating in the developing receptacle network. Furthermore, almost 25% of the F-box genes in the genome are transiently induced in developing anthers at the meiosis stage, indicating active protein degradation. Together, this work provides important insights into the molecular networks underlying strawberry's unique reproductive developmental processes. This extensive floral transcriptome data set is publicly available and can be readily queried at the project Web site, serving as an important genomic resource for the plant biology research community. [ABSTRACT FROM AUTHOR]

RV144 is the first phase 3 HIV vaccine clinical trial to demonstrate efficacy. This study consisted of more than 8,000 individuals in each arm of the trial, representing the four major regions of Thailand. Human leukocyte antigen (HLA) class I and killer cell immunoglobulin-like receptor (KIR) genes, as well as 96 genome-wide ancestry informative markers (AIMs) were genotyped in 450 placebo HIV-1-uninfected individuals to identify the immunogenetic diversity and population structure of this cohort. High-resolution genotyping identified the common HLA alleles as A*02:03, A*02:07, A*11:01, A*24:02, A*24:07, A*33:03, B*13:01, B*15:02, B*18:01, B*40:01, B*44:03, B*46:01, B*58:01, C*01:02, C*03:02, C*03:04, C*07:01, C*07:02, C*07:04, and C*08:01. The most frequent three-loci haplotype was B*46:01-C*01:02-A*02:07. Framework genes KIR2DL4, 3DL2, and 3DL3 were present in all samples, and KIR2DL1, 2DL3, 3DL1, 2DS4, and 2DP1 occurred at frequencies greater than 90 %. The combined HLA and KIR profile suggests admixture with neighboring Asian populations. Principal component and correspondence analyses comparing the RV144 samples to the phase 3 International HapMap Project (HapMap3) populations using AIMs corroborated these findings. Structure analyses identified a distinct profile in the Thai population that did not match the Asian or other HapMap3 samples. This shows genetic variability unique to Thais in RV144, making it essential to take into account population stratification while performing genetic association studies. The overall analyses from all three genetic markers indicate that the RV144 samples are representative of the Thai population. This will inform subsequent host genetic analyses in the RV144 cohort and provide insight for future genetic association studies in the Thai population. [ABSTRACT FROM AUTHOR]

Background Fragaria vesca, a diploid strawberry species commonly known as the alpine or woodland strawberry, is a versatile experimental plant system and an emerging model for the Rosaceae family. An ancestral F. vesca genome contributed to the genome of the octoploid dessert strawberry (F. ×ananassa), and the extant genome exhibits synteny with other commercially important members of the Rosaceae family such as apple and peach. To provide a molecular description of floral organ and fruit development at the resolution of specific tissues and cell types, RNAs from flowers and early developmental stage fruit tissues of the inbred F. vesca line YW5AF7 were extracted and the resulting cDNA libraries sequenced using an Illumina HiSeq2000. To enable easy access as well as mining of this two-dimensional (stage and tissue) transcriptome dataset, a web-based database, the Strawberry Genomic Resource (SGR), is developed. Description SGR is a web accessible database that contains sample description, sample statistics, gene annotation, and gene expression analysis. This information can be accessed publicly from a web-based interface at http://bioinformatics.towson.edu/strawberry/Default.aspx. The SGR website provides user friendly search and browse capabilities for all the data stored in the database. Users are able to search for genes using a gene ID or description or obtain differentially expressed genes by entering different comparison parameters. Search results can be downloaded in a tabular format compatible with Microsoft excel application. Aligned reads to individual genes and exon/intron structures are displayed using the genome browser, facilitating gene re-annotation by individual users. Conclusions The SGR database is developed to facilitate dissemination and data mining of extensive floral and fruit transcriptome data in the woodland strawberry. It enables users to mine the data in different ways to study different pathways or biological processes during reproductive development. [ABSTRACT FROM AUTHOR]