Matthew D. MacManes

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.

Comparative genomics explains the evolutionary success of reef-forming corals

Transcriptome and genome data from twenty stony coral species and a selection of reference bilaterians were studied to elucidate coral evolutionary history. We identified genes that encode the proteins responsible for the precipitation and aggregation of the aragonite skeleton on which the organisms live, and revealed a network of environmental sensors that coordinate responses of the host animals to temperature, light, and pH. Furthermore, we describe a variety of stress-related pathways, including apoptotic pathways that allow the host animals to detoxify reactive oxygen and nitrogen species that are generated by their intracellular photosynthetic symbionts, and determine the fate of corals under environmental stress. Some of these genes arose through horizontal gene transfer and comprise at least 0.2% of the animal gene inventory. Our analysis elucidates the evolutionary strategies that have allowed symbiotic corals to adapt and thrive for hundreds of millions of years. DOI: http://dx.doi.org/10.7554/eLife.13288.001

Establishing evidenced-based best practice for the de novo assembly and evaluation of transcriptomes from non-model organisms

Characterizing transcriptomes in both model and non-model organisms has resulted in a massive increase in our understanding of biological phenomena. This boon, largely made possible via high-throughput sequencing, means that studies of functional, evolutionary and population genomics are now being done by hundreds or even thousands of labs around the world. For many, these studies begin with a de novo transcriptome assembly, which is a technically complicated process involving several discrete steps. Each step may be accomplished in one of several different ways, using different software packages, each producing different results. This analytical complexity begs the question – Which method(s) are optimal? Using reference and non-reference based evaluative methods, I propose a set of guidelines that aim to standardize and facilitate the process of transcriptome assembly. These recommendations include the generation of between 20 million and 40 million sequencing reads from single individual where possible, error correction of reads, gentle quality trimming, assembly filtering using Transrate and/or gene expression, annotation using dammit, and appropriate reporting. These recommendations have been extensively benchmarked and applied to publicly available transcriptomes, resulting in improvements in both content and contiguity. To facilitate the implementation of the proposed standardized methods, I have released a set of version controlled open-sourced code, The Oyster River Protocol for Transcriptome Assembly, available at http://oyster-river-protocol.rtfd.org/.Characterizing transcriptomes in both model and non-model organisms has resulted in a massive increase in our understanding of biological phenomena. This boon, largely made possible via high-throughput sequencing, means that studies of functional, evolutionary and population genomics are now being done by hundreds or even thousands of labs around the world. For many, these studies begin with a de novo transcriptome assembly, which is a technically complicated process involving several discrete steps. Each step may be accomplished in one of several different ways, using different software packages, each producing different results. This analytical complexity begs the question -- Which method(s) are optimal? Using reference and non-reference based evaluative methods, I propose a set of guidelines that aim to standardize and facilitate the process of transcriptome assembly. These recommendations include the generation of between 20 million and 40 million sequencing reads from single individual where possible, error correction of reads, gentle quality trimming, assembly filtering using Transrate and/or gene expression, annotation using dammit, and appropriate reporting. These recommendations have been extensively benchmarked and applied to publicly available transcriptomes, resulting in improvements in both content and contiguity. To facilitate the implementation of the proposed standardized methods, I have released a set of version controlled open-sourced code, The Oyster River Protocol for Transcriptome Assembly, available at http://oyster-river-protocol.rtfd.org/.

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.

Physiological and biochemical changes associated with acute experimental dehydration in the desert adapted mouse, Peromyscus eremicus

Characterizing traits critical for adaptation to a given environment is an important first step in understanding how phenotypes evolve. How animals adapt to the extreme heat and aridity commonplace to deserts is an exceptionally interesting example of these processes, and has been the focus of study for decades. In contrast to those studies, where experiments are conducted on either wild animals or captive animals held in non‐desert conditions, the study described here leverages a unique environmental chamber that replicates desert conditions for captive Peromyscus eremicus (cactus mouse). Here, we establish baseline values for daily water intake and for serum electrolytes, as well as the response of these variables to acute experimental dehydration. In brief, P. eremicus daily water intake is very low. Its serum electrolytes are distinct from many previously studied animals, and its response to acute dehydration is profound, though not suggestive of renal impairment, which is atypical of mammals.

Severe acute dehydration in a desert rodent elicits a transcriptional response that effectively prevents kidney injury

Animals living in desert environments are forced to survive despite severe heat, intense solar radiation, and both acute and chronic dehydration. These animals have evolved phenotypes that effectively address these environmental stressors. To begin to understand the ways in which the desert-adapted rodent Peromyscus eremicus survives, reproductively mature adults were subjected to 72 h of water deprivation, during which they lost, on average, 23% of their body weight. The animals reacted via a series of changes in the kidney, which included modulating expression of genes responsible for reducing the rate of transcription and maintaining water and salt balance. Extracellular matrix turnover appeared to be decreased, and apoptosis was limited. In contrast to the canonical human response, serum creatinine and other biomarkers of kidney injury were not elevated, suggesting that changes in gene expression related to acute dehydration may effectively prohibit widespread kidney damage in the cactus mouse.

