Sunday, September 23, 2018

Is the NorCal Rapist Arrest Another Success Story for Forensic Genomic Genealogy?

On the morning of September 20, police arrested Roy Charles Waller, 58, of Benicia, California, for 10 rapes over a 15-year period starting in 1991. A Washington Post article described the crucial lead in the case as "what's called a familial DNA search." According to the Post,
    The arrest was reminiscent of that of Joseph James DeAngelo, who was arrested in April on the suspicion that he was the so-called “Golden State Killer,” who was wanted for raping dozens of women and killing at least 12 people in a bloody swath of crime that spanned decades in the state.
    Like DeAngelo, Waller was arrested after police searched the genealogy site GEDmatch for leads. In DeAngelo’s case, officials did what’s called a familial DNA search of GEDmatch, in which they sought to find someone who was closely genetically related to and worked backward to find a suspect. Familial DNA searching, particularly as it relates to government-run DNA databases, has come into wider use around the country, but it raises complicated questions about whether it means that the privacy rights of people are forfeited, in effect, by the decisions made by their relatives.
    It was not immediately clear if police did a familial search in Waller’s case or found his profile in GEDmatch.
In forensic genetics, "familial searching" normally refers to trawling a law-enforcement database of DNA profiles after failing to find a match between a profile in the database and a profile from a crime-scene sample at a mere 20 or so STRs (short tandem repeats of a short string of base pairs) in the genome of billions of base pairs. But a near miss to a profile in the database might be the result of kinship -- the database inhabitant does not match, but a very close relative outside the database does. A more precise description of such "familial searching" is outer-directed kinship trawling of law-enforcement databases (Kaye 2013).

The Post's suggestion that the NorCal Rapist may have been identified through "familial DNA searching" and the statement that "[f]amilial DNA searching, particularly as it relates to government-run DNA databases ... raises complicated questions about whether it means that the privacy rights of people are forfeited, in effect, by the decisions made by their relatives" is badly confused. The inhabitants, so to speak, of the law-enforcement databases, have not elected to be there. Their relatives, near or distant, have not chosen to enroll in the database. DNA profiles are in these databases because of a conviction or, sometimes, an arrest for certain crimes. The profiles are not useful for much of anything other than personal identification and kinship testing for close relatives (lineal and siblings).

The "genealogy site GEDmatch" is very different. It is populated by individuals who have elected to make their DNA searchable for kinship by other people who are looking for possible biological relatives. These trawls do not use the limited STR profiles developed for both inner- and outer-directed trawls of offender databases. They use the much more extensive set of SNPs (single-nucleotide polymorphisms) -- hundreds of  thousands of them -- that genetic testing companies such as 23-and-me provide to consumers who are curious about their traits and ancestry. Those data can be used to infer more attenuated kinship. Knowing that long, shared stretches of DNA in the crime-scene sample and the sample in a private (that is, nongovernmental) genealogy database such as GEDmatch could reflect a common ancestor several generations back enables genealogists using public information about families to walk up the family tree to a putative common ancestor and then down again to living relatives.

It is clear enough that California officials used the latter kind of kinship trawling in the publicly accessible genealogy database GEDmatch to find Mr. Waller. At a half-an-hour long press conference, Sacramento County  District Attorney Anne Marie Schubert refused to give direct answers to questions such as "Did a relative upload something?" and did you have to go through a lot of family members to get to him? However, she did state that "[t]he link was made through genetic genealogy through the use of GEDmatch." Eschewing details, she repeated "It's genetic genealogy, that's what I'll say." If Mr. Waller had put his own genome-wide scan on GEDmatch, why would Ms. Schubert add "a lot of kudos to the folks in our office ... that I call the experts in tree building"? The district attorney's desire not to say anything about specific family members reflects a commendable concern about unnecessarily divulging information about the family, but she could have been more transparent about the investigative procedure without revealing any particularly sensitive information on the family or any material that would interfere with the prosecution.

REFERENCES

Friday, September 21, 2018

