Serratia marcescens Outbreak at a Correctional Facility: Environmental Sampling, Laboratory Analyses and Genomic Characterization to Assess Sources and Persistence.

Serratia marcescens is an environmental bacterium and clinical pathogen that can cause an array of infections. We describe an environmental sampling and comparative genomics approach used to investigate a multi-year outbreak of S. marcescens at a correctional facility. Whole genome sequencing analysis revealed a predominant cluster of clonally related S. marcescens from nine patient cases and items associated with illicit drug use. Closely related strains found among items associated with case-patient cells and diluted Cell Block 64 (CB64), a quaternary ammonium disinfectant, and Break Out (BO), a multipurpose cleaner, highlighted their role as environmental reservoirs for S. marcescens in this outbreak. Comparative genomic analysis suggested outbreak strains were both persistent (identical strains found over long periods and in multiple locations of the correctional facility) and diverse (strains clustered with multiple global samples from NCBI database). No correlation was found between antimicrobial resistance (AMR) genes of outbreak strains; NCBI strains have more AMR genes. Principal component analysis (PCA) of virulence factors associated with persistence and infectivity indicated variation based on phylogroups, including the predominant cluster; identifiable variations among environmental versus clinical strains were not observed. Identification of multiple distinct genetic groups highlights the importance of putting epidemiological genomic studies in a proper genetic context.

Trophic ontogeny of a generalist predator is conserved across space.

Consumers can influence ecological patterns and processes through their trophic roles and contributions to the flow of energy through ecosystems. However, the diet and associated trophic roles of consumers commonly change during ontogeny. Despite the prevalence of ontogenetic variation in trophic roles of most animals, we lack an understanding of whether they change consistently across local populations and broad geographic gradients. We examined how the diet and trophic position of a generalist marine predator varied with ontogeny across seven broadly separated locations (~ 750 km). We observed a high degree of heterogeneity in prey consumed without evidence of spatial structuring in this variability. However, compound-specific isotope analysis of amino acids revealed remarkably consistent patterns of increasing trophic position through ontogeny across local populations, suggesting that the roles of this generalist predator scaled with its body size across space. Given the high degree of diet heterogeneity we observed, this finding suggests that even though the dietary patterns differed, the underlying food web architecture transcended variation in prey species across locations for this generalist consumer. Our research addresses a gap in empirical field work regarding the interplay between stage-structured populations and food webs, and suggests ontogenetic changes in trophic position can be consistent in generalist consumers.

Genetic pathways disrupted by ENPP1 deficiency provide insight into mechanisms of osteoporosis, osteomalacia, and paradoxical mineralization.

Ectonucleotide phosphatase/phosphodiesterase 1 (ENPP1) deficiency results in either lethal arterial calcifications ('Generalized Arterial Calcification of Infancy' - GACI), phosphate wasting rickets ('Autosomal Recessive Hypophosphatemic Rickets type 2' - ARHR2), early onset osteoporosis, or progressive spinal rigidity ('Ossification of the Posterior Longitudinal Ligament' - OPLL). As ENPP1 generates a strong endogenous mineralization inhibitor - extracellular pyrophosphate (PPi) - ENPP1 deficiency should not result in reduced bone volume, and therefore the mechanism ENPP1 associated osteoporosis is not apparent given current understanding of the enzyme's function. To investigate genetic pathways driving the skeletal phenotype of ENPP1 deficiency we compared gene expression in Enpp1asj/asj mice and WT sibling pairs by RNAseq and qPCR in whole bones, and in the liver and kidney by qPCR, directly correlating gene expression with measures of bone microarchitectural and biomechanical phenotypes. Unbiased analysis of the differentially expressed genes compared to relevant human disease phenotypes revealed that Enpp1asj/asj mice exhibit strong signatures of osteoporosis, ARHR2 and OPLL. We found that ENPP1 deficient mice exhibited reduced gene transcription of Wnt ligands in whole bone and increased transcription of soluble Wnt inhibitors in the liver and kidney, suggestive of multiorgan inhibition of Wnt activity. Consistent with Wnt suppression in bone, Collagen gene pathways in bone were significantly decreased and Fgf23 was significantly increased, all of which directly correlated with bone microarchitectural defects and fracture risk in Enpp1asj/asj mice. Moreover, the bone findings in 10-week old mice correlated with Enpp1 transcript counts but not plasma [PPi], suggesting that the skeletal phenotype at 10 weeks is driven by catalytically independent ENPP1 function. In contrast, the bone findings in 23-week Enpp1asj/asj mice strongly correlated with plasma PPi, suggesting that chronically low PPi drives the skeletal phenotype in older mice. Finally, correlation between Enpp1 and Fgf23 transcription suggested ENPP1 regulation of Fgf23, which we confirmed by dosing Enpp1asj/asj mice with soluble ENPP1-Fc and observing suppression of intact plasma FGF23 and ALP. In summary, our findings suggest that osteoporosis associated with ENPP1 deficiency involves the suppression of Wnt via catalytically independent Enpp1 pathways, and validates Enpp1asj/asj mice as tools to better understand OPLL and Paradoxical Mineralization Disorders.

