Haplotying of HLA loci with ultra-deep shotgun sequencing

Patent metadata

Field	Value
Publication number	US-9920370-B2
Application number	US-201514696278-A
Country	US
Kind code	B2
Filing date	Apr 24, 2015
Priority date	Jan 22, 2013
Publication date	Mar 20, 2018
Grant date	Mar 20, 2018

Abstract

Methods are provided to determine the entire genomic region of a particular HLA locus including both intron and exons. The resultant consensus sequences provides linkage information between different exons, and produces the unique sequence from each of the two genes from the individual sample being typed. The sequence information in intron regions along with the exon sequences provides an accurate HLA haplotype.

First claim

What is claimed is:

1. A method for determining a haplotype of an HLA locus, the method comprising:
amplifying an HLA gene to produce an amplified HLA gene from a sample comprising a source of genomic DNA and having a genotype, wherein said amplifying step comprises performing a long-range polymerase chain reaction (PCR); and wherein the amplified HLA gene is fragmented to produce a fragmented amplified HLA gene and ligated to barcode primers comprising: (i) a target specific identifier for the source of the genomic DNA; (ii) a target specific identifier for the HLA gene; and (iii) a sequencing adaptor;
sequencing the amplified HLA gene to produce a plurality of sequences; and
performing deconvolution analysis to resolve the haplotype of the locus, wherein deconvolution analysis comprises the steps of: (a) mapping the plurality of sequences to a reference chromatid generating a tiling pattern, wherein a continuous tiling pattern indicates a correct mapping of the sequences to the reference chromatid and a staggered tiling pattern indicates an incorrect mapping of the sequences to a reference chromatid; (b) filtering out reference chromatids with incorrect mapping;
and (c) assigning the reference chromatid with a correct mapping as the haplotype of the locus.

2. The method of claim 1 , wherein the HLA locus is an HLA Class I locus.

3. The method of claim 1 , wherein the fragmented amplified HLA gene is sequenced to a depth of at least 100 reads per sequence.

4. The method of claim 1 , wherein the fragmented amplified HLA gene is sequenced to a depth of at least 1000 reads per sequence.

5. The method of claim 4 , wherein the at least 1000 reads per sequence are mapped to a reference chromatid.

6. The method of claim 5 , wherein the step of mapping to a reference chromatid comprises the steps of:
counting a number of “central reads” for any given point, where central reads are empirically defined as mapped reads for which a ratio between a length of a left arm and that of a right arm related to a particular point is between 0.5 and 2;
computing the minimum coverage of overall reads (MOOR) and the minimum coverage of central reads (MCCR) for each reference;
eliminating references with an MOOR less than 20 and an MCCR less than 10;
enumerating all possible combinations of either one reference (homozygous allele) or two references (heterozygous alleles) of the same locus from the remaining references;
counting a number of distinct reads that map to each combination; wherein the combination with a maximum number of distinct reads is assigned as the genotype of that sample.

7. The method of claim 1 , wherein the amplified HLA gene comprises at least two exons and at least two introns.

Assignees

• Univ Leland Stanford Junior

Inventors

• Wang Chunlin
• Mindrinos Michael N
• Davis Mark M
• Davis Ronald W
• Krishnakumar Sujatha

Classifications

• C12Q2600/156

  Polymorphic or mutational markers · CPC title

• C12Q2600/172

  Haplotypes · CPC title

• C12Q1/6881 Primary

  for tissue or cell typing, e.g. human leukocyte antigen [HLA] probes · CPC title

• C12Q2600/118

  Prognosis of disease development · CPC title

• G06F19/18

  Physics · mapped topic