Widespread patterns of sexually dimorphic gene expression in an avian hypothalamic–pituitary–gonadal (HPG) axis

The hypothalamic-pituitary-gonadal (HPG) axis is a key biological system required for reproduction and associated sexual behaviors to occur. In the avian reproductive model of the rock dove (Columba livia), we characterized the transcript community of each tissue of the HPG axis in both sexes, thereby significantly expanding our mechanistic insight into HPG activity. We report greater sex-biased differential expression in the pituitary as compared to the hypothalamus, with multiple genes more highly expressed in the male pituitary being related to secretory function, and multiple genes more highly expressed in the female pituitary being related to reproduction, growth, and development. We report tissue-specific and sex-biased expression in genes commonly investigated when studying reproduction, highlighting the need for sex parity in future studies. In addition, we uncover new targets of investigation in both sexes, which could potentially change our understanding of HPG function.

PALADIN: protein alignment for functional profiling whole metagenome shotgun data

Motivation: Whole metagenome shotgun sequencing is a powerful approach for assaying the functional potential of microbial communities. We currently lack tools that efficiently and accurately align DNA reads against protein references, the technique necessary for constructing a functional profile. Here, we present PALADIN—a novel modification of the Burrows‐Wheeler Aligner that provides accurate alignment, robust reporting capabilities and orders‐of‐magnitude improved efficiency by directly mapping in protein space. Results: We compared the accuracy and efficiency of PALADIN against existing tools that employ nucleotide or protein alignment algorithms. Using simulated reads, PALADIN consistently outperformed the popular DNA read mappers BWA and NovoAlign in detected proteins, percentage of reads mapped and ontological similarity. We also compared PALADIN against four existing protein alignment tools: BLASTX, RAPSearch2, DIAMOND and Lambda, using empirically obtained reads. PALADIN yielded results seven times faster than the best performing alternative, DIAMOND and nearly 8000 times faster than BLASTX. PALADINs accuracy was comparable to all tested solutions. Availability and Implementation: PALADIN was implemented in C, and its source code and documentation are available at https://github.com/twestbrookunh/paladin Contact: anthonyw@wildcats.unh.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Characterization of a male reproductive transcriptome for Peromyscus eremicus (Cactus mouse)

Rodents of the genus Peromyscus have become increasingly utilized models for investigations into adaptive biology. This genus is particularly powerful for research linking genetics with adaptive physiology or behaviors, and recent research has capitalized on the unique opportunities afforded by the ecological diversity of these rodents. Well characterized genomic and transcriptomic data is intrinsic to explorations of the genetic architecture responsible for ecological adaptations. Therefore, this study characterizes the transcriptome of three male reproductive tissues (testes, epididymis and vas deferens) of Peromyscus eremicus (Cactus mouse), a desert specialist. The transcriptome assembly process was optimized in order to produce a high quality and substantially complete annotated transcriptome. This composite transcriptome was generated to characterize the expressed transcripts in the male reproductive tract of P. eremicus, which will serve as a crucial resource for future research investigating our hypothesis that the male Cactus mouse possesses an adaptive reproductive phenotype to mitigate water-loss from ejaculate. This study reports genes under positive selection in the male Cactus mouse reproductive transcriptome relative to transcriptomes from Peromyscus maniculatus (deer mouse) and Mus musculus. Thus, this study expands upon existing genetic research in this species, and we provide a high quality transcriptome to enable further explorations of our proposed hypothesis for male Cactus mouse reproductive adaptations to minimize seminal fluid loss.

Brain Transcriptional Profiles of Male Alternative Reproductive Tactics and Females in Bluegill Sunfish

Bluegill sunfish (Lepomis macrochirus) are one of the classic systems for studying male alternative reproductive tactics (ARTs) in teleost fishes. In this species, there are two distinct life histories: parental and cuckolder, encompassing three reproductive tactics, parental, satellite, and sneaker. The parental life history is fixed, whereas individuals who enter the cuckolder life history transition from sneaker to satellite tactic as they grow. For this study, we used RNAseq to characterize the brain transcriptome of the three male tactics and females during spawning to identify gene ontology (GO) categories and potential candidate genes associated with each tactic. We found that sneaker males had higher levels of gene expression differentiation compared to the other two male tactics. Sneaker males also had higher expression in ionotropic glutamate receptor genes, specifically AMPA receptors, compared to other males, which may be important for increased spatial working memory while attempting to cuckold parental males at their nests. Larger differences in gene expression also occurred among male tactics than between males and females. We found significant expression differences in several candidate genes that were previously identified in other species with ARTs and suggest a previously undescribed role for cAMP-responsive element modulator (crem) in influencing parental male behaviors during spawning.