Canine Oral Melanoma Genomic and Transcriptomic Study Defines Two Molecular Subgroups with Different Therapeutical Targets

Simple Summary

In humans, mucosal melanoma (MM) is a rare and aggressive cancer. The canine model is frequently and spontaneously affected by MM, thus facilitating the collection of samples and the study of its genetic bases. Thanks to an integrative genomic and transcriptomic analysis of 32 canine MM samples, we identified two molecular subgroups of MM with a different microenvironment and structural variant (SV) content. We demonstrated that SVs are associated with recurrently amplified regions, and identified new candidate oncogenes (TRPM7, GABPB1, and SPPL2A) for MM. Our findings suggest the existence of two MM molecular subgroups that could benefit from dedicated therapies, such as immune checkpoint inhibitors or targeted therapies, for both human and veterinary medicine.

Abstract

Mucosal melanoma (MM) is a rare, aggressive clinical cancer. Despite recent advances in genetics and treatment, the prognosis of MM remains poor. Canine MM offers a relevant spontaneous and immunocompetent model to decipher the genetic bases and explore treatments for MM. We performed an integrative genomic and transcriptomic analysis of 32 canine MM samples, which identified two molecular subgroups with a different microenvironment and structural variant (SV) content. The overexpression of genes related to the microenvironment and T-cell response was associated with tumors harboring a lower content of SVs, whereas the overexpression of pigmentation-related pathways and oncogenes, such as TERT, was associated with a high SV burden. Using whole-genome sequencing, we showed that focal amplifications characterized complex chromosomal rearrangements targeting oncogenes, such as MDM2 or CDK4, and a recurrently amplified region on canine chromosome 30. We also demonstrated that the genes TRPM7, GABPB1, and SPPL2A, located in this CFA30 region, play a role in cell proliferation, and thus, may be considered as new candidate oncogenes for human MM. Our findings suggest the existence of two MM molecular subgroups that may benefit from dedicated therapies, such as immune checkpoint inhibitors or targeted therapies, for both human and veterinary medicine.

Circulating tumor DNA is detectable in canine histiocytic sarcoma, oral malignant melanoma, and multicentric lymphoma

Circulating tumor DNA (ctDNA) has become an attractive biomarker in human oncology, and its use may be informative in canine cancer. Thus, we used droplet digital PCR or PCR for antigen receptor rearrangement, to explore tumor-specific point mutations, copy number alterations, and chromosomal rearrangements in the plasma of cancer-affected dogs. We detected ctDNA in 21/23 (91.3%) of histiocytic sarcoma (HS), 2/8 (25%) of oral melanoma, and 12/13 (92.3%) of lymphoma cases. The utility of ctDNA in diagnosing HS was explored in 133 dogs, including 49 with HS, and the screening of recurrent PTPN11 mutations in plasma had a specificity of 98.8% and a sensitivity between 42.8 and 77% according to the clinical presentation of HS. Sensitivity was greater in visceral forms and especially related to pulmonary location. Follow-up of four dogs by targeting lymphoma-specific antigen receptor rearrangement in plasma showed that minimal residual disease detection was concordant with clinical evaluation and treatment response. Thus, our study shows that ctDNA is detectable in the plasma of cancer-affected dogs and is a promising biomarker for diagnosis and clinical follow-up. ctDNA detection appears to be useful in comparative oncology research due to growing interest in the study of natural canine tumors and exploration of new therapies.

Cross-species analysis of enhancer logic using deep learning

Deciphering the genomic regulatory code of enhancers is a key challenge in biology because this code underlies cellular identity. A better understanding of how enhancers work will improve the interpretation of noncoding genome variation and empower the generation of cell type-specific drivers for gene therapy. Here, we explore the combination of deep learning and cross-species chromatin accessibility profiling to build explainable enhancer models. We apply this strategy to decipher the enhancer code in melanoma, a relevant case study owing to the presence of distinct melanoma cell states. We trained and validated a deep learning model, called DeepMEL, using chromatin accessibility data of 26 melanoma samples across six different species. We show the accuracy of DeepMEL predictions on the CAGI5 challenge, where it significantly outperforms existing models on the melanoma enhancer of IRF4 Next, we exploit DeepMEL to analyze enhancer architectures and identify accurate transcription factor binding sites for the core regulatory complexes in the two different melanoma states, with distinct roles for each transcription factor, in terms of nucleosome displacement or enhancer activation. Finally, DeepMEL identifies orthologous enhancers across distantly related species, where sequence alignment fails, and the model highlights specific nucleotide substitutions that underlie enhancer turnover. DeepMEL can be used from the Kipoi database to predict and optimize candidate enhancers and to prioritize enhancer mutations. In addition, our computational strategy can be applied to other cancer or normal cell types.

