The Global Invertebrate Genomics Alliance (GIGA): developing community resources to study diverse invertebrate genomes

Brief History of Placozoa

Placozoans are morphologically the simplest free-living animals. They represent a unique window of opportunities to understand both the origin of the animal organization and the rules of life for the system and synthetic biology of the future. However, despite more than 100 years of their investigations, we know little about their organization, natural habitats, and life strategies. Here, we introduce this unique animal phylum and highlight some directions vital to broadening the frontiers of the biomedical sciences. In particular, understanding the genomic bases of placozoan biodiversity, cell identity, connectivity, reproduction, and cellular bases of behavior are critical hot spots for future studies.

Creating, curating and evaluating a mitogenomic reference database to improve regional species identification using environmental DNA

Species detection using eDNA is revolutionizing global capacity to monitor biodiversity. However, the lack of regional, vouchered, genomic sequence information-especially sequence information that includes intraspecific variation-creates a bottleneck for management agencies wanting to harness the complete power of eDNA to monitor taxa and implement eDNA analyses. eDNA studies depend upon regional databases of mitogenomic sequence information to evaluate the effectiveness of such data to detect and identify taxa. We created the Oregon Biodiversity Genome Project to create a database of complete, nearly error-free mitogenomic sequences for all of Oregon's fishes. We have successfully assembled the complete mitogenomes of 313 specimens of freshwater, anadromous and estuarine fishes representing 24 families, 55 genera and 129 species and lineages. Comparative analyses of these sequences illustrate that many regions of the mitogenome are taxonomically informative, that the short (~150 bp) mitochondrial 'barcode' regions typically used for eDNA assays do not consistently diagnose for species and that complete single or multiple genes of the mitogenome are preferable for identifying Oregon's fishes. This project provides a blueprint for other researchers to follow as they build regional databases, illustrates the taxonomic value and limits of complete mitogenomic sequences and offers clues as to how current eDNA assays and environmental genomics methods of the future can best leverage this information.

Why do hybrids turn down sex?

Asexual reproduction is ancestral in prokaryotes; the switch to sexuality in eukaryotes is one of the major transitions in the history of life. The study of the maintenance of sex in eukaryotes has raised considerable interest for decades and is still one of evolutionary biology's most prominent question. The observation that many asexual species are of hybrid origin has led some to propose that asexuality in hybrids results from sexual processes being disturbed because of incompatibilities between the two parental species' genomes. However, in some cases, failure to produce asexual F1s in the lab may indicate that this mechanism is not the only road to asexuality in hybrid species. Here, we present a mathematical model and propose an alternative, adaptive route for the evolution of asexuality from previously sexual hybrids. Under some reproductive alterations, we show that asexuality can evolve to rescue hybrids' reproduction. Importantly, we highlight that when incompatibilities only affect the fusion of sperm and egg's genomes, the two traits that characterize asexuality, namely unreduced meiosis and the initiation of embryogenesis without the incorporation of the sperm's pronucleus, can evolve separately, greatly facilitating the overall evolutionary route. Taken together, our results provide an alternative, potentially complementary explanation for the link between asexuality and hybridization.

A review of the endangered mollusks transcriptome under the threatened species initiative of Korea

Transcriptome studies for conservation of endangered mollusks is a proactive approach towards managing threats and uncertainties facing these species in natural environments. The population of these species is declining due to habitat destruction, illicit wildlife trade, and global climate change. These activities risk the free movement of species across the wild landscape, loss of breeding grounds, and restrictions in displaying the physiological attributes so crucial for faunal welfare. Gastropods face the most negative ecological effects and have been enlisted under Korea's protective species consortium based on their population dynamics in the last few years. Moreover, with the genetic resources restricted for such species, conservation by informed planning is not possible. This review provides insights into the activities under the threatened species initiative of Korea with special reference to the transcriptome assemblies of endangered mollusks. The gastropods such as Ellobium chinense, Aegista chejuensis, Aegista quelpartensis, Incilaria fruhstorferi, Koreanohadra kurodana, Satsuma myomphala, and Clithon retropictus have been represented. Moreover, the transcriptome summary of bivalve Cristaria plicata and Caenogastropoda Charonia lampas sauliae is also discussed. Sequencing, de novo assembly, and annotation identified transcripts or homologs for the species and, based on an understanding of the biochemical and molecular pathways, were ascribed to predictive gene function. Mining for simple sequence repeats from the transcriptome have successfully assisted genetic polymorphism studies. A comparison of the transcriptome scheme of Korean endangered mollusks with the genomic resources of other endangered mollusks have been discussed with homologies and analogies for dictating future research.

Chromosome-level genome assembly for the Aldabra giant tortoise enables insights into the genetic health of a threatened population

BACKGROUND

The Aldabra giant tortoise (Aldabrachelys gigantea) is one of only two giant tortoise species left in the world. The species is endemic to Aldabra Atoll in Seychelles and is listed as Vulnerable on the International Union for Conservation of Nature Red List (v2.3) due to its limited distribution and threats posed by climate change. Genomic resources for A. gigantea are lacking, hampering conservation efforts for both wild and ex situpopulations. A high-quality genome would also open avenues to investigate the genetic basis of the species' exceptionally long life span.

FINDINGS

We produced the first chromosome-level de novo genome assembly of A. gigantea using PacBio High-Fidelity sequencing and high-throughput chromosome conformation capture. We produced a 2.37-Gbp assembly with a scaffold N50 of 148.6 Mbp and a resolution into 26 chromosomes. RNA sequencing-assisted gene model prediction identified 23,953 protein-coding genes and 1.1 Gbp of repetitive sequences. Synteny analyses among turtle genomes revealed high levels of chromosomal collinearity even among distantly related taxa. To assess the utility of the high-quality assembly for species conservation, we performed a low-coverage resequencing of 30 individuals from wild populations and two zoo individuals. Our genome-wide population structure analyses detected genetic population structure in the wild and identified the most likely origin of the zoo-housed individuals. We further identified putatively deleterious mutations to be monitored.

CONCLUSIONS

We establish a high-quality chromosome-level reference genome for A. gigantea and one of the most complete turtle genomes available. We show that low-coverage whole-genome resequencing, for which alignment to the reference genome is a necessity, is a powerful tool to assess the population structure of the wild population and reveal the geographic origins of ex situ individuals relevant for genetic diversity management and rewilding efforts.

An Outsider's Perspective on Why We Climb Mountains and Why Projects Like the i5k Matter

Initiatives like the i5k are creating evermore genome assemblies. These initiatives are resource heavy, and their justifications and economics deserve attention. Scientifically, these initiatives are important, paving the way for cross-species analysis, requiring the building of new computational analysis and tools, and creating other new resources. However, an open question remains of how we quantitively measure the impact of genomes, and by extension these initiatives. This forum article discusses one such method which is to look at the publications about a species over time, however, this method does not show any signal from a published genome, leaving an open question of how to measure impact.

