Identification of common predisposing loci to hematopoietic cancers in four dog breeds.

Histiocytic sarcoma (HS) is a rare but aggressive cancer in both humans and dogs. The spontaneous canine model, which has clinical, epidemiological, and histological similarities with human HS and specific breed predispositions, provides a unique opportunity to unravel the genetic basis of this cancer. In this study, we aimed to identify germline risk factors associated with the development of HS in canine-predisposed breeds. We used a methodology that combined several genome-wide association studies in a multi-breed and multi-cancer approach as well as targeted next-generation sequencing, and imputation We combined several dog breeds (Bernese mountain dogs, Rottweilers, flat-coated retrievers, and golden retrievers), and three hematopoietic cancers (HS, lymphoma, and mast cell tumor). Results showed that we not only refined the previously identified HS risk CDKN2A locus, but also identified new loci on canine chromosomes 2, 5, 14, and 20. Capture and targeted sequencing of specific loci suggested the existence of regulatory variants in non-coding regions and methylation mechanisms linked to risk haplotypes, which lead to strong cancer predisposition in specific dog breeds. We also showed that these canine cancer predisposing loci appeared to be due to the additive effect of several risk haplotypes involved in other hematopoietic cancers such as lymphoma or mast cell tumors as well. This illustrates the pleiotropic nature of these canine cancer loci as observed in human oncology, thereby reinforcing the interest of predisposed dog breeds to study cancer initiation and progression.

Multitrait GWAS to connect disease variants and biological mechanisms.

Genome-wide association studies (GWASs) have uncovered a wealth of associations between common variants and human phenotypes. Here, we present an integrative analysis of GWAS summary statistics from 36 phenotypes to decipher multitrait genetic architecture and its link with biological mechanisms. Our framework incorporates multitrait association mapping along with an investigation of the breakdown of genetic associations into clusters of variants harboring similar multitrait association profiles. Focusing on two subsets of immunity and metabolism phenotypes, we then demonstrate how genetic variants within clusters can be mapped to biological pathways and disease mechanisms. Finally, for the metabolism set, we investigate the link between gene cluster assignment and the success of drug targets in randomized controlled trials.

IRF8 regulates efficacy of therapeutic anti-CD20 monoclonal antibodies.

Anti-CD20 monoclonal antibodies such as Rituximab, Ofatumumab, and Obinutuzumab are widely used to treat lymphomas and autoimmune diseases. They act by depleting B cells, mainly through Fc-dependent effectors functions. Some patients develop resistance to treatment but the underlying mechanisms are poorly understood. Here, we performed a genome-wide CRISPR/Cas9 screen to identify genes regulating the efficacy of anti-CD20 antibodies. We used as a model the killing of RAJI B cells by Rituximab through complement-dependent-cytotoxicity (CDC). As expected, the screen identified MS4A1, encoding CD20, the target of Rituximab. Among other identified genes, the role of Interferon Regulatory Factor 8 (IRF8) was validated in two B-cell lines. IRF8 KO also decreased the efficacy of antibody-dependent cellular cytotoxicity and phagocytosis (ADCC and ADCP) induced by anti-CD20 antibodies. We further show that IRF8 is necessary for efficient CD20 transcription. Levels of IRF8 and CD20 RNA or proteins correlated in normal B cells and in hundreds of malignant B cells. Therefore, IRF8 regulates CD20 expression and controls the depleting capacity of anti-CD20 antibodies. Our results bring novel insights into the pathways underlying resistance to CD20-targeting immunotherapies.

A novel role for ciliary function in atopy: ADGRV1 and DNAH5 interactions.

