HLA DQA1 genes generate multiple transcripts by alternative splicing and polyadenylation of the 3' untranslated region

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.

An atlas of alternative polyadenylation quantitative trait loci contributing to complex trait and disease heritability

Genome-wide association studies have identified thousands of noncoding variants associated with human traits and diseases. However, the functional interpretation of these variants is a major challenge. Here, we constructed a multi-tissue atlas of human 3'UTR alternative polyadenylation (APA) quantitative trait loci (3'aQTLs), containing approximately 0.4 million common genetic variants associated with the APA of target genes, identified in 46 tissues isolated from 467 individuals (Genotype-Tissue Expression Project). Mechanistically, 3'aQTLs can alter poly(A) motifs, RNA secondary structure and RNA-binding protein-binding sites, leading to thousands of APA changes. Our CRISPR-based experiments indicate that such 3'aQTLs can alter APA regulation. Furthermore, we demonstrate that mapping 3'aQTLs can identify APA regulators, such as La-related protein 4. Finally, 3'aQTLs are colocalized with approximately 16.1% of trait-associated variants and are largely distinct from other QTLs, such as expression QTLs. Together, our findings show that 3'aQTLs contribute substantially to the molecular mechanisms underlying human complex traits and diseases.

The changing landscape of HLA typing: Understanding how and when HLA typing data can be used with confidence from bench to bedside

Human leukocyte antigen (HLA) genes are extraordinary for their extreme diversity and widespread impact on human health and disease. More than 30,000 HLA alleles have been officially named and more alleles continue to be discovered at a rapid pace. HLA typing systems which have been developed to detect HLA diversity have advanced rapidly and are revolutionizing our understanding of HLA's clinical importance. However, continuous improvements in knowledge and technology have created challenges for clinicians and scientists. This review explains how differences in HLA typing systems can impact the HLA types that are assigned. The consequences of differences in laboratory testing methods and reference databases are described. The challenges of using HLA types that are not equivalent are illustrated. A fundamental understanding of the continual expansion of our understanding of HLA diversity and limitations in some of the typing data is essential for using typing data appropriately in clinical and research settings.

A long road/read to rapid high-resolution HLA typing: The nanopore perspective

Next-generation sequencing (NGS) has been widely adopted for clinical HLA typing and advanced immunogenetics researches. Current methodologies still face challenges in resolving cis-trans ambiguity involving distant variant positions, and the turnaround time is affected by testing volume and batching. Nanopore sequencing may become a promising addition to the existing options for HLA typing. The technology delivered by the MinION sequencer of Oxford Nanopore Technologies (ONT) can record the ionic current changes during the translocation of DNA/RNA strands through transmembrane pores and translate the signals to sequence reads. It features simple and flexible library preparations, long sequencing reads, portable and affordable sequencing devices, and rapid, real-time sequencing. However, the error rate of the sequencing reads is high and remains a hurdle for its broad application. This review article will provide a brief overview of this technology and then focus on the opportunities and challenges of using nanopore sequencing for high-resolution HLA typing and immunogenetics research.

The Length of the Expressed 3' UTR Is an Intermediate Molecular Phenotype Linking Genetic Variants to Complex Diseases

In the last decades, genome-wide association studies (GWAS) have uncovered tens of thousands of associations between common genetic variants and complex diseases. However, these statistical associations can rarely be interpreted functionally and mechanistically. As the majority of the disease-associated variants are located far from coding sequences, even the relevant gene is often unclear. A way to gain insight into the relevant mechanisms is to study the genetic determinants of intermediate molecular phenotypes, such as gene expression and transcript structure. We propose a computational strategy to discover genetic variants affecting the relative expression of alternative 3′ untranslated region (UTR) isoforms, generated through alternative polyadenylation, a widespread posttranscriptional regulatory mechanism known to have relevant functional consequences. When applied to a large dataset in which whole genome and RNA sequencing data are available for 373 European individuals, 2,530 genes with alternative polyadenylation quantitative trait loci (apaQTL) were identified. We analyze and discuss possible mechanisms of action of these variants, and we show that they are significantly enriched in GWAS hits, in particular those concerning immune-related and neurological disorders. Our results point to an important role for genetically determined alternative polyadenylation in affecting predisposition to complex diseases, and suggest new ways to extract functional information from GWAS data.

Factors affecting HLA expression: A review

