ArCH: improving the performance of clonal hematopoiesis variant calling and interpretation.

MOTIVATION: The acquisition of somatic mutations in hematopoietic stem and progenitor stem cells with resultant clonal expansion, termed clonal hematopoiesis (CH), is associated with increased risk of hematologic malignancies and other adverse outcomes. CH is generally present at low allelic fractions, but clonal expansion and acquisition of additional mutations leads to hematologic cancers in a small proportion of individuals. With high depth and high sensitivity sequencing, CH can be detected in most adults and its clonal trajectory mapped over time. However, accurate CH variant calling is challenging due to the difficulty in distinguishing low frequency CH mutations from sequencing artifacts. The lack of well-validated bioinformatic pipelines for CH calling may contribute to lack of reproducibility in studies of CH. RESULTS: Here, we developed ArCH, an Artifact filtering Clonal Hematopoiesis variant calling pipeline for detecting single nucleotide variants and short insertions/deletions by combining the output of four variant calling tools and filtering based on variant characteristics and sequencing error rate estimation. ArCH is an end-to-end cloud-based pipeline optimized to accept a variety of inputs with customizable parameters adaptable to multiple sequencing technologies, research questions, and datasets. Using deep targeted sequencing data generated from six acute myeloid leukemia patient tumor: normal dilutions, 31 blood samples with orthogonal validation, and 26 blood samples with technical replicates, we show that ArCH improves the sensitivity and positive predictive value of CH variant detection at low allele frequencies compared to standard application of commonly used variant calling approaches. AVAILABILITY AND IMPLEMENTATION: The code for this workflow is available at: https://github.com/kbolton-lab/ArCH.

Biallelic BRCA2 Mutations Shape the Somatic Mutational Landscape of Aggressive Prostate Tumors.

To identify clinically important molecular subtypes of prostate cancer (PCa), we characterized the somatic landscape of aggressive tumors via deep, whole-genome sequencing. In our discovery set of ten tumor/normal subject pairs with Gleason scores of 8-10 at diagnosis, coordinated analysis of germline and somatic variants, including single-nucleotide variants, indels, and structural variants, revealed biallelic BRCA2 disruptions in a subset of samples. Compared to the other samples, the PCa BRCA2-deficient tumors exhibited a complex and highly specific mutation signature, featuring a 2.88-fold increased somatic mutation rate, depletion of context-specific C>T substitutions, and an enrichment for deletions, especially those longer than 10 bp. We next performed a BRCA2 deficiency-targeted reanalysis of 150 metastatic PCa tumors, and each of the 18 BRCA2-mutated samples recapitulated the BRCA2 deficiency-associated mutation signature, underscoring the potent influence of these lesions on somatic mutagenesis and tumor evolution. Among all 21 individuals with BRCA2-deficient tumors, only about half carried deleterious germline alleles. Importantly, the somatic mutation signature in tumors with one germline and one somatic risk allele was indistinguishable from those with purely somatic mutations. Our observations clearly demonstrate that BRCA2-disrupted tumors represent a unique and clinically relevant molecular subtype of aggressive PCa, highlighting both the promise and utility of this mutation signature as a prognostic and treatment-selection biomarker. Further, any test designed to leverage BRCA2 status as a biomarker for PCa must consider both germline and somatic mutations and all types of deleterious mutations.

Comparison against 186 canid whole-genome sequences reveals survival strategies of an ancient clonally transmissible canine tumor.

Canine transmissible venereal tumor (CTVT) is a parasitic cancer clone that has propagated for thousands of years via sexual transfer of malignant cells. Little is understood about the mechanisms that converted an ancient tumor into the world's oldest known continuously propagating somatic cell lineage. We created the largest existing catalog of canine genome-wide variation and compared it against two CTVT genome sequences, thereby separating alleles derived from the founder's genome from somatic mutations that must drive clonal transmissibility. We show that CTVT has undergone continuous adaptation to its transmissible allograft niche, with overlapping mutations at every step of immunosurveillance, particularly self-antigen presentation and apoptosis. We also identified chronologically early somatic mutations in oncogenesis- and immune-related genes that may represent key initiators of clonal transmissibility. Thus, we provide the first insights into the specific genomic aberrations that underlie CTVT's dogged perseverance in canids around the world.

A copy number variant at the KITLG locus likely confers risk for canine squamous cell carcinoma of the digit.