PTPN11 mutations in canine and human disseminated histiocytic sarcoma

In humans, histiocytic sarcoma (HS) is an aggressive cancer involving histiocytes. Its rarity and heterogeneity explain that treatment remains a challenge. Sharing high clinical and histopathological similarities with human HS, the canine HS is conversely frequent in specific breeds and thus constitutes a unique spontaneous model for human HS to decipher the genetic bases and to explore therapeutic options. We identified sequence alterations in the MAPK pathway in at least 63.9% (71/111) of HS cases with mutually exclusive BRAF (0.9%; 1/111), KRAS (7.2%; 8/111) and PTPN11 (56.75%; 63/111) mutations concentrated at hotspots common to human cancers. Recurrent PTPN11 mutations are associated to visceral disseminated HS subtype in dogs, the most aggressive clinical presentation. We then identified PTPN11 mutations in 3/19 (15.7%) human HS patients. Thus, we propose PTPN11 mutations as key events for a specific subset of human and canine HS: the visceral disseminated form. Finally, by testing drugs targeting the MAPK pathway in eight canine HS cell lines, we identified a better anti-proliferation activity of MEK inhibitors than PTPN11 inhibitors in canine HS neoplastic cells. In combination, these results illustrate the relevance of naturally affected dogs in deciphering genetic mechanisms and selecting efficient targeted therapies for such rare and aggressive cancers in humans.

Prognostic value of somatic focal amplifications on chromosome 30 in canine oral melanoma

Canine oral melanoma is the first malignancy of the oral cavity in dogs and is characterized by a local invasiveness and a high metastatic propensity. A better knowledge of genetic alterations is expected to improve management of this tumour. Copy number alterations are known characteristics of mucosal melanomas both in dogs and humans. The goal of this study was to explore the prognostic value of somatic focal amplifications on chromosomes (Canis Familiaris [CFA]) 10 and 30 in canine oral melanoma. The cohort included 73 dogs with oral melanoma confirmed by histology, removed surgically without adjuvant therapy and with a minimal follow-up of 6 months. Epidemiological, clinical and histological data were collected and quantitative-PCR were performed on formalin-fixed paraffin-embedded (FFPE) samples to identify specific focal amplifications. The 73 dogs included in the study had a median survival time of 220 days. Focal amplifications on CFA 10 and 30 were recurrent (49.3% and 50.7% of cases, respectively) and CFA 30 amplification was significantly associated with the amelanotic phenotype (P = .046) and high mitotic index (MI; P = .0039). CFA 30 amplification was also linked to poor prognosis (P = .0005). Other negative prognostic factors included gingiva location (P = .003), lymphadenomegaly (P = .026), tumour ulceration at diagnosis (P = .003), MI superior to 6 mitoses over 10 fields (P = .001) and amelanotic tumour (P = .029). In multivariate analyses using Cox proportional hazards regression, CFA 30 amplification (Hazard ratio [HR] = 2.08; P = .011), tumour location (HR = 2.20; P = .005) and histological pigmentation (HR = 1.87; P = .036) were significantly associated with shorter survival time. Focal amplification of CFA 30 is linked to an aggressive subset and constitutes a new prognostic factor.

Genome-Wide Analysis of Long Non-Coding RNA Profiles in Canine Oral Melanomas

Mucosal melanomas (MM) are rare aggressive cancers in humans, and one of the most common forms of oral cancers in dogs. Similar biological and histological features are shared between MM in both species, making dogs a powerful model for comparative oncology studies of melanomas. Although exome sequencing recently identified recurrent coding mutations in canine MM, little is known about changes in non-coding gene expression, and more particularly, in canine long non-coding RNAs (lncRNAs), which are commonly dysregulated in human cancers. Here, we sampled a large cohort (n = 52) of canine normal/tumor oral MM from three predisposed breeds (poodles, Labrador retrievers, and golden retrievers), and used deep transcriptome sequencing to identify more than 400 differentially expressed (DE) lncRNAs. We further prioritized candidate lncRNAs by comparative genomic analysis to pinpoint 26 dog–human conserved DE lncRNAs, including SOX21-AS, ZEB2-AS, and CASC15 lncRNAs. Using unsupervised co-expression network analysis with coding genes, we inferred the potential functions of the DE lncRNAs, suggesting associations with cancer-related genes, cell cycle, and carbohydrate metabolism Gene Ontology (GO) terms. Finally, we exploited our multi-breed design to identify DE lncRNAs within breeds. This study provides a unique transcriptomic resource for studying oral melanoma in dogs, and highlights lncRNAs that may potentially be diagnostic or therapeutic targets for human and veterinary medicine.

Identification of a Missense Variant in MFSD12 Involved in Dilution of Phaeomelanin Leading to White or Cream Coat Color in Dogs