The state of Medusozoa genomics: current evidence and future challenges

Medusozoa is a widely distributed ancient lineage that harbors one-third of Cnidaria diversity divided into 4 classes. This clade is characterized by the succession of stages and modes of reproduction during metagenic lifecycles, and includes some of the most plastic body plans and life cycles among animals. The characterization of traditional genomic features, such as chromosome numbers and genome sizes, was rather overlooked in Medusozoa and many evolutionary questions still remain unanswered. Modern genomic DNA sequencing in this group started in 2010 with the publication of the Hydra vulgaris genome and has experienced an exponential increase in the past 3 years. Therefore, an update of the state of Medusozoa genomics is warranted. We reviewed different sources of evidence, including cytogenetic records and high-throughput sequencing projects. We focused on 4 main topics that would be relevant for the broad Cnidaria research community: (i) taxonomic coverage of genomic information; (ii) continuity, quality, and completeness of high-throughput sequencing datasets; (iii) overview of the Medusozoa specific research questions approached with genomics; and (iv) the accessibility of data and metadata. We highlight a lack of standardization in genomic projects and their reports, and reinforce a series of recommendations to enhance future collaborative research.

Chromosome-level genome assembly, annotation, and phylogenomics of the gooseneck barnacle Pollicipes pollicipes

BACKGROUND

The barnacles are a group of >2,000 species that have fascinated biologists, including Darwin, for centuries. Their lifestyles are extremely diverse, from free-swimming larvae to sessile adults, and even root-like endoparasites. Barnacles also cause hundreds of millions of dollars of losses annually due to biofouling. However, genomic resources for crustaceans, and barnacles in particular, are lacking.

RESULTS

Using 62× Pacific Biosciences coverage, 189× Illumina whole-genome sequencing coverage, 203× HiC coverage, and 69× CHi-C coverage, we produced a chromosome-level genome assembly of the gooseneck barnacle Pollicipes pollicipes. The P. pollicipes genome is 770 Mb long and its assembly is one of the most contiguous and complete crustacean genomes available, with a scaffold N50 of 47 Mb and 90.5% of the BUSCO Arthropoda gene set. Using the genome annotation produced here along with transcriptomes of 13 other barnacle species, we completed phylogenomic analyses on a nearly 2 million amino acid alignment. Contrary to previous studies, our phylogenies suggest that the Pollicipedomorpha is monophyletic and sister to the Balanomorpha, which alters our understanding of barnacle larval evolution and suggests homoplasy in a number of naupliar characters. We also compared transcriptomes of P. pollicipes nauplius larvae and adults and found that nearly one-half of the genes in the genome are differentially expressed, highlighting the vastly different transcriptomes of larvae and adult gooseneck barnacles. Annotation of the genes with KEGG and GO terms reveals that these stages exhibit many differences including cuticle binding, chitin binding, microtubule motor activity, and membrane adhesion.

CONCLUSION

This study provides high-quality genomic resources for a key group of crustaceans. This is especially valuable given the roles P. pollicipes plays in European fisheries, as a sentinel species for coastal ecosystems, and as a model for studying barnacle adhesion as well as its key position in the barnacle tree of life. A combination of genomic, phylogenetic, and transcriptomic analyses here provides valuable insights into the evolution and development of barnacles.

Threatened Species Initiative: Empowering conservation action using genomic resources

Globally, 15,521 animal species are listed as threatened by the International Union for the Conservation of Nature, and of these less than 3% have genomic resources that can inform conservation management. To combat this, global genome initiatives are developing genomic resources, yet production of a reference genome alone does not conserve a species. The reference genome allows us to develop a suite of tools to understand both genome-wide and functional diversity within and between species. Conservation practitioners can use these tools to inform their decision-making. But, at present there is an implementation gap between the release of genome information and the use of genomic data in applied conservation by conservation practitioners. In May 2020, we launched the Threatened Species Initiative and brought a consortium of genome biologists, population biologists, bioinformaticians, population geneticists, and ecologists together with conservation agencies across Australia, including government, zoos, and nongovernment organizations. Our objective is to create a foundation of genomic data to advance our understanding of key Australian threatened species, and ultimately empower conservation practitioners to access and apply genomic data to their decision-making processes through a web-based portal. Currently, we are developing genomic resources for 61 threatened species from a range of taxa, across Australia, with more than 130 collaborators from government, academia, and conservation organizations. Developed in direct consultation with government threatened-species managers and other conservation practitioners, herein we present our framework for meeting their needs and our systematic approach to integrating genomics into threatened species recovery.

Predictability and parallelism in the contemporary evolution of hybrid genomes

Hybridization between species is widespread across the tree of life. As a result, many species, including our own, harbor regions of their genome derived from hybridization. Despite the recognition that this process is widespread, we understand little about how the genome stabilizes following hybridization, and whether the mechanisms driving this stabilization tend to be shared across species. Here, we dissect the drivers of variation in local ancestry across the genome in replicated hybridization events between two species pairs of swordtail fish: Xiphophorus birchmanni × X. cortezi and X. birchmanni × X. malinche. We find unexpectedly high levels of repeatability in local ancestry across the two types of hybrid populations. This repeatability is attributable in part to the fact that the recombination landscape and locations of functionally important elements play a major role in driving variation in local ancestry in both types of hybrid populations. Beyond these broad scale patterns, we identify dozens of regions of the genome where minor parent ancestry is unusually low or high across species pairs. Analysis of these regions points to shared sites under selection across species pairs, and in some cases, shared mechanisms of selection. We show that one such region is a previously unknown hybrid incompatibility that is shared across X. birchmanni × X. cortezi and X. birchmanni × X. malinche hybrid populations.

Evolution of glutamatergic signaling and synapses

Glutamate (Glu) is the primary excitatory transmitter in the mammalian brain. But, we know little about the evolutionary history of this adaptation, including the selection of l-glutamate as a signaling molecule in the first place. Here, we used comparative metabolomics and genomic data to reconstruct the genealogy of glutamatergic signaling. The origin of Glu-mediated communications might be traced to primordial nitrogen and carbon metabolic pathways. The versatile chemistry of L-Glu placed this molecule at the crossroad of cellular biochemistry as one of the most abundant metabolites. From there, innovations multiplied. Many stress factors or injuries could increase extracellular glutamate concentration, which led to the development of modular molecular systems for its rapid sensing in bacteria and archaea. More than 20 evolutionarily distinct families of ionotropic glutamate receptors (iGluRs) have been identified in eukaryotes. The domain compositions of iGluRs correlate with the origins of multicellularity in eukaryotes. Although L-Glu was recruited as a neuro-muscular transmitter in the early-branching metazoans, it was predominantly a non-neuronal messenger, with a possibility that glutamatergic synapses evolved more than once. Furthermore, the molecular secretory complexity of glutamatergic synapses in invertebrates (e.g., Aplysia) can exceed their vertebrate counterparts. Comparative genomics also revealed 15+ subfamilies of iGluRs across Metazoa. However, most of this ancestral diversity had been lost in the vertebrate lineage, preserving AMPA, Kainate, Delta, and NMDA receptors. The widespread expansion of glutamate synapses in the cortical areas might be associated with the enhanced metabolic demands of the complex brain and compartmentalization of Glu signaling within modular neuronal ensembles.

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.