The domestic dog is a robust model for studying the genetics of complex disease susceptibility. The strategies used to develop and propagate modern breeds have resulted in an elevated risk for specific diseases in particular breeds. One example is that of Standard Poodles (STPOs), who have increased risk for squamous cell carcinoma of the digit (SCCD), a locally aggressive cancer that causes lytic bone lesions, sometimes with multiple toe recurrence. However, only STPOs of dark coat color are at high risk; light colored STPOs are almost entirely unaffected, suggesting that interactions between multiple pathways are necessary for oncogenesis. We performed a genome-wide association study (GWAS) on STPOs, comparing 31 SCCD cases to 34 unrelated black STPO controls. The peak SNP on canine chromosome 15 was statistically significant at the genome-wide level (P(raw) = 1.60 × 10(-7); P(genome) = 0.0066). Additional mapping resolved the region to the KIT Ligand (KITLG) locus. Comparison of STPO cases to other at-risk breeds narrowed the locus to a 144.9-Kb region. Haplotype mapping among 84 STPO cases identified a minimal region of 28.3 Kb. A copy number variant (CNV) containing predicted enhancer elements was found to be strongly associated with SCCD in STPOs (P = 1.72 × 10(-8)). Light colored STPOs carry the CNV risk alleles at the same frequency as black STPOs, but are not susceptible to SCCD. A GWAS comparing 24 black and 24 light colored STPOs highlighted only the MC1R locus as significantly different between the two datasets, suggesting that a compensatory mutation within the MC1R locus likely protects light colored STPOs from disease. Our findings highlight a role for KITLG in SCCD susceptibility, as well as demonstrate that interactions between the KITLG and MC1R loci are potentially required for SCCD oncogenesis. These findings highlight how studies of breed-limited diseases are useful for disentangling multigene disorders.

Confirmation of genetic variants associated with lethal prostate cancer in a cohort of men from hereditary prostate cancer families.

Germline genetic variants have been suggested as prognostic biomarkers for identifying patients at high risk for lethal prostate cancer (PCa). Validation studies have confirmed the association of several single nucleotide polymorphisms (SNPs) with fatal PCa, but whether these variants affect PCa-specific mortality (PCSM) in patients with an inherited predisposition to PCa, based on familial history, is unknown. For this study, a cohort of 957 PCa patients from 270 hereditary prostate cancer families of European ancestry was genotyped for a panel of 22 PCSM-associated SNPs. Death certificates were reviewed to confirm cause of death. Mixed-effect Cox proportional hazards models were used to assess survival according to genotypes, accounting for relatedness and clinicopathological factors. Within this cohort, 98 PCa deaths were confirmed over an average follow-up period of 12.7 years after diagnosis. Variant allele carriers for three SNPs had significantly altered risk for PCSM [rs635261 at RNASEL, hazard ratio (HR), 0.35, 95% CI, 0.18-0.66; p = 0.002; rs915927 in XRCC1, HR, 1.91, 95% CI, 1.21-3.02; p = 0.009; and rs2494750 at AKT1, HR, 0.45, 95% CI, 0.23-0.90; p = 0.016). These results confirm the association of genetic variation in three genes with PCa lethality in a cohort of men with an inherited susceptibility to the disease and provide validation evidence that germline SNPs provide prognostic information for PCa patients. Development of a panel of germline biomarkers with clinical utility for distinguishing patients at detection who have an increased risk for fatal PCa is warranted.

Polygenic risk scores, radiation treatment exposures and subsequent cancer risk in childhood cancer survivors.

Survivors of childhood cancer are at increased risk for subsequent cancers attributable to the late effects of radiotherapy and other treatment exposures; thus, further understanding of the impact of genetic predisposition on risk is needed. Combining genotype data for 11,220 5-year survivors from the Childhood Cancer Survivor Study and the St Jude Lifetime Cohort, we found that cancer-specific polygenic risk scores (PRSs) derived from general population, genome-wide association study, cancer loci identified survivors of European ancestry at increased risk of subsequent basal cell carcinoma (odds ratio per s.d. of the PRS: OR = 1.37, 95% confidence interval (CI) = 1.29-1.46), female breast cancer (OR = 1.42, 95% CI = 1.27-1.58), thyroid cancer (OR = 1.48, 95% CI = 1.31-1.67), squamous cell carcinoma (OR = 1.20, 95% CI = 1.00-1.44) and melanoma (OR = 1.60, 95% CI = 1.31-1.96); however, the association for colorectal cancer was not significant (OR = 1.19, 95% CI = 0.94-1.52). An investigation of joint associations between PRSs and radiotherapy found more than additive increased risks of basal cell carcinoma, and breast and thyroid cancers. For survivors with radiotherapy exposure, the cumulative incidence of subsequent cancer by age 50 years was increased for those with high versus low PRS. These findings suggest a degree of shared genetic etiology for these malignancy types in the general population and survivors, which remains evident in the context of strong radiotherapy-related risk.

