The tuatara genome reveals ancient features of amniote evolution.

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.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

PHI-base in 2022: a multi-species phenotype database for Pathogen-Host Interactions.

Since 2005, the Pathogen-Host Interactions Database (PHI-base) has manually curated experimentally verified pathogenicity, virulence and effector genes from fungal, bacterial and protist pathogens, which infect animal, plant, fish, insect and/or fungal hosts. PHI-base (www.phi-base.org) is devoted to the identification and presentation of phenotype information on pathogenicity and effector genes and their host interactions. Specific gene alterations that did not alter the in host interaction phenotype are also presented. PHI-base is invaluable for comparative analyses and for the discovery of candidate targets in medically and agronomically important species for intervention. Version 4.12 (September 2021) contains 4387 references, and provides information on 8411 genes from 279 pathogens, tested on 228 hosts in 18, 190 interactions. This provides a 24% increase in gene content since Version 4.8 (September 2019). Bacterial and fungal pathogens represent the majority of the interaction data, with a 54:46 split of entries, whilst protists, protozoa, nematodes and insects represent 3.6% of entries. Host species consist of approximately 54% plants and 46% others of medical, veterinary and/or environmental importance. PHI-base data is disseminated to UniProtKB, FungiDB and Ensembl Genomes. PHI-base will migrate to a new gene-centric version (version 5.0) in early 2022. This major development is briefly described.

PHI-base: the pathogen-host interactions database.

The pathogen-host interactions database (PHI-base) is available at www.phi-base.org. PHI-base contains expertly curated molecular and biological information on genes proven to affect the outcome of pathogen-host interactions reported in peer reviewed research articles. PHI-base also curates literature describing specific gene alterations that did not affect the disease interaction phenotype, in order to provide complete datasets for comparative purposes. Viruses are not included, due to their extensive coverage in other databases. In this article, we describe the increased data content of PHI-base, plus new database features and further integration with complementary databases. The release of PHI-base version 4.8 (September 2019) contains 3454 manually curated references, and provides information on 6780 genes from 268 pathogens, tested on 210 hosts in 13,801 interactions. Prokaryotic and eukaryotic pathogens are represented in almost equal numbers. Host species consist of approximately 60% plants (split 50:50 between cereal and non-cereal plants), and 40% other species of medical and/or environmental importance. The information available on pathogen effectors has risen by more than a third, and the entries for pathogens that infect crop species of global importance has dramatically increased in this release. We also briefly describe the future direction of the PHI-base project, and some existing problems with the PHI-base curation process.

A genome-wide investigation of adaptive signatures in protein-coding genes related to tool behaviour in New Caledonian and Hawaiian crows.

Very few animals habitually manufacture and use tools. It has been suggested that advanced tool behaviour co-evolves with a suite of behavioural, morphological and life history traits. In fact, there are indications for such an adaptive complex in tool-using crows (genus Corvus species). Here, we sequenced the genomes of two habitually tool-using and ten non-tool-using crow species to search for genomic signatures associated with a tool-using lifestyle. Using comparative genomic and population genetic approaches, we screened for signals of selection in protein-coding genes in the tool-using New Caledonian and Hawaiian crows. While we detected signals of recent selection in New Caledonian crows near genes associated with bill morphology, our data indicate that genetic changes in these two lineages are surprisingly subtle, with little evidence at present for convergence. We explore the biological explanations for these findings, such as the relative roles of gene regulation and protein-coding changes, as well as the possibility that statistical power to detect selection in recently diverged lineages may have been insufficient. Our study contributes to a growing body of literature aiming to decipher the genetic basis of recently evolved complex behaviour.

PomBase 2018: user-driven reimplementation of the fission yeast database provides rapid and intuitive access to diverse, interconnected information.