Signals of Positive Selection in Sea Slug Transcriptomes

AbstractUnderstanding how species may respond to climate change is of paramount importance. Species that occupy highly heterogenous environments, such as intertidal zone estuarine habitats, provide an ideal test case for examining phenotypic and genomic adaptations to different environmental conditions, which may influence their response to rapidly shifting climatic conditions. The California coast is projected to experience changes in both temperature and salinity, which currently vary seasonally and latitudinally. Using comparative transcriptomics, I documented patterns of positive selection between the northern-dwelling planktotrophic sacoglossan sea slug Alderia modesta, which is remarkably tolerant of low temperatures and low salinities, and its southern congener Alderia willowi, which exhibits a striking flexibility for larval type in response to seasonally shifting changes in temperature and salinity. Out of over 4000 1-to-1 orthologous genes, I found a signal of positive selection between A. willowi and A. modesta for genes involved in cell membrane and cell transport, particularly ion homeostasis (aquaporin), cell-cell signal transduction, and phosphorylation (reduced nicotinamide adenine dinucleotide [NADH] dehydrogenase). Positive selection for ion homeostasis in A. modesta has implications for its ability to tolerate the lower salinity of its northern range, and in A. willowi substitutions in NADH may assist in the high temperature tolerance of its southern California habitats. Identifying these candidate genes enables future studies of their functionalization as we seek to understand the relationship between phenotype and genotype in species whose phenotypes are influenced by environmental conditions.

The balance of crystalline and amorphous regions in the fibroin structure underpins the tensile strength of bagworm silk

Protein-based materials are considered versatile biomaterials, and their biodegradability is an advantage for sustainable development. Bagworm produces strong silk for use in unique situations throughout its life stages. Rigorous molecular analyses of Eumeta variegata suggested that the particular mechanical properties of its silk are due to the coexistence of poly-A and GA motifs. However, little molecular information on closely related species is available, and it is not understood how these properties were acquired evolutionarily or whether the motif combination is a conserved trait in other bagworms. Here, we performed a transcriptome analysis of two other bagworm species (Canephora pungelerii and Bambalina sp.) belonging to the family Psychidae to elucidate the relationship between the fibroin gene and silk properties. The obtained transcriptome assemblies and tensile tests indicated that the motif combination and silk properties were conserved among the bagworms. Furthermore, our analysis showed that C. pungelerii produces extraordinarily strong silk (breaking strength of 1.4 GPa) and indicated that the cause may be the C. pungelerii -specific balance of crystalline/amorphous regions in the H-fibroin repetitive domain. This particular H-fibroin architecture may have been evolutionarily acquired to produce strong thread to maintain bag stability during the relatively long development period of Canephora species relative to other bagworms.

AniProtDB: A Collection of Consistently Generated Metazoan Proteomes for Comparative Genomics Studies

To address the void in the availability of high-quality proteomic data traversing the animal tree, we have implemented a pipeline for generating de novo assemblies based on publicly available data from the NCBI Sequence Read Archive, yielding a comprehensive collection of proteomes from 100 species spanning 21 animal phyla. We have also created the Animal Proteome Database (AniProtDB), a resource providing open access to this collection of high-quality metazoan proteomes, along with information on predicted proteins and protein domains for each taxonomic classification and the ability to perform sequence similarity searches against all proteomes generated using this pipeline. This solution vastly increases the utility of these data by removing the barrier to access for research groups who do not have the expertise or resources to generate these data themselves and enables the use of data from non-traditional research organisms that have the potential to address key questions in biomedicine.

An Overview of Studies on Meiofaunal Traits of the Littoral Zone of Lakes

We carried out an overview of the studies on the traits of the meiofauna of the littoral zone of lakes to investigate the question relating to the Raunkiaeran shortfall (lack of knowledge on biological traits). For this purpose, we selected a series of keywords associated with response and effect traits (feeding habits, locomotion and substrate relation, body size, shape and mass, life history, reproductive strategy, respiration and thermal tolerance) and we counted the relative frequency of occurrence in a set of scientific papers retrieved from Web of Science. The results showed that, except for the traits related to diet and feeding habits, the Raunkiaeran shortfall is very pronounced for all meiofaunal taxa of the littoral zone of lakes, especially for those related to soft-bodied organisms. The reason behind this deficiency concerns many aspects ranging from the high taxonomic expertise required to the intrinsic difficulties of observing organisms of such a small size. The relationship with temperature has not been sufficiently explored and formalized in any of the examined traits; this research aspect needs to be rapidly addressed since the prospects of climate change impacts on lake littorals are expected to be particularly severe.

Whole-genome resequencing reveals aberrant autosomal SNPs affect chicken feathering rate

Previous studies have shown that the feather growth rate of chicks is determined by two alleles located on the sex chromosome Z; however, in chicken production, feathering is usually not consistently controlled by the sex chromosome. To identify whether the feathering rate is related to autosomal inheritance, whole-genome resequencing was performed in eight chickens with slow- and fast-feathering rate. A total of 54,984 autosomal single nucleotide polymorphisms (SNPs) were identified, including 393 and 376 exonic SNPs in slow-feathering and fast-feathering chickens, respectively. Mutated genes were mainly involved in response to stimuli and growth and reproduction processes. Mutated genes related to slow-feathering rate were mainly involved in wingless-type MMTV integration site signaling pathway and mitogen-activated protein kinase signaling pathway, whereas mutated genes associated with fast-feathering rate were primarily enriched in autophagy, calcium signaling pathway, extracellular matrix-receptor interaction, and Focal adhesion processes. Importantly, two SNPs, involved in feather development, were found in the exonic regions of Wnt signaling genes. These results shed new light on the relationship between genetic mutation and feather growth rate from the perspective of autosomal inheritance and may have economic significance in chicken breeding.

Draft genome of Bugula neritina, a colonial animal packing powerful symbionts and potential medicines

Many animal phyla have no representatives within the catalog of whole metazoan genome sequences. This dataset fills in one gap in the genome knowledge of animal phyla with a draft genome of Bugula neritina (phylum Bryozoa). Interest in this species spans ecology and biomedical sciences because B. neritina is the natural source of bioactive compounds called bryostatins. Here we present a draft assembly of the B. neritina genome obtained from PacBio and Illumina HiSeq data, as well as genes and proteins predicted de novo and verified using transcriptome data, along with the functional annotation. These sequences will permit a better understanding of host-symbiont interactions at the genomic level, and also contribute additional phylogenomic markers to evaluate Lophophorate or Lophotrochozoa phylogenetic relationships. The effort also fits well with plans to ultimately sequence all orders of the Metazoa.

