Variable effects of captivity on microbiomes in populations of IUCN-endangered Blanding's turtles (Emydoidea blandingii).

AIMS: Microbiome composition is increasingly considered in species reintroduction efforts and may influence survival and reproductive success. Many turtle species are threatened by anthropogenic pressures and are frequently raised in captivity for reintroduction efforts, yet little is known about turtle microbiome composition in either wild or captive settings. Here, we investigated trends in microbiome composition of captive and wild IUCN-endangered Blanding's turtles (Emydoidea blandingii). METHODS AND RESULTS: We amplified and sequenced the V4 region of the 16S rDNA locus from plastron, cloaca, and water samples of wild E. blandingii adults and two populations of captive E. blandingii juveniles being raised for headstarting. Plastron, cloaca, and water-associated microbiomes differed strongly from each other and were highly variable among captive sites and between captive and wild sites. Across plastron, cloaca, and water-associated microbial communities, microbial diversity changed over time, but not in a predictable direction between captive sites. Plastron beta diversity correlated with growth rate in captive samples, indicating that external microbiomes may correlate with individual fitness. CONCLUSIONS: Our results indicate that external and internal microbiomes vary between captive and wild turtles and may reflect differences in fitness of captive-raised individuals.

Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time.

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

Turtle species and ecology drive carapace microbiome diversity in three seasonally interconnected wetland habitats.

Different species of freshwater turtles exhibit primary behaviours ranging from aerial basking to benthic bottom-walking, cycle between wet and dry conditions at different time intervals, and undertake short-distance overland movements between aquatic habitats. These behaviours in turn may impact the accumulation of microbes on external shell surfaces of turtles and provide novel niches for differentiation of microbial communities. We assessed microbial diversity using 16S and 18S rRNA metabarcoding on carapace surfaces of six species of freshwater turtles residing in three adjacent and seasonally interconnected wetland habitats in southeast Oklahoma (United States). Communities were highly diverse, with nearly 4200 prokaryotic and 500 micro-eukaryotic amplicon sequence variants recovered, and included taxa previously reported as common or differentially abundant on turtle shells. The 16S rRNA alpha diversity tended to be highest for two species of benthic turtles, while 18S rRNA alpha diversity was highest for two basking and one shallow-water benthic species. Beta diversity of communities was more strongly differentiated by turtle species than by collection site, and ordination patterns were largely reflective of turtle species' primary habits (i.e. benthic, basking, or benthic-basking). Our data support that freshwater turtles could play a role in microbial ecology and evolution in freshwater habitats and warrant additional exploration including in areas with high native turtle diversity and inter-habitat turtle movements.

The salivary metatranscriptome as an accurate diagnostic indicator of oral cancer.

Despite advances in cancer treatment, the 5-year mortality rate for oral cancers (OC) is 40%, mainly due to the lack of early diagnostics. To advance early diagnostics for high-risk and average-risk populations, we developed and evaluated machine-learning (ML) classifiers using metatranscriptomic data from saliva samples (n = 433) collected from oral premalignant disorders (OPMD), OC patients (n = 71) and normal controls (n = 171). Our diagnostic classifiers yielded a receiver operating characteristics (ROC) area under the curve (AUC) up to 0.9, sensitivity up to 83% (92.3% for stage 1 cancer) and specificity up to 97.9%. Our metatranscriptomic signature incorporates both taxonomic and functional microbiome features, and reveals a number of taxa and functional pathways associated with OC. We demonstrate the potential clinical utility of an AI/ML model for diagnosing OC early, opening a new era of non-invasive diagnostics, enabling early intervention and improved patient outcomes.

Consistent patterns in 16S and 18S microbial diversity from the shells of the common and widespread red-eared slider turtle (Trachemys scripta).

