Food deprivation reduces sensitivity of liver Igf1 synthesis pathways to growth hormone in juvenile gopher rockfish (Sebastes carnatus).

Growth hormone (Gh) regulates growth in part by stimulating the liver to synthesize and release insulin-like growth factor-1 (Igf1), which then promotes somatic growth. However, for fish experiencing food limitation, elevated blood Gh can occur even with low circulating Igf1 and slow growth, suggesting that nutritional stress can alter the sensitivity of liver Igf1 synthesis pathways to Gh. Here, we examined how recent feeding experience affected Gh regulation of liver Igf1 synthesis pathways in juvenile gopher rockfish (Sebastes carnatus) to illuminate mechanisms underlying the nutritional modulation of Igf1 production. Juvenile gopher rockfish were maintained under conditions of feeding or complete food deprivation (fasting) for 14 d and then treated with recombinant sea bream (Sparus aurata) Gh or saline control. Gh upregulated hepatic igf1 mRNA levels in fed fish but not in fasted fish. The liver of fasted rockfish also showed a lower relative abundance of gene transcripts encoding teleost Gh receptors 1 (ghr1) and 2 (ghr2), as well as reduced protein levels of phosphorylated janus tyrosine kinase 2 (pJak2) and signal transducer and activator of transcription 5 (pStat5), which function to induce igf1 gene transcription following Gh binding to Gh receptors. Relative hepatic mRNA levels for suppressors of cytokine signaling (Socs) genes socs2, socs3a, and socs3b were also lower in fasted rockfish. Socs2 can suppress Gh activation of Jak2/Stat5, and fasting-related variation in socs expression may reflect modulated inhibitory control of igf1 gene transcription. Fasted rockfish also had elevated liver mRNA abundances for lipolytic hormone-sensitive lipase 1 (hsl1) and Igf binding proteins igfbp1a, -1b and -3a, reduced liver mRNAs encoding igfbp2b and an Igfbp acid labile subunit-like (igfals) gene, and higher transcript abundances for Igf1 receptors igf1ra and igf1rb in skeletal muscle. Together, these findings suggest that food deprivation impacts liver Igf1 responsiveness to Gh via multiple mechanisms that include a downregulation of hepatic Gh receptors, modulation of the intracellular Jak2/Stat5 transduction pathway, and possible shifts in Socs-inhibitory control of igf1 gene transcription, while also demonstrating that these changes occur in concert with shifts in liver Igfbp expression and muscle Gh/Igf1 signaling pathway components.

Mitochondrial genome of Geomydoecus aurei, a pocket-gopher louse.

Parasitic lice demonstrate an unusual array of mitochondrial genome architectures and gene arrangements. We characterized the mitochondrial genome of Geomydoecus aurei, a chewing louse (Phthiraptera: Trichodectidae) found on pocket gophers (Rodentia: Geomyidae) using reads from both Illumina and Oxford Nanopore sequencing coupled with PCR, cloning, and Sanger sequencing to verify structure and arrangement for each chromosome. The genome consisted of 12 circular mitochondrial chromosomes ranging in size from 1,318 to 2,088 nucleotides (nt). Total genome size was 19,015 nt. All 37 genes typical of metazoans (2 rRNA genes, 22 tRNA genes, and 13 protein-coding genes) were present. An average of 26% of each chromosome was composed of non-gene sequences. Within the non-gene region of each chromosome, there was a 79-nt nucleotide sequence that was identical among chromosomes and a conserved sequence with secondary structure that was always followed by a poly-T region. We hypothesize that these regions may be important in the initiation of transcription and DNA replication, respectively. The G. aurei genome shares 8 derived gene clusters with other chewing lice of mammals, but in G. aurei, genes on several chromosomes are not contiguous.

Locomotory adaptations in entoptychine gophers (Rodentia: Geomyidae) and the mosaic evolution of fossoriality.

