Oryza sativa or Asian cultivated rice is one of the major cereal grass species domesticated for human food use during the Neolithic. Domestication of this species from the wild grass Oryza rufipogon was accompanied by changes in several traits, including seed shattering, percent seed set, tillering, grain weight, and flowering time. Quantitative trait locus (QTL) mapping has identified three genomic regions in chromosome 3 that appear to be associated with these traits. We would like to study whether these regions show signatures of selection and whether the same genetic basis underlies the domestication of different rice varieties. Fragments of 88 genes spanning these three genomic regions were sequenced from multiple accessions of two major varietal groups in O. sativa--indica and tropical japonica--as well as the ancestral wild rice species O. rufipogon. In tropical japonica, the levels of nucleotide variation in these three QTL regions are significantly lower compared to genome-wide levels, and coalescent simulations based on a complex demographic model of rice domestication indicate that these patterns are consistent with selection. In contrast, there is no significant reduction in nucleotide diversity in the homologous regions in indica rice. These results suggest that there are differences in the genetic and selective basis for domestication between these two Asian rice varietal groups.

Analysis of polymorphism and divergence in the non-coding portion of the human genome yields crucial information about factors driving the evolution of gene regulation. Candidate cis-regulatory regions spanning more than 15,000 genes in 15 African Americans and 20 European Americans were re-sequenced and aligned to the chimpanzee genome in order to identify potentially functional polymorphism and to characterize and quantify departures from neutral evolution. Distortions of the site frequency spectra suggest a general pattern of selective constraint on conserved non-coding sites in the flanking regions of genes (CNCs). Moreover, there is an excess of fixed differences that cannot be explained by a Gamma model of deleterious fitness effects, suggesting the presence of positive selection on CNCs. Extensions of the McDonald-Kreitman test identified candidate cis-regulatory regions with high probabilities of positive and negative selection near many known human genes, the biological characteristics of which exhibit genome-wide trends that differ from patterns observed in protein-coding regions. Notably, there is a higher probability of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, suggesting that a larger portion of adaptive regulatory changes has occurred in genes expressed during brain development. Overall we find that natural selection has played an important role in the evolution of candidate cis-regulatory regions throughout hominid evolution.

Quantifying the distribution of fitness effects among newly arising mutations in the human genome is key to resolving important debates in medical and evolutionary genetics. Here, we present a method for inferring this distribution using Single Nucleotide Polymorphism (SNP) data from a population with non-stationary demographic history (such as that of modern humans). Application of our method to 47,576 coding SNPs found by direct resequencing of 11,404 protein coding-genes in 35 individuals (20 European Americans and 15 African Americans) allows us to assess the relative contribution of demographic and selective effects to patterning amino acid variation in the human genome. We find evidence of an ancient population expansion in the sample with African ancestry and a relatively recent bottleneck in the sample with European ancestry. After accounting for these demographic effects, we find strong evidence for great variability in the selective effects of new amino acid replacing mutations. In both populations, the patterns of variation are consistent with a leptokurtic distribution of selection coefficients (e.g., gamma or log-normal) peaked near neutrality. Specifically, we predict 27-29% of amino acid changing (nonsynonymous) mutations are neutral or nearly neutral (|s|<0.01%), 30-42% are moderately deleterious (0.01%<|s|<1%), and nearly all the remainder are highly deleterious or lethal (|s|>1%). Our results are consistent with 10-20% of amino acid differences between humans and chimpanzees having been fixed by positive selection with the remainder of differences being neutral or nearly neutral. Our analysis also predicts that many of the alleles identified via whole-genome association mapping may be selectively neutral or (formerly) positively selected, implying that deleterious genetic variation affecting disease phenotype may be missed by this widely used approach for mapping genes underlying complex traits.

State-of-the-art health care includes genome sequencing of the patient to identify genetic variants that contribute to either the cause of their malady or variants that can be targeted to improve treatment. The goal was to introduce state-of-the-art health care to cats using genomics and a precision medicine approach. To test the feasibility of a precision medicine approach in domestic cats, a single cat that presented to the University of Missouri, Veterinary Health Center with an undiagnosed neurologic disease was whole-genome sequenced. The DNA variants from the cat were compared to the DNA variant database produced by the 99 Lives Cat Genome Sequencing Consortium. Approximately 25× genomic coverage was produced for the cat. A predicted p.H441P missense mutation was identified in NPC1, the gene causing Niemann-Pick type C1 on cat chromosome D3.47456793 caused by an adenine-to-cytosine transversion, c.1322A>C. The cat was homozygous for the variant. The variant was not identified in any other 73 domestic and 9 wild felids in the sequence database or 190 additionally genotyped cats of various breeds. The successful effort suggested precision medicine is feasible for cats and other undiagnosed cats may benefit from a genomic analysis approach. The 99 Lives DNA variant database was sufficient but would benefit from additional cat sequences. Other cats with the mutation may be identified and could be introduced as a new biomedical model for NPC1. A genetic test could eliminate the disease variant from the population.

Tabby patterns of fur coats are defining characteristics in wild and domestic felids. Historically, three autosomal alleles at one locus (Tabby): Abyssinian (T^a ; a.k.a. ticked), mackerel (T^m ; a.k.a. striped) and blotched (t^b ; a.k.a. classic, blotched) were thought to control these patterns in domestic cats and their breeds. Currently, at least three loci influence cat tabby markings, two of which are designated Tabby and Ticked. The Tabby locus is laeverin (LVRN) and affects the mackerel and blotched patterns. The unidentified gene for the Ticked locus on cat chromosome B1 was suggested to control the presence or absence of the ticked pattern (Tabby - Abyssinian (T^a ; a.k.a. ticked). The cat reference genome (Cinnamon, the Abyssinian) has the ticked phenotype and the variant dataset and coat phenotypes from the 99 Lives Cat Genome Consortium (195 cats) were used to identify candidate genes and variants associated with the Ticked locus. Two strategies were used to find the Ticked allele(s), one considered Cinnamon with the reference allele or heterozygous (Strategy A) and the other considered Cinnamon as having the variant allele or heterozygous (Strategy B). For Strategy A, two variants in Dickkopf Wnt Signaling Pathway Inhibitor 4 (DKK4), a p.Cys63Tyr (B1:41621481, c.188G>A) and a less common p.Ala18Val (B1:42620835, c.53C>T) variant are suggested as two alleles influencing the Ticked phenotype. Bioinformatic and molecular modeling analysis suggests that these changes disrupt a key disulfide bond in the Dkk4 cysteine-rich domain 1 or Dkk4 signal peptide cleavage respectively. All coding variants were excluded as Ticked alleles using Strategy B.