Answers.Body: <p>To add to all the other great answers, I would mention that the question is somewhat misleading. If the reference genome is for a single individual, then it should be diploid. However, it's a reference for all humans. It should really contain billions of copies to fully account for all the diversity. Since that is not realistic, the reference serves as a simple approximation.</p>
  
  <p>This point was addressed by the recent <a href="https://www.nature.com/articles/ncomms13637" rel="nofollow noreferrer">Korean genome paper</a>:</p>
  
  <blockquote>
    <p>Human genomes are routinely compared against a universal reference.
    However, this strategy could miss population-specific and personal
    genomic variations, which may be detected more efficiently using an
    ethnically relevant or personal reference. ... Systematic comparison
    of human assemblies shows the importance of assembly quality,
    suggesting the necessity of new technologies to comprehensively map
    ethnic and personal genomic structure variations.</p>
  </blockquote>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-26T20:07:16.807Z
Answers.LastActivityDate: 2017-06-26T20:07:16.807Z
Answers.UserDisplayName: burger
Answers.Body: <blockquote>
    <p>At the moment, the standard reference genomes (e.g. hg19, hg38) are haploid genomes. We know that the human genome is diploid. Naturally, the latter would be the respectively correct representation of the human genome.</p>
  </blockquote>
  
  <p>The premise of the OP's question is false.  The natural reference representation of the human genome is not diploid.</p>
  
  <h3>Think of a reference genome as a map, and not as a specific example of a human being's DNA.</h3>
  
  <p>Not only is the human genome reference haploid, but <a href="https://en.wikipedia.org/wiki/Human_genome#Human_reference_genome" rel="nofollow noreferrer">it is also a composite genome</a>. This means that the human genome reference sequence is composed of sequences from multiple individuals.  In other words, the human reference does not correspond to any one human sequence. </p>
  
  <p>Any particular read from a DNA sequencer will be a read from a human genome that diverges from the reference genome. So an algorithm that tries to match the read to the reference genome will always need to handle potential discrepancies. Adding a second map against which to match a read would not change that fact.  Therefore, there is little value in providing a second haploid reference genome.</p>
  
  <p><strong>Side note</strong>:
  There are parts of the human genome that "are too complex to be represented by a single path", and the <a href="https://www.ncbi.nlm.nih.gov/grc/human" rel="nofollow noreferrer">Genome Reference Consortium</a> provides "alternate loci" for such regions of the genome.</p>
  
  <h3>Answer to the OP's second question</h3>
  
  <blockquote>
    <p>More importantly, how have large-scale genomic studies dealt with the fact that the haploid reference genome is a consensus-based "half" of a human genome?</p>
  </blockquote>
  
  <p>I interpret this question as "How do large genomic studies represent genome diversity that cannot be represented by a haploid reference genome?" </p>
  
  <p>One standard approach that has been very popular in genome-wide association studies is the use of a haplotype reference panel such as for example from the <a href="http://www.internationalgenome.org/about" rel="nofollow noreferrer">1000 Genomes Project</a>.  Modern phasing programs exploit the shared ancestry of the samples in a study and reference haplotypes to phase a study's samples. One popular phasing program is <a href="https://mathgen.stats.ox.ac.uk/impute/impute.html" rel="nofollow noreferrer">Impute2</a></p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-26T20:13:20.157Z