PomBase (www.pombase.org), the model organism database for the fission yeast Schizosaccharomyces pombe, has undergone a complete redevelopment, resulting in a more fully integrated, better-performing service. The new infrastructure supports daily data updates as well as fast, efficient querying and smoother navigation within and between pages. New pages for publications and genotypes provide routes to all data curated from a single source and to all phenotypes associated with a specific genotype, respectively. For ontology-based annotations, improved displays balance comprehensive data coverage with ease of use. The default view now uses ontology structure to provide a concise, non-redundant summary that can be expanded to reveal underlying details and metadata. The phenotype annotation display also offers filtering options to allow users to focus on specific areas of interest. An instance of the JBrowse genome browser has been integrated, facilitating loading of and intuitive access to, genome-scale datasets. Taken together, the new data and pages, along with improvements in annotation display and querying, allow users to probe connections among different types of data to form a comprehensive view of fission yeast biology. The new PomBase implementation also provides a rich set of modular, reusable tools that can be deployed to create new, or enhance existing, organism-specific databases.

Molecular structure of sauropsid ß-keratins from tuatara (Sphenodon punctatus).

The birds and reptiles, collectively known as the sauropsids, can be subdivided phylogenetically into the archosaurs (birds, crocodiles), the testudines (turtles), the squamates (lizards, snakes) and the rhynchocephalia (tuatara). The structural framework of the epidermal appendages from the sauropsids, which include feathers, claws and scales, has previously been characterised by electron microscopy, infrared spectroscopy and X-ray diffraction analyses, as well as by studies of the amino acid sequences of the constituent ß-keratin proteins (also referred to as the corneous ß-proteins). An important omission in this work, however, was the lack of sequence and structural data relating to the epidermal appendages of the rhynchocephalia (tuatara), one of the two branches of the lepidosaurs. Considerable effort has gone into sequencing the tuatara genome and while this is not yet complete, there are now sufficient sequence data for conclusions to be drawn on the similarity of the ß-keratins from the tuatara to those of other members of the sauropsids. These results, together with a comparison of the X-ray diffraction pattern of tuatara claw with those from seagull feather and goanna claw, confirm that there is a common structural plan in the ß-keratins of all of the sauropsids, and not just those that comprise the archosaurs (birds and crocodiles), the testudines (turtles) and the squamates (lizards and snakes).

Female Mimicry by Sneaker Males Has a Transcriptomic Signature in Both the Brain and the Gonad in a Sex-Changing Fish.

Phenotypic plasticity represents an elegant adaptive response of individuals to a change in their environment. Bluehead wrasses (Thalassoma bifasciatum) exhibit astonishing sexual plasticity, including female-to-male sex change and discrete male morphs that differ strikingly in behavior, morphology, and gonadal investment. Using RNA-seq transcriptome profiling, we examined the genes and physiological pathways underlying flexible behavioral and gonadal differences among female, dominant (bourgeois) male, and female-mimic (sneaker) male blueheads. For the first time in any organism, we find that female mimicry by sneaker males has a transcriptional signature in both the brain and the gonad. Sneaker males shared striking similarity in neural gene expression with females, supporting the idea that males with alternative reproductive phenotypes have "female-like brains." Sneaker males also overexpressed neuroplasticity genes, suggesting that their opportunistic reproductive strategy requires a heightened capacity for neuroplasticity. Bourgeois males overexpressed genes associated with socio-sexual behaviors (e.g., isotocin), but also neuroprotective genes and biomarkers of oxidative stress and aging, indicating a hitherto unexplored cost to these males of attaining the reproductively privileged position at the top of the social hierarchy. Our novel comparison of testicular transcriptomes in a fish with male sexual polymorphism associates greater gonadal investment by sneaker males with overexpression of genes involved in cell proliferation and sperm quality control. We propose that morphological female-mimicry by sneaker male teleosts entails pervasive downregulation of androgenesis genes, consistent with low androgen production in males lacking well-developed secondary sexual characters.

PHI-base: a new interface and further additions for the multi-species pathogen-host interactions database.