Genomic characterization of cervical lymph node metastases in papillary thyroid carcinoma following the Chornobyl accident.

Childhood radioactive iodine exposure from the Chornobyl accident increased papillary thyroid carcinoma (PTC) risk. While cervical lymph node metastases (cLNM) are well-recognized in pediatric PTC, the PTC metastatic process and potential radiation association are poorly understood. Here, we analyze cLNM occurrence among 428 PTC with genomic landscape analyses and known drivers (131I-exposed = 349, unexposed = 79; mean age = 27.9 years). We show that cLNM are more frequent in PTC with fusion (55%) versus mutation (30%) drivers, although the proportion varies by specific driver gene (RET-fusion = 71%, BRAF-mutation = 38%, RAS-mutation = 5%). cLNM frequency is not associated with other characteristics, including radiation dose. cLNM molecular profiling (N = 47) demonstrates 100% driver concordance with matched primary PTCs and highly concordant mutational spectra. Transcriptome analysis reveals 17 differentially expressed genes, particularly in the HOXC cluster and BRINP3; the strongest differentially expressed microRNA also is near HOXC10. Our findings underscore the critical role of driver alterations and provide promising candidates for elucidating the biological underpinnings of PTC cLNM.

HOXB13 mutations in a population-based, case-control study of prostate cancer.

BACKGROUND: Prostate cancer (PC) is the most frequently diagnosed non-skin malignancy in men in the Western world, yet few disease-associated mutations have been found. Recently, a low frequency recurring mutation in the HOXB13 gene was reported among both hereditary PC families and men from the general population. MATERIALS AND METHODS: We determined the distribution and frequency of the G84E HOXB13 variant in 1,310 incipient PC cases and 1,259 age-mated controls from a population-based, case-control study of PC. RESULTS: The G84E mutation was more frequent in cases than controls (1.3% vs. 0.4%, respectively), and men with the HOXB13 G84E variant had a 3.3-fold higher relative risk of PC compared with noncarriers (95% CI, 1.21-8.96). There was a stronger association between the G84E variant and PC among men with no first-degree relative with PC (OR, 4.04; 95% CI, 1.12-14.51) compared to men with a family history of PC (OR, 1.49; 95% CI, 0.30-7.50; P = 0.36 for interaction). We observed some evidence of higher risk estimates associated with the variant for men with higher versus lower Gleason score (OR, 4.13; 95% CI, 1.38-12.38 vs. OR, 2.71; 95% CI, 0.88-8.30), and advanced versus local stage (OR, 4.47; 95% CI, 1.28-15.57 vs. OR, 2.98; 95% CI, 1.04-8.49), however these differences were not statistically different. CONCLUSIONS: These results confirm the association of a rare HOXB13 mutation with PC in the general population and suggest that this variant may be associated with features of more aggressive disease.

Inflated expectations: Rare-variant association analysis using public controls.

The use of publicly available sequencing datasets as controls (hereafter, "public controls") in studies of rare variant disease associations has great promise but can increase the risk of false-positive discovery. The specific factors that could contribute to inflated distribution of test statistics have not been systematically examined. Here, we leveraged both public controls, gnomAD v2.1 and several datasets sequenced in our laboratory to systematically investigate factors that could contribute to the false-positive discovery, as measured by λΔ95, a measure to quantify the degree of inflation in statistical significance. Analyses of datasets in this investigation found that 1) the significantly inflated distribution of test statistics decreased substantially when the same variant caller and filtering pipelines were employed, 2) differences in library prep kits and sequencers did not affect the false-positive discovery rate and, 3) joint vs. separate variant-calling of cases and controls did not contribute to the inflation of test statistics. Currently available methods do not adequately adjust for the high false-positive discovery. These results, especially if replicated, emphasize the risks of using public controls for rare-variant association tests in which individual-level data and the computational pipeline are not readily accessible, which prevents the use of the same variant-calling and filtering pipelines on both cases and controls. A plausible solution exists with the emergence of cloud-based computing, which can make it possible to bring containerized analytical pipelines to the data (rather than the data to the pipeline) and could avert or minimize these issues. It is suggested that future reports account for this issue and provide this as a limitation in reporting new findings based on studies that cannot practically analyze all data on a single pipeline.