Microbial communities associated with freshwater aquatic habitats and resident species are both critical to and indicative of ecosystem status and organismal health. External surfaces of turtle shells readily accumulate microbial growth and could carry representation of habitat-wide microbial diversity, since they are in regular contact with multiple elements of freshwater environments. Yet, microbial diversity residing on freshwater turtle shells is poorly understood. We applied 16S and 18S metabarcoding to characterize microbiota associated with external shell surfaces of 20 red-eared slider (Trachemys scripta) turtles collected from varied habitats in central and western Oklahoma, and ranging to southeast Iowa. Shell-associated microbial communities were highly diverse, with samples dominated by Bacteroidia and alpha-/gamma-proteobacteria, and ciliophoran alveolates. Alpha diversity was lower on turtle shells compared to shallow-water-associated environmental samples, likely resulting from basking-drying behavior and seasonal scute shedding, while alpha diversity was higher on carapace than plastron surfaces. Beta diversity of turtle shells was similarly differentiated from environmental samples, although sampling site was consistently a significant factor. Deinococcus-Thermus bacteria and ciliophoran alveolates were recovered with significantly higher abundance on turtle shells versus environmental samples, while bacterial taxa known to include human-pathogenic species were variably more abundant between shell and environmental samples. Microbial communities from a single, shared-site collection of the ecologically similar river cooter (P. concinna) largely overlapped with those of T. scripta. These data add to a foundation for further characterization of turtle shell microbial communities across species and habitats, with implications for freshwater habitat assessment, microbial ecology and wildlife conservation efforts.

A clinically validated human capillary blood transcriptome test for global systems biology studies.

To prevent and treat chronic diseases, including cancer, a global application of systems biology is needed. We report here a whole blood transcriptome test that needs only 50 µl of capillary (fingerprick) blood. This test is suitable for global applications because the samples are preserved at ambient temperature for up to 4 weeks and the RNA preservative inactivates all pathogens, enabling safe transportation. Both the laboratory and bioinformatic steps are automated and performed in a clinical lab, which minimizes batch effects and creates unbiased datasets. Given its clinical testing performance and accessibility to traditionally underrepresented and diverse populations, this test offers a unique ability to reveal molecular mechanisms of disease and enable longitudinal, population-scale studies.

Correction to: Implications of sequence variation on the evolution of rRNA.

The article Implications of sequence variation on the evolution of rRNA, written by Matthew M. Parks, Chad M. Kurylo, Jake E. Batchelder, C. Theresa Vincent and Scott C. Blanchard, was originally published electronically on the publisher's internet portal (currently SpringerLink).

Ribosome biogenesis during cell cycle arrest fuels EMT in development and disease.

Ribosome biogenesis is a canonical hallmark of cell growth and proliferation. Here we show that execution of Epithelial-to-Mesenchymal Transition (EMT), a migratory cellular program associated with development and tumor metastasis, is fueled by upregulation of ribosome biogenesis during G1/S arrest. This unexpected EMT feature is independent of species and initiating signal, and is accompanied by release of the repressive nucleolar chromatin remodeling complex (NoRC) from rDNA, together with recruitment of the EMT-driving transcription factor Snai1 (Snail1), RNA Polymerase I (Pol I) and the Upstream Binding Factor (UBF). EMT-associated ribosome biogenesis is also coincident with increased nucleolar recruitment of Rictor, an essential component of the EMT-promoting mammalian target of rapamycin complex 2 (mTORC2). Inhibition of rRNA synthesis in vivo differentiates primary tumors to a benign, Estrogen Receptor-alpha (ERα) positive, Rictor-negative phenotype and reduces metastasis. These findings implicate the EMT-associated ribosome biogenesis program with cellular plasticity, de-differentiation, cancer progression and metastatic disease.

Implications of sequence variation on the evolution of rRNA.

The evolution of the multi-copy family of ribosomal RNA (rRNA) genes is unique in regard to its genetics and genome evolution. Paradoxically, rRNA genes are highly homogenized within and between individuals, yet they are globally distinct between species. Here, we discuss the implications for models of rRNA gene evolution in light of our recent discoveries that ribosomes bearing rRNA sequence variants can affect gene expression and physiology and that intra-individual rRNA alleles exhibit both context- and tissue-specific expression.

ADAM10 controls the differentiation of the coronary arterial endothelium.

The coronary vasculature is crucial for normal heart function, yet much remains to be learned about its development, especially the maturation of coronary arterial endothelium. Here, we show that endothelial inactivation of ADAM10, a key regulator of Notch signaling, leads to defects in coronary arterial differentiation, as evidenced by dysregulated genes related to Notch signaling and arterial identity. Moreover, transcriptome analysis indicated reduced EGFR signaling in A10ΔEC coronary endothelium. Further analysis revealed that A10ΔEC mice have enlarged dysfunctional hearts with abnormal myocardial compaction, and increased expression of venous and immature endothelium markers. These findings provide the first evidence for a potential role for endothelial ADAM10 in cardioprotective homeostatic EGFR signaling and implicate ADAM10/Notch signaling in coronary arterial cell specification, which is vital for normal heart development and function. The ADAM10/Notch signaling pathway thus emerges as a potential therapeutic target for improving the regenerative capacity and maturation of the coronary vasculature.

