Evolution of the ancestral mammalian karyotype and syntenic regions.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

The genomic substrate for adaptive radiation in African cichlid fish.

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Low guanylyl cyclase activity in Weddell seals: implications for peripheral vasoconstriction and perfusion of the brain during diving.

Nitric oxide (NO) is a potent vasodilator, which improves perfusion and oxygen delivery during tissue hypoxia in terrestrial animals. The vertebrate dive response involves vasoconstriction in select tissues, which persists despite profound hypoxia. Using tissues collected from Weddell seals at necropsy, we investigated whether vasoconstriction is aided by downregulation of local hypoxia signaling mechanisms. We focused on NO-soluble guanylyl cyclase (GC)-cGMP signaling, a well-known vasodilatory transduction pathway. Seals have a lower GC protein abundance, activity, and capacity to respond to NO stimulation than do terrestrial mammals. In seal lung homogenates, GC produced less cGMP (20.1 ± 3.7 pmol·mg protein-1·min-1) than the lungs of dogs (-80 ± 144 pmol·mg protein-1·min-1 less than seals), sheep (-472 ± 96), rats (-664 ± 104) or mice (-1,160 ± 104, P < 0.0001). Amino acid sequences of the GC enzyme α-subunits differed between seals and terrestrial mammals, potentially affecting their structure and function. Vasoconstriction in diving Weddell seals is not consistent across tissues; perfusion is maintained in the brain and heart but decreased in other organs such as the kidney. A NO donor increased median GC activity 49.5-fold in the seal brain but only 27.4-fold in the kidney, consistent with the priority of cerebral perfusion during diving. Nos3 expression was high in the seal brain, which could improve NO production and vasodilatory potential. Conversely, Pde5a expression was high in the seal renal artery, which may increase cGMP breakdown and vasoconstriction in the kidney. Taken together, the results of this study suggest that alterations in the NO-cGMP pathway facilitate the diving response.

Exome sequencing of lymphomas from three dog breeds reveals somatic mutation patterns reflecting genetic background.

Lymphoma is the most common hematological malignancy in developed countries. Outcome is strongly determined by molecular subtype, reflecting a need for new and improved treatment options. Dogs spontaneously develop lymphoma, and the predisposition of certain breeds indicates genetic risk factors. Using the dog breed structure, we selected three lymphoma predisposed breeds developing primarily T-cell (boxer), primarily B-cell (cocker spaniel), and with equal distribution of B- and T-cell lymphoma (golden retriever), respectively. We investigated the somatic mutations in B- and T-cell lymphomas from these breeds by exome sequencing of tumor and normal pairs. Strong similarities were evident between B-cell lymphomas from golden retrievers and cocker spaniels, with recurrent mutations in TRAF3-MAP3K14 (28% of all cases), FBXW7 (25%), and POT1 (17%). The FBXW7 mutations recurrently occur in a specific codon; the corresponding codon is recurrently mutated in human cancer. In contrast, T-cell lymphomas from the predisposed breeds, boxers and golden retrievers, show little overlap in their mutation pattern, sharing only one of their 15 most recurrently mutated genes. Boxers, which develop aggressive T-cell lymphomas, are typically mutated in the PTEN-mTOR pathway. T-cell lymphomas in golden retrievers are often less aggressive, and their tumors typically showed mutations in genes involved in cellular metabolism. We identify genes with known involvement in human lymphoma and leukemia, genes implicated in other human cancers, as well as novel genes that could allow new therapeutic options.

Genome-wide association study identifies shared risk loci common to two malignancies in golden retrievers.