White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability.

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans.

Dosage Compensation of the X Chromosomes in Bovine Germline, Early Embryos, and Somatic Tissues.

Dosage compensation of the mammalian X chromosome (X) was proposed by Susumu Ohno as a mechanism wherein the inactivation of one X in females would lead to doubling the expression of the other. This would resolve the dosage imbalance between eutherian females (XX) versus male (XY) and between a single active X versus autosome pairs (A). Expression ratio of X- and A-linked genes has been relatively well studied in humans and mice, despite controversial results over the existence of upregulation of X-linked genes. Here we report the first comprehensive test of Ohno's hypothesis in bovine preattachment embryos, germline, and somatic tissues. Overall an incomplete dosage compensation (0.5 < X:A < 1) of expressed genes and an excess X dosage compensation (X:A > 1) of ubiquitously expressed "dosage-sensitive" genes were seen. No significant differences in X:A ratios were observed between bovine female and male somatic tissues, further supporting Ohno's hypothesis. Interestingly, preimplantation embryos manifested a unique pattern of X dosage compensation dynamics. Specifically, X dosage decreased after fertilization, indicating that the sperm brings in an inactive X to the matured oocyte. Subsequently, the activation of the bovine embryonic genome enhanced expression of X-linked genes and increased the X dosage. As a result, an excess compensation was exhibited from the 8-cell stage to the compact morula stage. The X dosage peaked at the 16-cell stage and stabilized after the blastocyst stage. Together, our findings confirm Ohno's hypothesis of X dosage compensation in the bovine and extend it by showing incomplete and over-compensation for expressed and "dosage-sensitive" genes, respectively.

Tissue-Specific Transcriptome for Poeciliopsis prolifica Reveals Evidence for Genetic Adaptation Related to the Evolution of a Placental Fish.

The evolution of the placenta is an excellent model to examine the evolutionary processes underlying adaptive complexity due to the recent, independent derivation of placentation in divergent animal lineages. In fishes, the family Poeciliidae offers the opportunity to study placental evolution with respect to variation in degree of post-fertilization maternal provisioning among closely related sister species. In this study, we present a detailed examination of a new reference transcriptome sequence for the live-bearing, matrotrophic fish, Poeciliopsis prolifica, from multiple-tissue RNA-seq data. We describe the genetic components active in liver, brain, late-stage embryo, and the maternal placental/ovarian complex, as well as associated patterns of positive selection in a suite of orthologous genes found in fishes. Results indicate the expression of many signaling transcripts, "non-coding" sequences and repetitive elements in the maternal placental/ovarian complex. Moreover, patterns of positive selection in protein sequence evolution were found associated with live-bearing fishes, generally, and the placental P. prolifica, specifically, that appear independent of the general live-bearer lifestyle. Much of the observed patterns of gene expression and positive selection are congruent with the evolution of placentation in fish functionally converging with mammalian placental evolution and with the patterns of rapid evolution facilitated by the teleost-specific whole genome duplication event.