Pocket gophers (family Geomyidae) are the dominant burrowing rodents in North America today. Their fossil record is also incredibly rich; in particular, entoptychine gophers, a diverse extinct subfamily of the Geomyidae, are known from countless teeth and jaws from Oligocene and Miocene-aged deposits of the western United States and Mexico. Their postcranial remains, however, are much rarer and little studied. Yet, they offer the opportunity to investigate the locomotion of fossil gophers, shed light on the evolution of fossoriality, and enable ecomorphological comparisons with contemporaneous rodents. We present herein a quantitative study of the cranial and postcranial remains of eight different species of entoptychine gophers as well as many contemporary rodent species. We find a range of burrowing capabilities within Entoptychinae, including semifossorial scratch-digging animals and fossorial taxa with cranial adaptations to burrowing. Our results suggest the repeated evolution of chisel-tooth digging across genera. Comparisons between entoptychine gophers and contemporaneous rodent taxa show little ecomorphological overlap and suggest that the succession of burrowing rodent taxa on the landscape may have had more to do with habitat partitioning than competition.

Getting a head in hard soils: Convergent skull evolution and divergent allometric patterns explain shape variation in a highly diverse genus of pocket gophers (Thomomys).

BACKGROUND: High morphological diversity can occur in closely related animals when selection favors morphologies that are subject to intrinsic biological constraints. A good example is subterranean rodents of the genus Thomomys, one of the most taxonomically and morphologically diverse mammalian genera. Highly procumbent, tooth-digging rodent skull shapes are often geometric consequences of increased body size. Indeed, larger-bodied Thomomys species tend to inhabit harder soils. We used geometric morphometric analyses to investigate the interplay between soil hardness (the main extrinsic selection pressure on fossorial mammals) and allometry (i.e. shape change due to size change; generally considered the main intrinsic factor) on crania and humeri in this fast-evolving mammalian clade. RESULTS: Larger Thomomys species/subspecies tend to have more procumbent cranial shapes with some exceptions, including a small-bodied species inhabiting hard soils. Counter to earlier suggestions, cranial shape within Thomomys does not follow a genus-wide allometric pattern as even regional subpopulations differ in allometric slopes. In contrast, humeral shape varies less with body size and with soil hardness. Soft-soil taxa have larger humeral muscle attachment sites but retain an orthodont (non-procumbent) cranial morphology. In intermediate soils, two pairs of sister taxa diverge through differential modifications on either the humerus or the cranium. In the hardest soils, both humeral and cranial morphology are derived through large muscle attachment sites and a high degree of procumbency. CONCLUSIONS: Our results show that conflict between morphological function and intrinsic allometric patterning can quickly and differentially alter the rodent skeleton, especially the skull. In addition, we found a new case of convergent evolution of incisor procumbency among large-, medium-, and small-sized species inhabiting hard soils. This occurs through different combinations of allometric and non-allometric changes, contributing to shape diversity within the genus. The strong influence of allometry on cranial shape appears to confirm suggestions that developmental change underlies mammalian cranial shape divergences, but this requires confirmation from ontogenetic studies. Our findings illustrate how a variety of intrinsic processes, resulting in species-level convergence, could sustain a genus-level range across a variety of extrinsic environments. This might represent a mechanism for observations of genus-level niche conservation despite species extinctions in mammals.

Cophylogenetic Reconciliation with ILP.

In this paper, we present an integer linear programming (ILP) approach, called CoRe-ILP, for finding an optimal time consistent cophylogenetic host-parasite reconciliation under the cophylogenetic event model with the events cospeciation, duplication, sorting, host switch, and failure to diverge. Instead of assuming event costs, a simplified model is used, maximizing primarily for cospeciations and secondarily minimizing host switching events. Duplications, sortings, and failure to diverge events are not explicitly scored. Different from existing event based reconciliation methods, CoRe-ILP can use (approximate) phylogenetic branch lengths for filtering possible ancestral host-parasite interactions. Experimentally, it is shown that CoRe-ILP can successfully use branch length information and performs well for biological and simulated data sets. The results of CoRe-ILP are compared with the results of the reconciliation tools Jane 4, Treemap 3b, NOTUNG 2.8 Beta, and Ranger-DTL. Algorithm CoRe-ILP is implemented using IBM ILOG CPLEX Optimizer 12.6 and is freely available from http://pacosy.informatik.uni-leipzig.de/core-ilp.

DISSECT: an assignment-free Bayesian discovery method for species delimitation under the multispecies coalescent.