Family-based association analysis of 42 hereditary prostate cancer families identifies the Apolipoprotein L3 region on chromosome 22q12 as a risk locus.

Multiple genome-wide scans for hereditary prostate cancer (HPC) have identified susceptibility loci on nearly every chromosome. However, few results have been replicated with statistical significance. One exception is chromosome 22q, for which five independent linkage studies yielded strong evidence for a susceptibility locus in HPC families. Previously, we refined this region to a 2.53 Mb interval, using recombination mapping in 42 linked pedigrees. We now refine this locus to a 15 kb interval, spanning Apolipoprotein L3 (APOL3), using family-based association analyses of 150 total prostate cancer (PC) cases from two independent family collections with 506 unrelated population controls. Analysis of the two independent sets of PC cases highlighted single nucleotide polymorphisms (SNPs) within the APOL3 locus showing the strongest associations with HPC risk, with the most robust results observed when all 150 cases were combined. Analysis of 15 tagSNPs across the 5' end of the locus identified six SNPs with P-values < or =2 × 10(-4). The two independent sets of HPC cases highlight the same 15 kb interval at the 5' end of the APOL3 gene and provide strong evidence that SNPs within this 15 kb interval, or in strong linkage disequilibrium with it, contribute to HPC risk. Further analyses of this locus in an independent population-based, case-control study revealed an association between an SNP within the APOL3 locus and PC risk, which was not confirmed in the Cancer Genetic Markers of Susceptibility data set. This study further characterizes the 22q locus in HPC risk and suggests that the role of this region in sporadic PC warrants additional studies.

Analysis of Xq27-28 linkage in the international consortium for prostate cancer genetics (ICPCG) families.

BACKGROUND: Genetic variants are likely to contribute to a portion of prostate cancer risk. Full elucidation of the genetic etiology of prostate cancer is difficult because of incomplete penetrance and genetic and phenotypic heterogeneity. Current evidence suggests that genetic linkage to prostate cancer has been found on several chromosomes including the X; however, identification of causative genes has been elusive. METHODS: Parametric and non-parametric linkage analyses were performed using 26 microsatellite markers in each of 11 groups of multiple-case prostate cancer families from the International Consortium for Prostate Cancer Genetics (ICPCG). Meta-analyses of the resultant family-specific linkage statistics across the entire 1,323 families and in several predefined subsets were then performed. RESULTS: Meta-analyses of linkage statistics resulted in a maximum parametric heterogeneity lod score (HLOD) of 1.28, and an allele-sharing lod score (LOD) of 2.0 in favor of linkage to Xq27-q28 at 138 cM. In subset analyses, families with average age at onset less than 65 years exhibited a maximum HLOD of 1.8 (at 138 cM) versus a maximum regional HLOD of only 0.32 in families with average age at onset of 65 years or older. Surprisingly, the subset of families with only 2-3 affected men and some evidence of male-to-male transmission of prostate cancer gave the strongest evidence of linkage to the region (HLOD = 3.24, 134 cM). For this subset, the HLOD was slightly increased (HLOD = 3.47 at 134 cM) when families used in the original published report of linkage to Xq27-28 were excluded. CONCLUSIONS: Although there was not strong support for linkage to the Xq27-28 region in the complete set of families, the subset of families with earlier age at onset exhibited more evidence of linkage than families with later onset of disease. A subset of families with 2-3 affected individuals and with some evidence of male to male disease transmission showed stronger linkage signals. Our results suggest that the genetic basis for prostate cancer in our families is much more complex than a single susceptibility locus on the X chromosome, and that future explorations of the Xq27-28 region should focus on the subset of families identified here with the strongest evidence of linkage to this region.

In-utero exposure to zidovudine-containing antiretroviral therapy and clonal hematopoiesis in HIV-exposed uninfected newborns.