Rapid Evolutionary Rates and Unique Genomic Signatures Discovered in the First Reference Genome for the Southern Ocean Salp, Salpa thompsoni (Urochordata, Thaliacea).

A preliminary genome sequence has been assembled for the Southern Ocean salp, Salpa thompsoni (Urochordata, Thaliacea). Despite the ecological importance of this species in Antarctic pelagic food webs and its potential role as an indicator of changing Southern Ocean ecosystems in response to climate change, no genomic resources are available for S. thompsoni or any closely related urochordate species. Using a multiple-platform, multiple-individual approach, we have produced a 318,767,936-bp genome sequence, covering >50% of the estimated 602 Mb (±173 Mb) genome size for S. thompsoni Using a nonredundant set of predicted proteins, >50% (16,823) of sequences showed significant homology to known proteins and ∼38% (12,151) of the total protein predictions were associated with Gene Ontology functional information. We have generated 109,958 SNP variant and 9,782 indel predictions for this species, serving as a resource for future phylogenomic and population genetic studies. Comparing the salp genome to available assemblies for four other urochordates, Botryllus schlosseri, Ciona intestinalis, Ciona savignyi and Oikopleura dioica, we found that S. thompsoni shares the previously estimated rapid rates of evolution for these species. High mutation rates are thus independent of genome size, suggesting that rates of evolution >1.5 times that observed for vertebrates are a broad taxonomic characteristic of urochordates. Tests for positive selection implemented in PAML revealed a small number of genes with sites undergoing rapid evolution, including genes involved in ribosome biogenesis and metabolic and immune process that may be reflective of both adaptation to polar, planktonic environments as well as the complex life history of the salps. Finally, we performed an initial survey of small RNAs, revealing the presence of known, conserved miRNAs, as well as novel miRNA genes; unique piRNAs; and mature miRNA signatures for varying developmental stages. Collectively, these resources provide a genomic foundation supporting S. thompsoni as a model species for further examination of the exceptional rates and patterns of genomic evolution shown by urochordates. Additionally, genomic data will allow for the development of molecular indicators of key life history events and processes and afford new understandings and predictions of impacts of climate change on this key species of Antarctic pelagic ecosystems.

From shelf to shelf: assessing historical and contemporary genetic differentiation and connectivity across the Gulf of Mexico in Gag, Mycteroperca microlepis.

Describing patterns of connectivity among populations of species with widespread distributions is particularly important in understanding the ecology and evolution of marine species. In this study, we examined patterns of population differentiation, migration, and historical population dynamics using microsatellite and mitochondrial loci to test whether populations of the epinephelid fish, Gag, Mycteroperca microlepis, an important fishery species, are genetically connected across the Gulf of Mexico and if so, whether that connectivity is attributable to either contemporary or historical processes. Populations of Gag on the Campeche Bank and the West Florida Shelf show significant, but low magnitude, differentiation. Time since divergence/expansion estimates associated with historical population dynamics indicate that any population or spatial expansions indicated by population genetics would have likely occurred in the late Pleistocene. Using coalescent-based approaches, we find that the best model for explaining observed spatial patterns of contemporary genetic variation is one of asymmetric gene flow, with movement from Campeche Bank to the West Florida Shelf. Both estimated migration rates and ecological data support the hypothesis that Gag populations throughout the Gulf of Mexico are connected via present day larval dispersal. Demonstrating this greatly expanded scale of connectivity for Gag highlights the influence of "ghost" populations (sensu Beerli) on genetic patterns and presents a critical consideration for both fisheries management and conservation of this and other species with similar genetic patterns.

Determination of dosage compensation of the mammalian X chromosome by RNA-seq is dependent on analytical approach.