Measurement(s)	DNA • genome • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Bugula neritina

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12988355

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

Multi-model seascape genomics identifies distinct environmental drivers of selection among sympatric marine species

BACKGROUND

As global change and anthropogenic pressures continue to increase, conservation and management increasingly needs to consider species' potential to adapt to novel environmental conditions. Therefore, it is imperative to characterise the main selective forces acting on ecosystems, and how these may influence the evolutionary potential of populations and species. Using a multi-model seascape genomics approach, we compare putative environmental drivers of selection in three sympatric southern African marine invertebrates with contrasting ecology and life histories: Cape urchin (Parechinus angulosus), Common shore crab (Cyclograpsus punctatus), and Granular limpet (Scutellastra granularis).

RESULTS

Using pooled (Pool-seq), restriction-site associated DNA sequencing (RAD-seq), and seven outlier detection methods, we characterise genomic variation between populations along a strong biogeographical gradient. Of the three species, only S. granularis showed significant isolation-by-distance, and isolation-by-environment driven by sea surface temperatures (SST). In contrast, sea surface salinity (SSS) and range in air temperature correlated more strongly with genomic variation in C. punctatus and P. angulosus. Differences were also found in genomic structuring between the three species, with outlier loci contributing to two clusters in the East and West Coasts for S. granularis and P. angulosus, but not for C. punctatus.

CONCLUSION

The findings illustrate distinct evolutionary potential across species, suggesting that species-specific habitat requirements and responses to environmental stresses may be better predictors of evolutionary patterns than the strong environmental gradients within the region. We also found large discrepancies between outlier detection methodologies, and thus offer a novel multi-model approach to identifying the principal environmental selection forces acting on species. Overall, this work highlights how adding a comparative approach to seascape genomics (both with multiple models and species) can elucidate the intricate evolutionary responses of ecosystems to global change.

The Ecobiomics project: Advancing metagenomics assessment of soil health and freshwater quality in Canada

Transformative advances in metagenomics are providing an unprecedented ability to characterize the enormous diversity of microorganisms and invertebrates sustaining soil health and water quality. These advances are enabling a better recognition of the ecological linkages between soil and water, and the biodiversity exchanges between these two reservoirs. They are also providing new perspectives for understanding microorganisms and invertebrates as part of interacting communities (i.e. microbiomes and zoobiomes), and considering plants, animals, and humans as holobionts comprised of their own cells as well as diverse microorganisms and invertebrates often acquired from soil and water. The Government of Canada's Genomics Research and Development Initiative (GRDI) launched the Ecobiomics Project to coordinate metagenomics capacity building across federal departments, and to apply metagenomics to better characterize microbial and invertebrate biodiversity for advancing environmental assessment, monitoring, and remediation activities. The Project has adopted standard methods for soil, water, and invertebrate sampling, collection and provenance of metadata, and nucleic acid extraction. High-throughput sequencing is located at a centralized sequencing facility. A centralized Bioinformatics Platform was established to enable a novel government-wide approach to harmonize metagenomics data collection, storage and bioinformatics analyses. Sixteen research projects were initiated under Soil Microbiome, Aquatic Microbiome, and Invertebrate Zoobiome Themes. Genomic observatories were established at long-term environmental monitoring sites for providing more comprehensive biodiversity reference points to assess environmental change.

Estimation of the relative abundance of species in artificial mixtures of insects using low-coverage shotgun metagenomics

Amplicon metabarcoding is an established technique to analyse the taxonomic composition of communities of organisms using high-throughput DNA sequencing, but there are doubts about its ability to quantify the relative proportions of the species, as opposed to the species list. Here, we bypass the enrichment step and avoid the PCR-bias, by directly sequencing the extracted DNA using shotgun metagenomics. This approach is common practice in prokaryotes, but not in eukaryotes, because of the low number of sequenced genomes of eukaryotic species. We tested the metagenomics approach using insect species whose genome is already sequenced and assembled to an advanced degree. We shotgun-sequenced, at low-coverage, 18 species of insects in 22 single-species and 6 mixed-species libraries and mapped the reads against 110 reference genomes of insects. We used the single-species libraries to calibrate the process of assignation of reads to species and the libraries created from species mixtures to evaluate the ability of the method to quantify the relative species abundance. Our results showed that the shotgun metagenomic method is easily able to set apart closely-related insect species, like four species of Drosophila included in the artificial libraries. However, to avoid the counting of rare misclassified reads in samples, it was necessary to use a rather stringent detection limit of 0.001, so species with a lower relative abundance are ignored. We also identified that approximately half the raw reads were informative for taxonomic purposes. Finally, using the mixed-species libraries, we showed that it was feasible to quantify with confidence the relative abundance of individual species in the mixtures.

Improving Cancer Drug Discovery by Studying Cancer across the Tree of Life

Despite a considerable expenditure of time and resources and significant advances in experimental models of disease, cancer research continues to suffer from extremely low success rates in translating preclinical discoveries into clinical practice. The continued failure of cancer drug development, particularly late in the course of human testing, not only impacts patient outcomes, but also drives up the cost for those therapies that do succeed. It is clear that a paradigm shift is necessary if improvements in this process are to occur. One promising direction for increasing translational success is comparative oncology-the study of cancer across species, often involving veterinary patients that develop naturally-occurring cancers. Comparative oncology leverages the power of cross-species analyses to understand the fundamental drivers of cancer protective mechanisms, as well as factors contributing to cancer initiation and progression. Clinical trials in veterinary patients with cancer provide an opportunity to evaluate novel therapeutics in a setting that recapitulates many of the key features of human cancers, including genomic aberrations that underly tumor development, response and resistance to treatment, and the presence of comorbidities that can affect outcomes. With a concerted effort from basic scientists, human physicians and veterinarians, comparative oncology has the potential to enhance the cost-effectiveness and efficiency of pipelines for cancer drug discovery and other cancer treatments.

The Value of Reference Genomes in the Conservation of Threatened Species

Conservation initiatives are now more crucial than ever-over a million plant and animal species are at risk of extinction over the coming decades. The genetic management of threatened species held in insurance programs is recommended; however, few are taking advantage of the full range of genomic technologies available today. Less than 1% of the 13505 species currently listed as threated by the International Union for Conservation of Nature (IUCN) have a published genome. While there has been much discussion in the literature about the importance of genomics for conservation, there are limited examples of how having a reference genome has changed conservation management practice. The Tasmanian devil (Sarcophilus harrisii), is an endangered Australian marsupial, threatened by an infectious clonal cancer devil facial tumor disease (DFTD). Populations have declined by 80% since the disease was first recorded in 1996. A reference genome for this species was published in 2012 and has been crucial for understanding DFTD and the management of the species in the wild. Here we use the Tasmanian devil as an example of how a reference genome has influenced management actions in the conservation of a species.