OBJECTIVE: Zidovudine (ZDV) has been extensively used in pregnant women to prevent vertical transmission of HIV but few studies have evaluated potential mutagenic effects of ZDV during fetal development. DESIGN: Our study investigated clonal hematopoiesis in HIV-exposed uninfected (HEU) newborns, 94 of whom were ZDV-exposed and 91 antiretroviral therapy (ART)-unexposed and matched for potential confounding factors. METHODS: Utilizing high depth sequencing and genotyping arrays, we comprehensively examined blood samples collected during the first week after birth for potential clonal hematopoiesis associated with fetal ZDV exposure, including clonal single nucleotide variants (SNVs), small insertions and deletions (indels), and large structural copy number or copy neutral alterations. RESULTS: We observed no statistically significant difference in the number of SNVs and indels per person in ZDV-exposed children (adjusted ratio [95% confidence interval, CI] for expected number of mutations = 0.79 [0.50--1.22], P = 0.3), and no difference in the number of large structural alterations. Mutations in common clonal hematopoiesis driver genes were not found in the study population. Mutational signature analyses on SNVs detected no novel signatures unique to the ZDV-exposed children and the mutational profiles were similar between the two groups. CONCLUSION: Our results suggest that clonal hematopoiesis at levels detectable in our study is not strongly influenced by in-utero ZDV exposure; however, additional follow-up studies are needed to further evaluate the safety and potential long-term impacts of in-utero ZDV exposure in HEU children as well as better investigate genomic aberrations occurring late in pregnancy.

Frequency of Pathogenic Germline Variants in Cancer-Susceptibility Genes in the Childhood Cancer Survivor Study.

Background: Pediatric cancers are the leading cause of death by disease in children despite improved survival rates overall. The contribution of germline genetic susceptibility to pediatric cancer survivors has not been extensively characterized. We assessed the frequency of pathogenic or likely pathogenic (P/LP) variants in 5451 long-term pediatric cancer survivors from the Childhood Cancer Survivor Study. Methods: Exome sequencing was conducted on germline DNA from 5451 pediatric cancer survivors (cases who survived ≥5 years from diagnosis; n = 5105 European) and 597 European cancer-free adults (controls). Analyses focused on comparing the frequency of rare P/LP variants in 237 cancer-susceptibility genes and a subset of 60 autosomal dominant high-to-moderate penetrance genes, for both case-case and case-control comparisons. Results: Of European cases, 4.1% harbored a P/LP variant in high-to-moderate penetrance autosomal dominant genes compared with 1.3% in controls (2-sided P = 3 × 10-4). The highest frequency of P/LP variants was in genes typically associated with adult onset rather than pediatric cancers, including BRCA1/2, FH, PALB2, PMS2, and CDKN2A. A statistically significant excess of P/LP variants, after correction for multiple tests, was detected in patients with central nervous system cancers (NF1, SUFU, TSC1, PTCH2), Wilms tumor (WT1, REST), non-Hodgkin lymphoma (PMS2), and soft tissue sarcomas (SDHB, DICER1, TP53, ERCC4, FGFR3) compared with other pediatric cancers. Conclusion: In long-term pediatric cancer survivors, we identified P/LP variants in cancer-susceptibility genes not previously associated with pediatric cancer as well as confirmed known associations. Further characterization of variants in these genes in pediatric cancer will be important to provide optimal genetic counseling for patients and their families.

Radiation-related genomic profile of papillary thyroid carcinoma after the Chernobyl accident.

The 1986 Chernobyl nuclear power plant accident increased papillary thyroid carcinoma (PTC) incidence in surrounding regions, particularly for radioactive iodine (131I)-exposed children. We analyzed genomic, transcriptomic, and epigenomic characteristics of 440 PTCs from Ukraine (from 359 individuals with estimated childhood 131I exposure and 81 unexposed children born after 1986). PTCs displayed radiation dose-dependent enrichment of fusion drivers, nearly all in the mitogen-activated protein kinase pathway, and increases in small deletions and simple/balanced structural variants that were clonal and bore hallmarks of nonhomologous end-joining repair. Radiation-related genomic alterations were more pronounced for individuals who were younger at exposure. Transcriptomic and epigenomic features were strongly associated with driver events but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth after environmental radiation exposure.

Subsequent Neoplasm Risk Associated With Rare Variants in DNA Damage Response and Clinical Radiation Sensitivity Syndrome Genes in the Childhood Cancer Survivor Study.