The detection and semiquantitative measurement of circulating human leucocyte antigen (HLA)-specific antibodies is essential for the management of patients before and after transplantation. In addition, the pretransplant cross-match to assess the reactivity of recipient HLA antibody against donor lymphocytes has long been the gold standard to prevent hyperacute rejection. Whilst both of these tests assume that recipient HLA-specific antibody is the only variable in the assessment of transplant risk, this is not the case. Transplant immunologists recognize that some HLA antigens are expressed at levels a magnitude lower than others (e.g., HLA-C, HLA-DQ), but within loci, and between different cell types there are many factors that influence HLA expression in both resting and activated cells. HLA is not usually expressed without the specific promoter proteins NLRC5, for HLA class I, and CIITA, for class II. The quantity of HLA protein production is then affected by factors including promoter region polymorphisms, alternative exon splice sites, methylation and microRNA-directed degradation. Different loci are influenced by multiple combinations of these control mechanisms making prediction of HLA regulation difficult, but an ability to measure the cellular expression of each HLA antigen, in conjunction with knowledge of circulating HLA-specific antibody, would lead to a more informed algorithm to assess transplant risk.

The immunogenetics of neurological disease

Genes encoding antigen-presenting molecules within the human major histocompatibility complex (MHC) account for the highest component of genetic risk for many neurological diseases, such as multiple sclerosis, neuromyelitis optica, Parkinson's disease, Alzheimer's disease, schizophrenia, myasthenia gravis and amyotrophic lateral sclerosis. Myriad genetic, immunological and environmental factors may contribute to an individual's susceptibility to neurological disease. Here, we review and discuss the decades long research on the influence of genetic variation at the MHC locus and the role of immunogenetic killer cell immunoglobulin-like receptor (KIR) loci in neurological diseases, including multiple sclerosis, neuromyelitis optica, Parkinson's disease, Alzheimer's disease, schizophrenia, myasthenia gravis and amyotrophic lateral sclerosis. The findings of immunogenetic association studies are consistent with a polygenic model of inheritance in the heterogeneous and multifactorial nature of complex traits in various neurological diseases. Future investigation is highly recommended to evaluate both coding and non-coding variation in immunogenetic loci using high-throughput high-resolution next-generation sequencing technologies in diverse ethnic groups to fully appreciate their role in neurological diseases.

The immunogenetics of multiple sclerosis: A comprehensive review

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system and common cause of non-traumatic neurological disability in young adults. The likelihood for an individual to develop MS is strongly influenced by her or his ethnic background and family history of disease, suggesting that genetic susceptibility is a key determinant of risk. Over 100 loci have been firmly associated with susceptibility, whereas the main signal genome-wide maps to the class II region of the human leukocyte antigen (HLA) gene cluster and explains up to 10.5% of the genetic variance underlying risk. HLA-DRB1*15:01 has the strongest effect with an average odds ratio of 3.08. However, complex allelic hierarchical lineages, cis/trans haplotypic effects, and independent protective signals in the class I region of the locus have been described as well. Despite the remarkable molecular dissection of the HLA region in MS, further studies are needed to generate unifying models to account for the role of the MHC in disease pathogenesis. Driven by the discovery of combinatorial associations of Killer-cell Immunoglobulin-like Receptor (KIR) and HLA alleles with infectious, autoimmune diseases, transplantation outcome and pregnancy, multi-locus immunogenomic research is now thriving. Central to immunity and critically important for human health, KIR molecules and their HLA ligands are encoded by complex genetic systems with extraordinarily high levels of sequence and structural variation and complex expression patterns. However, studies to-date of KIR in MS have been few and limited to very low resolution genotyping. Application of modern sequencing methodologies coupled with state of the art bioinformatics and analytical approaches will permit us to fully appreciate the impact of HLA and KIR variation in MS.

The use of haplotype-specific transcripts improves sample annotation consistency

Background

Exact sample annotation in expression microarray datasets is essential for any type of pharmacogenomics research.

Results

Candidate markers were explored through the application of Hartigans’ dip test statistics to a publically available human whole genome microarray dataset. The marker performance was tested on 188 serial samples from 53 donors and of variable tissue origin from five public microarray datasets. A qualified transcript marker panel consisting of three probe sets for human leukocyte antigens HLA-DQA1 (2 probe sets) and HLA-DRB4 identified sample donor identifier inconsistencies in six of the 188 test samples. About 3% of the test samples require root-cause analysis due to unresolvable inaccuracies.

Conclusions

The transcript marker panel consisting of HLA-DQA1 and HLA-DRB4 represents a robust, tissue-independent composite marker to assist control donor annotation concordance at the transcript level. Allele-selectivity of HLA genes renders them good candidates for “fingerprinting” with donor specific expression pattern.

EBP1 and DRBP76/NF90 binding proteins are included in the major histocompatibility complex class II RNA operon

