Lan Zhu

Simultaneous inference of selection and population growth from patterns of variation in the human genome

Natural selection and demographic forces can have similar effects on patterns of DNA polymorphism. Therefore, to infer selection from samples of DNA sequences, one must simultaneously account for demographic effects. Here we take a model-based approach to this problem by developing predictions for patterns of polymorphism in the presence of both population size change and natural selection. If data are available from different functional classes of variation, and a priori information suggests that mutations in one of those classes are selectively neutral, then the putatively neutral class can be used to infer demographic parameters, and inferences regarding selection on other classes can be performed given demographic parameter estimates. This procedure is more robust to assumptions regarding the true underlying demography than previous approaches to detecting and analyzing selection. We apply this method to a large polymorphism data set from 301 human genes and find (i) widespread negative selection acting on standing nonsynonymous variation, (ii) that the fitness effects of nonsynonymous mutations are well predicted by several measures of amino acid exchangeability, especially site-specific methods, and (iii) strong evidence for very recent population growth.

Linkage Disequilibrium and Demographic History of Wild and Domestic Canids

Assessing the extent of linkage disequilibrium (LD) in natural populations of a nonmodel species has been difficult due to the lack of available genomic markers. However, with advances in genotyping and genome sequencing, genomic characterization of natural populations has become feasible. Using sequence data and SNP genotypes, we measured LD and modeled the demographic history of wild canid populations and domestic dog breeds. In 11 gray wolf populations and one coyote population, we find that the extent of LD as measured by the distance at which r2 = 0.2 extends <10 kb in outbred populations to >1.7 Mb in populations that have experienced significant founder events and bottlenecks. This large range in the extent of LD parallels that observed in 18 dog breeds where the r2 value varies from ∼20 kb to >5 Mb. Furthermore, in modeling demographic history under a composite-likelihood framework, we find that two of five wild canid populations exhibit evidence of a historical population contraction. Five domestic dog breeds display evidence for a minor population contraction during domestication and a more severe contraction during breed formation. Only a 5% reduction in nucleotide diversity was observed as a result of domestication, whereas the loss of nucleotide diversity with breed formation averaged 35%.

Estimation of heritabilities, genetic correlations, and breeding values of four traits that collectively define hip dysplasia in dogs

OBJECTIVE-To estimate heritabilities and genetic correlations among 4 traits of hip joints (distraction index [DI], dorsolateral subluxation [DLS] score, Norberg angle [NA], and extended-hip joint radiograph [EHR] score) and to derive the breeding values for these traits in dogs. ANIMALS-2,716 dogs of 17 breeds (1,551 dogs in which at least 1 hip joint trait was measured). PROCEDURES-The NA was measured, and an EHR score was assigned. Hip joint radiographs were obtained from some dogs to allow calculation of the DI and DLS score. Heritabilities, genetic correlations, and breeding values among the DI, DLS score, NA, and EHR score were calculated by use of a set of multiple-trait, derivative-free, restricted maximum likelihood computer programs. RESULTS-Among 2,716 dogs, 1,411 (52%) had an estimated inbreeding coefficient of 0%; the remaining dogs had a mean inbreeding coefficient of 6.21%. Estimated heritabilities were 0.61, 0.54, 0.73, and 0.76 for the DI, DLS score, NA, and EHR score, respectively. The EHR score was highly genetically correlated with the NA (r = -0.89) and was moderately genetically correlated with the DI (r = 0.69) and DLS score (r = -0.70). The NA was moderately genetically correlated with the DI (r = -0.69) and DLS score (r = 0.58). Genetic correlation between the DI and DLS score was high (r = -0.91). CONCLUSIONS AND CLINICAL RELEVANCE-Establishment of a selection index that makes use of breeding values jointly estimated from the DI, DLS score, NA, and EHR score should enhance breeding programs to reduce the incidence of hip dysplasia in dogs.

Retrospective analysis for genetic improvement of hip joints of cohort labrador retrievers in the United States: 1970-2007.