BACKGROUND: An enduring question surrounding sex chromosome evolution is whether effective hemizygosity in the heterogametic sex leads inevitably to dosage compensation of sex-linked genes, and whether this compensation has been observed in a variety of organisms. Incongruence in the conclusions reached in some recent reports has been attributed to different high-throughput approaches to transcriptome analysis. However, recent reports each utilizing RNA-seq to gauge X-linked gene expression relative to autosomal gene expression also arrived at diametrically opposed conclusions regarding X chromosome dosage compensation in mammals. RESULTS: Here we analyze RNA-seq data from X-monosomic female human and mouse tissues, which are uncomplicated by genes that escape X-inactivation, as well as published RNA-seq data to describe relative X expression (RXE). We find that the determination of RXE is highly dependent upon a variety of computational, statistical and biological assumptions underlying RNA-seq analysis. Parameters implemented in short-read mapping programs, choice of reference genome annotation, expression data distribution, tissue source for RNA and RNA-seq library construction method have profound effects on comparing expression levels across chromosomes. CONCLUSIONS: Our analysis shows that the high number of paralogous gene families on the mammalian X chromosome relative to autosomes contributes to the ambiguity in RXE calculations, RNA-seq analysis that takes into account that single- and multi-copy genes are compensated differently supports the conclusion that, in many somatic tissues, the mammalian X is up-regulated compared to the autosomes.

Large-scale natural disturbance alters genetic population structure of the sailfin molly, Poecilia latipinna.

Many inferences about contemporary rates of gene flow are based on the assumption that the observed genetic structure among populations is stable. Recent studies have uncovered several cases in which this assumption is tenuous. Most of those studies have focused on the effects that regular environmental fluctuations can have on genetic structure and gene flow patterns. Occasional catastrophic disturbances could also alter either the distribution of habitat or the spatial distribution of organisms in a way that affects population structure. However, evidence of such effects is sparse in the literature because it is difficult to obtain. Hurricanes, in particular, have the potential to exert dramatic effects on population structure of organisms found on islands or coral reefs or in near shore and coastal habitats. Here we draw on a historic genetic data set and new data to suggest that the genetic structure of sailfin molly (Poecilia latipinna) populations in north Florida was altered dramatically by an unusually large and uncommon type of storm surge associated with Hurricane Dennis in 2005. We compare the spatial pattern of genetic variation in these populations after Hurricane Dennis to the patterns described in an earlier study in this same area. We use comparable genetic data from another region of Florida, collected in the same two periods, to estimate the amount of change expected from typical temporal variation in population structure. The comparative natural history of sailfin mollies in these two regions indicates that the change in population structure produced by the storm surge is not the result of many local extinctions with recolonization from a few refugia but emerged from a pattern of mixing and redistribution.

Two different high throughput sequencing approaches identify thousands of de novo genomic markers for the genetically depleted Bornean elephant.

High throughput sequencing technologies are being applied to an increasing number of model species with a high-quality reference genome. The application and analyses of whole-genome sequence data in non-model species with no prior genomic information are currently under way. Recent sequencing technologies provide new opportunities for gathering genomic data in natural populations, laying the empirical foundation for future research in the field of conservation and population genomics. Here we present the case study of the Bornean elephant, which is the most endangered subspecies of Asian elephant and exhibits very low genetic diversity. We used two different sequencing platforms, the Roche 454 FLX (shotgun) and Illumina, GAIIx (Restriction site associated DNA, RAD) to evaluate the feasibility of the two methodologies for the discovery of de novo markers (single nucleotide polymorphism, SNPs and microsatellites) using low coverage data. Approximately, 6,683 (shotgun) and 14,724 (RAD) SNPs were detected within our elephant sequence dataset. Genotyping of a representative sample of 194 SNPs resulted in a SNP validation rate of ~83 to 94% and 17% of the loci were polymorphic with a low diversity (H(o)=0.057). Different numbers of microsatellites were identified through shotgun (27,226) and RAD (868) techniques. Out of all di-, tri-, and tetra-microsatellite loci, 1,706 loci had sufficient flanking regions (shotgun) while only 7 were found with RAD. All microsatellites were monomorphic in the Bornean but polymorphic in another elephant subspecies. Despite using different sample sizes, and the well known differences in the two platforms used regarding sequence length and throughput, the two approaches showed high validation rate. The approaches used here for marker development in a threatened species demonstrate the utility of high throughput sequencing technologies as a starting point for the development of genomic tools in a non-model species and in particular for a species with low genetic diversity.