Forensic Genomic Genealogy and Comparative Justice

Thirty years ago, the introduction of law-enforcement DNA databases for locating the sources of DNA samples recovered from crime-scenes and victims was greeted with unbridled enthusiasm from some quarters and deep distrust from others.  So too, reactions to the spate of recent arrests in cold cases made possible by forays into DNA databases created for genealogy research have ranged from visions of a golden era for police investigations to glimpses into a dark and dystopian future.

Particularly in the law-enforcement database context, one concern has been the disproportionate impact of confining DNA databases to profiles of individuals who have been arrested for or convicted of crimes. For a variety of reasons, racial minorities tend to be overrepresented in law-enforcement databases (as compared to their percentage of the general population). 1/ But it seems most unlikely that nonwhites are similarly concentrated in the private genealogy databases that have resulted from the growth of recreational genetics. It may be, as Peter Neufeld observed when interviewed about the Golden State Killer arrest, that "[t]here is a whole generation that says, ‘I don’t really care about privacy,’" 2/ but it seems odd to speak of minorities as being "disproportionately affected [by] the unintended consequences of this genetic data" in the newly exploited databases. 3/

If anything, to quote Professor Erin Murphy, "the racial composition of recreational DNA sites -- which heavily skew white -- may end up complementing and balancing that of government databases, which disproportionately contain profiles from persons of color." 4/ That is not much of an argument for widespread forensic genealogy (and Professor Murphy did not rely on it for that purpose). Given how labor intensive forensic genomic genealogy is for genealogists and police, it seems unlikely that the technique for developing investigative leads to distant relatives will be used often enough to produce or correct massive disparities in who is subject to arrest or conviction.

Still, if police routinely were able to obtain complete results on crime-scene DNA with the DNA chips used in genome-wide association studies and recreational genetics, they could easily check whether any of the DNA records on those databases are immediate matches (or indicative of close relatives who might be tracked down without too much effort). The database used in the Golden State Killer case, GEDmatch, has data from a million or so curious individuals in it. That is considerably less than the 16 or 17 million profiles in the FBI's national DNA database (NDIS), but it is far from insignificant.

NOTES
  1. David H. Kaye & Michael Smith, DNA Identification Databases: Legality, Legitimacy, and the Case for Population-Wide Coverage, 2003 Wisc. L. Rev. 41.
  2. Gina Kolata & Heather Murphy, The Golden State Killer Is Tracked Through a Thicket of DNA, and Experts Shudder, N.Y. Times, Apr. 27, 2018 (quoting Peter Neufeld).
  3. Id.
  4. Erin Murphy, Law and Policy Oversight of Familial Searches in Recreational Genealogy Databases, Forensic Sci. Int’l (2018) (in press).

Monday, September 17, 2018

P-values in "the home for big data in biology"

A p-value indicates the significance of the difference in frequency of the allele tested between cases and controls i.e. the probability that the allele is likely to be associated with the trait. 1/
Such is the explanation of "p-value" provided by "the EMBL-EBI Training Programme [which] is celebrating 10 amazing years of providing onsite, offsite and online training in bioinformatics" for "scientists at all levels." 2/ The explanation appears in the answer to the question "What are genome wide association studies (GWAS)?" 3/ It comes from people in the know -- "The home for big data in biology" 4/ at "Europe's flagship laboratory for the life sciences." 5/

Above the description is a Manhattan plot of the p-values for the differences between the frequencies of the single nucleotide alleles comprising a genome-wide SNP array in samples of "cases" (individuals with a trait) and "controls" (individuals without the trait). Sadly, the p-values in the plot cannot be equated to "the probability that the allele is likely to to be associated with the trait."

I. Transposition

To begin with, a p-value is the probability of a difference at least as large as the one observed in the sample of "cases" and the sample of "controls" -- when the probability that a random person has the allele is the same for both cases and controls in the entire population. For a single allele, the expected value of the difference in the sample proportions is zero, but the observed value will vary from one pair of samples to the next. Because most pairs will have small differences, a big difference in a single study is evidence against the hypothesis that there is no difference at the population level. Differences that rarely would occur by chance alone are indicative of a true association. They are not usually false positives. At least, that is the theory behind the p-value.