Major histocompatibility complex class II mRNAs encode heterodimeric proteins involved in the presentation of exogenous antigens during an immune response. Their 3'UTRs bind a protein complex in which we identified two factors: EBP1, an ErbB3 receptor-binding protein and DRBP76, a double-stranded RNA binding nuclear protein, also known as nuclear factor 90 (NF90). Both are well-characterized regulatory factors of several mRNA molecules processing. Using either EBP1 or DRBP76/NF90-specific knockdown experiments, we established that the two proteins play a role in regulating the expression of HLA-DRA, HLA-DRB1 and HLA-DQA1 mRNAs levels. Our study represents the first indication of the existence of a functional unit that includes different transcripts involved in the adaptive immune response. We propose that the concept of 'RNA operon' may be suitable for our system in which MHCII mRNAs are modulated via interaction of their 3'UTR with same proteins.

Variation analysis and gene annotation of eight MHC haplotypes: the MHC Haplotype Project

The human major histocompatibility complex (MHC) is contained within about 4 Mb on the short arm of chromosome 6 and is recognised as the most variable region in the human genome. The primary aim of the MHC Haplotype Project was to provide a comprehensively annotated reference sequence of a single, human leukocyte antigen-homozygous MHC haplotype and to use it as a basis against which variations could be assessed from seven other similarly homozygous cell lines, representative of the most common MHC haplotypes in the European population. Comparison of the haplotype sequences, including four haplotypes not previously analysed, resulted in the identification of >44,000 variations, both substitutions and indels (insertions and deletions), which have been submitted to the dbSNP database. The gene annotation uncovered haplotype-specific differences and confirmed the presence of more than 300 loci, including over 160 protein-coding genes. Combined analysis of the variation and annotation datasets revealed 122 gene loci with coding substitutions of which 97 were non-synonymous. The haplotype (A3-B7-DR15; PGF cell line) designated as the new MHC reference sequence, has been incorporated into the human genome assembly (NCBI35 and subsequent builds), and constitutes the largest single-haplotype sequence of the human genome to date. The extensive variation and annotation data derived from the analysis of seven further haplotypes have been made publicly available and provide a framework and resource for future association studies of all MHC-associated diseases and transplant medicine.

Cloning and expression of a novel serine protease from Japanese flounder, Paralichthys olivaceus

Two different cDNA clones of Japanese flounder (types I-1 and I-2) with lengths of 1096 and 1572bp, respectively, were found to encode the same serine protease consisting of 244 identical amino acid residues with three putative N-glycosylation sites, an 18-amino acid signal peptide and a 2-amino acid activation peptide. The amino acid sequence of the Japanese flounder serine protease shares 39-44% identity to known hematopoietic serine proteases. Genomic analysis showed that two different clones were alternatively spliced from the same gene. A phylogenetic tree analysis showed that the Japanese flounder serine protease clustered with a hypothetical fugu protein and this cluster belonged to the neutrophil serine protease family cluster, which includes myeloblastin, N-elastase, and azurocidin. Expression of the Japanese flounder serine protease gene was observed to be up-regulated in head kidney cells after infection with Hirame rhabdovirus and LPS induction. In situ hybridization indicated that cells expressing Japanese flounder serine protease are different from CD8(+) and immunoglobulin(+) cells.

Position-dependent repression and promotion of DQB1 intron 3 splicing by GGGG motifs

Alternative splicing of HLA-DQB1 exon 4 is allele-dependent and results in variable expression of soluble DQbeta. We have recently shown that differential inclusion of this exon in mature transcripts is largely due to intron 3 variants in the branch point sequence (BPS) and polypyrimidine tract. To identify additional regulatory cis-elements that contribute to haplotype-specific splicing of DQB1, we systematically examined the effect of guanosine (G) repeats on intron 3 removal. We found that the GGG or GGGG repeats generally improved splicing of DQB1 intron 3, except for those that were adjacent to the 5' splice site where they had the opposite effect. The most prominent splicing enhancement was conferred by GGGG motifs arranged in tandem upstream of the BPS. Replacement of a G-rich segment just 5' of the BPS with a series of random sequences markedly repressed splicing, whereas substitutions of a segment further upstream that lacked the G-rich elements and had the same size did not result in comparable splicing inhibition. Systematic mutagenesis of both suprabranch guanosine quadruplets (G(4)) revealed a key role of central G residues in splicing enhancement, whereas cytosines in these positions had the most prominent repressive effects. Together, these results show a significant role of tandem G(4)NG(4) structures in splicing of both complete and truncated DQB1 intron 3, support position dependency of G repeats in splicing promotion and inhibition, and identify positively and negatively acting sequences that contribute to the haplotype-specific DQB1 expression.