White coat color in mammals has been selected several times during the domestication process. Numerous dog breeds are fixed for one form of white coat color that involves darkly pigmented skin. The genetic basis of this color, due to the absence of pigment in the hairs, was suggested to correspond to extreme dilution of the phaeomelanin, by both the expression of only phaeomelanin (locus E) and its extreme dilution (locus I). To go further, we performed genome-wide association studies (GWAS) using a multiple breed approach. The first GWAS, using 34 white dogs and 128 non-white dogs, including White Shepherds, Poodles, Cotons de Tulear and Bichons allowed us to identify two significantly associated loci on the locus E and a novel locus on chromosome 20. A second GWAS using 15 other breeds presenting extreme phaeomelanin dilution confirmed the position of locus I on the chromosome 20 (position 55 Mb pcorrected = 6 × 10-13). Using whole-genome sequencing, we identified a missense variant in the first exon of MFSD12, a gene recently identified to be involved in human, mouse and horse pigmentation. We confirmed the role of this variant in phaeomelanin dilution of numerous canine breeds, and the conserved role of MFSD12 in mammalian pigmentation.

Characterisation and functional predictions of canine long non-coding RNAs

Long non-coding RNAs (lncRNAs) are a family of heterogeneous RNAs that play major roles in multiple biological processes. We recently identified an extended repertoire of more than 10,000 lncRNAs of the domestic dog however, predicting their biological functionality remains challenging. In this study, we have characterised the expression profiles of 10,444 canine lncRNAs in 26 distinct tissue types, representing various anatomical systems. We showed that lncRNA expressions are mainly clustered by tissue type and we highlighted that 44% of canine lncRNAs are expressed in a tissue-specific manner. We further demonstrated that tissue-specificity correlates with specific families of canine transposable elements. In addition, we identified more than 900 conserved dog-human lncRNAs for which we show their overall reproducible expression patterns between dog and human through comparative transcriptomics. Finally, co-expression analyses of lncRNA and neighbouring protein-coding genes identified more than 3,400 canine lncRNAs, suggesting that functional roles of these lncRNAs act as regulatory elements. Altogether, this genomic and transcriptomic integrative study of lncRNAs constitutes a major resource to investigate genotype to phenotype relationships and biomedical research in the dog species.

MKLN1 splicing defect in dogs with lethal acrodermatitis

Lethal acrodermatitis (LAD) is a genodermatosis with monogenic autosomal recessive inheritance in Bull Terriers and Miniature Bull Terriers. The LAD phenotype is characterized by poor growth, immune deficiency, and skin lesions, especially at the paws. Utilizing a combination of genome wide association study and haplotype analysis, we mapped the LAD locus to a critical interval of ~1.11 Mb on chromosome 14. Whole genome sequencing of an LAD affected dog revealed a splice region variant in the MKLN1 gene that was not present in 191 control genomes (chr14:5,731,405T>G or MKLN1:c.400+3A>C). This variant showed perfect association in a larger combined Bull Terrier/Miniature Bull Terrier cohort of 46 cases and 294 controls. The variant was absent from 462 genetically diverse control dogs of 62 other dog breeds. RT-PCR analysis of skin RNA from an affected and a control dog demonstrated skipping of exon 4 in the MKLN1 transcripts of the LAD affected dog, which leads to a shift in the MKLN1 reading frame. MKLN1 encodes the widely expressed intracellular protein muskelin 1, for which diverse functions in cell adhesion, morphology, spreading, and intracellular transport processes are discussed. While the pathogenesis of LAD remains unclear, our data facilitate genetic testing of Bull Terriers and Miniature Bull Terriers to prevent the unintentional production of LAD affected dogs. This study may provide a starting point to further clarify the elusive physiological role of muskelin 1 in vivo.

Lethal acrodermatitis (LAD) is an autosomal recessive hereditary disease in dogs. It is characterized by poor growth, immune deficiency and characteristic skin lesions of the paws and of the face. We mapped the LAD locus to a ~1.11 Mb segment on canine chromosome 14. Whole genome sequence data of an LAD affected dog and 191 controls revealed a candidate causative variant in the MKLN1 gene, encoding muskelin 1. The identified variant, a single nucleotide substitution, MKLN1:c.400+3A>C, altered the 5’-splice site at the beginning of intron 4. We experimentally confirmed that this variant leads to complete skipping of exon 4 in the MKLN1 mRNA in skin. Various cellular functions have been postulated for muskelin 1 including roles in intracellular transport processes, cell morphology, cell spreading, and cell adhesion. Our data from dogs reveal a novel in vivo role for muskelin 1 that is related to the immune system and skin. MKLN1 thus represents a novel candidate gene for human patients with unsolved acrodermatitis and/or immune deficiency phenotypes. LAD affected dogs may serve as models to gain more insights into the function of muskelin 1.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