Evidence for Faster X Chromosome Evolution in Spiders

In species with chromosomal sex determination, X chromosomes are predicted to evolve faster than autosomes because of positive selection on recessive alleles or weak purifying selection. We investigated X chromosome evolution in Stegodyphus spiders that differ in mating system, sex ratio, and population dynamics. We assigned scaffolds to X chromosomes and autosomes using a novel method based on flow cytometry of sperm cells and reduced representation sequencing. We estimated coding substitution patterns (dN/dS) in a subsocial outcrossing species (S. africanus) and its social inbreeding and female-biased sister species (S. mimosarum), and found evidence for faster-X evolution in both species. X chromosome-to-autosome diversity (piX/piA) ratios were estimated in multiple populations. The average piX/piA estimates of S. africanus (0.57 [95% CI: 0.55-0.60]) was lower than the neutral expectation of 0.75, consistent with more hitchhiking events on X-linked loci and/or a lower X chromosome mutation rate, and we provide evidence in support of both. The social species S. mimosarum has a significantly higher piX/piA ratio (0.72 [95% CI: 0.65-0.79]) in agreement with its female-biased sex ratio. Stegodyphus mimosarum also have different piX/piA estimates among populations, which we interpret as evidence for recurrent founder events. Simulations show that recurrent founder events are expected to decrease the piX/piA estimates in S. mimosarum, thus underestimating the true effect of female-biased sex ratios. Finally, we found lower synonymous divergence on X chromosomes in both species, and the male-to-female substitution ratio to be higher than 1, indicating a higher mutation rate in males.

Polymorphic microsatellite loci and partial mitogenome for the Chestnut-bellied Seed-finch Sporophila angolensis (Aves, Passeriformes) using next generation sequencing

Brazil is one of the major contributors to international trade in wildlife and species of the bird genus Sporophila are currently under threat due to illegal trade. Microsatellite loci and mitochondrial DNA constitute important molecular markers for population genetics studies and parentage analyses, and hold great potential to help authorities manage illegal trafficking and control commercial breeders. We describe and characterize 19 polymorphic microsatellite loci and recover part of the mitochondrial genome for Sporophila angolensis using massive parallel sequencing with the Illumina platform. DNA sequencing resulted in a dataset with 2,379,295 paired reads, of which 392 were mapped to the mitogenome of S. maximiliani, resulting in a partial mitogenome of 16,785 bp for S. angolensis. The microsatellite search identified a total of 4737 loci, from which 27 primer pairs were tested on 24 individuals of unknown geographic origin. Nineteen of the 27 loci were successfully amplified and exhibited high levels of genetic variation, with a mean of 11.2 alleles per locus, a mean observed heterozygosity of 0.588 and a mean expected heterozygosity of 0.852. About half of the loci showed significant evidence for the presence of a null allele and significant deviation from Hardy-Weinberg equilibrium; the remaining eight loci had high paternity exclusion probabilities and low identity probabilities. The high levels of polymorphism for these loci, as well as their high paternity exclusion probability and low identity probability, indicate that they hold potential for parentage analyses and population genetics studies of S. angolensis.

SQUAT: a Sequencing Quality Assessment Tool for data quality assessments of genome assemblies

Background

With the rapid increase in genome sequencing projects for non-model organisms, numerous genome assemblies are currently in progress or available as drafts, but not made available as satisfactory, usable genomes. Data quality assessment of genome assemblies is gaining importance not only for people who perform the assembly/re-assembly processes, but also for those who attempt to use assemblies as maps in downstream analyses. Recent studies of the quality control, quality evaluation/ assessment of genome assemblies have focused on either quality control of reads before assemblies or evaluation of the assemblies with respect to their contiguity and correctness. However, correctness assessment depends on a reference and is not applicable for de novo assembly projects. Hence, development of methods providing both post-assembly and pre-assembly quality assessment reports for examining the quality/correctness of de novo assemblies and the input reads is worth studying.

Results

We present SQUAT, an efficient tool for both pre-assembly and post-assembly quality assessment of de novo genome assemblies. The pre-assembly module of SQUAT computes quality statistics of reads and presents the analysis in a well-designed interface to visualize the distribution of high- and poor-quality reads in a portable HTML report. The post-assembly module of SQUAT provides read mapping analytics in an HTML format. We categorized reads into several groups including uniquely mapped reads, multiply mapped, unmapped reads; for uniquely mapped reads, we further categorized them into perfectly matched, with substitutions, containing clips, and the others. We carefully defined the poorly mapped (PM) reads into several groups to prevent the underestimation of unmapped reads; indeed, a high PM% would be a sign of a poor assembly that requires researchers’ attention for further examination or improvements before using the assembly. Finally, we evaluate SQUAT with six datasets, including the genome assemblies for eel, worm, mushroom, and three bacteria. The results show that SQUAT reports provide useful information with details for assessing the quality of assemblies and reads.

Availability

The SQUAT software with links to both its docker image and the on-line manual is freely available at https://github.com/luke831215/SQUAT.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5445-3) contains supplementary material, which is available to authorized users.

Multiple Facets of Marine Invertebrate Conservation Genomics

Conservation genomics aims to preserve the viability of populations and the biodiversity of living organisms. Invertebrate organisms represent 95% of animal biodiversity; however, few genomic resources currently exist for the group. The subset of marine invertebrates includes the most ancient metazoan lineages and possesses codes for unique gene products and possible keys to adaptation. The benefits of supporting invertebrate conservation genomics research (e.g., likely discovery of novel genes, protein regulatory mechanisms, genomic innovations, and transposable elements) outweigh the various hurdles (rare, small, or polymorphic starting materials). Here we review best conservation genomics practices in the laboratory and in silico when applied to marine invertebrates and also showcase unique features in several case studies of acroporid corals, crown-of-thorns starfish, apple snails, and abalone. Marine conservation genomics should also address how diversity can lead to unique marine innovations, the impact of deleterious variation, and how genomic monitoring and profiling could positively affect broader conservation goals (e.g., value of baseline data for in situ/ex situ genomic stocks).

Biodiversity estimates and ecological interpretations of meiofaunal communities are biased by the taxonomic approach