The pathogen-host interactions database (PHI-base) is available at www.phi-base.org PHI-base contains expertly curated molecular and biological information on genes proven to affect the outcome of pathogen-host interactions reported in peer reviewed research articles. In addition, literature that indicates specific gene alterations that did not affect the disease interaction phenotype are curated to provide complete datasets for comparative purposes. Viruses are not included. Here we describe a revised PHI-base Version 4 data platform with improved search, filtering and extended data display functions. A PHIB-BLAST search function is provided and a link to PHI-Canto, a tool for authors to directly curate their own published data into PHI-base. The new release of PHI-base Version 4.2 (October 2016) has an increased data content containing information from 2219 manually curated references. The data provide information on 4460 genes from 264 pathogens tested on 176 hosts in 8046 interactions. Prokaryotic and eukaryotic pathogens are represented in almost equal numbers. Host species belong ∼70% to plants and 30% to other species of medical and/or environmental importance. Additional data types included into PHI-base 4 are the direct targets of pathogen effector proteins in experimental and natural host organisms. The curation problems encountered and the future directions of the PHI-base project are briefly discussed.

Adipose transcriptome analysis provides novel insights into molecular regulation of prolonged fasting in northern elephant seal pups.

The physiological and cellular adaptations to extreme fasting in northern elephant seals ( Mirounga angustirostris, NES) are remarkable and may help to elucidate endocrine mechanisms that regulate lipid metabolism and energy homeostasis in mammals. Recent studies have highlighted the importance of thyroid hormones in the maintenance of a lipid-based metabolism during prolonged fasting in weaned NES pups. To identify additional molecular regulators of fasting, we used a transcriptomics approach to examine changes in global gene expression profiles before and after 6-8 wk of fasting in weaned NES pups. We produced a de novo assembly and identified 98 unique protein-coding genes that were differentially expressed between early and late fasting. Most of the downregulated genes were associated with lipid, carbohydrate, and protein metabolism. A number of downregulated genes were also associated with maintenance of the extracellular matrix, consistent with tissue remodeling during weight loss and the multifunctional nature of blubber tissue, which plays both metabolic and structural roles in marine mammals. Using this data set, we predict potential mechanisms by which NES pups sustain metabolism and regulate adipose stores throughout the fast, and provide a valuable resource for additional studies of extreme metabolic adaptations in mammals.

Identification of sex differences in zebrafish (Danio rerio) brains during early sexual differentiation and masculinization using 17α-methyltestoterone.

Sexual behavior in teleost fish is highly plastic. It can be attributed to the relatively few sex differences found in adult brain transcriptomes. Environmental and hormonal factors can influence sex-specific behavior. Androgen treatment stimulates behavioral masculinization. Sex dimorphic gene expression in developing teleost brains and the molecular basis for androgen-induced behavioral masculinization are poorly understood. In this study, juvenile zebrafish (Danio rerio) were treated with 100 ng/L of 17 alpha-methyltestosterone (MT) during sexual development from 20 days post fertilization to 40 days and 60 days post fertilization. We compared brain gene expression patterns in MT-treated zebrafish with control males and females using RNA-Seq to shed light on the dynamic changes in brain gene expression during sexual development and how androgens affect brain gene expression leading to behavior masculinization. We found modest differences in gene expression between juvenile male and female zebrafish brains. Brain aromatase (cyp19a1b), prostaglandin 3a synthase (ptges3a), and prostaglandin reductase 1 (ptgr1) were among the genes with sexually dimorphic expression patterns. MT treatment significantly altered gene expression relative to both male and female brains. Fewer differences were found among MT-treated brains and male brains compared to female brains, particularly at 60 dpf. MT treatment upregulated the expression of hydroxysteroid 11-beta dehydrogenase 2 (hsd11b2), deiodinase, iodothyronine, type II (dio2), and gonadotrophin releasing hormones (GnRH) 2 and 3 (gnrh2 and gnrh3) suggesting local synthesis of 11-ketotestosterone, triiodothyronine, and GnRHs in zebrafish brains which are influenced by androgens. Androgen, estrogen, prostaglandin, thyroid hormone, and GnRH signaling pathways likely interact to modulate teleost sexual behavior.

Histological and transcriptomic effects of 17α-methyltestosterone on zebrafish gonad development.