Abstract IA08: Analysis of a histiocytic carcinoma in the dog: The canine's utility in fetching disease genes

Histiocytic sarcoma (HS) is a rare human hematopoietic neoplasm that accounts for less than 1% of hematolymphoid cancers and occurs either concurrently with or subsequently to B- or T-lymphoblastic lymphoma and other B-cell cancers. As such, its genetic etiology has been difficult to establish. HS in dogs, by comparison, shows strong breed specificity with 20% of flat-coated retrievers and 25% of Bernese mountain dogs (BMD) estimated to develop the disease in their lifetime. These two affected breeds share no recent common ancestors and display differences in disease presentation, offering a unique system in which to investigate the genetic predisposition and progression of HS. Our lab has previously undertaken genome wide association study (GWAS) on both BMD and FCR collected in North America and Europe and identified four significantly associated loci, two in BMD and two in FCR. In each breed, only one locus was segregating with disease in the American population while both were present in the European populations. There were no loci in common between the two breeds. Fine mapping of the common BMD HS locus on canine chromosome (CFA) CFA11 revealed a single haplotype spanning the MTAP (methylthioadenosine phosphorylase) gene and part of the cyclin-dependent kinase inhibitor 2A (CDKN2A) gene, which was present in 96% of BMD with HS. The European BMD also showed a significant locus on CFA14 that may be approaching fixation in the American dogs based on allele frequency calculations. We have done targeted resequencing of the CFA14 locus in a subset of European BMD and are in the process of filtering the identified mutations. Promising candidate genes in this region include GRM8, which inhibits the cyclic AMP cascade in the central nervous system and has shown to be involved in several of cancers; POT1, which is involved in telomere length maintenance and protection; and HYAL4, a hyaluronidase potentially involved in extracellular matrix degradation. In addition, the laboratory has generated 30x sequence from two HS cases and a control BMD, two FCR cases, and a control and low pass sequence from 16 dogs (8 cases and controls) in order to identify the likely causative mutations in each of the associated loci as well as breed specific mutations that may predispose the dogs to aggressive cancers.

Citation Format: Heidi G. Parker, Edouard Cadieu, Abigail Shearin, Gerard Rutteman, Elaine A. Ostrander. Analysis of a histiocytic carcinoma in the dog: The canine's utility in fetching disease genes [abstract]. In: Proceedings of the AACR Conference on Advances in Sarcomas: From Basic Science to Clinical Translation; May 16-19, 2017; Philadelphia, PA. Philadelphia (PA): AACR; Clin Cancer Res 2018;24(2_Suppl):Abstract nr IA08.

An exome sequencing based approach for genome-wide association studies in the dog

Genome-wide association studies (GWAS) are widely used to identify loci associated with phenotypic traits in the domestic dog that has emerged as a model for Mendelian and complex traits. However, a disadvantage of GWAS is that it always requires subsequent fine-mapping or sequencing to pinpoint causal mutations. Here, we performed whole exome sequencing (WES) and canine high-density (cHD) SNP genotyping of 28 dogs from 3 breeds to compare the SNP and linkage disequilibrium characteristics together with the power and mapping precision of exome-guided GWAS (EG-GWAS) versus cHD-based GWAS. Using simulated phenotypes, we showed that EG-GWAS has a higher power than cHD to detect associations within target regions and less power outside target regions, with power being influenced further by sample size and SNP density. We analyzed two real phenotypes (hair length and furnishing), that are fixed in certain breeds to characterize mapping precision of the known causal mutations. EG-GWAS identified the associated exonic and 3'UTR variants within the FGF5 and RSPO2 genes, respectively, with only a few samples per breed. In conclusion, we demonstrated that EG-GWAS can identify loci associated with Mendelian phenotypes both within and across breeds.

FEELnc: a tool for long non-coding RNA annotation and its application to the dog transcriptome

Whole transcriptome sequencing (RNA-seq) has become a standard for cataloguing and monitoring RNA populations. One of the main bottlenecks, however, is to correctly identify the different classes of RNAs among the plethora of reconstructed transcripts, particularly those that will be translated (mRNAs) from the class of long non-coding RNAs (lncRNAs). Here, we present FEELnc (FlExible Extraction of LncRNAs), an alignment-free program that accurately annotates lncRNAs based on a Random Forest model trained with general features such as multi k-mer frequencies and relaxed open reading frames. Benchmarking versus five state-of-the-art tools shows that FEELnc achieves similar or better classification performance on GENCODE and NONCODE data sets. The program also provides specific modules that enable the user to fine-tune classification accuracy, to formalize the annotation of lncRNA classes and to identify lncRNAs even in the absence of a training set of non-coding RNAs. We used FEELnc on a real data set comprising 20 canine RNA-seq samples produced by the European LUPA consortium to substantially expand the canine genome annotation to include 10 374 novel lncRNAs and 58 640 mRNA transcripts. FEELnc moves beyond conventional coding potential classifiers by providing a standardized and complete solution for annotating lncRNAs and is freely available at https://github.com/tderrien/FEELnc.