Accurate assessments of biodiversity are crucial to advising ecosystem-monitoring programs and understanding ecosystem function. Nevertheless, a standard operating procedure to assess biodiversity accurately and consistently has not been established. This is especially true for meiofauna, a diverse community (>20 phyla) of small benthic invertebrates that have fundamental ecological roles. Recent studies show that metabarcoding is a cost-effective and time-effective method to estimate meiofauna biodiversity, in contrast to morphological-based taxonomy. Here, we compare biodiversity assessments of a diverse meiofaunal community derived by applying multiple taxonomic methods based on comparative morphology, molecular phylogenetic analysis, DNA barcoding of individual specimens, and metabarcoding of environmental DNA. We show that biodiversity estimates are strongly biased across taxonomic methods and phyla. Such biases affect understanding of community structures and ecological interpretations. This study supports the urgency of improving aspects of environmental high-throughput sequencing and the value of taxonomists in correctly understanding biodiversity estimates.

Earth BioGenome Project: Sequencing life for the future of life

Increasing our understanding of Earth's biodiversity and responsibly stewarding its resources are among the most crucial scientific and social challenges of the new millennium. These challenges require fundamental new knowledge of the organization, evolution, functions, and interactions among millions of the planet's organisms. Herein, we present a perspective on the Earth BioGenome Project (EBP), a moonshot for biology that aims to sequence, catalog, and characterize the genomes of all of Earth's eukaryotic biodiversity over a period of 10 years. The outcomes of the EBP will inform a broad range of major issues facing humanity, such as the impact of climate change on biodiversity, the conservation of endangered species and ecosystems, and the preservation and enhancement of ecosystem services. We describe hurdles that the project faces, including data-sharing policies that ensure a permanent, freely available resource for future scientific discovery while respecting access and benefit sharing guidelines of the Nagoya Protocol. We also describe scientific and organizational challenges in executing such an ambitious project, and the structure proposed to achieve the project's goals. The far-reaching potential benefits of creating an open digital repository of genomic information for life on Earth can be realized only by a coordinated international effort.

Understanding Animal Evolution: The Added Value of Sponge Transcriptomics and Genomics: The disconnect between gene content and body plan evolution

Sponges are important but often-neglected organisms. The absence of classical animal traits (nerves, digestive tract, and muscles) makes sponges challenging for non-specialists to work with and has delayed getting high quality genomic data compared to other invertebrates. Yet analyses of sponge genomes and transcriptomes currently available have radically changed our understanding of animal evolution. Sponges are of prime evolutionary importance as one of the best candidates to form the sister group of all other animals, and genomic data are essential to understand the mechanisms that control animal evolution and diversity. Here we review the most significant outcomes of current genomic and transcriptomic analyses of sponges, and discuss limitations and future directions of sponge transcriptomic and genomic studies.

Ancient Origin of the CARD-Coiled Coil/Bcl10/MALT1-Like Paracaspase Signaling Complex Indicates Unknown Critical Functions

The CARD–coiled coil (CC)/Bcl10/MALT1-like paracaspase (CBM) signaling complexes composed of a CARD–CC family member (CARD-9, -10, -11, or -14), Bcl10, and the type 1 paracaspase MALT1 (PCASP1) play a pivotal role in immunity, inflammation, and cancer. Targeting MALT1 proteolytic activity is of potential therapeutic interest. However, little is known about the evolutionary origin and the original functions of the CBM complex. Type 1 paracaspases originated before the last common ancestor of planulozoa (bilaterians and cnidarians). Notably in bilaterians, Ecdysozoa (e.g., nematodes and insects) lacks Bcl10, whereas other lineages have a Bcl10 homolog. A survey of invertebrate CARD–CC homologs revealed such homologs only in species with Bcl10, indicating an ancient common origin of the entire CBM complex. Furthermore, vertebrate-like Syk/Zap70 tyrosine kinase homologs with the ITAM-binding SH2 domain were only found in invertebrate organisms with CARD–CC/Bcl10, indicating that this pathway might be related to the original function of the CBM complex. Moreover, the type 1 paracaspase sequences from invertebrate organisms that have CARD–CC/Bcl10 are more similar to vertebrate paracaspases. Functional analysis of protein–protein interactions, NF-κB signaling, and CYLD cleavage for selected invertebrate type 1 paracaspase and Bcl10 homologs supports this scenario and indicates an ancient origin of the CARD–CC/Bcl10/paracaspase signaling complex. By contrast, many of the known MALT1-associated activities evolved fairly recently, indicating that unknown functions are at the basis of the protein conservation. As a proof-of-concept, we provide initial evidence for a CBM- and NF-κB-independent neuronal function of the Caenorhabditis elegans type 1 paracaspase malt-1. In conclusion, this study shows how evolutionary insights may point at alternative functions of MALT1.

Extreme Mitogenomic Variation in Natural Populations of Chaetognaths

The extent of within-species genetic variation across the diversity of animal life is an underexplored problem in ecology and evolution. Although neutral genetic variation should scale positively with population size, mitochondrial diversity levels are believed to show little variation across animal species. Here, we report an unprecedented case of extreme mitochondrial diversity within natural populations of two morphospecies of chaetognaths (arrow worms). We determine that this diversity is composed of deep sympatric mitochondrial lineages, which are in some cases as divergent as human and platypus. Additionally, based on 54 complete mitogenomes, we observed mitochondrial gene order differences between several of these lineages. We examined nuclear divergence patterns (18S, 28S, and an intron) to determine the possible origin of these lineages, but did not find congruent patterns between mitochondrial and nuclear markers. We also show that extreme mitochondrial divergence in chaetognaths is not driven by positive selection. Hence, we propose that the extreme levels of mitochondrial variation could be the result of either a complex scenario of reproductive isolation, or a combination of large population size and accelerated mitochondrial mutation rate. These findings emphasize the importance of characterizing genome-wide levels of nuclear variation in these species and promote chaetognaths as a remarkable model to study mitochondrial evolution.

Scaling up: A guide to high-throughput genomic approaches for biodiversity analysis

The purpose of this review is to present the most common and emerging DNA-based methods used to generate data for biodiversity and biomonitoring studies. As environmental assessment and monitoring programmes may require biodiversity information at multiple levels, we pay particular attention to the DNA metabarcoding method and discuss a number of bioinformatic tools and considerations for producing DNA-based indicators using operational taxonomic units (OTUs), taxa at a variety of ranks and community composition. By developing the capacity to harness the advantages provided by the newest technologies, investigators can "scale up" by increasing the number of samples and replicates processed, the frequency of sampling over time and space, and even the depth of sampling such as by sequencing more reads per sample or more markers per sample. The ability to scale up is made possible by the reduced hands-on time and cost per sample provided by the newest kits, platforms and software tools. Results gleaned from broad-scale monitoring will provide opportunities to address key scientific questions linked to biodiversity and its dynamics across time and space as well as being more relevant for policymakers, enabling science-based decision-making, and provide a greater socio-economic impact. As genomic approaches are continually evolving, we provide this guide to methods used in biodiversity genomics.

Comparative immunogenomics of molluscs

