Benefits and Limits of Phasing Alleles for Network Inference of Allopolyploid Complexes.

Accurately reconstructing the reticulate histories of polyploids remains a central challenge for understanding plant evolution. Although phylogenetic networks can provide insights into relationships among polyploid lineages, inferring networks may be hindered by the complexities of homology determination in polyploid taxa. We use simulations to show that phasing alleles from allopolyploid individuals can improve phylogenetic network inference under the multispecies coalescent by obtaining the true network with fewer loci compared to haplotype consensus sequences or sequences with heterozygous bases represented as ambiguity codes. Phased allelic data can also improve divergence time estimates for networks, which is helpful for evaluating allopolyploid speciation hypotheses and proposing mechanisms of speciation. To achieve these outcomes in empirical data, we present a novel pipeline that leverages a recently developed phasing algorithm to reliably phase alleles from polyploids. This pipeline is especially appropriate for target enrichment data, where depth of coverage is typically high enough to phase entire loci. We provide an empirical example in the North American Dryopteris fern complex that demonstrates insights from phased data as well as the challenges of network inference. We establish that our pipeline (PATÉ: Phased Alleles from Target Enrichment data) is capable of recovering a high proportion of phased loci from both diploids and polyploids. These data may improve network estimates compared to using haplotype consensus assemblies by accurately inferring the direction of gene flow, but statistical non-identifiability of phylogenetic networks poses a barrier to inferring the evolutionary history of reticulate complexes.

Interspecies Papillomavirus Type Infection and a Novel Papillomavirus Type in Red Ruffed Lemurs (Varecia rubra).

The Papillomaviridae are a family of vertebrate-infecting viruses of oncogenic potential generally thought to be host species- and tissue-specific. Despite their phylogenetic relatedness to humans, there is a scarcity of data on papillomaviruses (PVs) in speciose non-human primate lineages, particularly the lemuriform primates. Varecia variegata (black-and-white ruffed lemurs) and Varecia rubra (red ruffed lemurs), two closely related species comprising the Varecia genus, are critically endangered with large global captive populations. Varecia variegata papillomavirus (VavPV) types -1 and -2, the first PVs in lemurs with a fully identified genome, were previously characterized from captive V. variegata saliva. To build upon this discovery, saliva samples were collected from captive V. rubra with the following aims: (1) to identify PVs shared between V. variegata and V. rubra and (2) to characterize novel PVs in V. rubra to better understand PV diversity in the lemuriform primates. Three complete PV genomes were determined from V. rubra samples. Two of these PV genomes share 98% L1 nucleotide identity with VavPV2, denoting interspecies infection of V. rubra by VavPV2. This work represents the first reported case of interspecies PV infection amongst the strepsirrhine primates. The third PV genome shares <68% L1 nucleotide identity with that of all PVs. Thus, it represents a new PV species and has been named Varecia rubra papillomavirus 1 (VarPV1). VavPV1, VavPV2, and VarPV1 form a new clade within the Papillomaviridae family, likely representing a novel genus. Future work diversifying sample collection (i.e., lemur host species from multiple genera, sample type, geographic location, and wild populations) is likely to uncover a world of diverse lemur PVs.

Identification of diverse papillomaviruses in captive black-and-white ruffed lemurs (Varecia variegata).

Papillomaviruses (PVs) are host-species-specific and tissue-specific viruses that infect a diverse array of vertebrate hosts, including humans and non-human primates, with varying pathogenic outcomes. Although primate PVs have been studied extensively, no complete genome sequences of PVs from lemurs have been determined to date. Saliva samples from three critically endangered, captive black-and-white ruffed lemurs (Varecia variegata variegata) at the Duke Lemur Center (USA) were analyzed, using high-throughput sequencing, for the presence of oral papillomaviruses. We identified three PVs from two individuals, one of which had a coinfection with two different PVs. Two of the three PVs share 99.6% nucleotide sequence identity, and we have named these isolates "Varecia variegata papillomavirus 1" (VavPV1). The third PV shares ~63% nucleotide sequence identity with VavPV1, and thus, we have named it "Varecia variegata papillomavirus 2" (VavPV2). Based on their E1 + E2 + L1 protein sequence phylogeny, the VavPVs form a distinct clade. This clade likely represents a novel genus, with VavPV1 and VavPV2 belonging to two distinct species. Our findings represent the first complete genome sequences of PVs found in lemuriform primates, with their presence suggesting the potential existence of diverse PVs across the over 100 species of lemurs.

Applications of 3D printing in small animal magnetic resonance imaging.

Three-dimensional (3D) printing has significantly impacted the quality, efficiency, and reproducibility of preclinical magnetic resonance imaging. It has vastly expanded the ability to produce MR-compatible parts that readily permit customization of animal handling, achieve consistent positioning of anatomy and RF coils promptly, and accelerate throughput. It permits the rapid and cost-effective creation of parts customized to a specific imaging study, animal species, animal weight, or even one unique animal, not routinely used in preclinical research. We illustrate the power of this technology by describing five preclinical studies and specific solutions enabled by different 3D printing processes and materials. We describe fixtures, assemblies, and devices that were created to ensure the safety of anesthetized lemurs during an MR examination of their brain or to facilitate localized, contrast-enhanced measurements of white blood cell concentration in a mouse model of pancreatitis. We illustrate expansive use of 3D printing to build a customized birdcage coil and components of a ventilator to enable imaging of pulmonary gas exchange in rats using hyperpolarized Xe 129 . Finally, we present applications of 3D printing to create high-quality, dual RF coils to accelerate brain connectivity mapping in mouse brain specimens and to increase the throughput of brain tumor examinations in a mouse model of pituitary adenoma.