Background Canine Hip Dysplasia (CHD) is a common inherited disease that affects dog wellbeing and causes a heavy financial and emotional burden to dog owners and breeders due to secondary hip osteoarthritis. The Orthopedic Foundation for Animals (OFA) initiated a program in the 1960s to radiograph hip and elbow joints and release the OFA scores to the public for breeding dogs against CHD. Over last four decades, more than one million radiographic scores have been released. Methodology/Principal Findings The pedigrees in the OFA database consisted of 258,851 Labrador retrievers, the major breed scored by the OFA (25% of total records). Of these, 154,352 dogs had an OFA hip score reported between 1970 and 2007. The rest of the dogs (104,499) were the ancestors of the 154,352 dogs to link the pedigree relationships. The OFA hip score is based on a 7-point scale with the best ranked as 1 (excellent) and the worst hip dysplasia as 7. A mixed linear model was used to estimate the effects of age, sex, and test year period and to predict the breeding value for each dog. Additive genetic and residual variances were estimated using the average information restricted maximum likelihood procedure. The analysis also provided an inbreeding coefficient for each dog. The hip scores averaged 1.93 (±SD = 0.59) and the heritability was 0.21. A steady genetic improvement has accrued over the four decades. The breeding values decreased (improved) linearly. By the end of 2005, the total genetic improvement was 0.1 units, which is equivalent to 17% of the total phenotypic standard deviation. Conclusion/Significance A steady genetic improvement has been achieved through the selection based on the raw phenotype released by the OFA. As the heritability of the hip score was on the low end (0.21) of reported ranges, we propose that selection based on breeding values will result in more rapid genetic improvement than breeding based on phenotypic selection alone.

Differential Genetic Regulation of Canine Hip Dysplasia and Osteoarthritis

Background Canine hip dysplasia (HD) is a common polygenic trait characterized by hip malformation that results in osteoarthritis (OA). The condition in dogs is very similar to developmental dysplasia of the human hip which also leads to OA. Methodology/Principal Findings A total of 721 dogs, including both an association and linkage population, were genotyped. The association population included 8 pure breeds (Labrador retriever, Greyhounds, German Shepherd, Newfoundland, Golden retriever, Rottweiler, Border Collie and Bernese Mountain Dog). The linkage population included Labrador retrievers, Greyhounds, and their crosses. Of these, 366 dogs were genotyped at ∼22,000 single nucleotide polymorphism (SNP) loci and a targeted screen across 8 chromosomes with ∼3,300 SNPs was performed on 551 dogs (196 dogs were common to both sets). A mixed linear model approach was used to perform an association study on this combined association and linkage population. The study identified 4 susceptibility SNPs associated with HD and 2 SNPs associated with hip OA. Conclusion/Significance The identified SNPs included those near known genes (PTPRD, PARD3B, and COL15A1) reported to be associated with, or expressed in, OA in humans. This suggested that the canine model could provide a unique opportunity to identify genes underlying natural HD and hip OA, which are common and debilitating conditions in both dogs and humans.

Canine hip dysplasia is predictable by genotyping

OBJECTIVE To establish a predictive method using whole genome genotyping for early intervention in canine hip dysplasia (CHD) risk management, for the prevention of the progression of secondary osteoarthritis (OA), and for selective breeding. DESIGN Two sets of dogs (six breeds) were genotyped with dense SNPs covering the entire canine genome. The first set contained 359 dogs upon which a predictive formula for genomic breeding value (GBV) was derived by using their estimated breeding value (EBV) of the Norberg angle (a measure of CHD) and their genotypes. To investigate how well the formula would work for an individual dog with genotype only (without using EBV), a cross validation was performed by masking the EBV of one dog at a time. The genomic data and the EBV of the remaining dogs were used to predict the GBV for the single dog that was left out. The second set of dogs included 38 new Labrador retriever dogs, which had no pedigree relationship to the dogs in the first set. RESULTS The cross validation showed a strong correlation (R>0.7) between the EBV and the GBV. The independent validation showed a moderate correlation (R=0.5) between GBV for the Norberg angle and the observed Norberg angle (no EBV was available for the new 38 dogs). Sensitivity, specificity, positive and negative predictive values of the genomic data were all above 70%. CONCLUSIONS Prediction of CHD from genomic data is feasible, and can be applied for risk management of CHD and early selection for genetic improvement to reduce the prevalence of CHD in breeding programs. The prediction can be implemented before maturity, at which age current radiographic screening programs are traditionally applied, and as soon as DNA is available.

Evaluation of a fibrillin 2 gene haplotype associated with hip dysplasia and incipient osteoarthritis in dogs