Dogs, with their breed-determined limited genetic background, are great models of human disease including cancer. Canine B-cell lymphoma and hemangiosarcoma are both malignancies of the hematologic system that are clinically and histologically similar to human B-cell non-Hodgkin lymphoma and angiosarcoma, respectively. Golden retrievers in the US show significantly elevated lifetime risk for both B-cell lymphoma (6%) and hemangiosarcoma (20%). We conducted genome-wide association studies for hemangiosarcoma and B-cell lymphoma, identifying two shared predisposing loci. The two associated loci are located on chromosome 5, and together contribute ~20% of the risk of developing these cancers. Genome-wide p-values for the top SNP of each locus are 4.6×10-7 and 2.7×10-6, respectively. Whole genome resequencing of nine cases and controls followed by genotyping and detailed analysis identified three shared and one B-cell lymphoma specific risk haplotypes within the two loci, but no coding changes were associated with the risk haplotypes. Gene expression analysis of B-cell lymphoma tumors revealed that carrying the risk haplotypes at the first locus is associated with down-regulation of several nearby genes including the proximal gene TRPC6, a transient receptor Ca2+-channel involved in T-cell activation, among other functions. The shared risk haplotype in the second locus overlaps the vesicle transport and release gene STX8. Carrying the shared risk haplotype is associated with gene expression changes of 100 genes enriched for pathways involved in immune cell activation. Thus, the predisposing germ-line mutations in B-cell lymphoma and hemangiosarcoma appear to be regulatory, and affect pathways involved in T-cell mediated immune response in the tumor. This suggests that the interaction between the immune system and malignant cells plays a common role in the tumorigenesis of these relatively different cancers.

The genome of the green anole lizard and a comparative analysis with birds and mammals.

The evolution of the amniotic egg was one of the great evolutionary innovations in the history of life, freeing vertebrates from an obligatory connection to water and thus permitting the conquest of terrestrial environments. Among amniotes, genome sequences are available for mammals and birds, but not for non-avian reptiles. Here we report the genome sequence of the North American green anole lizard, Anolis carolinensis. We find that A. carolinensis microchromosomes are highly syntenic with chicken microchromosomes, yet do not exhibit the high GC and low repeat content that are characteristic of avian microchromosomes. Also, A. carolinensis mobile elements are very young and diverse-more so than in any other sequenced amniote genome. The GC content of this lizard genome is also unusual in its homogeneity, unlike the regionally variable GC content found in mammals and birds. We describe and assign sequence to the previously unknown A. carolinensis X chromosome. Comparative gene analysis shows that amniote egg proteins have evolved significantly more rapidly than other proteins. An anole phylogeny resolves basal branches to illuminate the history of their repeated adaptive radiations.

Genome-wide immunity studies in the rabbit: transcriptome variations in peripheral blood mononuclear cells after in vitro stimulation by LPS or PMA-Ionomycin.

BACKGROUND: Our purpose was to obtain genome-wide expression data for the rabbit species on the responses of peripheral blood mononuclear cells (PBMCs) after in vitro stimulation by lipopolysaccharide (LPS) or phorbol myristate acetate (PMA) and ionomycin. This transcriptome profiling was carried out using microarrays enriched with immunity-related genes, and annotated with the most recent data available for the rabbit genome. RESULTS: The LPS affected 15 to 20 times fewer genes than PMA-Ionomycin after both 4 hours (T4) and 24 hours (T24), of in vitro stimulation, in comparison with mock-stimulated PBMCs. LPS induced an inflammatory response as shown by a significant up-regulation of IL12A and CXCL11 at T4, followed by an increased transcription of IL6, IL1B, IL1A, IL36, IL37, TNF, and CCL4 at T24. Surprisingly, we could not find an up-regulation of IL8 either at T4 or at T24, and detected a down-regulation of DEFB1 and BPI at T24. A concerted up-regulation of SAA1, S100A12 and F3 was found upon stimulation by LPS. PMA-Ionomycin induced a very early expression of Th1, Th2, Treg, and Th17 responses by PBMCs at T4. The Th1 response increased at T24 as shown by the increase of the transcription of IFNG and by contrast to other cytokines which significantly decreased from T4 to T24 (IL2, IL4, IL10, IL13, IL17A, CD69) by comparison to mock-stimulation. The granulocyte-macrophage colony-stimulating factor (CSF2) was by far the most over-expressed gene at both T4 and T24 by comparison to mock-stimulated cells, confirming a major impact of PMA-Ionomycin on cell growth and proliferation. A significant down-regulation of IL16 was observed at T4 and T24, in agreement with a role of IL16 in PBMC apoptosis. CONCLUSIONS: We report new data on the responses of PBMCs to LPS and PMA-Ionomycin in the rabbit species, thus enlarging the set of mammalian species for which such reports exist. The availability of the rabbit genome assembly together with high throughput genomic tools should pave the way for more intense genomic studies for this species, which is known to be a very relevant biomedical model in immunology and physiology.

