The social brain: transcriptome assembly and characterization of the hippocampus from a social subterranean rodent, the colonial tuco-tuco (Ctenomys sociabilis)

Unusual occurrence of domestication syndrome amongst African mole-rats: Is the naked mole-rat a domestic animal?

The Naked mole-rat (NMR) is becoming a prominent model organism due to its peculiar traits, such as eusociality, extreme longevity, cancer resistance, and reduced pain sensitivity. It belongs to the African mole-rats (AMR), a family of subterranean rodents that includes solitary, cooperative breeding and eusocial species. We identified and quantified the domestication syndrome (DS) across AMR, a set of morphological and behavioural traits significantly more common and pronounced amongst domesticated animals than in their wild counterparts. Surprisingly, the NMR shows apparent DS traits when compared to the solitary AMR. Animals can self-domesticate when a reduction of the fear response is naturally selected, such as in islands with no predators, or to improve the group’s harmony in cooperative breeding species. The DS may be caused by alterations in the physiology of the neural crest cells (NCC), a transient population of cells that generate a full range of tissues during development. The NCC contribute to organs responsible for transmitting the fear response and various other tissues, including craniofacial bones. Therefore, mutations affecting the NCC can manifest as behavioural and morphological alterations in many structures across the body, as seen in neurocristopathies. We observed that all social AMRs are chisel-tooth diggers, an adaption to hard soils that requires the flattening of the skull. We hypothesise that chisel-tooth digging could impose a selective pressure on the NCC that triggered the DS’s appearance, possibly facilitating the evolution of sociality. Finally, we discuss how DS traits are neutral or beneficial for the subterranean niche, strategies to test this hypothesis and report well-studied mutations in the NMR that are associated with the NCC physiology or with the control of the fear response. In conclusion, we argue that many of the NMR’s unconventional traits are compatible with the DS and provide a hypothesis about its origins. Our model proposes a novel avenue to enhance the understanding of the extraordinary biology of the NMR.

A comparative genomics multitool for scientific discovery and conservation

The Zoonomia Project is investigating the genomics of shared and specialized traits in eutherian mammals. Here we provide genome assemblies for 131 species, of which all but 9 are previously uncharacterized, and describe a whole-genome alignment of 240 species of considerable phylogenetic diversity, comprising representatives from more than 80% of mammalian families. We find that regions of reduced genetic diversity are more abundant in species at a high risk of extinction, discern signals of evolutionary selection at high resolution and provide insights from individual reference genomes. By prioritizing phylogenetic diversity and making data available quickly and without restriction, the Zoonomia Project aims to support biological discovery, medical research and the conservation of biodiversity.

A whole-genome alignment of 240 phylogenetically diverse species of eutherian mammal—including 131 previously uncharacterized species—from the Zoonomia Project provides data that support biological discovery, medical research and conservation.

Giant Island Mice Exhibit Widespread Gene Expression Changes in Key Metabolic Organs

Island populations repeatedly evolve extreme body sizes, but the genomic basis of this pattern remains largely unknown. To understand how organisms on islands evolve gigantism, we compared genome-wide patterns of gene expression in Gough Island mice, the largest wild house mice in the world, and mainland mice from the WSB/EiJ wild-derived inbred strain. We used RNA-seq to quantify differential gene expression in three key metabolic organs: gonadal adipose depot, hypothalamus, and liver. Between 4,000 and 8,800 genes were significantly differentially expressed across the evaluated organs, representing between 20% and 50% of detected transcripts, with 20% or more of differentially expressed transcripts in each organ exhibiting expression fold changes of at least 2×. A minimum of 73 candidate genes for extreme size evolution, including Irs1 and Lrp1, were identified by considering differential expression jointly with other data sets: 1) genomic positions of published quantitative trait loci for body weight and growth rate, 2) whole-genome sequencing of 16 wild-caught Gough Island mice that revealed fixed single-nucleotide differences between the strains, and 3) publicly available tissue-specific regulatory elements. Additionally, patterns of differential expression across three time points in the liver revealed that Arid5b potentially regulates hundreds of genes. Functional enrichment analyses pointed to cell cycling, mitochondrial function, signaling pathways, inflammatory response, and nutrient metabolism as potential causes of weight accumulation in Gough Island mice. Collectively, our results indicate that extensive gene regulatory evolution in metabolic organs accompanied the rapid evolution of gigantism during the short time house mice have inhabited Gough Island.