Diversity of long-chain toxins in Tityus zulianus and Tityus discrepans venoms (Scorpiones, Buthidae): molecular, immunological, and mass spectral analyses

In Venezuela, stings by Tityus zulianus scorpions produce cardiorespiratory arrest, whereas envenoming by Tityus discrepans involves gastrointestinal/pancreatic complications, suggesting structural and/or functional differences. We sought to compare their toxin repertoires through immunological, molecular, and mass spectral analyses. First, in vivo tests showed that neutralization of T. zulianus venom toxicity by the anti-T. discrepans antivenom was not complete. To compare T. discrepans and T. zulianus long-chain (sodium channel-active) toxins, their most toxic Sephadex G-50 fractions, TdII and TzII, were subjected to acid-urea PAGE, which showed differences in composition. Amplification of toxin-encoding mRNAs using a leader peptide-based oligonucleotide rendered cDNAs representing twelve T. discrepans and two T. zulianus distinct toxin transcripts, including only one shared component, indicating divergence between T. zulianus and T. discrepans 5' region-encoded, toxin signal peptides. A 3'-UTR polymorphism was also noticed among the transcripts encoding shared components Tz1 and Td4. MALDI-TOF MS profiling of TdII and TzII produced species-specific spectra, with seven of the individual masses matching those predicted by cDNA sequencing. Phylogenetic analysis showed that the unique T. zulianus transcript-encoded sequence, Tz2, is structurally related to Tityus serrulatus and Centruroides toxins. Together with previous reports, this work indicates that T. zulianus and T. discrepans toxin repertoires differ structurally and functionally.

HLA-DQA1 Introns 2 and 3 Sequencing: DQA1 Sequencing-Based Typing and Characterization of a Highly Polymorphic Microsatellite at Intron 3 of DQA1*0505

DQA1 class II gene encodes the alpha-chain of the human leukocyte antigen (HLA)-DQ heterodimer. Sequencing-based typing (SBT) for HLA genes is the most powerful methodology described. However, most of the SBT procedures reported for HLA class II genes are not able to define complete exon 2 region. For that purpose, we have characterized introns 2 and 3 from most DQA1 alleles to design amplification procedures that were able to obtain complete exon 2 and 3 sequences from DQA1 genes. This coding information allowed us to reduce the number of ambiguities for DQA1 typing. DQA1 intron 2 and 3 characterization demonstrated the presence of two polymorphisms for alleles with the same exons 2 and 3 sequence from DQA1*05 group. Different samples including the DQA1*050101 alleles showed a single nucleotide polymorphism at position 53 of intron 2 (G53T). Additional haplotypic analysis showed the possible association of T53 allele with the Ax-Cw5-B18-DR17-DQ2 extended haplotype. On the other hand, DQA1*0505 sequencing from different control samples noticed the existence of a microsatellite (TTTC/AAAG)n located at position 126 of intron 3. Fragment length analysis demonstrated a high polymorphism for this short tandem repeat system (0505STR), defining alleles that ranged from 8 to 20 repetitions in our population.

A new splicing acceptor site and poly(A)+ sequence signal within DQA1*0401 and DQA1*0501 mRNA 3'UTR contribute to increase the extraordinary diversity of mRNA isoforms

In this paper, we have analysed the diversity of mRNA species generated by DQA1*0501 and DQA1*0401 alleles in homozygous B-lymphoblastoid cell lines. As we previously reported, six mRNA isoforms that differ in the 3'UTR have been identified in these cells. This diversity of mRNA species results both from the alternative use of two acceptor spliced sites and the differential selection of two poly(A)+ sequence signals by the processing machinery. In this report we describe a new acceptor sequence signal that allows generation of a new alternative spliced mRNA species. This acceptor sequence signal was also present in all of the seven DQA1 homozygous cell lines analysed. In addition, we have identified a previously undetected, non-conventional but functional, poly(A)+ sequence signal that lacks an identifiable AATAAA hexamer, one of the most important element of the core. This sequence signal allows the generation of two additional mRNA isoforms both in DQA1*0501 and DQA1*0401 homozygous cell lines but not in the others. We show that DQA1*0501 and DQA1*0401 primary transcripts can be processed into nine mRNA isoforms that differ in the 3'UTR. Finally, we summarized all the DQA1 mRNA species deriving from DQA1*0101, DQA1*0102, DQA1*0103, DQA1*0201, DQA1*0301, DQA1*0401 and DQA1*0501 alleles and shown in B-lymphoblastoid cell lines.