For example, a p-value of 1/100,000, is a rare occurrence when the two samples come from a population in which the probability of the trait is the same with or without the allele. Consequently, the difference that corresponds to this p-value is thought to be strong evidence that the allele really is more (or less) common among cases than controls.

But even if the reasoning that "we would not expect it if H is true, therefore H is likely to be false" is correct, the p-value of 1/100,000 is not "the probability that the allele is likely to to be associated with the trait." It is the probability of that sort of a discrepancy if the allele has absolutely no association with the trait. In contrast, the probability of "no association" is not even defined in the statistical framework that gives us p-values.

Another way to say it: The p-value is a statement about the evidence (the observed difference) given the hypothesis of no association. It does not represent the probability of the hypothesis of zero association given an observed association. Equating the probability associated with the evidence with the probability associated with the hypothesis is so common that it has a name -- the transposition fallacy.

II. Multiple Comparisons

A second defect of the EMBL-EBI's definition is that a p-value of, say, 1/100,000 is not a good measure of surprise for GWAS data. Suppose there are 500,000 SNPs in the array and none of their alleles has any true association with the trait. If all apparent associations are independent, the expected number with the p-value is 500,000 × (1/100,000) = 5. Because of the many opportunities for individually impressive differences to appear, it is no surprise that some alleles have this small p-value. The p-value would have to be much smaller than 1/100,000 for the apparent association to be as surprising as the reported p-value would suggest. P-values that would produce a reasonable false discovery rate could be very small indeed.

An oversimplified analogy 6/ is this: Flipping a fair coin 18 times and getting 18 heads or tails has a probability on the order of 1/100,000. 7/ Flipping 500,000 coins 18 times each and finding that some of these experiments yielded 18 heads or tails is not strong evidence against the proposition that all the coins are fair. It is just what we would expect to see if the coins are all fair.

Of course, the bioinformaticists at EMBL-EBI are acutely aware of the effect of multiple comparisons. Their catalog of findings only includes "variant-trait associations ... if they have a p-value <1.0 × 10-5 in the overall (initial GWAS + replication) population." 8/ But why insist on so small a number if this p-value is "the probability that the allele is likely to to be associated with the trait"? For multiple comparisons, the p-value of 1/100,000 is not the measure of surprise that it is supposed to be.

NOTES
  1. EMBL-EBI, What Are Genome Wide Association Studies (GWAS)?, https://www.ebi.ac.uk/training/online/course/gwas-catalog-exploring-snp-trait-associations/why-do-we-need-gwas-catalog/what-are-genom, last visited Sept. 16, 2018.
  2. EMBL-EBI, EMBLI-EBI Training, https://www.ebi.ac.uk/training/online/course/gwas-catalog-exploring-snp-trait-associations/why-do-we-need-gwas-catalog/what-are-genome, last visited Sept. 16, 2018.
  3. EMBL-EBI, supra note 1.
  4. EMBL-EBI, https://www.ebi.ac.uk/, last visited Sept. 16, 2017.
  5. Id. 
  6. It is oversimplified because not all associations are independent. Nearby SNPs tend to be inherited together, but methods that take account of dependencies and enable researchers to pick out the associations that are real have been studied. 
  7. The more exact probability is 1/131,072.
  8. EMBL-EBI, Where Does the Data Come From?, https://www.ebi.ac.uk/training/online/course/gwas-catalog-exploring-snp-trait-associations/what-gwas-catalog/where-does-data-come, last visited Sept. 17, 2018. The phrase "the overall (initial GWAS + replication) population" is puzzling. It sounds like the data from an exploratory study are combined with those from the replication study to give a p-value for a larger sample (not a population). If so, the p-values for each study could be more than 1/100,000.