Species discovery and validation in a cryptic radiation of endangered primates: coalescent-based species delimitation in Madagascar's mouse lemurs.

Implementation of the coalescent model in a Bayesian framework is an emerging strength in genetically based species delimitation studies. By providing an objective measure of species diagnosis, these methods represent a quantitative enhancement to the analysis of multilocus data, and complement more traditional methods based on phenotypic and ecological characteristics. Recognized as two species 20 years ago, mouse lemurs (genus Microcebus) now comprise more than 20 species, largely diagnosed from mtDNA sequence data. With each new species description, enthusiasm has been tempered with scientific scepticism. Here, we present a statistically justified and unbiased Bayesian approach towards mouse lemur species delimitation. We perform validation tests using multilocus sequence data and two methodologies: (i) reverse-jump Markov chain Monte Carlo sampling to assess the likelihood of different models defined a priori by a guide tree, and (ii) a Bayes factor delimitation test that compares different species-tree models without a guide tree. We assess the sensitivity of these methods using randomized individual assignments, which has been used in bpp studies, but not with Bayes factor delimitation tests. Our results validate previously diagnosed taxa, as well as new species hypotheses, resulting in support for three new mouse lemur species. As the challenge of multiple researchers using differing criteria to describe diversity is not unique to Microcebus, the methods used here have significant potential for clarifying diversity in other taxonomic groups. We echo previous studies in advocating that multiple lines of evidence, including use of the coalescent model, should be trusted to delimit new species.

Evolution of a tumorigenic property conferred by glycophosphatidyl-inositol membrane anchors of carcinoembryonic antigen gene family members during the primate radiation.

GPI membrane anchors of cell surface glycoproteins have been shown to confer functional properties that are different from their transmembrane (TM)-anchored counterparts. For the human carcinoembryonic antigen (CEA) family, a subfamily of the immunoglobulin superfamily, conversion of the mode of membrane linkage from TM to GPI confers radical changes in function: from tumor suppression or neutrality toward inhibition of differentiation and anoikis and distortion of tissue architecture, thereby contributing to tumorigenesis. We show here that GPI anchorage in the CEA family evolved twice independently in primates, very likely from more primitive TM anchors, by different packages of mutations. Both mutational packages, one package found in many primates, including humans, and a second, novel package found only in the Cebidae radiation of New World monkeys, give rise to efficiently processed GPI-linked proteins. Both types of GPI anchors mediate inhibition of cell differentiation. The estimated rate of nonsynonymous mutations (Ka) in the anchor-determining domain for conversion from TM to GPI anchorage in the CEA family that were fixed during evolution in these primates is 7 times higher than the average Ka in primates, indicating positive selection. These results suggest therefore that the functional changes mediated by CEA GPI anchors, including the inhibition of differentiation and anoikis, could be adaptive and advantageous.

Comparison of likelihood and Bayesian methods for estimating divergence times using multiple gene Loci and calibration points, with application to a radiation of cute-looking mouse lemur species.

Divergence time and substitution rate are seriously confounded in phylogenetic analysis, making it difficult to estimate divergence times when the molecular clock (rate constancy among lineages) is violated. This problem can be alleviated to some extent by analyzing multiple gene loci simultaneously and by using multiple calibration points. While different genes may have different patterns of evolutionary rate change, they share the same divergence times. Indeed, the fact that each gene may violate the molecular clock differently leads to the advantage of simultaneous analysis of multiple loci. Multiple calibration points provide the means for characterizing the local evolutionary rates on the phylogeny. In this paper, we extend previous likelihood models of local molecular clock for estimating species divergence times to accommodate multiple calibration points and multiple genes. Heterogeneity among different genes in evolutionary rate and in substitution process is accounted for by the models. We apply the likelihood models to analyze two mitochondrial protein-coding genes, cytochrome oxidase II and cytochrome b, to estimate divergence times of Malagasy mouse lemurs and related outgroups. The likelihood method is compared with the Bayes method of Thorne et al. (1998, Mol. Biol. Evol. 15:1647-1657), which uses a probabilistic model to describe the change in evolutionary rate over time and uses the Markov chain Monte Carlo procedure to derive the posterior distribution of rates and times. Our likelihood implementation has the drawbacks of failing to accommodate uncertainties in fossil calibrations and of requiring the researcher to classify branches on the tree into different rate groups. Both problems are avoided in the Bayes method. Despite the differences in the two methods, however, data partitions and model assumptions had the greatest impact on date estimation. The three codon positions have very different substitution rates and evolutionary dynamics, and assumptions in the substitution model affect date estimation in both likelihood and Bayes analyses. The results demonstrate that the separate analysis is unreliable, with dates variable among codon positions and between methods, and that the combined analysis is much more reliable. When the three codon positions were analyzed simultaneously under the most realistic models using all available calibration information, the two methods produced similar results. The divergence of the mouse lemurs is dated to be around 7-10 million years ago, indicating a surprisingly early species radiation for such a morphologically uniform group of primates.