Endogenous rRNA Sequence Variation Can Regulate Stress Response Gene Expression and Phenotype.

Prevailing dogma holds that ribosomes are uniform in composition and function. Here, we show that nutrient limitation-induced stress in E. coli changes the relative expression of rDNA operons to alter the rRNA composition within the actively translating ribosome pool. The most upregulated operon encodes the unique 16S rRNA, rrsH, distinguished by conserved sequence variation within the small ribosomal subunit. rrsH-bearing ribosomes affect the expression of functionally coherent gene sets and alter the levels of the RpoS sigma factor, the master regulator of the general stress response. These impacts are associated with phenotypic changes in antibiotic sensitivity, biofilm formation, and cell motility and are regulated by stress response proteins, RelA and RelE, as well as the metabolic enzyme and virulence-associated protein, AdhE. These findings establish that endogenously encoded, naturally occurring rRNA sequence variation can modulate ribosome function, central aspects of gene expression regulation, and cellular physiology.

Using controls to limit false discovery in the era of big data.

BACKGROUND: Procedures for controlling the false discovery rate (FDR) are widely applied as a solution to the multiple comparisons problem of high-dimensional statistics. Current FDR-controlling procedures require accurately calculated p-values and rely on extrapolation into the unknown and unobserved tails of the null distribution. Both of these intermediate steps are challenging and can compromise the reliability of the results. RESULTS: We present a general method for controlling the FDR that capitalizes on the large amount of control data often found in big data studies to avoid these frequently problematic intermediate steps. The method utilizes control data to empirically construct the distribution of the test statistic under the null hypothesis and directly compares this distribution to the empirical distribution of the test data. By not relying on p-values, our control data-based empirical FDR procedure more closely follows the foundational principles of the scientific method: that inference is drawn by comparing test data to control data. The method is demonstrated through application to a problem in structural genomics. CONCLUSIONS: The method described here provides a general statistical framework for controlling the FDR that is specifically tailored for the big data setting. By relying on empirically constructed distributions and control data, it forgoes potentially problematic modeling steps and extrapolation into the unknown tails of the null distribution. This procedure is broadly applicable insofar as controlled experiments or internal negative controls are available, as is increasingly common in the big data setting.

Recurrent Loss, Horizontal Transfer, and the Obscure Origins of Mitochondrial Introns in Diatoms (Bacillariophyta).

We sequenced mitochondrial genomes from five diverse diatoms (Toxarium undulatum, Psammoneis japonica, Eunotia naegelii, Cylindrotheca closterium, and Nitzschia sp.), chosen to fill important phylogenetic gaps and help us characterize broadscale patterns of mitochondrial genome evolution in diatoms. Although gene content was strongly conserved, intron content varied widely across species. The vast majority of introns were of group II type and were located in the cox1 or rnl genes. Although recurrent intron loss appears to be the principal underlying cause of the sporadic distributions of mitochondrial introns across diatoms, phylogenetic analyses showed that intron distributions superficially consistent with a recurrent-loss model were sometimes more complicated, implicating horizontal transfer as a likely mechanism of intron acquisition as well. It was not clear, however, whether diatoms were the donors or recipients of horizontally transferred introns, highlighting a general challenge in resolving the evolutionary histories of many diatom mitochondrial introns. Although some of these histories may become clearer as more genomes are sampled, high rates of intron loss suggest that the origins of many diatom mitochondrial introns are likely to remain unclear.

Phylogenomics reveals an extensive history of genome duplication in diatoms (Bacillariophyta).