MOTIVATION: The multispecies coalescent model provides a formal framework for the assignment of individual organisms to species, where the species are modeled as the branches of the sp tree. None of the available approaches so far have simultaneously co-estimated all the relevant parameters in the model, without restricting the parameter space by requiring a guide tree and/or prior assignment of individuals to clusters or species. RESULTS: We present DISSECT, which explores the full space of possible clusterings of individuals and species tree topologies in a Bayesian framework. It uses an approximation to avoid the need for reversible-jump Markov Chain Monte Carlo, in the form of a prior that is a modification of the birth-death prior for the species tree. It incorporates a spike near zero in the density for node heights. The model has two extra parameters: one controls the degree of approximation and the second controls the prior distribution on the numbers of species. It is implemented as part of BEAST and requires only a few changes from a standard *BEAST analysis. The method is evaluated on simulated data and demonstrated on an empirical dataset. The method is shown to be insensitive to the degree of approximation, but quite sensitive to the second parameter, suggesting that large numbers of sequences are needed to draw firm conclusions. AVAILABILITY AND IMPLEMENTATION: http://tree.bio.ed.ac.uk/software/beast/, http://www.indriid.com/dissectinbeast.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cophylogeny on a fine scale: Geomydoecus chewing lice and their pocket gopher hosts, Pappogeomys bulleri.

Many species of pocket gophers and their ectoparasitic chewing lice have broadly congruent phylogenies, indicating a history of frequent codivergence. For a variety of reasons, phylogenies of codiverging hosts and parasites are expected to be less congruent for more recently diverged taxa. This study is the first of its scale in the pocket gopher and chewing louse system, with its focus entirely on comparisons among populations within a single species of host and 3 chewing louse species in the Geomydoecus bulleri species complex. We examined mitochondrial DNA from a total of 46 specimens of Geomydoecus lice collected from 11 populations of the pocket gopher host, Pappogeomys bulleri. We also examined nuclear DNA from a subset of these chewing lice. Louse phylogenies were compared with a published pocket gopher phylogeny. Contrary to expectations, we observed a statistically significant degree of parallel cladogenesis in these closely related hosts and their parasites. We also observed a higher rate of evolution in chewing louse lineages than in their corresponding pocket gopher hosts. In addition, we found that 1 louse species (Geomydoecus burti) may not be a valid species, that subspecies within G. bulleri are not reciprocally monophyletic, and that morphological and genetic evidence support recognition of a new species of louse, Geomydoecus pricei.

Phylogenetic structure of the Thomomys bottae-umbrinus complex in North America.

The phylogeography of the Thomomys bottae-umbrinus complex in the United States and Mexico was assessed with sequences of the mitochondrial cytochrome b gene. These sequences were obtained from 225 individuals representing 108 locations over the range, including 56 sequences from GenBank. 110 (500bp) sequences were used for Bayesian inference and neighbor-joining analyses, and 34 (1140bp) specimens from the main clades obtained from the Bayesian inference were used in maximum-parsimony and maximum-likelihood analyses. The different analyses indicate significant variation within the species complex that averages 13% among major groups of genetic differences among Thomomys bottae-umbrinus. The overall pattern of geographic variation is not concordant with the current taxonomy. To the contrary, eight monophyletic groups are supported by all analyses and can be considered phylogenetic species. Overall divergence among these groups appears influenced by historical biogeographic events active during the Pliocene and Pleistocene.

A novel test for host-symbiont codivergence indicates ancient origin of fungal endophytes in grasses.

Significant phylogenetic codivergence between plant or animal hosts (H) and their symbionts or parasites (P) indicates the importance of their interactions on evolutionary time scales. However, valid and realistic methods to test for codivergence are not fully developed. One of the systems where possible codivergence has been of interest involves the large subfamily of temperate grasses (Pooideae) and their endophytic fungi (epichloae). These widespread symbioses often help protect host plants from herbivory and stresses and affect species diversity and food web structures. Here we introduce the MRCALink (most-recent-common-ancestor link) method and use it to investigate the possibility of grass-epichloë codivergence. MRCALink applied to ultrametric H and P trees identifies all corresponding nodes for pairwise comparisons of MRCA ages. The result is compared to the space of random H and P tree pairs estimated by a Monte Carlo method. Compared to tree reconciliation, the method is less dependent on tree topologies (which often can be misleading), and it crucially improves on phylogeny-independent methods such as ParaFit or the Mantel test by eliminating an extreme (but previously unrecognized) distortion of node-pair sampling. Analysis of 26 grass species-epichloë species symbioses did not reject random association of H and P MRCA ages. However, when five obvious host jumps were removed, the analysis significantly rejected random association and supported grass-endophyte codivergence. Interestingly, early cladogenesis events in the Pooideae corresponded to early cladogenesis events in epichloae, suggesting concomitant origins of this grass subfamily and its remarkable group of symbionts. We also applied our method to the well-known gopher-louse data set.

