happi: a hierarchical approach to pangenomics inference.

Recovering metagenome-assembled genomes (MAGs) from shotgun sequencing data is an increasingly common task in microbiome studies, as MAGs provide deeper insight into the functional potential of both culturable and non-culturable microorganisms. However, metagenome-assembled genomes vary in quality and may contain omissions and contamination. These errors present challenges for detecting genes and comparing gene enrichment across sample types. To address this, we propose happi, an approach to testing hypotheses about gene enrichment that accounts for genome quality. We illustrate the advantages of happi over existing approaches using published Saccharibacteria MAGs, Streptococcus thermophilus MAGs, and via simulation.

Evaluating replicability in microbiome data.

High-throughput sequencing is widely used to study microbial communities. However, choice of laboratory protocol is known to affect the resulting microbiome data, which has an unquantified impact on many comparisons between communities of scientific interest. We propose a novel approach to evaluating replicability in high-dimensional data and apply it to assess the cross-laboratory replicability of signals in microbiome data using the Microbiome Quality Control Project data set. We learn distinctions between samples as measured by a single laboratory and evaluate whether the same distinctions hold in data produced by other laboratories. While most sequencing laboratories can consistently distinguish between samples (median correct classification 87% on genus-level proportion data), these distinctions frequently fail to hold in data from other laboratories (median correct classification 55% across laboratory on genus-level proportion data). As identical samples processed by different laboratories generate substantively different quantitative results, we conclude that 16S sequencing does not reliably resolve differences in human microbiome samples. However, because we observe greater replicability under certain data transformations, our results inform the analysis of microbiome data.

Segment number threshold determines juvenile onset of germline cluster expansion in Platynereis dumerilii.

Development of sexual characters and generation of gametes are tightly coupled with growth. Platynereis dumerilii is a marine annelid that has been used to study germline development and gametogenesis. P. dumerilii has germ cell clusters found across the body in the juvenile worms, and the clusters eventually form the gametes. Like other segmented worms, P. dumerilii grows by adding new segments at its posterior end. The number of segments reflect the growth state of the worms and therefore is a useful and measurable growth state metric to study the growth-reproduction crosstalk. To understand how growth correlates with progression of gametogenesis, we investigated germline development across several developmental stages. We discovered a distinct transition period when worms increase the number of germline clusters at a particular segment number threshold. Additionally, we found that keeping worms short in segment number, by manipulating environmental conditions or via amputations, supported a segment number threshold requirement for germline development. Finally, we asked if these clusters in P. dumerilii play a role in regeneration (as similar free-roaming cells are observed in Hydra and planarian regeneration) and found that the clusters were not required for regeneration in P. dumerilii, suggesting a strictly germline nature. Overall, these molecular analyses suggest a previously unidentified developmental transition dependent on the growth state of juvenile P. dumerilii leading to substantially increased germline expansion.

A Genomic Analysis of Factors Driving lincRNA Diversification: Lessons from Plants.

Transcriptomic analyses from across eukaryotes indicate that most of the genome is transcribed at some point in the developmental trajectory of an organism. One class of these transcripts is termed long intergenic noncoding RNAs (lincRNAs). Recently, attention has focused on understanding the evolutionary dynamics of lincRNAs, particularly their conservation within genomes. Here, we take a comparative genomic and phylogenetic approach to uncover factors influencing lincRNA emergence and persistence in the plant family Brassicaceae, to which Arabidopsis thaliana belongs. We searched 10 genomes across the family for evidence of > 5000 lincRNA loci from A. thaliana From loci conserved in the genomes of multiple species, we built alignments and inferred phylogeny. We then used gene tree/species tree reconciliation to examine the duplication history and timing of emergence of these loci. Emergence of lincRNA loci appears to be linked to local duplication events, but, surprisingly, not whole genome duplication events (WGD), or transposable elements. Interestingly, WGD events are associated with the loss of loci for species having undergone relatively recent polyploidy. Lastly, we identify 1180 loci of the 6480 previously annotated A. thaliana lincRNAs (18%) with elevated levels of conservation. These conserved lincRNAs show higher expression, and are enriched for stress-responsiveness and cis-regulatory motifs known as conserved noncoding sequences (CNSs). These data highlight potential functional pathways and suggest that CNSs may regulate neighboring genes at both the genomic and transcriptomic level. In sum, we provide insight into processes that may influence lincRNA diversification by providing an evolutionary context for previously annotated lincRNAs.