BACKGROUND: Atopy, an endotype underlying allergic diseases, has a substantial genetic component. OBJECTIVE: Our goal was to identify novel genes associated with atopy in asthma-ascertained families. METHODS: We implemented a 3-step analysis strategy in 3 data sets: the Epidemiological Study on the Genetics and Environment of Asthma (EGEA) data set (1660 subjects), the Saguenay-Lac-Saint-Jean study data set (1138 subjects), and the Medical Research Council (MRC) data set (446 subjects). This strategy included a single nucleotide polymorphism (SNP) genome-wide association study (GWAS), the selection of related gene pairs based on statistical filtering of GWAS results, and text-mining filtering using Gene Relationships Across Implicated Loci and SNP-SNP interaction analysis of selected gene pairs. RESULTS: We identified the 5q14 locus, harboring the adhesion G protein-coupled receptor V1 (ADGRV1) gene, which showed genome-wide significant association with atopy (rs4916831, meta-analysis P value = 6.8 × 10-9). Statistical filtering of GWAS results followed by text-mining filtering revealed relationships between ADGRV1 and 3 genes showing suggestive association with atopy (P ≤ 10-4). SNP-SNP interaction analysis between ADGRV1 and these 3 genes showed significant interaction between ADGRV1 rs17554723 and 2 correlated SNPs (rs2134256 and rs1354187) within the dynein axonemal heavy chain 5 (DNAH5) gene (Pmeta-int = 3.6 × 10-5 and 6.1 × 10-5, which met the multiple-testing corrected threshold of 7.3 × 10-5). Further conditional analysis indicated that rs2134256 alone accounted for the interaction signal with rs17554723. CONCLUSION: Because both DNAH5 and ADGRV1 contribute to ciliary function, this study suggests that ciliary dysfunction might represent a novel mechanism underlying atopy. Combining GWAS and epistasis analysis driven by statistical and knowledge-based evidence represents a promising approach for identifying new genes involved in complex traits.

A comprehensive genome-wide analysis of melanoma Breslow thickness identifies interaction between CDC42 and SCIN genetic variants.

Breslow thickness (BT) is a major prognostic factor of cutaneous melanoma (CM), the most fatal skin cancer. The genetic component of BT has only been explored by candidate gene studies with inconsistent results. Our objective was to uncover the genetic factors underlying BT using an hypothesis-free genome-wide approach. Our analysis strategy integrated a genome-wide association study (GWAS) of single nucleotide polymorphisms (SNPs) for BT followed by pathway analysis of GWAS outcomes using the gene-set enrichment analysis (GSEA) method and epistasis analysis within BT-associated pathways. This strategy was applied to two large CM datasets with Hapmap3-imputed SNP data: the French MELARISK study for discovery (966 cases) and the MD Anderson Cancer Center study (1,546 cases) for replication. While no marginal effect of individual SNPs was revealed through GWAS, three pathways, defined by gene ontology (GO) categories were significantly enriched in genes associated with BT (false discovery rate ≤5% in both studies): hormone activity, cytokine activity and myeloid cell differentiation. Epistasis analysis, within each significant GO, identified a statistically significant interaction between CDC42 and SCIN SNPs (pmeta-int =2.2 × 10(-6) , which met the overall multiple-testing corrected threshold of 2.5 × 10(-6) ). These two SNPs (and proxies) are strongly associated with CDC42 and SCIN gene expression levels and map to regulatory elements in skin cells. This interaction has important biological relevance since CDC42 and SCIN proteins have opposite effects in actin cytoskeleton organization and dynamics, a key mechanism underlying melanoma cell migration and invasion.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk.

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Identification of genomic regions associated with phenotypic variation between dog breeds using selection mapping.

The extraordinary phenotypic diversity of dog breeds has been sculpted by a unique population history accompanied by selection for novel and desirable traits. Here we perform a comprehensive analysis using multiple test statistics to identify regions under selection in 509 dogs from 46 diverse breeds using a newly developed high-density genotyping array consisting of >170,000 evenly spaced SNPs. We first identify 44 genomic regions exhibiting extreme differentiation across multiple breeds. Genetic variation in these regions correlates with variation in several phenotypic traits that vary between breeds, and we identify novel associations with both morphological and behavioral traits. We next scan the genome for signatures of selective sweeps in single breeds, characterized by long regions of reduced heterozygosity and fixation of extended haplotypes. These scans identify hundreds of regions, including 22 blocks of homozygosity longer than one megabase in certain breeds. Candidate selection loci are strongly enriched for developmental genes. We chose one highly differentiated region, associated with body size and ear morphology, and characterized it using high-throughput sequencing to provide a list of variants that may directly affect these traits. This study provides a catalogue of genomic regions showing extreme reduction in genetic variation or population differentiation in dogs, including many linked to phenotypic variation. The many blocks of reduced haplotype diversity observed across the genome in dog breeds are the result of both selection and genetic drift, but extended blocks of homozygosity on a megabase scale appear to be best explained by selection. Further elucidation of the variants under selection will help to uncover the genetic basis of complex traits and disease.

The colibactin-producing Escherichia coli alters the tumor microenvironment to immunosuppressive lipid overload facilitating colorectal cancer progression and chemoresistance.