Multilocus phylogenetics of a rapid radiation in the genus Thomomys (Rodentia: Geomyidae).

Species complexes undergoing rapid radiation present a challenge in molecular systematics because of the possibility that ancestral polymorphism is retained in component gene trees. Coalescent theory has demonstrated that gene trees often fail to match lineage trees when taxon divergence times are less than the ancestral effective population sizes. Suggestions to increase the number of loci and the number of individuals per taxon have been proposed; however, phylogenetic methods to adequately analyze these data in a coalescent framework are scarce. We compare two approaches to estimating lineage (species) trees using multiple individuals and multiple loci: the commonly used partitioned Bayesian analysis of concatenated sequences and a modification of a newly developed hierarchical Bayesian method (BEST) that simultaneously estimates gene trees and species trees from multilocus data. We test these approaches on a phylogeny of rapidly radiating species wherein divergence times are likely to be smaller than effective population sizes, and incomplete lineage sorting is known, in the rodent genus, Thomomys. We use seven independent noncoding nuclear sequence loci (total approximately 4300 bp) and between 1 and 12 individuals per taxon to construct a phylogenetic hypothesis for eight Thomomys species. The majority-rule consensus tree from the partitioned concatenated analysis included 14 strongly supported bipartitions, corroborating monophyletic species status of five of the eight named species. The BEST tree strongly supported only the split between the two subgenera and showed very low support for any other clade. Comparison of both lineage trees to individual gene trees revealed that the concatenation method appears to ignore conflicting signals among gene trees, whereas the BEST tree considers conflicting signals and downweights support for those nodes. Bayes factor analysis of posterior tree distributions from both analyses strongly favor the model underlying the BEST analysis. This comparison underscores the risks of overreliance on results from concatenation, and ignoring the properties of coalescence, especially in cases of recent, rapid radiations.

Cophylogeny and disparate rates of evolution in sympatric lineages of chewing lice on pocket gophers.

Although molecular-based phylogenetic studies of hosts and parasites are increasingly common in the literature, no study to date has examined two congeneric lineages of parasites that live in sympatry on the same lineage of hosts. This study examines phylogenetic relationships among chewing lice (Phthiraptera: Trichodectidae) of the Geomydoecus coronadoi and Geomydoecus mexicanus species complexes and compares these to phylogenetic patterns in their hosts (pocket gophers of the rodent family Geomyidae). Sympatry of congeneric lice provides a natural experiment to test the hypothesis that closely related lineages of parasites will respond similarly to the same host. Sequence data from the mitochondrial COI and the nuclear EF-1alpha genes confirm that the two louse complexes are reciprocally monophyletic and that individual clades within each species complex parasitize a different species of pocket gopher. Phylogenetic comparisons reveal that both louse complexes show a significant pattern of cophylogeny with their hosts. Comparisons of rates of nucleotide substitution at 4-fold degenerate sites in the COI gene indicate that both groups of lice have significantly higher basal mutation rates than their hosts. The two groups of lice have similar basal rates of mutation, but lice of the G. coronadoi complex show significantly elevated rates of nucleotide substitution at all sites. These rate differences are hypothesized to result from population-level phenomena, such as effective population size, founder effects, and drift, that influence rates of nucleotide substitution.

The genetics of adaptive coat color in gophers: coding variation at Mc1r is not responsible for dorsal color differences.

The genetics of adaptation is a key problem in evolutionary biology. Pocket gophers of the species Thomomys bottae provide one of the most striking examples of coat color variation in mammals. Dorsal pelage color is strongly correlated with soil color across the range of the species, presumably reflecting the selective pressure exerted by predation. To investigate the genetic basis of coat color variation in T. bottae, we cloned and sequenced the melanocortin-1 receptor locus (Mc1r), a candidate pigmentation gene, in 5 dark and 5 light populations of the species. Our results show that, in contrast to many other species of mammals and other vertebrates, coding variation at Mc1r is not the main determinant of coat color variation in T. bottae. These results demonstrate that similar phenotypic variation may have a different genetic basis among different mammalian species.