Naturally occurring melanomas in dogs as models for non-UV pathways of human melanomas

Spontaneously occurring melanomas are frequent in dogs. They appear at the same localizations as in humans, i.e. skin, mucosal sites, nail matrix and eyes. They display variable behaviors: tumors at oral localizations are more frequent and aggressive than at other anatomical sites. Interestingly, dog melanomas are associated with strong breed predispositions and overrepresentation of black-coated dogs. Epidemiological analysis of 2350 affected dogs showed that poodles are at high risk of developing oral melanoma, while schnauzers or Beauce shepherds mostly developped cutaneous melanoma. Clinical and histopathological analyses were performed on a cohort of 153 cases with a 4-yr follow-up. Histopathological characterization showed that most canine tumors are intradermal and homologous to human rare morphological melanomas types - 'nevocytoid type' and 'animal type'-. Tumor cDNA sequencing data, obtained from 95 dogs for six genes, relevant to human melanoma classification, detected somatic mutations in oral melanoma, in NRAS and PTEN genes, at human hotspot sites, but not in BRAF. Altogether, these findings support the relevance of the dog model for comparative oncology of melanomas, especially for the elucidation of non-UV induced pathways.

The MTAP-CDKN2A locus confers susceptibility to a naturally occurring canine cancer

BACKGROUND

Advantages offered by canine population substructure, combined with clinical presentations similar to human disorders, makes the dog an attractive system for studies of cancer genetics. Cancers that have been difficult to study in human families or populations are of particular interest. Histiocytic sarcoma is a rare and poorly understood neoplasm in humans that occurs in 15% to 25% of Bernese Mountain Dogs (BMD).

METHODS

Genomic DNA was collected from affected and unaffected BMD in North America and Europe. Both independent and combined genome-wide association studies (GWAS) were used to identify cancer-associated loci. Fine mapping and sequencing narrowed the primary locus to a single gene region.

RESULTS

Both populations shared the same primary locus, which features a single haplotype spanning MTAP and part of CDKN2A and is present in 96% of affected BMD. The haplotype is within the region homologous to human chromosome 9p21, which has been implicated in several types of cancer.

CONCLUSIONS

We present the first GWAS for histiocytic sarcoma in any species. The data identify an associated haplotype in the highly cited tumor suppressor locus near CDKN2A. These data show the power of studying distinctive malignancies in highly predisposed dog breeds.

IMPACT

Here, we establish a naturally occurring model of cancer susceptibility due to CDKN2 dysregulation, thus providing insight about this cancer-associated, complex, and poorly understood genomic region.

An insertion in the RSPO2 gene correlates with improper coat in the Portuguese water dog

We recently showed that genes at 3 loci account for the majority of variation in canine fur. Allelic variation at genes controlling length of fur, texture, and curl is responsible for the striking phenotypic variety observed among purebred dogs in the United States today. In this paper, we investigate the phenomenon of "improper coat" (IC) or a coat that is not typical of the breed. IC is occasionally observed among specific breeds, such as the Portuguese Water Dog (PWD), and is characterized by short hair on the head, face, and lower legs, rather than a thick and even coat covering the whole body. The IC is reminiscent of that observed on the curly or flat-coated retriever, thus making such dogs unable to compete effectively in conformation events. We have found that the presence of the wild-type allele, rather than the expected variant allele at the R-spondin 2 (RSPO2) gene, accounts for this phenotype. The development of a genetic test that distinguishes these 2 allelic types would allow breeders to easily avoid producing PWD with ICs.

Coat variation in the domestic dog is governed by variants in three genes

Coat color and type are essential characteristics of domestic dog breeds. Although the genetic basis of coat color has been well characterized, relatively little is known about the genes influencing coat growth pattern, length, and curl. We performed genome-wide association studies of more than 1000 dogs from 80 domestic breeds to identify genes associated with canine fur phenotypes. Taking advantage of both inter- and intrabreed variability, we identified distinct mutations in three genes, RSPO2, FGF5, and KRT71 (encoding R-spondin-2, fibroblast growth factor-5, and keratin-71, respectively), that together account for most coat phenotypes in purebred dogs in the United States. Thus, an array of varied and seemingly complex phenotypes can be reduced to the combinatorial effects of only a few genes.

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.

Facilitating genome navigation: survey sequencing and dense radiation-hybrid gene mapping