The genomic landscape of canine osteosarcoma cell lines reveals conserved structural complexity and pathway alterations.

The characterization of immortalized canine osteosarcoma (OS) cell lines used for research has historically been based on phenotypic features such as cellular morphology and expression of bone specific markers. With the increasing use of these cell lines to investigate novel therapeutic approaches prior to in vivo translation, a much more detailed understanding regarding the genomic landscape of these lines is required to ensure accurate interpretation of findings. Here we report the first whole genome characterization of eight canine OS cell lines, including single nucleotide variants, copy number variants and other structural variants. Many alterations previously characterized in primary canine OS tissue were observed in these cell lines, including TP53 mutations, MYC copy number gains, loss of CDKN2A, PTEN, DLG2, MAGI2, and RB1 and structural variants involving SETD2, DLG2 and DMD. These data provide a new framework for understanding how best to incorporate in vitro findings generated using these cell lines into the design of future clinical studies involving dogs with spontaneous OS.

The Antarctic Weddell seal genome reveals evidence of selection on cardiovascular phenotype and lipid handling.

The Weddell seal (Leptonychotes weddellii) thrives in its extreme Antarctic environment. We generated the Weddell seal genome assembly and a high-quality annotation to investigate genome-wide evolutionary pressures that underlie its phenotype and to study genes implicated in hypoxia tolerance and a lipid-based metabolism. Genome-wide analyses included gene family expansion/contraction, positive selection, and diverged sequence (acceleration) compared to other placental mammals, identifying selection in coding and non-coding sequence in five pathways that may shape cardiovascular phenotype. Lipid metabolism as well as hypoxia genes contained more accelerated regions in the Weddell seal compared to genomic background. Top-significant genes were SUMO2 and EP300; both regulate hypoxia inducible factor signaling. Liver expression of four genes with the strongest acceleration signals differ between Weddell seals and a terrestrial mammal, sheep. We also report a high-density lipoprotein-like particle in Weddell seal serum not present in other mammals, including the shallow-diving harbor seal.

Ancestry-inclusive dog genomics challenges popular breed stereotypes.

Behavioral genetics in dogs has focused on modern breeds, which are isolated subgroups with distinctive physical and, purportedly, behavioral characteristics. We interrogated breed stereotypes by surveying owners of 18,385 purebred and mixed-breed dogs and genotyping 2155 dogs. Most behavioral traits are heritable [heritability (h2) > 25%], and admixture patterns in mixed-breed dogs reveal breed propensities. Breed explains just 9% of behavioral variation in individuals. Genome-wide association analyses identify 11 loci that are significantly associated with behavior, and characteristic breed behaviors exhibit genetic complexity. Behavioral loci are not unusually differentiated in breeds, but breed propensities align, albeit weakly, with ancestral function. We propose that behaviors perceived as characteristic of modern breeds derive from thousands of years of polygenic adaptation that predates breed formation, with modern breeds distinguished primarily by aesthetic traits.

Genomically Complex Human Angiosarcoma and Canine Hemangiosarcoma Establish Convergent Angiogenic Transcriptional Programs Driven by Novel Gene Fusions.

Sporadic angiosarcomas are aggressive vascular sarcomas whose rarity and genomic complexity present significant obstacles in deciphering the pathogenic significance of individual genetic alterations. Numerous fusion genes have been identified across multiple types of cancers, but their existence and significance remain unclear in sporadic angiosarcomas. In this study, we leveraged RNA-sequencing data from 13 human angiosarcomas and 76 spontaneous canine hemangiosarcomas to identify fusion genes associated with spontaneous vascular malignancies. Ten novel protein-coding fusion genes, including TEX2-PECAM1 and ATP8A2-FLT1, were identified in seven of the 13 human tumors, with two tumors showing mutations of TP53. HRAS and NRAS mutations were found in angiosarcomas without fusions or TP53 mutations. We found 15 novel protein-coding fusion genes including MYO16-PTK2, GABRA3-FLT1, and AKT3-XPNPEP1 in 11 of the 76 canine hemangiosarcomas; these fusion genes were seen exclusively in tumors of the angiogenic molecular subtype that contained recurrent mutations in TP53, PIK3CA, PIK3R1, and NRAS. In particular, fusion genes and mutations of TP53 cooccurred in tumors with higher frequency than expected by random chance, and they enriched gene signatures predicting activation of angiogenic pathways. Comparative transcriptomic analysis of human angiosarcomas and canine hemangiosarcomas identified shared molecular signatures associated with activation of PI3K/AKT/mTOR pathways. Our data suggest that genome instability induced by TP53 mutations might create a predisposition for fusion events that may contribute to tumor progression by promoting selection and/or enhancing fitness through activation of convergent angiogenic pathways in this vascular malignancy. IMPLICATIONS: This study shows that, while drive events of malignant vasoformative tumors of humans and dogs include diverse mutations and stochastic rearrangements that create novel fusion genes, convergent transcriptional programs govern the highly conserved morphologic organization and biological behavior of these tumors in both species.

