Search

Cancer is a disease that affects nearly all multicellular life, including birds. However, little is known about what factors explain the variance in cancer prevalence among species. Litter size is positively correlated with cancer prevalence in managed species of mammals, and larger body size, but not incubation or nestling period, is linked to tumor prevalence in wild birds. Also, birds that produce more elaborate sexual traits are expected to have fewer resources for cancer defenses and thus higher cancer prevalence. In this study, we examined whether cancer prevalence is associated with a wide variety of life history traits (clutch size, incubation length, body mass, lifespan, and the extent of sexual dimorphism) across 108 species of managed birds in 25 different zoological facilities, sanctuaries, and veterinary clinics. We found that clutch size was positively correlated with cancer and neoplasia (both benign and malignant) prevalence, even after controlling for body mass. Cancer prevalence was not associated with incubation length, body mass, lifespan, or sexual dimorphism. The positive correlations of clutch size with cancer prevalence and neoplasia prevalence suggest that there may be life-history trade-offs between reproductive investment and somatic maintenance (in the form of cancer prevention mechanisms) in managed birds.

Chimerism is a widespread phenomenon across the tree of life. It is defined as a multicellular organism composed of cells from other genetically distinct entities. This ability to ‘tolerate’ non-self cells may be linked to susceptibility to diseases like cancer. Here we test whether chimerism is associated with cancers across obligately multicellular organisms in the tree of life. We classified 12 obligately multicellular taxa from lowest to highest chimerism levels based on the existing literature on the presence of chimerism in these species. We then tested for associations of chimerism with tumour invasiveness, neoplasia (benign or malignant) prevalence and malignancy prevalence in 11 terrestrial mammalian species. We found that taxa with higher levels of chimerism have higher tumour invasiveness, though there was no association between malignancy or neoplasia and chimerism among mammals. This suggests that there may be an important biological relationship between chimerism and susceptibility to tissue invasion by cancerous cells. Studying chimerism might help us identify mechanisms underlying invasive cancers and also could provide insights into the detection and management of emerging transmissible cancers.

The ecology in which species live and evolve likely affects their health and vulnerability to diseases including cancer. Using 14,267 necropsy records across 244 vertebrate species, we tested if animals in low productivity habitats, with large habitat range, high body temperature and weight-inferred estimates of metabolic rates, and in high trophic levels (from lowest to highest: herbivores, invertivores, primary carnivores, and secondary carnivores) are linked with having increased prevalence of neoplasia. This study found that: (1) habitat productivity negatively correlated with the prevalence of malignancy and neoplasia across tissues, and malignancy and neoplasia in gastrointestinal tissues; (2) inferred metabolic rates negatively correlated with the prevalence of neoplasia; and (3) trophic levels positively correlated with malignancy and neoplasia prevalence in both mammals and non-mammals. However, only the correlations with trophic levels remained significant after Bonferroni corrections for multiple testing. There are several mechanisms that might explain these findings, including the biomagnification of carcinogens in higher trophic levels, as well as tradeoffs between cancer suppression versus reproduction and survival in low productivity environments.

Disease susceptibility and resistance are important factors for the conservation of endangered species, including elephants. We analyzed pathology data from 26 zoos and report that Asian elephants have increased neoplasia and malignancy prevalence compared with African bush elephants. This is consistent with observed higher susceptibility to tuberculosis and elephant endotheliotropic herpesvirus (EEHV) in Asian elephants. To investigate genetic mechanisms underlying disease resistance, including differential responses between species, among other elephant traits, we sequenced multiple elephant genomes. We report a draft assembly for an Asian elephant, and defined 862 and 1,017 conserved potential regulatory elements in Asian and African bush elephants, respectively. In the genomes of both elephant species, conserved elements were significantly enriched with genes differentially expressed between the species. In Asian elephants, these putative regulatory regions were involved in immunity pathways including tumor-necrosis factor, which plays an important role in EEHV response. Genomic sequences of African bush, forest, and Asian elephant genomes revealed extensive sequence conservation at TP53 retrogene loci across three species, which may be related to TP53 functionality in elephant cancer resistance. Positive selection scans revealed outlier genes related to additional elephant traits. Our study suggests that gene regulation plays an important role in the differential inflammatory response of Asian and African elephants, leading to increased infectious disease and cancer susceptibility in Asian elephants. These genomic discoveries can inform future functional and translational studies aimed at identifying effective treatment approaches for ill elephants, which may improve conservation.