Identifying similar transcripts in a related organism from de Bruijn graphs of RNA-Seq data, with applications to the study of salt and waterlogging tolerance in Melilotus

Background

A popular strategy to study alternative splicing in non-model organisms starts from sequencing the entire transcriptome, then assembling the reads by using de novo transcriptome assembly algorithms to obtain predicted transcripts. A similarity search algorithm is then applied to a related organism to infer possible function of these predicted transcripts. While some of these predictions may be inaccurate and transcripts with low coverage are often missed, we observe that it is possible to obtain a more complete set of transcripts to facilitate possible functional assignments by starting the search from the intermediate de Bruijn graph that contains all branching possibilities.

Results

We develop an algorithm to extract similar transcripts in a related organism by starting the search from the de Bruijn graph that represents the transcriptome instead of from predicted transcripts. We show that our algorithm is able to recover more similar transcripts than existing algorithms, with large improvements in obtaining longer transcripts and a finer resolution of isoforms. We apply our algorithm to study salt and waterlogging tolerance in two Melilotus species by constructing new RNA-Seq libraries.

Conclusions

We have developed an algorithm to identify paths in the de Bruijn graph that correspond to similar transcripts in a related organism directly. Our strategy bypasses the transcript prediction step in RNA-Seq data and makes use of support from evolutionary information.

Electronic supplementary material

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

De novo assembly, annotation, and characterization of the whole brain transcriptome of male and female Syrian hamsters

Hamsters are an ideal animal model for a variety of biomedical research areas such as cancer, virology, circadian rhythms, and behavioural neuroscience. The use of hamsters has declined, however, most likely due to the dearth of genetic tools available for these animals. Our laboratory uses hamsters to study acute social stress, and we are beginning to investigate the genetic mechanisms subserving defeat-induced behavioural change. We have been limited, however, by the lack of genetic resources available for hamsters. In this study, we sequenced the brain transcriptome of male and female Syrian hamsters to generate the necessary resources to continue our research. We completed a de novo assembly and after assembly optimization, there were 113,329 transcripts representing 14,530 unique genes. This study is the first to characterize transcript expression in both female and male hamster brains and offers invaluable information to promote understanding of a host of important biomedical research questions for which hamsters are an excellent model.

Peromyscus transcriptomics: Understanding adaptation and gene expression plasticity within and between species of deer mice

Deer mice in the genus Peromyscus occupy nearly every terrestrial habitat in North America, and have a long history as subjects of behavioral, ecological, evolutionary, and physiological study. Recent advances in transcriptomics, the study of the complete set of RNA transcripts produced by certain cell types or under certain conditions, have contributed to the development of Peromyscus as a model system. We review the recent use of transcriptomics to investigate how natural selection and gene expression plasticity contribute to the existence of deer mice in challenging environments such as highlands, deserts, and cities across North America. Transcriptomics also holds great promise for elucidating the genetic basis of mating systems and other behaviors in Peromyscus, but has to date been underutilized for developmental biology and disease studies. Future Peromyscus studies should apply robust comparative frameworks to analyze the transcriptomics of multiple populations of the same species across varying environmental conditions, as well as multiple species that vary in traits of interest.

Characterization of the transcriptome, nucleotide sequence polymorphism, and natural selection in the desert adapted mouse Peromyscus eremicus

As a direct result of intense heat and aridity, deserts are thought to be among the most harsh of environments, particularly for their mammalian inhabitants. Given that osmoregulation can be challenging for these animals, with failure resulting in death, strong selection should be observed on genes related to the maintenance of water and solute balance. One such animal, Peromyscus eremicus, is native to the desert regions of the southwest United States and may live its entire life without oral fluid intake. As a first step toward understanding the genetics that underlie this phenotype, we present a characterization of the P. eremicus transcriptome. We assay four tissues (kidney, liver, brain, testes) from a single individual and supplement this with population level renal transcriptome sequencing from 15 additional animals. We identified a set of transcripts undergoing both purifying and balancing selection based on estimates of Tajima's D. In addition, we used the branch-site test to identify a transcript-Slc2a9, likely related to desert osmoregulation-undergoing enhanced selection in P. eremicus relative to a set of related non-desert rodents.