Accurate and comprehensive sequence coverage for large genomes has been restricted to only a few species of specific interest. Lower sequence coverage (survey sequencing) of related species can yield a wealth of information about gene content and putative regulatory elements. But survey sequences lack long-range continuity and provide only a fragmented view of a genome. Here we show the usefulness of combining survey sequencing with dense radiation-hybrid (RH) maps for extracting maximum comparative genome information from model organisms. Based on results from the canine system, we propose that from now on all low-pass sequencing projects should be accompanied by a dense, gene-based RH map-construction effort to extract maximum information from the genome with a marginal extra cost.

An integrated 4249 marker FISH/RH map of the canine genome

BACKGROUND

The 156 breeds of dog recognized by the American Kennel Club offer a unique opportunity to map genes important in genetic variation. Each breed features a defining constellation of morphological and behavioral traits, often generated by deliberate crossing of closely related individuals, leading to a high rate of genetic disease in many breeds. Understanding the genetic basis of both phenotypic variation and disease susceptibility in the dog provides new ways in which to dissect the genetics of human health and biology.

RESULTS

To facilitate both genetic mapping and cloning efforts, we have constructed an integrated canine genome map that is both dense and accurate. The resulting resource encompasses 4249 markers, and was constructed using the RHDF5000-2 whole genome radiation hybrid panel. The radiation hybrid (RH) map features a density of one marker every 900 Kb and contains 1760 bacterial artificial chromosome clones (BACs) localized to 1423 unique positions, 851 of which have also been mapped by fluorescence in situ hybridization (FISH). The two data sets show excellent concordance. Excluding the Y chromosome, the map features an RH/FISH mapped BAC every 3.5 Mb and an RH mapped BAC-end, on average, every 2 Mb. For 2233 markers, the orthologous human genes have been established, allowing the identification of 79 conserved segments (CS) between the dog and human genomes, dramatically extending the length of most previously described CS.

CONCLUSIONS

These results provide a necessary resource for the canine genome mapping community to undertake positional cloning experiments and provide new insights into the comparative canine-human genome maps.

Comparison of the canine and human olfactory receptor gene repertoires

In this study, 817 novel canine olfactory receptor (OR) sequences were identified, and 640 have been characterized. Of the 661 characterized OR sequences, representing half of the canine repertoire, 18% are predicted to be pseudogenes, compared with 63% in human and 20% in mouse.

Background

Olfactory receptors (ORs), the first dedicated molecules with which odorants physically interact to arouse an olfactory sensation, constitute the largest gene family in vertebrates, including around 900 genes in human and 1,500 in the mouse. Whereas dogs, like many other mammals, have a much keener olfactory potential than humans, only 21 canine OR genes have been described to date.

Results

In this study, 817 novel canine OR sequences were identified, and 640 have been characterized. Of the 661 characterized OR sequences, representing half of the canine repertoire, 18% are predicted to be pseudogenes, compared with 63% in human and 20% in mouse. Phylogenetic analysis of 403 canine OR sequences identified 51 families, and radiation-hybrid mapping of 562 showed that they are distributed on 24 dog chromosomes, in 37 distinct regions. Most of these regions constitute clusters of 2 to 124 closely linked genes. The two largest clusters (124 and 109 OR genes) are located on canine chromosomes 18 and 21. They are orthologous to human clusters located on human chromosomes 11q11-q13 and HSA11p15, containing 174 and 115 ORs respectively.

Conclusions

This study shows a strongly conserved genomic distribution of OR genes between dog and human, suggesting that OR genes evolved from a common mammalian ancestral repertoire by successive duplications. In addition, the dog repertoire appears to have expanded relative to that of humans, leading to the emergence of specific canine OR genes.

A 1-Mb resolution radiation hybrid map of the canine genome

The purebred dog population consists of >300 partially inbred genetic isolates or breeds. Restriction of gene flow between breeds, together with strong selection for traits, has led to the establishment of a unique resource for dissecting the genetic basis of simple and complex mammalian traits. Toward this end, we present a comprehensive radiation hybrid map of the canine genome composed of 3,270 markers including 1,596 microsatellite-based markers, 900 cloned gene sequences and ESTs, 668 canine-specific bacterial artificial chromosome (BAC) ends, and 106 sequence-tagged sites. The map was constructed by using the RHDF5000-2 whole-genome radiation hybrid panel and computed by using MULTIMAP and TSP/CONCORDE. The 3,270 markers map to 3,021 unique positions and define an average intermarker distance corresponding to 1 Mb. We also define a minimal screening set of 325 highly informative well spaced markers, to be used in the initiation of genome-wide scans. The well defined synteny between the dog and human genomes, established in part as a function of this work by the identification of 85 conserved fragments, will allow follow-up of initial findings of linkage by selection of candidate genes from the human genome sequence. This work continues to define the canine system as the method of choice in the pursuit of the genes causing mammalian variation and disease.