Intratumoral bacteria flexibly contribute to cellular and molecular tumor heterogeneity for supporting cancer recurrence through poorly understood mechanisms. Using spatial metabolomic profiling technologies and 16SrRNA sequencing, we herein report that right-sided colorectal tumors are predominantly populated with Colibactin-producing Escherichia coli (CoPEC) that are locally establishing a high-glycerophospholipid microenvironment with lowered immunogenicity. It coincided with a reduced infiltration of CD8+ T lymphocytes that produce the cytotoxic cytokines IFN-γ where invading bacteria have been geolocated. Mechanistically, the accumulation of lipid droplets in infected cancer cells relied on the production of colibactin as a measure to limit genotoxic stress to some extent. Such heightened phosphatidylcholine remodeling by the enzyme of the Land's cycle supplied CoPEC-infected cancer cells with sufficient energy for sustaining cell survival in response to chemotherapies. This accords with the lowered overall survival of colorectal patients at stage III-IV who were colonized by CoPEC when compared to patients at stage I-II. Accordingly, the sensitivity of CoPEC-infected cancer cells to chemotherapies was restored upon treatment with an acyl-CoA synthetase inhibitor. By contrast, such metabolic dysregulation leading to chemoresistance was not observed in human colon cancer cells that were infected with the mutant strain that did not produce colibactin (11G5∆ClbQ). This work revealed that CoPEC locally supports an energy trade-off lipid overload within tumors for lowering tumor immunogenicity. This may pave the way for improving chemoresistance and subsequently outcome of CRC patients who are colonized by CoPEC.

Parvimonas micra, an oral pathobiont associated with colorectal cancer, epigenetically reprograms human colonocytes.

Recently, an intestinal dysbiotic microbiota with enrichment in oral cavity bacteria has been described in colorectal cancer (CRC) patients. Here, we characterize and investigate one of these oral pathobionts, the Gram-positive anaerobic coccus Parvimonas micra. We identified two phylotypes (A and B) exhibiting different phenotypes and adhesion capabilities. We observed a strong association of phylotype A with CRC, with its higher abundance in feces and in tumoral tissue compared with the normal homologous colonic mucosa, which was associated with a distinct methylation status of patients. By developing an in vitro hypoxic co-culture system of human primary colonic cells with anaerobic bacteria, we show that P. micra phylotype A alters the DNA methylation profile promoters of key tumor-suppressor genes, oncogenes, and genes involved in epithelial-mesenchymal transition. In colonic mucosa of CRC patients carrying P. micra phylotype A, we found similar DNA methylation alterations, together with significant enrichment of differentially expressed genes in pathways involved in inflammation, cell adhesion, and regulation of actin cytoskeleton, providing evidence of P. micra's possible role in the carcinogenic process.

Annotation of the domestic dog genome sequence: finding the missing genes.

There are over 350 genetically distinct breeds of domestic dog that present considerable variation in morphology, physiology, and disease susceptibility. The genome sequence of the domestic dog was assembled and released in 2005, providing an estimated 20,000 protein-coding genes that are a great asset to the scientific community that uses the dog system as a genetic biomedical model and for comparative and evolutionary studies. Although the canine gene set had been predicted using a combination of ab initio methods, homology studies, motif analysis, and similarity-based programs, it still requires a deep annotation of noncoding genes, alternative splicing, pseudogenes, regulatory regions, and gain and loss events. Such analyses could benefit from new sequencing technologies (RNA-Seq) to better exploit the advantages of the canine genetic system in tracking disease genes. Here, we review the catalog of canine protein-coding genes and the search for missing genes, and we propose rationales for an accurate identification of noncoding genes though next-generation sequencing.

Mixed-model admixture mapping identifies smoking-dependent loci of lung function in African Americans.

Admixture mapping has led to the discovery of many genes associated with differential disease risk by ancestry, highlighting the importance of ancestry-based approaches to association studies. However, the potential of admixture mapping in deciphering the interplay between genes and environment exposures has been seldom explored. Here we performed a genome-wide screening of local ancestry-smoking interactions for five spirometric lung function phenotypes in 3300 African Americans from the COPDGene study. To account for population structure and outcome heterogeneity across exposure groups, we developed a multi-component linear mixed model for mapping gene-environment interactions and empirically showed its robustness and increased power. When applied to the COPDGene study, our approach identified two 11p15.2-3 and 2q37 loci, exhibiting local ancestry-smoking interactions at genome-wide significant level, which would have been missed by standard single-nucleotide polymorphism analyses. These two loci harbor the PARVA and RAB17 genes previously recognized to be involved in smoking behavior. Overall, our study provides the first evidence for potential synergistic effects between African ancestry and smoking on pulmonary function, and underlines the importance of ethnic diversity in genetic studies.

RNA profiles of rat olfactory epithelia: individual and age related variations.