PURPOSE: Radiotherapy for childhood cancer is associated with elevated subsequent neoplasm (SN) risk, but the contribution of rare variants in DNA damage response and radiation sensitivity genes to SN risk is unknown. PATIENTS AND METHODS: We conducted whole-exome sequencing in a cohort of childhood cancer survivors originally diagnosed during 1970 to 1986 (mean follow-up, 32.7 years), with reconstruction of doses to body regions from radiotherapy records. We identified patients who developed SN types previously reported to be related to radiotherapy (RT-SNs; eg, basal cell carcinoma [BCC], breast cancer, meningioma, thyroid cancer, sarcoma) and matched controls (sex, childhood cancer type/diagnosis, age, SN location, radiation dose, survival). Conditional logistic regression assessed SN risk associated with potentially protein-damaging rare variants (SnpEff, ClinVar) in 476 DNA damage response or radiation sensitivity genes with exact permutation-based P values using a Bonferroni-corrected significance threshold of P < 8.06 × 10-5. RESULTS: Among 5,105 childhood cancer survivors of European descent, 1,108 (21.7%) developed at least 1 RT-SN. Out-of-field RT-SN risk, excluding BCC, was associated with homologous recombination repair (HRR) gene variants (patient cases, 23.2%; controls, 10.8%; odds ratio [OR], 2.6; 95% CI, 1.7 to 3.9; P = 4.79 × 10-5), most notably but nonsignificantly for FANCM (patient cases, 4.0%; matched controls, 0.6%; P = 9.64 × 10-5). HRR variants were not associated with likely in/near-field RT-SNs, excluding BCC (patient cases, 12.7%; matched controls, 12.9%; P = .92). Irrespective of radiation dose, risk for RT-SNs was also associated with EXO1 variants (patient cases, 1.8%; controls, 0.4%; P = 3.31 × 10-5), another gene implicated in DNA double-strand break repair. CONCLUSION: In this large-scale discovery study, we identified novel associations between RT-SN risk after childhood cancer and potentially protein-damaging rare variants in genes involved in DNA double-strand break repair, particularly HRR. With replication, these results could affect screening recommendations for childhood cancer survivors and risk-benefit assessments of treatment approaches.

Association of megalin genetic polymorphisms with prostate cancer risk and prognosis.

PURPOSE: Megalin, an endocytic receptor expressed by prostate epithelial cells, can internalize biologically active androgens bound to sex hormone binding globulin. Genetic variation within megalin could potentially influence levels of steroid hormone uptake. EXPERIMENTAL DESIGN: Forty haplotype-tagging single-nucleotide polymorphisms (htSNP) were analyzed in a population-based, case-control study of 553 Caucasian men who were diagnosed with prostate cancer between the ages of 40 and 64 years from the Seattle-Puget Sound region and 534 control men. Prostate cancer risk was estimated using adjusted unconditional logistic regression for both individual SNPs and haplotypes. Risks of disease recurrence/progression and prostate-specific cancer mortality were estimated using Cox proportional hazards regression. RESULTS: We found no strong evidence of altered risk of developing prostate cancer for any of the htSNPs when they were assessed individually or in haplotypes. However, three htSNPs were significantly associated with both disease recurrence/progression and mortality. Risk of recurrence/progression alone was also associated with five additional htSNPs, and six other htSNPS showed evidence of modification by primary androgen deprivation therapy. Two additional htSNPs were significantly associated with altered risk of death from prostate cancer. CONCLUSIONS: Preliminary results suggest that common genetic variation within the megalin gene could alter both risk of recurrence/progression and prostate-specific cancer mortality. In addition, androgen deprivation therapy effectiveness may be modified by the activity of this gene. To our knowledge, this is the first study that has examined polymorphisms within the megalin gene for associations with prostate cancer risk and outcomes.

Whole genome sequencing of canids reveals genomic regions under selection and variants influencing morphology.

Domestic dog breeds are characterized by an unrivaled diversity of morphologic traits and breed-associated behaviors resulting from human selective pressures. To identify the genetic underpinnings of such traits, we analyze 722 canine whole genome sequences (WGS), documenting over 91 million single nucleotide and small indels, creating a large catalog of genomic variation for a companion animal species. We undertake both selective sweep analyses and genome wide association studies (GWAS) inclusive of over 144 modern breeds, 54 wild canids and a hundred village dogs. Our results identify variants of strong impact associated with 16 phenotypes, including body weight variation which, when combined with existing data, explain greater than 90% of body size variation in dogs. We thus demonstrate that GWAS and selection scans performed with WGS are powerful complementary methods for expanding the utility of companion animal systems for the study of mammalian growth and biology.

