Widespread non-coding polymorphism in HLA class II genes of International HLA and Immunogenetics Workshop cell lines.

As the primary genetic determinant of immune recognition of self and non-self, the hyperpolymorphic HLA genes play key roles in disease association and transplantation. The large, variably sized HLA class II genes have historically been less well characterized than the shorter HLA class I genes. Here, we have used Pacific Biosciences Single Molecule Real-Time (SMRT®) DNA sequencing to perform four-field resolution HLA typing of HLA-DRB1/3/4/5, -DQA1, -DQB1, -DPA1 and -DPB1 from a panel of 181 B-lymphoblastoid cell lines from the International HLA and Immunogenetics Workshops. By interrogating all exons, introns, and the untranslated regions of these important reference cells, we have improved their HLA typing resolution on the IPD-IMGT/HLA database. We observed widespread non-coding polymorphism, with over twice as many unique genomic sequences identified compared with coding sequences (CDS). We submitted 263 unique sequences to the IPD-IMGT/HLA Database, often from multiple cell lines, including 114 confirmations of existing alleles, of which 30 were also extensions to full-length genomic sequences where only CDS was available previously. A total of 149 novel alleles were identified, largely differing from their closest reference allele sequences by a single nucleotide polymorphism (SNP). However, some highly divergent alleles were deemed to be recombinants, only detectable by full-length sequencing with long, phased reads. The fourth-field variation we observed allowed fine mapping of linkage disequilibrium patterns and haplotypes to particular ancestries. This study has highlighted the under-appreciated non-coding diversity in HLA class II genes, with potential implications for population genetic and clinical studies.

The HLA diversity of the Anthony Nolan register.

While the success of allogeneic stem cell transplantation depends on a high degree of HLA compatibility between donor and patient, finding a suitable donor remains challenging due to the hyperpolymorphic nature of HLA genes. We calculated high-resolution allele, haplotype and phenotype frequencies for HLA-A, -C, -B, -DRB1 and -DQB1 for 10 subpopulations of the Anthony Nolan (AN) register using an in-house expectation-maximisation (EM) algorithm run on mixed resolution HLA data, covering 676 155 individuals. Sample sizes range from 599 410 for British/Irish North West European (BINWE) individuals, the largest subpopulation in the United Kingdom to 1105 for the British Bangladeshi population. Calculation of genetic distance between the subpopulations based on haplotype frequencies shows three broad clusters, each following a major continental group: European, African and Asian. We further analysed the HLA haplotype and phenotype diversity of each subpopulation, and found that 35.52% of BINWE individuals ranging to 98.34% of Middle Eastern individuals on the register had a unique phenotype within their subpopulation. These analyses and the allele, haplotype and phenotype frequency data of the subpopulation on the AN register are a valuable resource in understanding the HLA diversity in the United Kingdom and can be used to improve the accuracy of match likelihoods and to inform future donor recruitment strategies.

Single molecule real-time DNA sequencing of the full HLA-E gene for 212 reference cell lines.

We have developed a genotyping assay that produces fully phased, unambiguous HLA-E genotyping using Pacific Biosciences' single molecule real-time DNA sequencing. In total 212 cell lines were genotyped, including the panel of 107 established at the 10th International Histocompatibility Workshop. Our results matched the previously known HLA-E genotype in 94 (44.3%) cell lines, in all cases either improving or equalling previous genotyping resolution. Three (1.4%) cells had discrepant HLA-E genotyping data and 115 (54.2%) had no previous HLA-E data. The HLA-E genotypes for four (1.9%) cell lines resulted in a change of zygosity by identifying two distinct haplotypes. We discovered eight novel HLA-E alleles, extended the known reference sequence of seven and confirmed the existence of a further 10.

MicroRNA expression profiles in Kaposi's sarcoma.

Kaposi's sarcoma (KS) is a mesenchymal tumor, caused by Human herpesvirus 8 (HHV8) with molecular and cytogenetic changes poorly understood. To gain further insight on the underlying molecular changes in KS, we performed microRNA (miRNA) microarray analysis of 17 Kaposi's sarcoma specimens. Three normal skin specimens were used as controls. The most significant differentially expressed miRNA were confirmed by quantitative reverse transcriptase polymerase chain reaction (RT-PCR). We detected in KS versus normal skin 185 differentially expressed miRNAs, 76 were upregulated and 109 were downregulated. The most significantly downregulated miRNAs were miR-99a, miR-200 family, miR-199b-5p, miR-100 and miR-335, whereas kshv-miR-K12-4-3p, kshv-miR-K12-1, kshv-miR-K12-2, kshv-miR-K12-4-5p and kshv-miR-K12-8 were significantly upregulated. High expression levels of kshv-miR-K12-1 (p = 0.004) and kshv-miR-K12-4-3p (p = 0.001) was confirmed by RT-PCR. The predicted target genes for differentially expressed miRNAs included genes which are involved in a variety of cellular processes such as angiogenesis (i.e. THBS1) and apoptosis (i.e. CASP3, MCL1), suggesting a role for these miRNAs in Kaposi's sarcoma pathogenesis.

An integrated analysis of miRNA and gene copy numbers in xenografts of Ewing's sarcoma.

BACKGROUND: Xenografts have been shown to provide a suitable source of tumor tissue for molecular analysis in the absence of primary tumor material. We utilized ES xenograft series for integrated microarray analyses to identify novel biomarkers. METHOD: Microarray technology (array comparative genomic hybridization (aCGH) and micro RNA arrays) was used to screen and identify copy number changes and differentially expressed miRNAs of 34 and 14 passages, respectively. Incubated cells used for xenografting (Passage 0) were considered to represent the primary tumor. Four important differentially expressed miRNAs (miR-31, miR-31*, miR-145, miR-106) were selected for further validation by real time polymerase chain reaction (RT-PCR). Integrated analysis of aCGH and miRNA data was performed on 14 xenograft passages by bioinformatic methods. RESULTS: The most frequent losses and gains of DNA copy number were detected at 9p21.3, 16q and at 8, 15, 17q21.32-qter, 1q21.1-qter, respectively. The presence of these alterations was consistent in all tumor passages. aCGH profiles of xenograft passages of each series resembled their corresponding primary tumors (passage 0). MiR-21, miR-31, miR-31*, miR-106b, miR-145, miR-150*, miR-371-5p, miR-557 and miR-598 showed recurrently altered expression. These miRNAS were predicted to regulate many ES-associated genes, such as genes of the IGF1 pathway, EWSR1, FLI1 and their fusion gene (EWS-FLI1). Twenty differentially expressed miRNAs were pinpointed in regions carrying altered copy numbers. CONCLUSION: In the present study, ES xenografts were successfully applied for integrated microarray analyses. Our findings showed expression changes of miRNAs that were predicted to regulate many ES associated genes, such as IGF1 pathway genes, FLI1, EWSR1, and the EWS-FLI1 fusion genes.