Brain transcriptome sequencing and assembly of three songbird model systems for the study of social behavior

Emberizid sparrows (emberizidae) have played a prominent role in the study of avian vocal communication and social behavior. We present here brain transcriptomes for three emberizid model systems, song sparrow Melospiza melodia, white-throated sparrow Zonotrichia albicollis, and Gambel's white-crowned sparrow Zonotrichia leucophrys gambelii. Each of the assemblies covered fully or in part, over 89% of the previously annotated protein coding genes in the zebra finch Taeniopygia guttata, with 16,846, 15,805, and 16,646 unique BLAST hits in song, white-throated and white-crowned sparrows, respectively. As in previous studies, we find tissue of origin (auditory forebrain versus hypothalamus and whole brain) as an important determinant of overall expression profile. We also demonstrate the successful isolation of RNA and RNA-sequencing from post-mortem samples from building strikes and suggest that such an approach could be useful when traditional sampling opportunities are limited. These transcriptomes will be an important resource for the study of social behavior in birds and for data driven annotation of forthcoming whole genome sequences for these and other bird species.

On the optimal trimming of high-throughput mRNA sequence data

The widespread and rapid adoption of high-throughput sequencing technologies has afforded researchers the opportunity to gain a deep understanding of genome level processes that underlie evolutionary change, and perhaps more importantly, the links between genotype and phenotype. In particular, researchers interested in functional biology and adaptation have used these technologies to sequence mRNA transcriptomes of specific tissues, which in turn are often compared to other tissues, or other individuals with different phenotypes. While these techniques are extremely powerful, careful attention to data quality is required. In particular, because high-throughput sequencing is more error-prone than traditional Sanger sequencing, quality trimming of sequence reads should be an important step in all data processing pipelines. While several software packages for quality trimming exist, no general guidelines for the specifics of trimming have been developed. Here, using empirically derived sequence data, I provide general recommendations regarding the optimal strength of trimming, specifically in mRNA-Seq studies. Although very aggressive quality trimming is common, this study suggests that a more gentle trimming, specifically of those nucleotides whose Phred score <2 or <5, is optimal for most studies across a wide variety of metrics.

Signatures of rapid evolution in urban and rural transcriptomes of white-footed mice (Peromyscus leucopus) in the New York metropolitan area

Urbanization is a major cause of ecological degradation around the world, and human settlement in large cities is accelerating. New York City (NYC) is one of the oldest and most urbanized cities in North America, but still maintains 20% vegetation cover and substantial populations of some native wildlife. The white-footed mouse, Peromyscusleucopus, is a common resident of NYC's forest fragments and an emerging model system for examining the evolutionary consequences of urbanization. In this study, we developed transcriptomic resources for urban P. leucopus to examine evolutionary changes in protein-coding regions for an exemplar "urban adapter." We used Roche 454 GS FLX+ high throughput sequencing to derive transcriptomes from multiple tissues from individuals across both urban and rural populations. From these data, we identified 31,015 SNPs and several candidate genes potentially experiencing positive selection in urban populations of P. leucopus. These candidate genes are involved in xenobiotic metabolism, innate immune response, demethylation activity, and other important biological phenomena in novel urban environments. This study is one of the first to report candidate genes exhibiting signatures of directional selection in divergent urban ecosystems.

Improving transcriptome assembly through error correction of high-throughput sequence reads

The study of functional genomics, particularly in non-model organisms, has been dramatically improved over the last few years by the use of transcriptomes and RNAseq. While these studies are potentially extremely powerful, a computationally intensive procedure, the de novo construction of a reference transcriptome must be completed as a prerequisite to further analyses. The accurate reference is critically important as all downstream steps, including estimating transcript abundance are critically dependent on the construction of an accurate reference. Though a substantial amount of research has been done on assembly, only recently have the pre-assembly procedures been studied in detail. Specifically, several stand-alone error correction modules have been reported on and, while they have shown to be effective in reducing errors at the level of sequencing reads, how error correction impacts assembly accuracy is largely unknown. Here, we show via use of a simulated and empiric dataset, that applying error correction to sequencing reads has significant positive effects on assembly accuracy, and should be applied to all datasets. A complete collection of commands which will allow for the production of Reptile corrected reads is available at https://github.com/macmanes/error_correction/tree/master/scripts and as File S1.