BACKGROUND: Mammalian genomes contain a large number (approximately 1000) of olfactory receptor (OR) genes, many of which (20 to 50%) are pseudogenes. OR gene transcription is not restricted to the olfactory epithelium, but is found in numerous tissues. Using microarray hybridization and RTqPCR, we analyzed the mRNA profiles of the olfactory epithelium of male and female Brown Norway rats of different origins and ages (newborn, adult and old). RESULTS: (1) We observed very little difference between males and females and between rats from two different suppliers. (2) Different OR genes were expressed at varying levels, rather than uniformly across the four endoturbinates. (3) A large proportion of the gene transcripts (2/3 of all probes) were detected in all three age groups. Adult and older rats expressed similar numbers of OR genes, both expressing more OR genes than newborns. (4) Comparisons of whole transcriptomes or transcription profiles of expressed OR genes only showed a clear clustering of the samples as a function of age. (5) Most OR genes were expressed at lower levels at birth than in older animals, but a small number of OR genes were expressed specifically or were overexpressed in newborns. CONCLUSION: Not all OR genes are expressed at a detectable level. Pups expressed fewer OR genes than adult rats, and generally at a lower level; however, a small subset of OR genes were more strongly expressed in these newborn rats. The reasons for these differences are not understood. However, the specific expression of some OR genes in newborn olfactory epithelia may be related to the blindness and deafness of pups at birth, when these pups are heavily reliant on olfaction and their mother.

Genetic diversity of canine olfactory receptors.

BACKGROUND: Evolution has resulted in large repertoires of olfactory receptor (OR) genes, forming the largest gene families in mammalian genomes. Knowledge of the genetic diversity of olfactory receptors is essential if we are to understand the differences in olfactory sensory capability between individuals. Canine breeds constitute an attractive model system for such investigations. RESULTS: We sequenced 109 OR genes considered representative of the whole OR canine repertoire, which consists of more than 800 genes, in a cohort of 48 dogs of six different breeds. SNP frequency showed the overall level of polymorphism to be high. However, the distribution of SNP was highly heterogeneous among OR genes. More than 50% of OR genes were found to harbour a large number of SNP, whereas the rest were devoid of SNP or only slightly polymorphic. Heterogeneity was also observed across breeds, with 25% of the SNP breed-specific. Linkage disequilibrium within OR genes and OR clusters suggested a gene conversion process, consistent with a mean level of polymorphism higher than that observed for introns and intergenic sequences. A large proportion (47%) of SNP induced amino-acid changes and the Ka/Ks ratio calculated for all alleles with a complete ORF indicated a low selective constraint with respect to the high level of redundancy of the olfactory combinatory code and an ongoing pseudogenisation process, which affects dog breeds differently. CONCLUSION: Our demonstration of a high overall level of polymorphism, likely to modify the ligand-binding capacity of receptors distributed differently within the six breeds tested, is the first step towards understanding why Labrador Retrievers and German Shepherd Dogs have a much greater potential for use as sniffer dogs than Pekingese dogs or Greyhounds. Furthermore, the heterogeneity in OR polymorphism observed raises questions as to why, in a context in which most OR genes are highly polymorphic, a subset of these genes is not? This phenomenon may be related to the nature of their ligands and their importance in everyday life.

Revisiting the missing protein-coding gene catalog of the domestic dog.

BACKGROUND: Among mammals for which there is a high sequence coverage, the whole genome assembly of the dog is unique in that it predicts a low number of protein-coding genes, ~19,000, compared to the over 20,000 reported for other mammalian species. Of particular interest are the more than 400 of genes annotated in primates and rodent genomes, but missing in dog. RESULTS: Using over 14,000 orthologous genes between human, chimpanzee, mouse rat and dog, we built multiple pairwise synteny maps to infer short orthologous intervals that were targeted for characterizing the canine missing genes. Based on gene prediction and a functionality test using the ratio of replacement to silent nucleotide substitution rates (d(N)/d(S)), we provide compelling structural and functional evidence for the identification of 232 new protein-coding genes in the canine genome and 69 gene losses, characterized as undetected gene or pseudogenes. Gene loss phyletic pattern analysis using ten species from chicken to human allowed us to characterize 28 canine-specific gene losses that have functional orthologs continuously from chicken or marsupials through human, and 10 genes that arose specifically in the evolutionary lineage leading to rodent and primates. CONCLUSION: This study demonstrates the central role of comparative genomics for refining gene catalogs and exploring the evolutionary history of gene repertoires, particularly as applied for the characterization of species-specific gene gains and losses.