BarkBase: Epigenomic Annotation of Canine Genomes.

Dogs are an unparalleled natural model for investigating the genetics of health and disease, particularly for complex diseases like cancer. Comprehensive genomic annotation of regulatory elements active in healthy canine tissues is crucial both for identifying candidate causal variants and for designing functional studies needed to translate genetic associations into disease insight. Currently, canine geneticists rely primarily on annotations of the human or mouse genome that have been remapped to dog, an approach that misses dog-specific features. Here, we describe BarkBase, a canine epigenomic resource available at barkbase.org. BarkBase hosts data for 27 adult tissue types, with biological replicates, and for one sample of up to five tissues sampled at each of four carefully staged embryonic time points. RNA sequencing is complemented with whole genome sequencing and with assay for transposase-accessible chromatin using sequencing (ATAC-seq), which identifies open chromatin regions. By including replicates, we can more confidently discern tissue-specific transcripts and assess differential gene expression between tissues and timepoints. By offering data in easy-to-use file formats, through a visual browser modeled on similar genomic resources for human, BarkBase introduces a powerful new resource to support comparative studies in dogs and humans.

Comparative Genomics Reveals Shared Mutational Landscape in Canine Hemangiosarcoma and Human Angiosarcoma.

Angiosarcoma is a highly aggressive cancer of blood vessel-forming cells with few effective treatment options and high patient mortality. It is both rare and heterogenous, making large, well-powered genomic studies nearly impossible. Dogs commonly suffer from a similar cancer, called hemangiosarcoma, with breeds like the golden retriever carrying heritable genetic factors that put them at high risk. If the clinical similarity of canine hemangiosarcoma and human angiosarcoma reflects shared genomic etiology, dogs could be a critically needed model for advancing angiosarcoma research. We assessed the genomic landscape of canine hemangiosarcoma via whole-exome sequencing (47 golden retriever hemangiosarcomas) and RNA sequencing (74 hemangiosarcomas from multiple breeds). Somatic coding mutations occurred most frequently in the tumor suppressor TP53 (59.6% of cases) as well as two genes in the PI3K pathway: the oncogene PIK3CA (29.8%) and its regulatory subunit PIK3R1 (8.5%). The predominant mutational signature was the age-associated deamination of cytosine to thymine. As reported in human angiosarcoma, CDKN2A/B was recurrently deleted and VEGFA, KDR, and KIT recurrently gained. We compared the canine data to human data recently released by The Angiosarcoma Project, and found many of the same genes and pathways significantly enriched for somatic mutations, particularly in breast and visceral angiosarcomas. Canine hemangiosarcoma closely models the genomic landscape of human angiosarcoma of the breast and viscera, and is a powerful tool for investigating the pathogenesis of this devastating disease. IMPLICATIONS: We characterize the genomic landscape of canine hemangiosarcoma and demonstrate its similarity to human angiosarcoma.

SETD2 Is Recurrently Mutated in Whole-Exome Sequenced Canine Osteosarcoma.