Answers.Text: Although technically composite, 70% of GRCh38 is from RP11, a single male donor (ref: http://genome.cshlp.org/content/27/5/849.full).
Answers.UserDisplayName: burger
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-26T07:23:26.087Z
Answers.LastActivityDate: 2018-01-25T10:57:23.753Z
Answers.LastEditDate: 2018-01-25T10:57:23.753Z
Answers.UserDisplayName: winni2k
Answers.Body: <p>There are some assemblers that produce assembly graphs that attempt to describe all the possible haploid paths within a set of reads. Such an assembly attempts to capture all the diploid variation (and/or population variation) in a sample at the expense of not having full-length chromosomes.</p>
  
  <p>Canu (for example) will produce contigs that are extended as long as consensus is maintained across different reads, but when there is a reliable break in coverage (i.e. an area where chromosomes are heterozygous) then the contigs will be broken up. Canu provides as output a GFA file (assembly graph) that can be used to determine which paths might combine together into a single chromosome.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-25T03:48:52.500Z
Answers.LastActivityDate: 2017-06-25T03:48:52.500Z
Answers.UserDisplayName: gringer
Answers.Body: <p>Since nobody has tackled you second question, I'll give it a go:</p>
  
  <p>The complementary strand to a genome sequence is the  5'-3' <em>reverse</em> complement, not the 3'-5' complement. This sequence is not captured in standard reference genomes, but the information is. Given a sequence, it's trivial to calculate it's reverse complement, and any tool that's designed to work with sequence data takes this feature into account.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2018-01-25T12:41:34.407Z
Answers.LastActivityDate: 2018-01-25T12:41:34.407Z
Answers.UserDisplayName: heathobrien
Answers.Body: <p>For calling small variants, the <em>standard way</em> is to simply call diploid genotypes. You can already do a variety of research with unphased genotypes. You may further phase genotypes with imputation, pedigree or with long reads/linked reads, but not many are doing this because phasing is more difficult, may add cost and may not always give you new insight into your data. For these analyses, we use a haploid genome. For human samples, the vast majority of "large-scale genomic studies" are done this way.</p>
  
  <p>A diploid reference actually doesn't help much with reference-based analysis; it only complicates algorithms. What could help a lot is a population reference, which may be represented by a graph or a compressed full-text index or both. In theory, if you have a comprehensive population reference and a capable mapping algorithm, you may call extra variants that would not be callable with short reads. In practice, however, there are quite a few technical challenges. Handling population references is a research topic. There are no "standards" yet.</p>
  
  <p>If the goal is to assemble a new reference genome from a diploid sample, we almost always prefer to produce a diploid assembly. Unfortunately, I believe there are no "standard" procedures, either. SuperNova from 10x genomics builds the diploid information into a graph. Falcon from PacBio uses "unzip". I don't think they have got widely used and evaluated so far.</p>
  
  <p>PS: saw your edit while writing the above. The fact that the genome only represents one strand does not mean we have to create the complement strand explicitly in analyses. We do most of reverse complement on the fly in algorithms as well as in mind.</p>
  
  <blockquote>
    <p>reference genomes are not truly haploid</p>
  </blockquote>
  
  <p>That depends on how the reference is assembled. If you sequence a haploid sample (e.g. bacteria), your assembly will be haploid. If you sequence an inbreed lab strain that is almost homozygous (e.g. mouse and fruit fly), your assembly will be nearly haploid. If you sequence a diploid sample, your assembly is very likely to be a mosaic of the two haplotypes. In case of the human reference genome, it is more complicated. It is largely a mosaic of several humans by stitching ~150kb haplotypes from these samples.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2018-01-25T10:32:59.113Z
Answers.Text: Sorry, but lots of groups phase their samples. That is because phasing of genotypes is a prerequisite for state of the art genotype imputation from a haplotype reference panel.  The paper of one such panel has 3437 citations according to the publisher's website (https://www.nature.com/articles/nature11632)
Answers.UserDisplayName: winni2k
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2018-01-25T13:09:48.443Z
Answers.Text: 1000g just *had to* do phasing. It didn't really have a choice. However, not many projects follow the 1000g design nowadays. I was referring to the new genome sequencing projects.
Answers.UserDisplayName: user172818
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2018-01-25T13:27:43.653Z
Answers.Text: I am not talking about 1000g. I am talking about the large subset of the 3437 studies that used the 1000g (phase1) reference panel. And then there are the 224 citations of the 2015 paper on 1000g phase3 https://www.nature.com/articles/nature15394, and the 63 citations of the 2016 Haplotype Rereferenc Consortium paper https://www.nature.com/articles/ng.3643. I think it really depends on what you are trying to achieve. For large GWAS, SNP chips + phasing + imputation are still the way to go. That may of course be different in your field.  Full disclosure: I am an author on the latter two papers.
Answers.UserDisplayName: winni2k
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2018-01-26T15:12:38.713Z
Answers.Text: Note the context of my sentence: "for calling small variants". It is just an example that you can work with genotypes directly. I of course know GWAS uses imputation all the time.
Answers.UserDisplayName: user172818
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-06-25T02:00:16.993Z
Answers.LastActivityDate: 2017-06-25T02:00:16.993Z
Answers.UserDisplayName: user172818
Body: <p>At the moment, the standard reference genomes (e.g. hg19, hg38) are haploid genomes. We know that the human genome is diploid. Naturally, the latter would be the respectively correct representation of the human genome. </p>
  
  <p>More and more biologists are using emerging technologies into order to capture the diploid nature of genetic information, e.g. phasing SNPs between mother and father chromosomes. </p>
  
  <p>What is the standard way that bioinformaticians generated a standard diploid reference genome? </p>
  
  <p>Actually, reference genomes are not truly haploid (by my understanding). Given that reference genomes are 5'-3' coordinated, in order to create a complementary strand, one would need to take the 3'-5' complement. In order to have a diploid genome, you need two reference genomes and two 3'-5' complements.</p>
  
  <p>More importantly, how have large-scale genomic studies dealt with the fact that the haploid reference genome is a consensus-based "half" of a human genome? </p>
  
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-06-25T06:59:02.577Z
Comments.Text: [This  BioStar](https://www.biostars.org/p/3423/) post might clarify bit the confusion about orientation.
Comments.UserDisplayName: Kamil S Jaron
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-06-26T08:40:36.780Z
Comments.Text: When you say "work with", are you only thinking of variant calling?
Comments.UserDisplayName: terdon
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2018-01-25T10:46:54.853Z
Comments.Text: This post is asking two questions. Perhaps it would be worth splitting the post into two?
Comments.UserDisplayName: winni2k
Tags: <human-genome><genome><reference-genome>
Title: What is the standard way to work with a diploid reference genome? Complementary strands?
ThreadText: # What is the standard way to work with a diploid reference genome? Complementary strands?
  
  (Asked by: username_0 on 2017-06-24)
  
  At the moment, the standard reference genomes (e.g. hg19, hg38) are haploid genomes. We know that the human genome is diploid. Naturally, the latter would be the respectively correct representation of the human genome.
  
  More and more biologists are using emerging technologies into order to capture the diploid nature of genetic information, e.g. phasing SNPs between mother and father chromosomes.
  
  What is the standard way that bioinformaticians generated a standard diploid reference genome?
  
  Actually, reference genomes are not truly haploid (by my understanding). Given that reference genomes are 5'-3' coordinated, in order to create a complementary strand, one would need to take the 3'-5' complement. In order to have a diploid genome, you need two reference genomes and two 3'-5' complements.
  
  More importantly, how have large-scale genomic studies dealt with the fact that the haploid reference genome is a consensus-based "half" of a human genome?
  
  **Tags:** human-genome, genome, reference-genome
  
  > **Comments:**
  > **username_1**: [This  BioStar](https://www.biostars.org/p/3423/) post might clarify bit the confusion about orientation.  
  > **username_2**: When you say "work with", are you only thinking of variant calling?  
  > **username_3**: This post is asking two questions. Perhaps it would be worth splitting the post into two?  
  
  ---
  
  ## 5 answers
  
  ### Answer by username_7 on 2017-06-25. Score: 5 | Accepted answer
  
  For calling small variants, the *standard way* is to simply call diploid genotypes. You can already do a variety of research with unphased genotypes. You may further phase genotypes with imputation, pedigree or with long reads/linked reads, but not many are doing this because phasing is more difficult, may add cost and may not always give you new insight into your data. For these analyses, we use a haploid genome. For human samples, the vast majority of "large-scale genomic studies" are done this way.
  
  A diploid reference actually doesn't help much with reference-based analysis; it only complicates algorithms. What could help a lot is a population reference, which may be represented by a graph or a compressed full-text index or both. In theory, if you have a comprehensive population reference and a capable mapping algorithm, you may call extra variants that would not be callable with short reads. In practice, however, there are quite a few technical challenges. Handling population references is a research topic. There are no "standards" yet.
  
  If the goal is to assemble a new reference genome from a diploid sample, we almost always prefer to produce a diploid assembly. Unfortunately, I believe there are no "standard" procedures, either. SuperNova from 10x genomics builds the diploid information into a graph. Falcon from PacBio uses "unzip". I don't think they have got widely used and evaluated so far.
  
  PS: saw your edit while writing the above. The fact that the genome only represents one strand does not mean we have to create the complement strand explicitly in analyses. We do most of reverse complement on the fly in algorithms as well as in mind.
  
  > reference genomes are not truly haploid
  
  That depends on how the reference is assembled. If you sequence a haploid sample (e.g. bacteria), your assembly will be haploid. If you sequence an inbreed lab strain that is almost homozygous (e.g. mouse and fruit fly), your assembly will be nearly haploid. If you sequence a diploid sample, your assembly is very likely to be a mosaic of the two haplotypes. In case of the human reference genome, it is more complicated. It is largely a mosaic of several humans by stitching ~150kb haplotypes from these samples.
  
  > **Comments:**
  > **username_3**: Sorry, but lots of groups phase their samples. That is because phasing of genotypes is a prerequisite for state of the art genotype imputation from a haplotype reference panel.  The paper of one such panel has 3437 citations according to the publisher's website (https://www.nature.com/articles/nature11632)  
  > **username_7**: 1000g just *had to* do phasing. It didn't really have a choice. However, not many projects follow the 1000g design nowadays. I was referring to the new genome sequencing projects.  
  > **username_3**: I am not talking about 1000g. I am talking about the large subset of the 3437 studies that used the 1000g (phase1) reference panel. And then there are the 224 citations of the 2015 paper on 1000g phase3 https://www.nature.com/articles/nature15394, and the 63 citations of the 2016 Haplotype Rereferenc Consortium paper https://www.nature.com/articles/ng.3643. I think it really depends on what you are trying to achieve. For large GWAS, SNP chips + phasing + imputation are still the way to go. That may of course be different in your field.  Full disclosure: I am an author on the latter two papers.  
  > **username_7**: Note the context of my sentence: "for calling small variants". It is just an example that you can work with genotypes directly. I of course know GWAS uses imputation all the time.  
  
  ---
  
  ### Answer by username_3 on 2017-06-26. Score: 6
  
  > At the moment, the standard reference genomes (e.g. hg19, hg38) are haploid genomes. We know that the human genome is diploid. Naturally, the latter would be the respectively correct representation of the human genome.
  
  The premise of the OP's question is false. The natural reference representation of the human genome is not diploid.
  
  ### Think of a reference genome as a map, and not as a specific example of a human being's DNA.
  
  Not only is the human genome reference haploid, but [it is also a composite genome](https://en.wikipedia.org/wiki/Human_genome#Human_reference_genome). This means that the human genome reference sequence is composed of sequences from multiple individuals. In other words, the human reference does not correspond to any one human sequence.
  
  Any particular read from a DNA sequencer will be a read from a human genome that diverges from the reference genome. So an algorithm that tries to match the read to the reference genome will always need to handle potential discrepancies. Adding a second map against which to match a read would not change that fact. Therefore, there is little value in providing a second haploid reference genome.
  
  **Side note**:
  There are parts of the human genome that "are too complex to be represented by a single path", and the [Genome Reference Consortium](https://www.ncbi.nlm.nih.gov/grc/human) provides "alternate loci" for such regions of the genome.
  
  ### Answer to the OP's second question
  
  > More importantly, how have large-scale genomic studies dealt with the fact that the haploid reference genome is a consensus-based "half" of a human genome?
  
  I interpret this question as "How do large genomic studies represent genome diversity that cannot be represented by a haploid reference genome?"
  
  One standard approach that has been very popular in genome-wide association studies is the use of a haplotype reference panel such as for example from the [1000 Genomes Project](http://www.internationalgenome.org/about). Modern phasing programs exploit the shared ancestry of the samples in a study and reference haplotypes to phase a study's samples. One popular phasing program is [Impute2](https://mathgen.stats.ox.ac.uk/impute/impute.html)
  
  > **Comments:**
  > **username_4**: Although technically composite, 70% of GRCh38 is from RP11, a single male donor (ref: http://genome.cshlp.org/content/27/5/849.full).  
  
  ---
  
  ### Answer by username_4 on 2017-06-26. Score: 3
  
  To add to all the other great answers, I would mention that the question is somewhat misleading. If the reference genome is for a single individual, then it should be diploid. However, it's a reference for all humans. It should really contain billions of copies to fully account for all the diversity. Since that is not realistic, the reference serves as a simple approximation.
  
  This point was addressed by the recent [Korean genome paper](https://www.nature.com/articles/ncomms13637):
  
  > Human genomes are routinely compared against a universal reference.
  > However, this strategy could miss population-specific and personal
  > genomic variations, which may be detected more efficiently using an
  > ethnically relevant or personal reference. ... Systematic comparison
  > of human assemblies shows the importance of assembly quality,
  > suggesting the necessity of new technologies to comprehensively map
  > ethnic and personal genomic structure variations.
  
  
  ---
  
  ### Answer by username_5 on 2017-06-25. Score: 2
  
  There are some assemblers that produce assembly graphs that attempt to describe all the possible haploid paths within a set of reads. Such an assembly attempts to capture all the diploid variation (and/or population variation) in a sample at the expense of not having full-length chromosomes.
  
  Canu (for example) will produce contigs that are extended as long as consensus is maintained across different reads, but when there is a reliable break in coverage (i.e. an area where chromosomes are heterozygous) then the contigs will be broken up. Canu provides as output a GFA file (assembly graph) that can be used to determine which paths might combine together into a single chromosome.
  
  
  ---
  
  ### Answer by username_6 on 2018-01-25. Score: 1
  
  Since nobody has tackled you second question, I'll give it a go:
  
  The complementary strand to a genome sequence is the 5'-3' *reverse* complement, not the 3'-5' complement. This sequence is not captured in standard reference genomes, but the information is. Given a sequence, it's trivial to calculate it's reverse complement, and any tool that's designed to work with sequence data takes this feature into account.
  
  
  ---
  
  