PREMISE OF THE STUDY: Diatoms are one of the most species-rich lineages of microbial eukaryotes. Similarities in clade age, species richness, and primary productivity motivate comparisons to angiosperms, whose genomes have been inordinately shaped by whole-genome duplication (WGD). WGDs have been linked to speciation, increased rates of lineage diversification, and identified as a principal driver of angiosperm evolution. We synthesized a large but scattered body of evidence that suggests polyploidy may be common in diatoms as well. METHODS: We used gene counts, gene trees, and distributions of synonymous divergence to carry out a phylogenomic analysis of WGD across a diverse set of 37 diatom species. KEY RESULTS: Several methods identified WGDs of varying age across diatoms. Determining the occurrence, exact number, and placement of events was greatly impacted by uncertainty in gene trees. WGDs inferred from synonymous divergence of paralogs varied depending on how redundancy in transcriptomes was assessed, gene families were assembled, and synonymous distances (Ks) were calculated. Our results highlighted a need for systematic evaluation of key methodological aspects of Ks-based approaches to WGD inference. Gene tree reconciliations supported allopolyploidy as the predominant mode of polyploid formation, with strong evidence for ancient allopolyploid events in the thalassiosiroid and pennate diatom clades. CONCLUSIONS: Our results suggest that WGD has played a major role in the evolution of diatom genomes. We outline challenges in reconstructing paleopolyploid events in diatoms that, together with these results, offer a framework for understanding the impact of genome duplication in a group that likely harbors substantial genomic diversity.

Variant ribosomal RNA alleles are conserved and exhibit tissue-specific expression.

The ribosome, the integration point for protein synthesis in the cell, is conventionally considered a homogeneous molecular assembly that only passively contributes to gene expression. Yet, epigenetic features of the ribosomal DNA (rDNA) operon and changes in the ribosome's molecular composition have been associated with disease phenotypes, suggesting that the ribosome itself may possess inherent regulatory capacity. Analyzing whole-genome sequencing data from the 1000 Genomes Project and the Mouse Genomes Project, we find that rDNA copy number varies widely across individuals, and we identify pervasive intra- and interindividual nucleotide variation in the 5S, 5.8S, 18S, and 28S ribosomal RNA (rRNA) genes of both human and mouse. Conserved rRNA sequence heterogeneities map to functional centers of the assembled ribosome, variant rRNA alleles exhibit tissue-specific expression, and ribosomes bearing variant rRNA alleles are present in the actively translating ribosome pool. These findings provide a critical framework for exploring the possibility that the expression of genomically encoded variant rRNA alleles gives rise to physically and functionally heterogeneous ribosomes that contribute to mammalian physiology and human disease.

Signal, Uncertainty, and Conflict in Phylogenomic Data for a Diverse Lineage of Microbial Eukaryotes (Diatoms, Bacillariophyta).

Diatoms (Bacillariophyta) are a species-rich group of eukaryotic microbes diverse in morphology, ecology, and metabolism. Previous reconstructions of the diatom phylogeny based on one or a few genes have resulted in inconsistent resolution or low support for critical nodes. We applied phylogenetic paralog pruning techniques to a data set of 94 diatom genomes and transcriptomes to infer perennially difficult species relationships, using concatenation and summary-coalescent methods to reconstruct species trees from data sets spanning a wide range of thresholds for taxon and column occupancy in gene alignments. Conflicts between gene and species trees decreased with both increasing taxon occupancy and bootstrap cutoffs applied to gene trees. Concordance between gene and species trees was lowest for short internodes and increased logarithmically with increasing edge length, suggesting that incomplete lineage sorting disproportionately affects species tree inference at short internodes, which are a common feature of the diatom phylogeny. Although species tree topologies were largely consistent across many data treatments, concatenation methods appeared to outperform summary-coalescent methods for sparse alignments. Our results underscore that approaches to species-tree inference based on few loci are likely to be misled by unrepresentative sampling of gene histories, particularly in lineages that may have diversified rapidly. In addition, phylogenomic studies of diatoms, and potentially other hyperdiverse groups, should maximize the number of gene trees with high taxon occupancy, though there is clearly a limit to how many of these genes will be available.

An exact test for comparing a fixed quantitative property between gene sets.