Comparison of MultiMap and TSP/CONCORDE for constructing radiation hybrid maps

Radiation hybrid (RH) map construction allows investigators to locate both type I and type II markers on a given genome map. The process is composed of two steps. The first consists of determining the pattern distribution of a set of markers within the different cell lines of an RH panel. This is mainly done by polymerase chain reaction (PCR) amplification and gel electrophoresis, and results in a series of numbers indicating the presence or the absence of each marker in each cell line. The second step consists of a comparison of these numbers, using various algorithms, to group and then order markers. Because different algorithms may provide (slightly) different orders, we have compared the merits of the MultiMap and TSP/CONCORDE packages using a data set of information currently under analysis for construction of the canine genome RH map.

The genome sequence of Schizosaccharomyces pombe

We have sequenced and annotated the genome of fission yeast (Schizosaccharomyces pombe), which contains the smallest number of protein-coding genes yet recorded for a eukaryote: 4,824. The centromeres are between 35 and 110 kilobases (kb) and contain related repeats including a highly conserved 1.8-kb element. Regions upstream of genes are longer than in budding yeast (Saccharomyces cerevisiae), possibly reflecting more-extended control regions. Some 43% of the genes contain introns, of which there are 4,730. Fifty genes have significant similarity with human disease genes; half of these are cancer related. We identify highly conserved genes important for eukaryotic cell organization including those required for the cytoskeleton, compartmentation, cell-cycle control, proteolysis, protein phosphorylation and RNA splicing. These genes may have originated with the appearance of eukaryotic life. Few similarly conserved genes that are important for multicellular organization were identified, suggesting that the transition from prokaryotes to eukaryotes required more new genes than did the transition from unicellular to multicellular organization.

Mapping and identification of essential gene functions on the X chromosome of Drosophila

The Drosophila melanogaster genome consists of four chromosomes that contain 165 Mb of DNA, 120 Mb of which are euchromatic. The two Drosophila Genome Projects, in collaboration with Celera Genomics Systems, have sequenced the genome, complementing the previously established physical and genetic maps. In addition, the Berkeley Drosophila Genome Project has undertaken large-scale functional analysis based on mutagenesis by transposable P element insertions into autosomes. Here, we present a large-scale P element insertion screen for vital gene functions and a BAC tiling map for the X chromosome. A collection of 501 X-chromosomal P element insertion lines was used to map essential genes cytogenetically and to establish short sequence tags (STSs) linking the insertion sites to the genome. The distribution of the P element integration sites, the identified genes and transcription units as well as the expression patterns of the P-element-tagged enhancers is described and discussed.

Analysis of the chromosome sequence of the legume symbiont Sinorhizobium meliloti strain 1021

Sinorhizobium meliloti is an alpha-proteobacterium that forms agronomically important N(2)-fixing root nodules in legumes. We report here the complete sequence of the largest constituent of its genome, a 62.7% GC-rich 3,654,135-bp circular chromosome. Annotation allowed assignment of a function to 59% of the 3,341 predicted protein-coding ORFs, the rest exhibiting partial, weak, or no similarity with any known sequence. Unexpectedly, the level of reiteration within this replicon is low, with only two genes duplicated with more than 90% nucleotide sequence identity, transposon elements accounting for 2.2% of the sequence, and a few hundred short repeated palindromic motifs (RIME1, RIME2, and C) widespread over the chromosome. Three regions with a significantly lower GC content are most likely of external origin. Detailed annotation revealed that this replicon contains all housekeeping genes except two essential genes that are located on pSymB. Amino acid/peptide transport and degradation and sugar metabolism appear as two major features of the S. meliloti chromosome. The presence in this replicon of a large number of nucleotide cyclases with a peculiar structure, as well as of genes homologous to virulence determinants of animal and plant pathogens, opens perspectives in the study of this bacterium both as a free-living soil microorganism and as a plant symbiont.

The composite genome of the legume symbiont Sinorhizobium meliloti

The scarcity of usable nitrogen frequently limits plant growth. A tight metabolic association with rhizobial bacteria allows legumes to obtain nitrogen compounds by bacterial reduction of dinitrogen (N2) to ammonium (NH4+). We present here the annotated DNA sequence of the alpha-proteobacterium Sinorhizobium meliloti, the symbiont of alfalfa. The tripartite 6.7-megabase (Mb) genome comprises a 3.65-Mb chromosome, and 1.35-Mb pSymA and 1.68-Mb pSymB megaplasmids. Genome sequence analysis indicates that all three elements contribute, in varying degrees, to symbiosis and reveals how this genome may have emerged during evolution. The genome sequence will be useful in understanding the dynamics of interkingdom associations and of life in soil environments.