Osteosarcoma is a debilitating bone cancer that affects humans, especially children and adolescents. A homologous form of osteosarcoma spontaneously occurs in dogs, and its differential incidence observed across breeds allows for the investigation of tumor mutations in the context of multiple genetic backgrounds. Using whole-exome sequencing and dogs from three susceptible breeds (22 golden retrievers, 21 Rottweilers, and 23 greyhounds), we found that osteosarcoma tumors show a high frequency of somatic copy-number alterations (SCNA), affecting key oncogenes and tumor-suppressor genes. The across-breed results are similar to what has been observed for human osteosarcoma, but the disease frequency and somatic mutation counts vary in the three breeds. For all breeds, three mutational signatures (one of which has not been previously reported) and 11 significantly mutated genes were identified. TP53 was the most frequently altered gene (83% of dogs have either mutations or SCNA in TP53), recapitulating observations in human osteosarcoma. The second most frequently mutated gene, histone methyltransferase SETD2, has known roles in multiple cancers, but has not previously been strongly implicated in osteosarcoma. This study points to the likely importance of histone modifications in osteosarcoma and highlights the strong genetic similarities between human and dog osteosarcoma, suggesting that canine osteosarcoma may serve as an excellent model for developing treatment strategies in both species.Significance: Canine osteosarcoma genomics identify SETD2 as a possible oncogenic driver of osteosarcoma, and findings establish the canine model as a useful comparative model for the corresponding human disease. Cancer Res; 78(13); 3421-31. ©2018 AACR.

Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication.

The genetic changes underlying the initial steps of animal domestication are still poorly understood. We generated a high-quality reference genome for the rabbit and compared it to resequencing data from populations of wild and domestic rabbits. We identified more than 100 selective sweeps specific to domestic rabbits but only a relatively small number of fixed (or nearly fixed) single-nucleotide polymorphisms (SNPs) for derived alleles. SNPs with marked allele frequency differences between wild and domestic rabbits were enriched for conserved noncoding sites. Enrichment analyses suggest that genes affecting brain and neuronal development have often been targeted during domestication. We propose that because of a truly complex genetic background, tame behavior in rabbits and other domestic animals evolved by shifts in allele frequencies at many loci, rather than by critical changes at only a few domestication loci.

An improved canine genome and a comprehensive catalogue of coding genes and non-coding transcripts.

The domestic dog, Canis familiaris, is a well-established model system for mapping trait and disease loci. While the original draft sequence was of good quality, gaps were abundant particularly in promoter regions of the genome, negatively impacting the annotation and study of candidate genes. Here, we present an improved genome build, canFam3.1, which includes 85 MB of novel sequence and now covers 99.8% of the euchromatic portion of the genome. We also present multiple RNA-Sequencing data sets from 10 different canine tissues to catalog ∼175,000 expressed loci. While about 90% of the coding genes previously annotated by EnsEMBL have measurable expression in at least one sample, the number of transcript isoforms detected by our data expands the EnsEMBL annotations by a factor of four. Syntenic comparison with the human genome revealed an additional ∼3,000 loci that are characterized as protein coding in human and were also expressed in the dog, suggesting that those were previously not annotated in the EnsEMBL canine gene set. In addition to ∼20,700 high-confidence protein coding loci, we found ∼4,600 antisense transcripts overlapping exons of protein coding genes, ∼7,200 intergenic multi-exon transcripts without coding potential, likely candidates for long intergenic non-coding RNAs (lincRNAs) and ∼11,000 transcripts were reported by two different library construction methods but did not fit any of the above categories. Of the lincRNAs, about 6,000 have no annotated orthologs in human or mouse. Functional analysis of two novel transcripts with shRNA in a mouse kidney cell line altered cell morphology and motility. All in all, we provide a much-improved annotation of the canine genome and suggest regulatory functions for several of the novel non-coding transcripts.

The draft genome sequence of the ferret (Mustela putorius furo) facilitates study of human respiratory disease.

The domestic ferret (Mustela putorius furo) is an important animal model for multiple human respiratory diseases. It is considered the 'gold standard' for modeling human influenza virus infection and transmission. Here we describe the 2.41 Gb draft genome assembly of the domestic ferret, constituting 2.28 Gb of sequence plus gaps. We annotated 19,910 protein-coding genes on this assembly using RNA-seq data from 21 ferret tissues. We characterized the ferret host response to two influenza virus infections by RNA-seq analysis of 42 ferret samples from influenza time-course data and showed distinct signatures in ferret trachea and lung tissues specific to 1918 or 2009 human pandemic influenza virus infections. Using microarray data from 16 ferret samples reflecting cystic fibrosis disease progression, we showed that transcriptional changes in the CFTR-knockout ferret lung reflect pathways of early disease that cannot be readily studied in human infants with cystic fibrosis disease.