BACKGROUND: Sex hormones play important roles in teleost ovarian and testicular development. In zebrafish, ovarian differentiation appears to be dictated by an oocyte-derived signal via Cyp19a1a aromatase-mediated estrogen production. Androgens and aromatase inhibitors can induce female-to-male sex reversal, however, the mechanisms underlying gonadal masculinisation are poorly understood. We used histological analyses together with RNA sequencing to characterise zebrafish gonadal transcriptomes and investigate the effects of 17α-methyltestosterone on gonadal differentiation. RESULTS: At a morphological level, 17α-methyltestosterone (MT) masculinised gonads and accelerated spermatogenesis, and these changes were paralleled in masculinisation and de-feminisation of gonadal transcriptomes. MT treatment upregulated expression of genes involved in male sex determination and differentiation (amh, dmrt1, gsdf and wt1a) and those involved in 11-oxygenated androgen production (cyp11c1 and hsd11b2). It also repressed expression of ovarian development and folliculogenesis genes (bmp15, gdf9, figla, zp2.1 and zp3b). Furthermore, MT treatment altered epigenetic modification of histones in zebrafish gonads. Contrary to expectations, higher levels of cyp19a1a or foxl2 expression in control ovaries compared to MT-treated testes and control testes were not statistically significant during early gonad development (40 dpf). CONCLUSION: Our study suggests that both androgen production and aromatase inhibition are important for androgen-induced gonadal masculinisation and natural testicular differentiation in zebrafish.

Analysis of the genome of the New Zealand giant collembolan (Holacanthella duospinosa) sheds light on hexapod evolution.

BACKGROUND: The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events. RESULTS: Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes. CONCLUSIONS: The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda.

Male-female relatedness at specific SNP-linkage groups influences cryptic female choice in Chinook salmon (Oncorhynchus tshawytscha).

In a range of taxa, the relatedness between mates influences both pre- and post-mating processes of sexual selection. However, relatively little is known about the genetic loci facilitating such a bias, with the exception of the major histocompatibility complex. Here, we performed tightly controlled replicated in vitro fertilization trials to explore the impact of relatedness on two possible mechanisms of cryptic female choice (CFC) in Chinook salmon (Oncorhynchus tshawytscha). We tested (i) whether relatedness of mates, assessed using 682 single nucleotide polymorphisms (SNPs) on 29 SNP-linkage groups (LGs), biases a male's sperm velocity in ovarian fluid (a parameter previously shown to predict male fertilization success), and (ii) whether relatedness of mates governs fertilization success via other mechanisms, probably via sperm-egg interactions. We found that relatedness on three LGs explained the variation in sperm velocity, and relatedness on two LGs explained fertilization success, which might indicate the presence of genes important in sperm-ovarian fluid and sperm-egg interactions in these genomic regions. Mapping of the SNPs on these LGs to the rainbow trout genome revealed two genes that affect fertility in humans and represent candidate genes for further studies. Our results thereby provide a novel contribution to the understanding of the mechanism of CFC.

PomBase 2015: updates to the fission yeast database.