Comparative immunology, studying both vertebrates and invertebrates, provided the earliest descriptions of phagocytosis as a general immune mechanism. However, the large scale of animal diversity challenges all-inclusive investigations and the field of immunology has developed by mostly emphasizing study of a few vertebrate species. In addressing the lack of comprehensive understanding of animal immunity, especially that of invertebrates, comparative immunology helps toward management of invertebrates that are food sources, agricultural pests, pathogens, or transmit diseases, and helps interpret the evolution of animal immunity. Initial studies showed that the Mollusca (second largest animal phylum), and invertebrates in general, possess innate defenses but lack the lymphocytic immune system that characterizes vertebrate immunology. Recognizing the reality of both common and taxon-specific immune features, and applying up-to-date cell and molecular research capabilities, in-depth studies of a select number of bivalve and gastropod species continue to reveal novel aspects of molluscan immunity. The genomics era heralded a new stage of comparative immunology; large-scale efforts yielded an initial set of full molluscan genome sequences that is available for analyses of full complements of immune genes and regulatory sequences. Next-generation sequencing (NGS), due to lower cost and effort required, allows individual researchers to generate large sequence datasets for growing numbers of molluscs. RNAseq provides expression profiles that enable discovery of immune genes and genome sequences reveal distribution and diversity of immune factors across molluscan phylogeny. Although computational de novo sequence assembly will benefit from continued development and automated annotation may require some experimental validation, NGS is a powerful tool for comparative immunology, especially increasing coverage of the extensive molluscan diversity. To date, immunogenomics revealed new levels of complexity of molluscan defense by indicating sequence heterogeneity in individual snails and bivalves, and members of expanded immune gene families are expressed differentially to generate pathogen-specific defense responses.

Whole-genome sequencing approaches for conservation biology: Advantages, limitations and practical recommendations

Whole-genome resequencing (WGR) is a powerful method for addressing fundamental evolutionary biology questions that have not been fully resolved using traditional methods. WGR includes four approaches: the sequencing of individuals to a high depth of coverage with either unresolved or resolved haplotypes, the sequencing of population genomes to a high depth by mixing equimolar amounts of unlabelled-individual DNA (Pool-seq) and the sequencing of multiple individuals from a population to a low depth (lcWGR). These techniques require the availability of a reference genome. This, along with the still high cost of shotgun sequencing and the large demand for computing resources and storage, has limited their implementation in nonmodel species with scarce genomic resources and in fields such as conservation biology. Our goal here is to describe the various WGR methods, their pros and cons and potential applications in conservation biology. WGR offers an unprecedented marker density and surveys a wide diversity of genetic variations not limited to single nucleotide polymorphisms (e.g., structural variants and mutations in regulatory elements), increasing their power for the detection of signatures of selection and local adaptation as well as for the identification of the genetic basis of phenotypic traits and diseases. Currently, though, no single WGR approach fulfils all requirements of conservation genetics, and each method has its own limitations and sources of potential bias. We discuss proposed ways to minimize such biases. We envision a not distant future where the analysis of whole genomes becomes a routine task in many nonmodel species and fields including conservation biology.

Mitochondrial RNA processing in absence of tRNA punctuations in octocorals

BACKGROUND

Mitogenome diversity is staggering among early branching animals with respect to size, gene density, content and order, and number of tRNA genes, especially in cnidarians. This last point is of special interest as tRNA cleavage drives the maturation of mitochondrial mRNAs and is a primary mechanism for mt-RNA processing in animals. Mitochondrial RNA processing in non-bilaterian metazoans, some of which possess a single tRNA gene in their mitogenomes, is essentially unstudied despite its importance in understanding the evolution of mitochondrial transcription in animals.

RESULTS

We characterized the mature mitochondrial mRNA transcripts in a species of the octocoral genus Sinularia (Alcyoniidae: Octocorallia), and defined precise boundaries of transcription units using different molecular methods. Most mt-mRNAs were polycistronic units containing two or three genes and 5' and/or 3' untranslated regions of varied length. The octocoral specific, mtDNA-encoded mismatch repair gene, the mtMutS, was found to undergo alternative polyadenylation, and exhibited differential expression of alternate transcripts suggesting a unique regulatory mechanism for this gene. In addition, a long noncoding RNA complementary to the ATP6 gene (lncATP6) potentially involved in antisense regulation was detected.

CONCLUSIONS

Mt-mRNA processing in octocorals possessing a single mt-tRNA is complex. Considering the variety of mitogenome arrangements known in cnidarians, and in general among non-bilaterian metazoans, our findings provide a first glimpse into the complex mtDNA transcription, mt-mRNA processing, and regulation among early branching animals and represent a first step towards understanding its functional and evolutionary implications.

Advancing Genomics through the Global Invertebrate Genomics Alliance (GIGA)

The Global Invertebrate Genomics Alliance (GIGA), a collaborative network of diverse scientists, marked its second anniversary with a workshop in Munich, Germany, where international attendees focused on discussing current progress, milestones and bioinformatics resources. The community determined the recruitment and training talented researchers as one of the most pressing future needs and identified opportunities for network funding. GIGA also promotes future research efforts to prioritize taxonomic diversity and create new synergies. Here, we announce the generation of a central and simple data repository portal with a wide coverage of available sequence data, via the compagen platform, in parallel with more focused and specialized organism databases to globally advance invertebrate genomics. Therefore this article serves the objectives of GIGA by disseminating current progress and future prospects in the science of invertebrate genomics with the aim of promotion and facilitation of interdisciplinary and international research.

An annotated draft genome for Radix auricularia (Gastropoda, Mollusca)

Molluscs are the second most species-rich phylum in the animal kingdom, yet only 11 genomes of this group have been published so far. Here, we present the draft genome sequence of the pulmonate freshwater snail Radix auricularia . Six whole genome shotgun libraries with different layouts were sequenced. The resulting assembly comprises 4,823 scaffolds with a cumulative length of 910 Mb and an overall read coverage of 72×. The assembly contains 94.6% of a metazoan core gene collection, indicating an almost complete coverage of the coding fraction. The discrepancy of ∼690 Mb compared with the estimated genome size of R. auricularia (1.6 Gb) results from a high repeat content of 70% mainly comprising DNA transposons. The annotation of 17,338 protein coding genes was supported by the use of publicly available transcriptome data. This draft will serve as starting point for further genomic and population genetic research in this scientifically important phylum.

Reefgenomics.Org - a repository for marine genomics data