Fine mapping of familial prostate cancer families narrows the interval for a susceptibility locus on chromosome 22q12.3 to 1.36 Mb.

Genetic studies suggest that hereditary prostate cancer is a genetically heterogeneous disease with multiple contributing loci. Studies of high-risk prostate cancer families selected for aggressive disease, analysis of large multigenerational families, and a meta-analysis from the International Consortium for Prostate Cancer Genetics (ICPCG), all highlight chromosome 22q12.3 as a susceptibility locus with strong statistical significance. Recently, two publications have narrowed the 22q12.3 locus to a 2.18 Mb interval using 54 high-risk families from the ICPCG collaboration, as defined by three recombination events on either side of the locus. In this paper, we present the results from fine mapping studies at 22q12.3 using both haplotype and recombination data from 42 high-risk families contributed from the Mayo Clinic and the Prostate Cancer Genetic Research Study (PROGRESS) mapping studies. No clear consensus interval is present when all families are used. However, in the subset of 14 families with >/=5 affected men per family, a 2.53-Mb shared consensus segment that overlaps with the previously published interval is identified. Combining these results with data from the earlier ICPCG study reduces the three-recombination interval at 22q12.3 to approximately 1.36 Mb.

Multiple independent genetic variants in the 8q24 region are associated with prostate cancer risk.

Recently, the 8q24 region has been identified as a prostate cancer susceptibility locus in a genome-wide scan of prostate cancer families in Iceland and an admixture scan of African Americans. Further investigations of variants at 8q24 have shown the existence of additional single nucleotide polymorphisms (SNPs) that enhance prostate cancer risk, suggesting the possibility of multiple regions harboring variants for the disease. In the present population-based study of Caucasians (1,308 cases and 1,266 controls) and African Americans (149 cases and 85 controls), we tested the association between prostate cancer and 23 SNPs in the 8q24 region. Fourteen SNPs in Caucasians and 5 SNPs in African Americans were significantly associated with risk of prostate cancer after adjusting for multiple comparisons; of these, 5 SNPs in Caucasians and 3 in African Americans were independently associated with risk. The strongest association was for rs6983561 (carriers of any C allele) with an odds ratio of 1.6 (95% confidence interval, 1.1-2.1) in Caucasians; variants in rs979200, rs1016343, rs7837328, and rs10090154 were also independently associated with risk. In African Americans, the strongest association was for rs7000448 (carriers of any T allele) with an odds ratio of 3.4 (95% confidence interval, 1.3-8.7). In addition, two SNPs that extend the boundaries of the 8q24 region were significantly associated with risk: rs979200 at the centromeric boundary and rs3891248, located in the first intron of the c-MYC gene (IVS1-355), which identifies a new telomeric boundary.

Multiple novel prostate cancer predisposition loci confirmed by an international study: the PRACTICAL Consortium.

A recent genome-wide association study found that genetic variants on chromosomes 3, 6, 7, 10, 11, 19 and X were associated with prostate cancer risk. We evaluated the most significant single-nucleotide polymorphisms (SNP) in these loci using a worldwide consortium of 13 groups (PRACTICAL). Blood DNA from 7,370 prostate cancer cases and 5,742 male controls was analyzed by genotyping assays. Odds ratios (OR) associated with each genotype were estimated using unconditional logistic regression. Six of the seven SNPs showed clear evidence of association with prostate cancer (P = 0.0007-P = 10(-17)). For each of these six SNPs, the estimated per-allele OR was similar to those previously reported and ranged from 1.12 to 1.29. One SNP on 3p12 (rs2660753) showed a weaker association than previously reported [per-allele OR, 1.08 (95% confidence interval, 1.00-1.16; P = 0.06) versus 1.18 (95% confidence interval, 1.06-1.31)]. The combined risks associated with each pair of SNPs were consistent with a multiplicative risk model. Under this model, and in combination with previously reported SNPs on 8q and 17q, these loci explain 16% of the familial risk of the disease, and men in the top 10% of the risk distribution have a 2.1-fold increased risk relative to general population rates. This study provides strong confirmation of these susceptibility loci in multiple populations and shows that they make an important contribution to prostate cancer risk prediction.