From first base: the sequence of the tip of the X chromosome of Drosophila melanogaster, a comparison of two sequencing strategies

We present the sequence of a contiguous 2.63 Mb of DNA extending from the tip of the X chromosome of Drosophila melanogaster. Within this sequence, we predict 277 protein coding genes, of which 94 had been sequenced already in the course of studying the biology of their gene products, and examples of 12 different transposable elements. We show that an interval between bands 3A2 and 3C2, believed in the 1970s to show a correlation between the number of bands on the polytene chromosomes and the 20 genes identified by conventional genetics, is predicted to contain 45 genes from its DNA sequence. We have determined the insertion sites of P-elements from 111 mutant lines, about half of which are in a position likely to affect the expression of novel predicted genes, thus representing a resource for subsequent functional genomic analysis. We compare the European Drosophila Genome Project sequence with the corresponding part of the independently assembled and annotated Joint Sequence determined through "shotgun" sequencing. Discounting differences in the distribution of known transposable elements between the strains sequenced in the two projects, we detected three major sequence differences, two of which are probably explained by errors in assembly; the origin of the third major difference is unclear. In addition there are eight sequence gaps within the Joint Sequence. At least six of these eight gaps are likely to be sites of transposable elements; the other two are complex. Of the 275 genes in common to both projects, 60% are identical within 1% of their predicted amino-acid sequence and 31% show minor differences such as in choice of translation initiation or termination codons; the remaining 9% show major differences in interpretation.

From sequence to chromosome: the tip of the X chromosome of D. melanogaster

One of the rewards of having a Drosophila melanogaster whole-genome sequence will be the potential to understand the molecular bases for structural features of chromosomes that have been a long-standing puzzle. Analysis of 2.6 megabases of sequence from the tip of the X chromosome of Drosophila identifies 273 genes. Cloned DNAs from the characteristic bulbous structure at the tip of the X chromosome in the region of the broad complex display an unusual pattern of in situ hybridization. Sequence analysis revealed that this region comprises 154 kilobases of DNA flanked by 1.2-kilobases of inverted repeats, each composed of a 350-base pair satellite related element. Thus, some aspects of chromosome structure appear to be revealed directly within the DNA sequence itself.

The genome sequence of Drosophila melanogaster

The fly Drosophila melanogaster is one of the most intensively studied organisms in biology and serves as a model system for the investigation of many developmental and cellular processes common to higher eukaryotes, including humans. We have determined the nucleotide sequence of nearly all of the approximately 120-megabase euchromatic portion of the Drosophila genome using a whole-genome shotgun sequencing strategy supported by extensive clone-based sequence and a high-quality bacterial artificial chromosome physical map. Efforts are under way to close the remaining gaps; however, the sequence is of sufficient accuracy and contiguity to be declared substantially complete and to support an initial analysis of genome structure and preliminary gene annotation and interpretation. The genome encodes approximately 13,600 genes, somewhat fewer than the smaller Caenorhabditis elegans genome, but with comparable functional diversity.

Disruption of six novel yeast genes reveals three genes essential for vegetative growth and one required for growth at low temperature

We describe here the construction of six deletion mutants and their basic phenotypic analysis. Six open reading frames (ORFs) from chromosome X, YJR039w, YJR041c, YJR043c, YJR046w, YJR053w and YJR065c, were disrupted by deletion cassettes with long (LFH) or short (SFH) flanking regions homologous to the target locus. The LFH deletion cassette was made by introducing into the kanMX4 marker module two polymerase chain reaction (PCR) fragments several hundred base pairs (bp) in size homologous to the promoter and terminator regions of a given ORF. The SFH gene disruption construct was obtained by PCR amplification of the kanMX4 marker with primers providing homology to the target gene. The region of homology to mediate homologous recombination was about 70 bp. Sporulation and tetrad analysis revealed that ORFs YJR041c, YJR046w and YJR065c are essential genes. Complementation tests by corresponding cognate gene clones confirmed this observation. The non-growing haploid segregants were observed under the microscope. The yjr041c delta haploid cells gave rise to microcolonies comprising about 20 to 50 cells. Most yjr046w delta cells were blocked after one or two cell cycles with heterogeneous bud sizes. The yjr065c delta cells displayed an unbudded spore or were arrested before completion of the first cell division cycle with a bud of variable size. The deduced protein of ORF YJR065c, that we named Act4, belongs to the Arp3 family of actin-related proteins. Three other ORFs, YJR039w, YJR043c and YJR053w are non-essential genes. The yjr043c delta cells hardly grew at 15 degrees C, indicating that this gene is required for growth at low temperature. Complementation tests confirmed that the disruption of YJR043c is responsible for this growth defect. In addition, the mating efficiency of yjr043c delta and yjr053w delta cells appear to be moderately affected.