Epidemiology, pathology, and genetics of histiocytic sarcoma in the Bernese mountain dog breed.

Histiocytic sarcoma (HS) refers to a highly aggressive and frequently disseminated neoplastic disease belonging to the class of canine histiocytic proliferative disorders. Disseminated HS (previously called malignant histiocytosis) is highly breed specific, with Bernese mountain dogs (BMDs), rottweilers, and retrievers having a high prevalence with a frequency of approximately 25% in the BMD breed. We collected DNA samples and clinical information from 800 BMDs, of which 200 are affected by HS. To better characterize the physiopathology and epidemiology, an in-depth analysis of 89 BMD cases has been performed. The mean age of onset was 6.5 years, males and females being equally affected. The clinical features, biochemical parameters, and pathological features have been determined. The life span after diagnosis has been estimated to be 49 days. A large BMD pedigree of 327 dogs, 121 of which are affected, was assembled. Using a subset of 160 BMDs, encompassing 21 complete sibships, we now propose an oligogenic transmission mode of the disease. Whole-genome linkage scans as well as association studies using a case/control analysis, in parallel with expression profiling of neoplastic versus normal histiocytes, are all underway. Altogether, these complementary approaches are expected to localize the genes for HS in the BMD, leading to advances in our knowledge of histiocyte diseases in dogs and humans.

A comprehensive study of metabolite genetics reveals strong pleiotropy and heterogeneity across time and context.

Genetic studies of metabolites have identified thousands of variants, many of which are associated with downstream metabolic and obesogenic disorders. However, these studies have relied on univariate analyses, reducing power and limiting context-specific understanding. Here we aim to provide an integrated perspective of the genetic basis of metabolites by leveraging the Finnish Metabolic Syndrome In Men (METSIM) cohort, a unique genetic resource which contains metabolic measurements, mostly lipids, across distinct time points as well as information on statin usage. We increase effective sample size by an average of two-fold by applying the Covariates for Multi-phenotype Studies (CMS) approach, identifying 588 significant SNP-metabolite associations, including 228 new associations. Our analysis pinpoints a small number of master metabolic regulator genes, balancing the relative proportion of dozens of metabolite levels. We further identify associations to changes in metabolic levels across time as well as genetic interactions with statin at both the master metabolic regulator and genome-wide level.

New susceptibility loci for cutaneous melanoma risk and progression revealed using a porcine model.

Despite major advances, it is estimated that a large part of melanoma predisposing genes remains to be discovered. Animal models of spontaneous diseases are valuable tools and experimental crosses can be used to identify and fine-map new susceptibility loci associated with melanoma. We performed a Genome-Wide Association Study (GWAS) of melanoma occurrence and progression (clinical ulceration and presence of metastasis) in a porcine model of spontaneous melanoma, the MeLiM pig. Five loci on chromosomes 2, 5, 7, 8 and 16 showed genome-wide significant associations (p < 5 × 10-6) with either one of these phenotypes. Suggestive associations (p < 5 × 10-5) were also found at 16 additional loci. Moreover, comparison of the porcine results to those reported by human melanoma GWAS indicated shared association signals notably at CDKAL1 and TERT loci but also nearby CCND1, FTO, PLA2G6 and TMEM38B-RAD23B loci. Extensive search of the literature revealed a potential key role of genes at the identified porcine loci in tumor invasion (DST, PLEKHA5, CBY1, LIMK2 and ETV5) and immune response modulation (ETV5, HERC3 and DICER1) of the progression phenotypes. These biological processes are consistent with the clinico-pathological features of MeLiM tumors and can open new routes for future melanoma research in humans.

Rare missense variants in POT1 predispose to familial cutaneous malignant melanoma.

Although CDKN2A is the most frequent high-risk melanoma susceptibility gene, the underlying genetic factors for most melanoma-prone families remain unknown. Using whole-exome sequencing, we identified a rare variant that arose as a founder mutation in the telomere shelterin gene POT1 (chromosome 7, g.124493086C>T; p.Ser270Asn) in five unrelated melanoma-prone families from Romagna, Italy. Carriers of this variant had increased telomere lengths and numbers of fragile telomeres, suggesting that this variant perturbs telomere maintenance. Two additional rare POT1 variants were identified in all cases sequenced in two separate Italian families, one variant per family, yielding a frequency for POT1 variants comparable to that for CDKN2A mutations in this population. These variants were not found in public databases or in 2,038 genotyped Italian controls. We also identified two rare recurrent POT1 variants in US and French familial melanoma cases. Our findings suggest that POT1 is a major susceptibility gene for familial melanoma in several populations.