PomBase (http://www.pombase.org) is the model organism database for the fission yeast Schizosaccharomyces pombe. PomBase provides a central hub for the fission yeast community, supporting both exploratory and hypothesis-driven research. It provides users easy access to data ranging from the sequence level, to molecular and phenotypic annotations, through to the display of genome-wide high-throughput studies. Recent improvements to the site extend annotation specificity, improve usability and allow for monthly data updates. Both in-house curators and community researchers provide manually curated data to PomBase. The genome browser provides access to published high-throughput data sets and the genomes of three additional Schizosaccharomyces species (Schizosaccharomyces cryophilus, Schizosaccharomyces japonicus and Schizosaccharomyces octosporus).

Uncovering the pathways underlying whole body regeneration in a chordate model, Botrylloides leachi using de novo transcriptome analysis.

BACKGROUND: Regenerative capacity differs greatly between animals. In vertebrates regenerative abilities are highly limited and tissue or organ specific. However the closest related chordate to the vertebrate clade, Botrylloides leachi, can undergo whole body regeneration (WBR). Therefore, research on WBR in B. leachi has focused on pathways known to be important for regeneration in vertebrates. To obtain a comprehensive vision of this unique process we have carried out the first de novo transcriptome sequencing for multiple stages of WBR occurring in B. leachi. The identified changes in gene expression during B. leachi WBR offer novel insights into this remarkable ability to regenerate. RESULTS: The transcriptome of B. leachi tissue undergoing WBR were analysed using differential gene expression, gene ontology and pathway analyses. We observed up-regulation in the expression of genes involved in wound healing and known developmental pathways including WNT, TGF-ß and Notch, during the earliest stages of WBR. Later in WBR, the expression patterns in several pathways required for protein synthesis, biogenesis and the organisation of cellular components were up-regulated. CONCLUSIONS: While the genes expressed early on are characteristic of a necessary wound healing response to an otherwise lethal injury, the subsequent vast increase in protein synthesis conceivably sustains the reestablishment of the tissue complexity and body axis polarity within the regenerating zooid. We have, for the first time, provided a global overview of the genes and their corresponding pathways that are modulated during WBR in B. leachi.

Canto: an online tool for community literature curation.

MOTIVATION: Detailed curation of published molecular data is essential for any model organism database. Community curation enables researchers to contribute data from their papers directly to databases, supplementing the activity of professional curators and improving coverage of a growing body of literature. We have developed Canto, a web-based tool that provides an intuitive curation interface for both curators and researchers, to support community curation in the fission yeast database, PomBase. Canto supports curation using OBO ontologies, and can be easily configured for use with any species. AVAILABILITY: Canto code and documentation are available under an Open Source license from http://curation.pombase.org/. Canto is a component of the Generic Model Organism Database (GMOD) project (http://www.gmod.org/).

Molecular evolution of Dmrt1 accompanies change of sex-determining mechanisms in reptilia.

In reptiles, sex-determining mechanisms have evolved repeatedly and reversibly between genotypic and temperature-dependent sex determination. The gene Dmrt1 directs male determination in chicken (and presumably other birds), and regulates sex differentiation in animals as distantly related as fruit flies, nematodes and humans. Here, we show a consistent molecular difference in Dmrt1 between reptiles with genotypic and temperature-dependent sex determination. Among 34 non-avian reptiles, a convergently evolved pair of amino acids encoded by sequence within exon 2 near the DM-binding domain of Dmrt1 distinguishes species with either type of sex determination. We suggest that this amino acid shift accompanied the evolution of genotypic sex determination from an ancestral condition of temperature-dependent sex determination at least three times among reptiles, as evident in turtles, birds and squamates. This novel hypothesis describes the evolution of sex-determining mechanisms as turnover events accompanied by one or two small mutations.

PomBase: a comprehensive online resource for fission yeast.

PomBase (www.pombase.org) is a new model organism database established to provide access to comprehensive, accurate, and up-to-date molecular data and biological information for the fission yeast Schizosaccharomyces pombe to effectively support both exploratory and hypothesis-driven research. PomBase encompasses annotation of genomic sequence and features, comprehensive manual literature curation and genome-wide data sets, and supports sophisticated user-defined queries. The implementation of PomBase integrates a Chado relational database that houses manually curated data with Ensembl software that supports sequence-based annotation and web access. PomBase will provide user-friendly tools to promote curation by experts within the fission yeast community. This will make a key contribution to shaping its content and ensuring its comprehensiveness and long-term relevance.

modMine: flexible access to modENCODE data.

In an effort to comprehensively characterize the functional elements within the genomes of the important model organisms Drosophila melanogaster and Caenorhabditis elegans, the NHGRI model organism Encyclopaedia of DNA Elements (modENCODE) consortium has generated an enormous library of genomic data along with detailed, structured information on all aspects of the experiments. The modMine database (http://intermine.modencode.org) described here has been built by the modENCODE Data Coordination Center to allow the broader research community to (i) search for and download data sets of interest among the thousands generated by modENCODE; (ii) access the data in an integrated form together with non-modENCODE data sets; and (iii) facilitate fine-grained analysis of the above data. The sophisticated search features are possible because of the collection of extensive experimental metadata by the consortium. Interfaces are provided to allow both biologists and bioinformaticians to exploit these rich modENCODE data sets now available via modMine.