White Shark Genome Reveals Ancient Elasmobranch Adaptations Associated with Wound Healing and the Maintenance of Genome Stability

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans

Pangolin genomes and the evolution of mammalian scales and immunity

Pangolins, unique mammals with scales over most of their body, no teeth, poor vision, and an acute olfactory system, comprise the only placental order (Pholidota) without a whole-genome map. To investigate pangolin biology and evolution, we developed genome assemblies of the Malayan (Manis javanica) and Chinese (M. pentadactyla) pangolins. Strikingly, we found that interferon epsilon (IFNE), exclusively expressed in epithelial cells and important in skin and mucosal immunity, is pseudogenized in all African and Asian pangolin species that we examined, perhaps impacting resistance to infection. We propose that scale development was an innovation that provided protection against injuries or stress and reduced pangolin vulnerability to infection. Further evidence of specialized adaptations was evident from positively selected genes involving immunity-related pathways, inflammation, energy storage and metabolism, muscular and nervous systems, and scale/hair development. Olfactory receptor gene families are significantly expanded in pangolins, reflecting their well-developed olfaction system. This study provides insights into mammalian adaptation and functional diversification, new research tools and questions, and perhaps a new natural IFNE-deficient animal model for studying mammalian immunity.University of Malaya and Ministry of Education, Malaysia [UM.C/HIR/MOHE/08]; UMRG grant from the University of Malaya and Ministry of Education, Malaysia [RG541-13HTM]; Russian Ministry of Science [11.G34.31.0068]; NIH-NHGRI grant [5U54HG00307907]SCI(E)[email protected]

Varidnaviruses in the human gut: A major expansion of the order Vinavirales

International audienceBacteriophages play key roles in the dynamics of the human microbiome. By far the most abundant components of the human gut virome are tailed bacteriophages of the realm Duplodnaviria, in particular, crAss-like phages. However, apart from duplodnaviruses, the gut virome has not been dissected in detail. Here we report a comprehensive census of a minor component of the gut virome, the tailless bacteriophages of the realm Varidnaviria. Tailless phages are primarily represented in the gut by prophages, that are mostly integrated in genomes of Alphaproteobacteria and Verrucomicrobia and belong to the order Vinavirales, which currently consists of the families Corticoviridae and Autolykiviridae. Phylogenetic analysis of the major capsid proteins (MCP) suggests that at least three new families should be established within Vinavirales to accommodate the diversity of prophages from the human gut virome. Previously, only the MCP and packaging ATPase genes were reported as conserved core genes of Vinavirales. Here we report an extended core set of 12 proteins, including MCP, packaging ATPase, and previously undetected lysis enzymes, that are shared by most of these viruses. We further demonstrate that replication system components are frequently replaced in the genomes of Vinavirales, suggestive of selective pressure for escape from yet unknown host defenses or avoidance of incompatibility with coinfecting related viruses. The results of this analysis show that, in a sharp contrast to marine viromes, varidnaviruses are a minor component of the human gut virome. Moreover, they are primarily represented by prophages, as indicated by the analysis of the flanking genes, suggesting that there are few, if any, lytic varidnavirus infections in the gut at any given time. These findings complement the existing knowledge of the human gut virome by exploring a group of viruses that has been virtually overlooked in previous work

Phylogenetic and Expression Analysis of CENH3 and APOLLO Genes in Sexual and Apomictic Boechera Species