Over the last decade, technological advancements have substantially decreased the cost and time of obtaining large amounts of sequencing data. Paired with the exponentially increased computing power, individual labs are now able to sequence genomes or transcriptomes to investigate biological questions of interest. This has led to a significant increase in available sequence data. Although the bulk of data published in articles are stored in public sequence databases, very often, only raw sequencing data are available; miscellaneous data such as assembled transcriptomes, genome annotations etc. are not easily obtainable through the same means. Here, we introduce our website (http://reefgenomics.org) that aims to centralize genomic and transcriptomic data from marine organisms. Besides providing convenient means to download sequences, we provide (where applicable) a genome browser to explore available genomic features, and a BLAST interface to search through the hosted sequences. Through the interface, multiple datasets can be queried simultaneously, allowing for the retrieval of matching sequences from organisms of interest. The minimalistic, no-frills interface reduces visual clutter, making it convenient for end-users to search and explore processed sequence data. DATABASE URL: http://reefgenomics.org.

Greater than X kb: a quantitative assessment of preservation conditions on genomic DNA quality, and a proposed standard for genome-quality DNA

Advances in biodiversity genomic sequencing will increasingly depend on the availability of DNA samples-and their quantifiable metadata-preserved in large institutional biorepositories that are discoverable to the scientific community. Improvements in sequencing technology constantly provide longer reads, such that longer fragment length, higher molecular weight, and overall "genome-quality" DNA (gDNA) will be desirable. Ideally, biorepositories should publish numerical scale measurements of DNA quality useful to the user community. However, the most widely used technique to evaluate DNA quality, the classic agarose gel, has yet to be quantified. Here we propose a simple and economical method using open source image analysis software to make gDNA gel images quantifiable, and propose percentage of gDNA "greater than X kb" as a standard of comparison, where X is a band from any widely used DNA ladder with desirably large band sizes. We employ two metadata standards ("DNA Threshold" and "Percent above Threshold") introduced as part of the Global Genome Biodiversity Network (GGBN) Darwin Core extension. We illustrate the method using the traditionally used HindIII ladder and the 9,416 base-pair (bp) band as a standard. We also present data, for two taxa, a vertebrate (fish) and an invertebrate (crab), on how gDNA quality varies with seven tissue preservation methods, time since death, preservation method (i.e. buffers vs. cold temperatures), and storage temperature of various buffers over time. Our results suggest that putting tissue into a buffer prior to freezing may be better than directly into ultra-cold conditions.

Greater than X kb: a quantitative assessment of preservation conditions on genomic DNA quality, and a proposed standard for genome-quality DNA

Advances in biodiversity genomic sequencing will increasingly depend on the availability of DNA samples—and their quantifiable metadata—preserved in large institutional biorepositories that are discoverable to the scientific community. Improvements in sequencing technology constantly provide longer reads, such that longer fragment length, higher molecular weight, and overall “genome-quality” DNA (gDNA) will be desirable. Ideally, biorepositories should publish numerical scale measurements of DNA quality useful to the user community. However, the most widely used technique to evaluate DNA quality, the classic agarose gel, has yet to be quantified. Here we propose a simple and economical method using open source image analysis software to make gDNA gel images quantifiable, and propose percentage of gDNA “greater than X kb” as a standard of comparison, where X is a band from any widely used DNA ladder with desirably large band sizes. We employ two metadata standards (“DNA Threshold” and “Percent above Threshold”) introduced as part of the Global Genome Biodiversity Network (GGBN) Darwin Core extension. We illustrate the method using the traditionally used HindIII ladder and the 9,416 base-pair (bp) band as a standard. We also present data, for two taxa, a vertebrate (fish) and an invertebrate (crab), on how gDNA quality varies with seven tissue preservation methods, time since death, preservation method (i.e. buffers vs. cold temperatures), and storage temperature of various buffers over time. Our results suggest that putting tissue into a buffer prior to freezing may be better than directly into ultra-cold conditions.

Animal Mitochondrial DNA as We Do Not Know It: mt-Genome Organization and Evolution in Nonbilaterian Lineages

Animal mitochondrial DNA (mtDNA) is commonly described as a small, circular molecule that is conserved in size, gene content, and organization. Data collected in the last decade have challenged this view by revealing considerable diversity in animal mitochondrial genome organization. Much of this diversity has been found in nonbilaterian animals (phyla Cnidaria, Ctenophora, Placozoa, and Porifera), which, from a phylogenetic perspective, form the main branches of the animal tree along with Bilateria. Within these groups, mt-genomes are characterized by varying numbers of both linear and circular chromosomes, extra genes (e.g. atp9, polB, tatC), large variation in the number of encoded mitochondrial transfer RNAs (tRNAs) (0-25), at least seven different genetic codes, presence/absence of introns, tRNA and mRNA editing, fragmented ribosomal RNA genes, translational frameshifting, highly variable substitution rates, and a large range of genome sizes. This newly discovered diversity allows a better understanding of the evolutionary plasticity and conservation of animal mtDNA and provides insights into the molecular and evolutionary mechanisms shaping mitochondrial genomes.

A New Orthology Assessment Method for Phylogenomic Data: Unrooted Phylogenetic Orthology

Current sequencing technologies are making available unprecedented amounts of genetic data for a large variety of species including nonmodel organisms. Although many phylogenomic surveys spend considerable time finding orthologs from the wealth of sequence data, these results do not transcend the original study and after being processed for specific phylogenetic purposes these orthologs do not become stable orthology hypotheses. We describe a procedure to detect and document the phylogenetic distribution of orthologs allowing researchers to use this information to guide selection of loci best suited to test specific evolutionary questions. At the core of this pipeline is a new phylogenetic orthology method that is neither affected by the position of the root nor requires explicit assignment of outgroups. We discuss the properties of this new orthology assessment method and exemplify its utility for phylogenomics using a small insects dataset. In addition, we exemplify the pipeline to identify and document stable orthologs for the group of orb-weaving spiders (Araneoidea) using RNAseq data. The scripts used in this study, along with sample files and additional documentation, are available at https://github.com/ballesterus/UPhO.

The miRNA biogenesis in marine bivalves

Small non-coding RNAs include powerful regulators of gene expression, transposon mobility and virus activity. Among the various categories, mature microRNAs (miRNAs) guide the translational repression and decay of several targeted mRNAs. The biogenesis of miRNAs depends on few gene products, essentially conserved from basal to higher metazoans, whose protein domains allow specific interactions with dsRNA. Here, we report the identification of key genes responsible of the miRNA biogenesis in 32 bivalves, with particular attention to the aquaculture species Mytilus galloprovincialis and Crassostrea gigas. In detail, we have identified and phylogenetically compared eight evolutionary conserved proteins: DROSHA, DGCR8, EXP5, RAN, DICER TARBP2, AGO and PIWI. In mussels, we recognized several other proteins participating in the miRNA biogenesis or in the subsequent RNA silencing. According to digital expression analysis, these genes display low and not inducible expression levels in adult mussels and oysters whereas they are considerably expressed during development. As miRNAs play an important role also in the antiviral responses, knowledge on their production and regulative effects can shed light on essential molecular processes and provide new hints for disease prevention in bivalves.