OBJECTIVE To determine whether a mutation in the fibrillin 2 gene (FBN2) is associated with canine hip dysplasia (CHD) and osteoarthritis in dogs. ANIMALS 1,551 dogs. Procedures-Hip conformation was measured radiographically. The FBN2 was sequenced from genomic DNA of 21 Labrador Retrievers and 2 Greyhounds, and a haplotype in intron 30 of FBN2 was sequenced in 90 additional Labrador Retrievers and 143 dogs of 6 other breeds. Steady-state values of FBN2 mRNA and control genes were measured in hip joint tissues of fourteen 8-month-old Labrador Retriever-Greyhound crossbreeds. RESULTS The Labrador Retrievers homozygous for a 10-bp deletion haplotype in intron 30 of FBN2 had significantly worse CHD as measured via higher distraction index and extended-hip joint radiograph score and a lower Norberg angle and dorsolateral subluxation score. Among 143 dogs of 6 other breeds, those homozygous for the same deletion haplotype also had significantly worse radiographic CHD. Among the 14 crossbred dogs, as the dorsolateral subluxation score decreased, the capsular FBN2 mRNA increased significantly. Those dogs with incipient hip joint osteoarthritis had significantly increased capsular FBN2 mRNA, compared with those dogs without osteoarthritis. Dogs homozygous for the FBN2 deletion haplotype had significantly less FBN2 mRNA in their femoral head articular cartilage. CONCLUSIONS AND CLINICAL RELEVANCE The FBN2 deletion haplotype was associated with CHD. Capsular gene expression of FBN2 was confounded by incipient secondary osteoarthritis in dysplastic hip joints. Genes influencing complex traits in dogs can be identified by genome-wide screening, fine mapping, and candidate gene screening.

Single nucleotide polymorphisms refine QTL intervals for hip joint laxity in dogs

Hip laxity is one characteristic of canine hip dysplasia (CHD), an inheritable disease that leads to hip osteoarthritis. Using a genome-wide screen with 250 microsatellites in a crossbreed pedigree of 159 dysplastic Labrador retrievers and unaffected greyhounds, we previously identified putative (P < 0.01) QTL on canine chromosomes 11 and 29 (CFA11 and CFA29). To refine these QTL locations, we have genotyped 257 dogs including 105 Labrador retrievers, seven greyhounds, four generations of their crossbreed offspring and three German shepherds for 111 and 171 SNPs on CFA11 and CFA29 respectively. The distraction index (DI, a measure of maximum hip laxity) was used as an intermediate phenotype that predicts whether a hip joint will or will not develop osteoarthritis. Using a multipoint linkage analysis, significant evidence (95% posterior probability) was found for QTL contributing to hip laxity in the 16.2-21 cM region on CFA11 that explained 15-18% of the total variance in DI. Evidence for an independent QTL on CFA29 was weaker than that on CFA11. Identification of the causative mutation(s) will lead to better understanding of biochemical pathways in both dogs and humans with hip laxity and dysplasia.

Disomic Inheritance and Segregation Distortion of SSR Markers in Two Populations of Cynodon dactylon (L.) Pers. var. dactylon.

Common bermudagrass [C. dactylon (L.) Pers. var. dactylon] is economically and environmentally the most important member among Cynodon species because of its extensive use for turf, forage and soil erosion control in the world. However, information regarding the inheritance within the taxon is limited. Accordingly, the objective of this study was to determine qualitative inheritance mode in common bermudagrass. Two tetraploid (2n = 4x = 36), first-generation selfed (S1) populations, 228 progenies of ‘Zebra’ and 273 from A12359, were analyzed for segregation with 21 and 12 simple sequence repeat (SSR) markers, respectively. It is concluded that the inheritance mode of tetraploid bermudagrass was complete or near complete disomic. It is evident that the two bermudagrass parents had an allotetraploid genome with two distinct subgenomes since 33 SSR primer pairs amplified 34 loci, each having two alleles. Severe transmission ratio distortions occurred in the Zebra population while less so in the A12359 population. The findings of disomic inheritance and segregation ratio distortion in common bermudagrass is significant in subsequent linkage map construction, quantitative trait locus mapping and marker-assisted selection in the species.

Identification of quantitative trait loci for canine hip dysplasia by two sequential multipoint linkage analyses

Canine hip dysplasia (CHD) is characterized by hip laxity and subluxation that can lead to hip osteoarthritis. Studies have shown the involvement of multiple genetic regions in the expression of CHD. Although we have associated some variants in the region of fibrillin 2 with CHD in a subset of dogs, no major disease-associated gene has been identified. The focus of this study is to identify quantitative trait loci (QTL) associated with CHD. Two sequential multipoint linkage analyses based on a reversible jump Markov chain Monte Carlo approach were applied on a cross-breed pedigree of 366 dogs. Hip radiographic trait (Norberg Angle, NA) on both hips of each dog was tested for linkage to 21,455 single nucleotide polymorphisms across 39 chromosomes. Putative QTL for the NA was found on 11 chromosomes (1, 2, 3, 4, 7, 14, 19, 21, 32, 36, and 39). Identification of genes in the QTL region(s) can assist in identification of the aberrant genes and biochemical pathways involving hip dysplasia in both dogs and humans.