Apomictic plants (reproducing via asexual seeds), unlike sexual individuals, avoid meiosis and egg cell fertilization. Consequently, apomixis is very important for fixing maternal genotypes in the next plant generations. Despite the progress in the study of apomixis, molecular and genetic regulation of the latter remains poorly understood. So far APOLLO gene encoding aspartate glutamate aspartate aspartate histidine exonuclease is one of the very few described genes associated with apomixis in Boechera species. The centromere-specific histone H3 variant encoded by CENH3 gene is essential for cell division. Mutations in CENH3 disrupt chromosome segregation during mitosis and meiosis since the attachment of spindle microtubules to a mutated form of the CENH3 histone fails. This paper presents in silico characteristic of APOLLO and CENH3 genes, which may affect apomixis. Furthermore, we characterize the structure of CENH3 by bioinformatic tools, study expression levels of APOLLO and CENH3 transcripts by Real-Time Polymerase Chain Reaction RT-PCR in gynoecium/siliques of the natural diploid apomictic and sexual Boechera species at the stages of meiosis and before and after fertilization. While CENH3 was a single copy gene in all Boechera species, the APOLLO gene have several polymorphic alleles associated with sexual and apomictic reproduction in the Boechera genera. Expression of the APOLLO apo-allele during meiosis was upregulated in gynoecium of apomict B. divaricarpa downregulating after meiosis until the 4th day after pollination (DAP). On the 5th DAP, expression in apomictic siliques increased again. In sexual B. stricta gynoecium and siliques APOLLO apo-allele did not express. Expression of the APOLLO sex-allele during and after meiosis in gynoecium of sexual plants was several times higher than that in apomictic gynoecium. However, after pollination the sex-allele was downregulated in sexual siliques to the level of apomicts and increased sharply on the 5th DAP, while in apomictic siliques it almost did not express. At the meiotic stage, the expression level of CENH3 in the gynoecium of apomicts was two times lower than that of the sexual Boechera, decreasing in both species after meiosis and keep remaining very low in siliques of both species for several days after artificial pollination until the 4th DAP, when the expression level raised in sexual B. stricta siliques exceeding 5 times the level in apomictic B. divaricarpa siliques. We also discuss polymorphism and phylogeny of the APOLLO and CENH3 genes. The results obtained may indicate to a role of the CENH3 and APOLLO genes in the development of apomixis in species of the genus Boechera

Assembly of the Boechera retrofracta Genome and Evolutionary Analysis of Apomixis-Associated Genes.

Closely related to the model plant Arabidopsis thaliana, the genus Boechera is known to contain both sexual and apomictic species or accessions. Boechera retrofracta is a diploid sexually reproducing species and is thought to be an ancestral parent species of apomictic species. Here we report the de novo assembly of the B. retrofracta genome using short Illumina and Roche reads from 1 paired-end and 3 mate pair libraries. The distribution of 23-mers from the paired end library has indicated a low level of heterozygosity and the presence of detectable duplications and triplications. The genome size was estimated to be equal 227 Mb. N50 of the assembled scaffolds was 2.3 Mb. Using a hybrid approach that combines homology-based and de novo methods 27,048 protein-coding genes were predicted. Also repeats, transfer RNA (tRNA) and ribosomal RNA (rRNA) genes were annotated. Finally, genes of B. retrofracta and 6 other Brassicaceae species were used for phylogenetic tree reconstruction. In addition, we explored the histidine exonuclease APOLLO locus, related to apomixis in Boechera, and proposed model of its evolution through the series of duplications. An assembled genome of B. retrofracta will help in the challenging assembly of the highly heterozygous genomes of hybrid apomictic species

Assembly of the Boechera retrofracta Genome and Evolutionary Analysis of Apomixis-Associated Genes

Closely related to the model plant Arabidopsis thaliana, the genus Boechera is known to contain both sexual and apomictic species or accessions. Boechera retrofracta is a diploid sexually reproducing species and is thought to be an ancestral parent species of apomictic species. Here we report the de novo assembly of the B. retrofracta genome using short Illumina and Roche reads from 1 paired-end and 3 mate pair libraries. The distribution of 23-mers from the paired end library has indicated a low level of heterozygosity and the presence of detectable duplications and triplications. The genome size was estimated to be equal 227 Mb. N50 of the assembled scaffolds was 2.3 Mb. Using a hybrid approach that combines homology-based and de novo methods 27,048 protein-coding genes were predicted. Also repeats, transfer RNA (tRNA) and ribosomal RNA (rRNA) genes were annotated. Finally, genes of B. retrofracta and 6 other Brassicaceae species were used for phylogenetic tree reconstruction. In addition, we explored the histidine exonuclease APOLLO locus, related to apomixis in Boechera, and proposed model of its evolution through the series of duplications. An assembled genome of B. retrofracta will help in the challenging assembly of the highly heterozygous genomes of hybrid apomictic species

GFF annotation file from maker annotation of genome including all scaffolds

This is the maker generated annotation file containing all repeat, gene, and protein annotations for the maker assembly that is uploaded on NCBI associated with the genome sequence for the manuscript "White shark genome: ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability.

proteins_from_10kbp_scaffolds

These are the predicted protein sequences from maker for scaffolds greater than 10 kbp in lengt

blast2go_blastx_results_proteins

This file is a file containing the blast descriptions for all of the predicted proteins (regardless of scaffold size) for the genome