Motivation: A significant difference in the distribution of a feature between two gene sets can provide insight into function or regulation. This statistical setting differs from much of hypothesis testing theory because the genome is often considered to be effectively fixed, finite and entirely known in commonly studied organisms, such as human. The Mann-Whitney U test is commonly employed in this scenario despite the assumptions of the test not being met, leading to unreliable and generally underpowered results. Permutation tests are also commonly employed for this purpose, but are computationally burdensome and are not tractable for obtaining small P values or for multiple comparisons. Results: We present an exact test for the null hypothesis that gene set membership is independent of the quantitative gene feature of interest. We derive an analytic expression for the randomization distribution of the median of the quantitative feature under the null hypothesis. Efficient implementation permits calculation of precise P values of arbitrary magnitude and makes thousands of simultaneous tests of transcriptome-sized gene sets computationally tractable. The flexibility of the hypothesis testing framework presented permits extension to a variety of related tests commonly found in genomics. The exact test is used to identify signatures of translation control and protein function in the human genome. Availability and implementation: The exact test presented here is implemented in R in the package kpmt available on CRAN. Contact: map2085@med.cornell.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Wnt5a Signals through DVL1 to Repress Ribosomal DNA Transcription by RNA Polymerase I.

Ribosome biogenesis is essential for cell growth and proliferation and is commonly elevated in cancer. Accordingly, numerous oncogene and tumor suppressor signaling pathways target rRNA synthesis. In breast cancer, non-canonical Wnt signaling by Wnt5a has been reported to antagonize tumor growth. Here, we show that Wnt5a rapidly represses rDNA gene transcription in breast cancer cells and generates a chromatin state with reduced transcription of rDNA by RNA polymerase I (Pol I). These effects were specifically dependent on Dishevelled1 (DVL1), which accumulates in nucleolar organizer regions (NORs) and binds to rDNA regions of the chromosome. Upon DVL1 binding, the Pol I transcription activator and deacetylase Sirtuin 7 (SIRT7) releases from rDNA loci, concomitant with disassembly of Pol I transcription machinery at the rDNA promoter. These findings reveal that Wnt5a signals through DVL1 to suppress rRNA transcription. This provides a novel mechanism for how Wnt5a exerts tumor suppressive effects and why disruption of Wnt5a signaling enhances mammary tumor growth in vivo.

Genome Sequence and Analysis of Escherichia coli MRE600, a Colicinogenic, Nonmotile Strain that Lacks RNase I and the Type I Methyltransferase, EcoKI.

Escherichia coli strain MRE600 was originally identified for its low RNase I activity and has therefore been widely adopted by the biomedical research community as a preferred source for the expression and purification of transfer RNAs and ribosomes. Despite its widespread use, surprisingly little information about its genome or genetic content exists. Here, we present the first de novo assembly and description of the MRE600 genome and epigenome. To provide context to these studies of MRE600, we include comparative analyses with E. coli K-12 MG1655 (K12). Pacific Biosciences Single Molecule, Real-Time sequencing reads were assembled into one large chromosome (4.83 Mb) and three smaller plasmids (89.1, 56.9, and 7.1 kb). Interestingly, the 7.1-kb plasmid possesses genes encoding a colicin E1 protein and its associated immunity protein. The MRE600 genome has a G + C content of 50.8% and contains a total of 5,181 genes, including 4,913 protein-encoding genes and 268 RNA genes. We identified 41,469 modified DNA bases (0.83% of total) and found that MRE600 lacks the gene for type I methyltransferase, EcoKI. Phylogenetic, taxonomic, and genetic analyses demonstrate that MRE600 is a divergent E. coli strain that displays features of the closely related genus, Shigella. Nevertheless, comparative analyses between MRE600 and E. coli K12 show that these two strains exhibit nearly identical ribosomal proteins, ribosomal RNAs, and highly homologous tRNA species. Substantiating prior suggestions that MRE600 lacks RNase I activity, the RNase I-encoding gene, rna, contains a single premature stop codon early in its open-reading frame.

Functional Dynamics within the Human Ribosome Regulate the Rate of Active Protein Synthesis.

The regulation of protein synthesis contributes to gene expression in both normal physiology and disease, yet kinetic investigations of the human translation mechanism are currently lacking. Using single-molecule fluorescence imaging methods, we have quantified the nature and timing of structural processes in human ribosomes during single-turnover and processive translation reactions. These measurements reveal that functional complexes exhibit dynamic behaviors and thermodynamic stabilities distinct from those observed for bacterial systems. Structurally defined sub-states of pre- and post-translocation complexes were sensitive to specific inhibitors of the eukaryotic ribosome, demonstrating the utility of this platform to probe drug mechanism. The application of three-color single-molecule fluorescence resonance energy transfer (smFRET) methods further revealed a long-distance allosteric coupling between distal tRNA binding sites within ribosomes bearing three tRNAs, which contributed to the rate of processive translation.