Most tissue collections of neoplasms are composed of formalin-fixed and paraffin-embedded (FFPE) excised tumor samples used for routine diagnostics. DNA sequencing is becoming increasingly important in cancer research and clinical management; however, it is difficult to accurately sequence DNA from FFPE samples. We developed and validated a new bioinformatic algorithm to robustly identify somatic single nucleotide variants (SNVs) from whole exome sequencing using small amounts of DNA extracted from archival FFPE samples of breast cancers. We optimized this strategy using 28 pairs of technical replicates. After optimization, the mean similarity between replicates increased 5-fold, reaching 88% (range 0-100%), with a mean of 21.4 SNVs (range 1-68) per sample, representing a markedly superior performance to existing algorithms. We found that the SNV-identification accuracy declined when there was less than 40ng of DNA available and that insertion-deletion variant calls are less reliable than single base substitutions. As the first application of the new algorithm, we compared samples of ductal carcinoma in situ (DCIS) of the breast to their adjacent invasive ductal carcinoma (IDC) samples. We observed an increased number of mutations (paired-samples sign test, pKey PointsThe sequencing of reduced quantities of DNA extracted from FFPE samples leads to substantial sequencing errors that require correction in order to obtain accurate detection of somatic mutations.We developed and validated a new bioinformatic algorithm to robustly identify somatic single nucleotide variants using small amounts of DNA extracted from archival FFPE samples of breast cancers.Variant calling software packages need to be optimized to reduce the impact of sequencing errors. Our bioinformatics pipeline represents a significant methodological advance compared to the currently available bioinformatic tools used for the analysis of small FFPE samples.

Disease susceptibility and resistance comprise important factors in conservation, particularly in elephants. To determine genetic mechanisms underlying disease resistance and other unique elephant traits, we estimated 862 and 1,017 potential regulatory elements in Asian and African elephants, respectively. These elements are significantly enriched in both species with differentially expressed genes involved in immunity pathways, including tumor-necrosis factor which plays a role in the response to elephant endotheliotropic herpesvirus (EEHV). Population genomics analyses indicate that amplifiedTP53retrogenes are maintained by purifying selection and may contribute to cancer resistance in elephants, including less malignancies in African vs. Asian elephants. Positive selection scans across elephant genomes revealed genes that may control iconic elephant traits such as tusk development, memory, and somatic maintenance. Our study supports the hypothesis that interspecies variation in gene regulation contributes to differential inflammatory responses leading to increased infectious disease and cancer susceptibility in Asian versus African elephants. Genomics can inform functional immunological studies which may improve both conservation for elephants and human therapies.

Ngatiwai Trust Board. et al., The tuatara (Sphenodon punctatus)—the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2—is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which—at approximately 5 Gb—is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing., N.J.G. acknowledges the support of Ngatiwai iwi, Allan Wilson Centre, University of Otago, New Zealand Department of Conservation, New Zealand Genomics and Illumina. J.I.A. was supported by CONICYT National Doctoral Scholarship No. 21130515. M.W. was supported by NIH grant R35 GM124827. Ensembl annotation was supported by the Wellcome Trust (WT108749/Z/15/Z) and the European Molecular Biology Laboratory.

The tuatara (Sphenodon punctatus), the only living member of the archaic reptilian order Rhynchocephalia (Sphenodontia) once widespread across Gondwana, is an iconic and enigmatic terrestrial vertebrate endemic to New Zealand. A key link to the now extinct stem reptiles from which dinosaurs, modern reptiles, birds and mammals evolved, the tuatara provides exclusive insights into the ancestral amniotes. The tuatara genome, at ∼5 Gbp, is among the largest vertebrate genomes assembled. Analysis of this genome and comparisons to other vertebrates reinforces the uniqueness of the tuatara. Phylogenetic analyses indicate tuatara diverged from the snakes and lizards ∼250 MYA. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Genome sequence analysis identifies expansions of protein, non-protein-coding RNA families, and repeat elements, the latter of which show an extraordinary amalgam of reptilian and mammalian features. Sequencing of this genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. It also provides important insights into both the technical challenges and the cultural obligations associated with genome sequencing.

Squamates include all lizards and snakes, and display some of the most diverse and extreme morphological adaptations among vertebrates. However, compared with birds and mammals, relatively few resources exist for comparative genomic analyses of squamates, hampering efforts to understand the molecular bases of phenotypic diversification in such a speciose clade. In particular, the ∼400 species of anole lizard represent an extensive squamate radiation. Here, we sequence and assemble the draft genomes of three anole species—Anolis frenatus, Anolis auratus, and Anolis apletophallus—for comparison with the available reference genome of Anolis carolinensis. Comparative analyses reveal a rapid background rate of molecular evolution consistent with a model of punctuated equilibrium, and strong purifying selection on functional genomic elements in anoles. We find evidence for accelerated evolution in genes involved in behavior, sensory perception, and reproduction, as well as in genes regulating limb bud development and hindlimb specification. Morphometric analyses of anole fore and hindlimbs corroborated these findings. We detect signatures of positive selection across several genes related to the development and regulation of the forebrain, hormones, and the iguanian lizard dewlap, suggesting molecular changes underlying behavioral adaptations known to reinforce species boundaries were a key component in the diversification of anole lizards.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana

Gamete size dimorphism between sexes (anisogamy) is predicted to have evolved from an isogamous system in which sexes have equal-sized, monomorphic gametes. Although adaptive explanations for the evolution of anisogamy abound, we lack comparable insights into molecular changes that bring about the transition from monomorphism to dimorphism. The basal fungal cladeAllomycesprovides unique opportunities to investigate genomic changes that are associated with this transition in closely related species that show either isogamous or anisogamous mating systems. The anisogamous species show sexual dimorphism in gamete size, number, pigmentation and motility. We sequenced transcriptomes of fiveAllomycesisolates representing the two mating systems, including both male and female phenotypes in the anisogamous species. Maximum likelihood ancestral character state reconstruction performed in MESQUITE using the de-novo assembled transcriptomes indicated that anisogamy likely evolved once inAllomyces, and is a derived character as predicted in theory. We found that sexual stages ofAllomycesexpress homologs of several genes known to be involved in sex determination in model organisms includingDrosophilaand humans. Furthermore, expression ofCatSperhomologs in male- and female-biased samples in our analysis support the hypothesis that gamete interaction in the anisogamous species ofAllomycesmay involve similar molecular events as the egg-sperm interaction in animals, including humans. Although the strains representing either mating system shared much of the transcriptome, supporting recent common ancestry, the analysis of rate of evolution using individual gene trees indicates high substitution rates and divergence between the strains. In summary, we find that anisogamy likely evolved once inAllomyces, using convergent mechanisms to those in other taxa.

In species with highly heteromorphic sex chromosomes, the degradation of one of the sex chromosomes will result in unequal gene expression between the sexes (e.g. between XX females and XY males) and between the sex chromosomes and the autosomes. Dosage compensation is a process whereby genes on the sex chromosomes achieve equal gene expression. We compared genome-wide levels of transcription between males and females, and between the X chromosome and the autosomes in the green anole, Anolis carolinensis. We present evidence for dosage compensation between the sexes, and between the sex chromosomes and the autosomes. When dividing the X chromosome into regions based on linkage groups, we discovered that genes in the first reported X-linked region, anole linkage group b (LGb), exhibit complete dosage compensation, although the rest of the X-linked genes exhibit incomplete dosage compensation. Our data further suggest that the mechanism of this dosage compensation is upregulation of the X chromosome in males. We report that approximately 10% of coding genes, most of which are on the autosomes, are differentially expressed between males and females. In addition, genes on the X chromosome exhibited higher ratios of nonsynonymous to synonymous substitution than autosomal genes, consistent with the fast-X effect. Our results from the green anole add an additional observation of dosage compensation in a species with XX/XY sex determination.