Scaling read aligners to hundreds of threads on general-purpose processors

The Abundant and Unique Transcripts and Alternative Splicing of the Artificially Autododecaploid London Plane (Platanus × acerifolia)

Transcription and alternative splicing (AS) are now appreciated in plants, but few studies have examined the effects of changing ploidy on transcription and AS. In this study, we showed that artificially autododecaploid plants of London plane (Platanus × acerifolia (Aiton) Willd) had few flowers relative to their hexaploid progenitors. Transcriptome analysis based on full-length Oxford Nanopore Technologies (ONTs) and next-generation sequencing (NGS) revealed that the increased ploidy level in P. × acerifolia led to more transcribed isoforms, accompanied by an increase in the number of isoforms per gene. The functional enrichment of genes indicated that novel genes transcribed specifically in the dodecaploids may have been highly correlated with the ability to maintain genome stability. The dodecaploids showed a higher number of genes with upregulated differentially expressed genes (DEGs) compared with the hexaploid counterpart. The genome duplication of P. × acerifolia resulted mainly in the DEGs involved in basic biological pathways. It was noted that there was a greater abundance of alternative splicing (AS) events and AS genes in the dodecaploids compared with the hexaploids in P. × acerifolia. In addition, a significant difference between the structure and expression of AS events between the hexaploids and dodecaploids of Platanus was found. Of note, some DEGs and differentially spliced genes (DSGs) related to floral transition and flower development were consistent with the few flower traits in the dodecaploids of P. × acerifolia. Collectively, our findings explored the difference in transcription and AS regulation between the hexaploids and dodecaploids of P. × acerifolia and gained new insight into the molecular mechanisms underlying the few-flower phenotype of P. × acerifolia. These results contribute to uncovering the regulatory role of transcription and AS in polyploids and breeding few-flower germplasms.

Structural and functional characterization of MrpR, the master repressor of the Bacillus subtilis prophage SPβ

Prophages control their lifestyle to either be maintained within the host genome or enter the lytic cycle. Bacillus subtilis contains the SPβ prophage whose lysogenic state depends on the MrpR (YopR) protein, a key component of the lysis-lysogeny decision system. Using a historic B. subtilis strain harboring the heat-sensitive SPβ c2 mutant, we demonstrate that the lytic cycle of SPβ c2 can be induced by heat due to a single nucleotide exchange in the mrpR gene, rendering the encoded MrpRG136E protein temperature-sensitive. Structural characterization revealed that MrpR is a DNA-binding protein resembling the overall fold of tyrosine recombinases. MrpR has lost its recombinase function and the G136E exchange impairs its higher-order structure and DNA binding activity. Genome-wide profiling of MrpR binding revealed its association with the previously identified SPbeta repeated element (SPBRE) in the SPβ genome. MrpR functions as a master repressor of SPβ that binds to this conserved element to maintain lysogeny. The heat-inducible excision of the SPβ c2 mutant remains reliant on the serine recombinase SprA. A suppressor mutant analysis identified a previously unknown component of the lysis-lysogeny management system that is crucial for the induction of the lytic cycle of SPβ.

Metagenomic sequencing reveals altered gut microbial compositions and gene functions in patients with non-segmental vitiligo

BACKGROUND

Vitiligo has been correlated with an abnormal gut microbiota. We aimed to systematically identify characteristics of the gut microbial compositions, genetic functions, and potential metabolic features in patients with non-segmental vitiligo.

METHODS

Twenty-five patients with non-segmental vitiligo and 25 matched healthy controls (HCs) were enrolled. Metagenomic sequencing and bioinformatic analysis were performed to determine the gut microbiota profiles. Differences in gut microbiota diversity and composition between patients with vitiligo and HCs were analyzed. Gene functions and gut metabolic modules were predicted with the Kyoto Encyclopedia of Gene and Genomes (KEGG) and MetaCyc databases.

RESULTS

Compared with HCs, alpha diversity of intestinal microbiome in vitiligo patients was significantly reduced. At the species level, the relative abundance of Staphylococcus thermophiles was decreased, and that of Bacteroides fragilis was increased in patients with vitiligo compared with those of the HCs. Linear discriminant analysis (LDA) effect size (LEfSe) analysis revealed representative microbial markers of Lachnospiraceae_bacterium_BX3, Massilioclostridium_coli, TM7_phylum_sp_oral_taxon_348 and Bacteroides_fragilis for patients with vitiligo. KEGG gene function analysis showed that the NOD-like receptor signaling pathway was significantly enriched in patients with vitiligo. Gut metabolic modules (GMMs) analysis showed that cysteine degradation was significantly down-regulated, and galactose degradation was up-regulated in patients with vitiligo. A panel of 28 microbial features was constructed to distinguish patients with vitiligo from HCs.

CONCLUSIONS

The gut microbial profiles and genetic functions of patients with vitiligo were distinct from those of the HCs. The identified gut microbial markers may potentially be used for earlier diagnosis and treatment targets.

Isoenzymes of the Flavonoid and Phenylpropanoid Pathways Show Organ-Specific Regulation during Apple Fruit Development

Elucidating the molecular mechanisms controlling fruit development is a primary target for the improvement of new apple (Malus × domestica Borkh.) cultivars. The first two weeks of development following pollination are crucial to determine fruit characteristics. During this period, a lot of changes take place in apple fruit, going from rapid cell division to the production of important metabolites. In this work, attention was focused on the phenylpropanoid and flavonoid pathways responsible for the production of numerous compounds contributing to fruit quality, such as flavonols, catechins, dihydrochalcones and anthocyanins. A total of 17 isoenzymes were identified, belonging to seven classes of the phenylpropanoid and flavonoid pathways that, despite showing more than 80% sequence identity, showed differential expression regulation during the first two weeks of apple fruit development. This feature seems to be quite common for most of the enzymes of both pathways. Differential regulation of isoenzymes was shown to be present in both 'Golden Delicious' and a wild relative (Malus mandshurica), even though differences were also present. Each isoenzyme showed a specific pattern of expression in the flower and fruit organs, suggesting that genes coding for enzymes with the same function may control different aspects of plant biology. Finally, promoter analysis was performed in order to highlight differences in the number and type of regulatory motifs. Overall, our results indicate that the control of the expression of genes involved in the phenylpropanoid and flavonoid pathways may be very complex as not only enzymes belonging to the same class, but even putative isoenzymes, can have different roles for the plant. Such genes may represent an important regulatory mechanism, as they would allow the plant to fine-tune the processing of metabolic intermediates towards different branches of the pathway, for example, in an organ-specific way.

Disentangling the interactions between nasopharyngeal and gut microbiome and their involvement in the modulation of COVID-19 infection

The human organism is inhabited by trillions of microorganisms, known as microbiota, which are considered to exploit a pivotal role in the regulation of host health and immunity. Recent investigations have suggested a relationship between the composition of the human microbiota and COVID-19 infection, highlighting a possible role of bacterial communities in the modulation of the disease severity. In this study, we performed a shotgun metagenomics analysis to explore and compare the nasopharyngeal microbiota of 38 hospitalized Italian patients with and without COVID-19 infection during the third and fourth pandemic waves. In detail, the metagenomic analysis combined with specific correlation analyses suggested a positive association of several microbial species, such as S. parasanguinis and P. melaninogenica, with the severity of COVID-19 infection. Furthermore, the comparison of the microbiota composition between the nasopharyngeal and their respective fecal samples highlighted an association between these different compartments represented by a sharing of several bacterial species. Additionally, lipidomic and deep-shotgun functional analyses of the fecal samples suggested a metabolic impact of the microbiome on the host's immune response, indicating the presence of key metabolic compounds in COVID-19 patients, such as lipid oxidation end products, potentially related to the inflammatory state. Conversely, the patients without COVID-19 displayed enzymatic patterns associated with the biosynthesis and degradation of specific compounds like lysine (synthesis) and phenylalanine (degradation) that could positively impact disease severity and contribute to modulating COVID-19 infection. IMPORTANCE The human microbiota is reported to play a major role in the regulation of host health and immunity, suggesting a possible impact on the severity of COVID-19 disease. This preliminary study investigated the possible correlation between nasopharyngeal microbiota and COVID-19 infection. In detail, the analysis of the nasopharyngeal microbiota of hospitalized Italian patients with and without COVID-19 infection suggested a positive association of several microbial species with the severity of the disease and highlighted a sharing of several bacteria species with the respective fecal samples. Moreover, the metabolic analyses suggested a possible impact of the microbiome on the host's immune response and the disease severity.

Measuring the transcriptome-wide effects of aging on murine adipocytes using RNAseq

Adipose tissue plays a central role in age-related diseases. While RNAseq protocols exist for many tissues, few data have been generated with this technology to explore gene expression in adipocytes, particularly during aging. Here, we present a protocol to analyze the transcriptional changes that occur in adipose tissue during normal and accelerated aging in mouse models. We describe steps for genotyping, diet control, euthanasia, and dissection. We then detail RNA purification and genome-wide data generation and analysis. For complete details on the use and execution of this protocol, please refer to De Cauwer et al. (2022) iScience. Sep 16;25(10):105149.

Comparison of the three-dimensional chromatin structures of adolescent and adult peripheral blood B cells: implications for the study of pediatric autoimmune diseases

Background/Purpose

Knowledge of the 3D genome is essential to elucidate genetic mechanisms driving autoimmune diseases. The 3D genome is distinct for each cell type, and it is uncertain whether cell lines faithfully recapitulate the 3D architecture of primary human cells or whether developmental aspects of the pediatric immune system require use of pediatric samples. We undertook a systematic analysis of B cells and B cell lines to compare 3D genomic features encompassing risk loci for juvenile idiopathic arthritis (JIA), systemic lupus (SLE), and type 1 diabetes (T1D).

Methods

We isolated B cells from healthy individuals, ages 9-17. HiChIP was performed using CTCF antibody, and CTCF peaks were identified. CTCF loops within the pediatric were compared to three datasets: 1) self-called CTCF consensus peaks called within the pediatric samples, 2) ENCODE's publicly available GM12878 CTCF ChIP-seq peaks, and 3) ENCODE's primary B cell CTCF ChIPseq peaks from two adult females. Differential looping was assessed within the pediatric samples and each of the three peak datasets.

Results

The number of consensus peaks called in the pediatric samples was similar to that identified in ENCODE's GM12878 and primary B cell datasets. We observed <1% of loops that demonstrated significantly differential looping between peaks called within the pediatric samples themselves and when called using ENCODE GM12878 peaks . Significant looping differences were even less when comparing loops of the pediatric called peaks to those of the ENCODE primary B cell peaks. When querying loops found in juvenile idiopathic arthritis, type 1 diabetes, or systemic lupus erythematosus risk haplotypes, we observed significant differences in only 2.2%, 1.0%, and 1.3% loops, respectively, when comparing peaks called within the pediatric samples and ENCODE GM12878 dataset. The differences were even less apparent when comparing loops called with the pediatric vs ENCODE adult primary B cell peak datasets.The 3D chromatin architecture in B cells is similar across pediatric, adult, and EBVtransformed cell lines. This conservation of 3D structure includes regions encompassing autoimmune risk haplotypes.

Conclusion

Thus, even for pediatric autoimmune diseases, publicly available adult B cell and cell line datasets may be sufficient for assessing effects exerted in the 3D genomic space.

Growth phase estimation for abundant bacterial populations sampled longitudinally from human stool metagenomes

Longitudinal sampling of the stool has yielded important insights into the ecological dynamics of the human gut microbiome. However, human stool samples are available approximately once per day, while commensal population doubling times are likely on the order of minutes-to-hours. Despite this mismatch in timescales, much of the prior work on human gut microbiome time series modeling has assumed that day-to-day fluctuations in taxon abundances are related to population growth or death rates, which is likely not the case. Here, we propose an alternative model of the human gut as a stationary system, where population dynamics occur internally and the bacterial population sizes measured in a bolus of stool represent a steady-state endpoint of these dynamics. We formalize this idea as stochastic logistic growth. We show how this model provides a path toward estimating the growth phases of gut bacterial populations in situ. We validate our model predictions using an in vitro Escherichia coli growth experiment. Finally, we show how this method can be applied to densely-sampled human stool metagenomic time series data. We discuss how these growth phase estimates may be used to better inform metabolic modeling in flow-through ecosystems, like animal guts or industrial bioreactors.

Comparison of the Whole-Genome Sequence of the African Swine Fever Virus from a Mongolian Wild Boar with Genotype II Viruses from Asia and Europe

African swine fever (ASF) is a highly contagious and severe viral hemorrhagic disease in domestic and wild pigs. ASF seriously affects the global swine industry as the mortality rate can reach 100% with highly virulent strains. In 2007, ASF was introduced into the Caucasus and spread to Russia and later into other European and Asian countries. This study reported the first whole-genome sequence (WGS) of the ASF virus (ASFV) that was detected in a Mongolian wild boar. This sequence was then compared to other WGS samples from Asia and Europe. Results show that the ASFV Genotype II from Mongolia is similar to the Asian Genotype II WGS. However, there were three nucleotide differences found between the Asian and European genome sequences, two of which were non-synonymous. It was also observed that the European Genotype II ASFV WGS was more diverse than that of the Asian counterparts. The study demonstrates that the ASFV Genotype II variants found in wild boars and domestic pigs are highly similar, suggesting these animals might have had direct or indirect contact, potentially through outdoor animal breeding. In conclusion, this study provides a WGS and mutation spectrum of the ASFV Genotype II WGS in Asia and Europe and thus provides important insights into the origin and spread of ASFV in Mongolia.

Gammaretroviruses, novel viruses and pathogenic bacteria in Australian bats with neurological signs, pneumonia and skin lesions

More than 70 bat species are found in mainland Australia. While most studies of bat viromes focus on sampling seemingly healthy individuals, little is known about the viruses and bacteria associated with diseased bats. We performed traditional diagnostic techniques and metatranscriptomic sequencing on tissue samples from 43 Australian bats, comprising three flying fox (Pteropodidae) and two microbat species experiencing a range of disease syndromes, including mass mortality, neurological signs, pneumonia and skin lesions. Of note, we identified the recently discovered Hervey pteropid gammaretrovirus in a bat with lymphoid leukemia, with evidence of replication consistent with an exogenous virus. The possible association of Hervey pteropid gammaretrovirus with lymphoid leukemia clearly merits additional investigation. One novel picornavirus and at least three new astroviruses and bat pegiviruses were also identified in a variety of tissue types, as well as a number of likely bacterial pathogens or opportunistic infections, most notably Pseudomonas aeruginosa.

A multi-omic Nicotiana benthamiana resource for fundamental research and biotechnology

Nicotiana benthamiana is an invaluable model plant and biotechnology platform with a ~3 Gb allotetraploid genome. To further improve its usefulness and versatility, we have produced high-quality chromosome-level genome assemblies, coupled with transcriptome, epigenome, microRNA and transposable element datasets, for the ubiquitously used LAB strain and a related wild accession, QLD. In addition, single nucleotide polymorphism maps have been produced for a further two laboratory strains and four wild accessions. Despite the loss of five chromosomes from the ancestral tetraploid, expansion of intergenic regions, widespread segmental allopolyploidy, advanced diploidization and evidence of recent bursts of Copia pseudovirus (Copia) mobility not seen in other Nicotiana genomes, the two subgenomes of N. benthamiana show large regions of synteny across the Solanaceae. LAB and QLD have many genetic, metabolic and phenotypic differences, including disparate RNA interference responses, but are highly interfertile and amenable to genome editing and both transient and stable transformation. The LAB/QLD combination has the potential to be as useful as the Columbia-0/Landsberg errecta partnership, utilized from the early pioneering days of Arabidopsis genomics to today.

Characterization of the Environmental Plasmidome of the Red Sea

Plasmids contribute to microbial diversity and adaptation, providing microorganisms with the ability to thrive in a wide range of conditions in extreme environments. However, while the number of marine microbiome studies is constantly increasing, very little is known about marine plasmids, and they are very poorly represented in public databases. To extend the repertoire of environmental marine plasmids, we established a pipeline for the de novo assembly of plasmids in the marine environment by analyzing available microbiome metagenomic sequencing data. By applying the pipeline to data from the Red Sea, we identified 362 plasmid candidates. We showed that the distribution of plasmids corresponds to environmental conditions, particularly, depth, temperature, and physical location. At least 7 of the 362 candidates are most probably real plasmids, based on a functional analysis of their open reading frames (ORFs). Only one of the seven has been described previously. Three plasmids were identified in other public marine metagenomic data from different locations all over the world; these plasmids contained different cassettes of functional genes at each location. Analysis of antibiotic and metal resistance genes revealed that the same positions that were enriched with genes encoding resistance to antibiotics were also enriched with resistance to metals, suggesting that plasmids contribute site-dependent phenotypic modules to their ecological niches. Finally, half of the ORFs (50.8%) could not be assigned to a function, emphasizing the untapped potential of the unique marine plasmids to provide proteins with multiple novel functions. IMPORTANCE Marine plasmids are understudied and hence underrepresented in databases. Plasmid functional annotation and characterization is complicated but, if successful, may provide a pool of novel genes and unknown functions. Newly discovered plasmids and their functional repertoire are potentially valuable tools for predicting the dissemination of antimicrobial resistance, providing vectors for molecular cloning and an understanding of plasmid-bacterial interactions in various environments.

Genetic screening in a Brazilian cohort with inborn errors of immunity

BACKGROUND

Inherited genetic defects in immune system-related genes can result in Inborn Errors of Immunity (IEI), also known as Primary Immunodeficiencies (PID). Diagnosis of IEI disorders is challenging due to overlapping clinical manifestations. Accurate identification of disease-causing germline variants is crucial for appropriate treatment, prognosis, and genetic counseling. However, genetic sequencing is challenging in low-income countries like Brazil. This study aimed to perform genetic screening on patients treated within Brazil's public Unified Health System to identify candidate genetic variants associated with the patient's phenotype.

METHODS

Thirteen singleton unrelated patients from three hospitals in Rio de Janeiro were enrolled in this study. Genomic DNA was extracted from the peripheral blood lymphocytes of each patient, and whole exome sequencing (WES) analyses were conducted using Illumina NextSeq. Germline genetic variants in IEI-related genes were prioritized using a computational framework considering their molecular consequence in coding regions; minor allele frequency ≤ 0.01; pathogenicity classification based on American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) guidelines gathered from the VarSome clinical database; and IEI-related phenotype using the Franklin tool. The genes classification into IEI categories follows internationally recognized guidelines informed by the International Union of Immunological Societies Expert Committee. Additional methods for confirmation of the variant included Sanger sequencing, phasing analysis, and splice site prediction.

RESULTS

A total of 16 disease-causing variants in nine genes, encompassing six different IEI categories, were identified. X-Linked Agammaglobulinemia, caused by BTK variations, emerged as the most prevalent IEI disorder in the cohort. However, pathogenic and likely pathogenic variants were also reported in other known IEI-related genes, namely CD40LG, CARD11, WAS, CYBB, C6, and LRBA. Interestingly, two patients with suspected IEI exhibited pathogenic variants in non-IEI-related genes, ABCA12 and SLC25A13, potentially explaining their phenotypes.

CONCLUSIONS

Genetic screening through WES enabled the detection of potentially harmful variants associated with IEI disorders. These findings contribute to a better understanding of patients' clinical manifestations by elucidating the genetic basis underlying their phenotypes.

MicroRNA-mediated translational pathways are regulated in the orbitofrontal cortex and peripheral blood samples during acute abstinence from heroin self-administration

Opioid misuse in the United States contributes to >70% of annual overdose deaths. To develop additional therapeutics that may prevent opioid misuse, further studies on the neurobiological consequences of opioid exposure are needed. Here we sought to characterize molecular neuroadaptations involving microRNA (miRNA) pathways in the brain and blood of adult male rats that self-administered the opioid heroin. miRNAs are ∼18-24 nucleotide RNAs that regulate protein expression by preventing mRNA translation into proteins. Manipulation of miRNAs and their downstream pathways can critically regulate drug seeking behavior. We performed small-RNA sequencing of miRNAs and proteomics profiling on tissue from the orbitofrontal cortex (OFC), a brain region associated with heroin seeking, following 2 days of forced abstinence from self-administration of 0.03 mg/kg/infusion heroin or sucrose. Heroin self-administration resulted in a robust shift of the OFC miRNA profile, regulating 77 miRNAs, while sucrose self-administration only regulated 9 miRNAs that did not overlap with the heroin-induced profile. Conversely, proteomics revealed dual regulation of seven proteins by both heroin and sucrose in the OFC. Pathway analysis determined that heroin-associated miRNA pathways are predicted to target genes associated with the term "prion disease," a term that was also enriched in the heroin-induced protein expression dataset. Lastly, we confirmed that a subset of heroin-induced miRNA expression changes in the OFC are regulated in peripheral serum and correlate with heroin infusions. These findings demonstrate that peripheral blood samples may have biomarker utility for assessment of drug-induced miRNA pathway alterations that occur in the brain following chronic drug exposure.

Systems-ecology designed bacterial consortium protects from severe Clostridioides difficile infection

Fecal Microbiota Transplant (FMT) is an emerging therapy that has had remarkable success in treatment and prevention of recurrent Clostridioides difficile infection (rCDI). FMT has recently been associated with adverse outcomes such as inadvertent transfer of antimicrobial resistance, necessitating development of more targeted bacteriotherapies. To address this challenge, we developed a novel systems biology pipeline to identify candidate probiotic strains that would be predicted to interrupt C. difficile pathogenesis. Utilizing metagenomic characterization of human FMT donor samples, we identified those metabolic pathways most associated with successful FMTs and reconstructed the metabolism of encoding species to simulate interactions with C. difficile . This analysis resulted in predictions of high levels of cross-feeding for amino acids in species most associated with FMT success. Guided by these in silico models, we assembled consortia of bacteria with increased amino acid cross-feeding which were then validated in vitro . We subsequently tested the consortia in a murine model of CDI, demonstrating total protection from severe CDI through decreased toxin levels, recovered gut microbiota, and increased intestinal eosinophils. These results support the novel framework that amino acid cross-feeding is likely a critical mechanism in the initial resolution of CDI by FMT. Importantly, we conclude that our predictive platform based on predicted and testable metabolic interactions between the microbiota and C. difficile led to a rationally designed biotherapeutic framework that may be extended to other enteric infections.

Rho enhancers play unexpectedly minor roles in Rhodopsin transcription and rod cell integrity

Enhancers function with a basal promoter to control the transcription of target genes. Enhancer regulatory activity is often studied using reporter-based transgene assays. However, unmatched results have been reported when selected enhancers are silenced in situ. In this study, using genomic deletion analysis in mice, we investigated the roles of two previously identified enhancers and the promoter of the Rho gene that codes for the visual pigment rhodopsin. The Rho gene is robustly expressed by rod photoreceptors of the retina, and essential for the subcellular structure and visual function of rod photoreceptors. Mutations in RHO cause severe vision loss in humans. We found that each Rho regulatory region can independently mediate local epigenomic changes, but only the promoter is absolutely required for establishing active Rho chromatin configuration and transcription and maintaining the cell integrity and function of rod photoreceptors. To our surprise, two Rho enhancers that enable strong promoter activation in reporter assays are largely dispensable for Rho expression in vivo. Only small and age-dependent impact is detectable when both enhancers are deleted. Our results demonstrate context-dependent roles of enhancers and highlight the importance of studying functions of cis-regulatory regions in the native genomic context.

Assessment of prokaryotic communities in Southwestern Atlantic deep-sea sediments reveals prevalent methanol-oxidising Methylomirabilales

Continental slopes can play a significant contribution to marine productivity and carbon cycling. These regions can harbour distinct geological features, such as salt diapirs and pockmarks, in which their depressions may serve as natural sediment traps where different compounds can accumulate. We investigated the prokaryotic communities in surface (0-2 cm) and subsurface (18-20 or 22-24 cm) sediments from a salt diapir and pockmark field in Santos Basin, Southwest Atlantic Ocean. Metabarcoding of 16 samples revealed that surface sediments were dominated by the archaeal class Nitrososphaeria, while the bacterial class Dehalococcoidia was the most prevalent in subsurface samples. Sediment strata were found to be a significant factor explaining 27% of the variability in community composition. However, no significant difference was observed among geomorphological features. We also performed a metagenomic analysis of three surface samples and analysed the highest quality metagenome-assembled genome retrieved, which belonged to the family CSP1-5, phylum Methylomirabilota. This non-methanotrophic methylotroph contains genes encoding for methanol oxidation and Calvin Cycle pathways, along with diverse functions that may contribute to its adaptation to deep-sea habitats and to oscillating environmental conditions. By integrating metabarcoding and metagenomic approaches, we reported that CSP1-5 is prevalent in the sediment samples from Santos Basin slope, indicating the potential importance of methanol metabolism in this region. Finally, using a phylogenetic approach integrating 16S rRNA sequences assigned to Methylomirabilota in this study with those from a public database, we argued that CSP1-5 public sequences might be misclassified as Methylomirabilaceae (the methanotrophic clade) and, therefore, the role of these organisms and the methanol cycling could also be neglected in other environments.

High-resolution mapping reveals the mechanism and contribution of genome insertions and deletions to RNA virus evolution

RNA viruses rapidly adapt to selective conditions due to the high intrinsic mutation rates of their RNA-dependent RNA polymerases (RdRps). Insertions and deletions (indels) in viral genomes are major contributors to both deleterious mutational load and evolutionary novelty, but remain understudied. To characterize the mechanistic details of their formation and evolutionary dynamics during infection, we developed a hybrid experimental-bioinformatic approach. This approach, called MultiMatch, extracts insertions and deletions from ultradeep sequencing experiments, including those occurring at extremely low frequencies, allowing us to map their genomic distribution and quantify the rates at which they occur. Mapping indel mutations in adapting poliovirus and dengue virus populations, we determine the rates of indel generation and identify mechanistic and functional constraints shaping indel diversity. Using poliovirus RdRp variants of distinct fidelity and genome recombination rates, we demonstrate tradeoffs between fidelity and Indel generation. Additionally, we show that maintaining translation frame and viral RNA structures constrain the Indel landscape and that, due to these significant fitness effects, Indels exert a significant deleterious load on adapting viral populations. Conversely, we uncover positively selected Indels that modulate RNA structure, generate protein variants, and produce defective interfering genomes in viral populations. Together, our analyses establish the kinetic and mechanistic tradeoffs between misincorporation, recombination, and Indel rates and reveal functional principles defining the central role of Indels in virus evolution, emergence, and the regulation of viral infection.

Ezh2 emerges as an epigenetic checkpoint regulator during monocyte differentiation limiting cardiac dysfunction post-MI

Epigenetic regulation of histone H3K27 methylation has recently emerged as a key step during alternative immunoregulatory M2-like macrophage polarization; known to impact cardiac repair after Myocardial Infarction (MI). We hypothesized that EZH2, responsible for H3K27 methylation, could act as an epigenetic checkpoint regulator during this process. We demonstrate for the first time an ectopic EZH2, and putative, cytoplasmic inactive localization of the epigenetic enzyme, during monocyte differentiation into M2 macrophages in vitro as well as in immunomodulatory cardiac macrophages in vivo in the post-MI acute inflammatory phase. Moreover, we show that pharmacological EZH2 inhibition, with GSK-343, resolves H3K27 methylation of bivalent gene promoters, thus enhancing their expression to promote human monocyte repair functions. In line with this protective effect, GSK-343 treatment accelerated cardiac inflammatory resolution preventing infarct expansion and subsequent cardiac dysfunction in female mice post-MI in vivo. In conclusion, our study reveals that pharmacological epigenetic modulation of cardiac-infiltrating immune cells may hold promise to limit adverse cardiac remodeling after MI.

Transcription factors interact with RNA to regulate genes

Transcription factors (TFs) orchestrate the gene expression programs that define each cell's identity. The canonical TF accomplishes this with two domains, one that binds specific DNA sequences and the other that binds protein coactivators or corepressors. We find that at least half of TFs also bind RNA, doing so through a previously unrecognized domain with sequence and functional features analogous to the arginine-rich motif of the HIV transcriptional activator Tat. RNA binding contributes to TF function by promoting the dynamic association between DNA, RNA, and TF on chromatin. TF-RNA interactions are a conserved feature important for vertebrate development and disrupted in disease. We propose that the ability to bind DNA, RNA, and protein is a general property of many TFs and is fundamental to their gene regulatory function.

Novel complete methanogenic pathways in longitudinal genomic study of monogastric age-associated archaea

BACKGROUND

Archaea perform critical roles in the microbiome system, including utilizing hydrogen to allow for enhanced microbiome member growth and influencing overall host health. With the majority of microbiome research focusing on bacteria, the functions of archaea are largely still under investigation. Understanding methanogenic functions during the host lifetime will add to the limited knowledge on archaeal influence on gut and host health. In our study, we determined lifelong archaea dynamics, including detection and methanogenic functions, while assessing global, temporal and host distribution of our novel archaeal metagenome-assembled genomes (MAGs). We followed 7 monogastric swine throughout their life, from birth to adult (1-156 days of age), and collected feces at 22 time points. The samples underwent gDNA extraction, Illumina sequencing, bioinformatic quality and assembly processes, MAG taxonomic assignment and functional annotation. MAGs were utilized in downstream phylogenetic analysis for global, temporal and host distribution in addition to methanogenic functional potential determination.

RESULTS

We generated 1130 non-redundant MAGs, representing 588 unique taxa at the species level, with 8 classified as methanogenic archaea. The taxonomic classifications were as follows: orders Methanomassiliicoccales (5) and Methanobacteriales (3); genera UBA71 (3), Methanomethylophilus (1), MX-02 (1), and Methanobrevibacter (3). We recovered the first US swine Methanobrevibacter UBA71 sp006954425 and Methanobrevibacter gottschalkii MAGs. The Methanobacteriales MAGs were identified primarily during the young, preweaned host whereas Methanomassiliicoccales primarily in the adult host. Moreover, we identified our methanogens in metagenomic sequences from Chinese swine, US adult humans, Mexican adult humans, Swedish adult humans, and paleontological humans, indicating that methanogens span different hosts, geography and time. We determined complete metabolic pathways for all three methanogenic pathways: hydrogenotrophic, methylotrophic, and acetoclastic. This study provided the first evidence of acetoclastic methanogenesis in archaea of monogastric hosts which indicated a previously unknown capability for acetate utilization in methanogenesis for monogastric methanogens. Overall, we hypothesized that the age-associated detection patterns were due to differential substrate availability via the host diet and microbial metabolism, and that these methanogenic functions are likely crucial to methanogens across hosts. This study provided a comprehensive, genome-centric investigation of monogastric-associated methanogens which will further improve our understanding of microbiome development and functions.

Gut dysbiosis in Thai intrahepatic cholangiocarcinoma and hepatocellular carcinoma

Primary liver cancer (PLC), which includes intrahepatic cholangiocarcinoma (iCCA) and hepatocellular carcinoma (HCC), has the highest incidence of all cancer types in Thailand. Known etiological factors, such as viral hepatitis and chronic liver disease do not fully account for the country's unusually high incidence. However, the gut-liver axis, which contributes to carcinogenesis and disease progression, is influenced by the gut microbiome. To investigate this relationship, fecal matter from 44 Thai PLC patients and 76 healthy controls were subjected to whole-genome metagenomic shotgun sequencing and then analyzed by marker gene-based and assembly based methods. Results revealed greater gut microbiome heterogeneity in iCCA compared to HCC and healthy controls. Two Veillonella species were found to be more abundant in iCCA samples and could distinguish iCCA from HCC and healthy controls. Conversely, Ruminococcus gnavus was depleted in iCCA patients and could distinguish HCC from iCCA samples. High Veillonella genus counts in the iCCA group were associated with enriched amino acid biosynthesis and glycolysis pathways, while enriched phospholipid and thiamine metabolism pathways characterized the HCC group with high Blautia genus counts. These findings reveal distinct landscapes of gut dysbiosis among Thai iCCA and HCC patients and warrant further investigation as potential biomarkers.

T-BET and EOMES sustain mature human NK cell identity and antitumor function

Since the T-box transcription factors (TFs) T-BET and EOMES are necessary for initiation of NK cell development, their ongoing requirement for mature NK cell homeostasis, function, and molecular programming remains unclear. To address this, T-BET and EOMES were deleted in unexpanded primary human NK cells using CRISPR/Cas9. Deleting these TFs compromised in vivo antitumor response of human NK cells. Mechanistically, T-BET and EOMES were required for normal NK cell proliferation and persistence in vivo. NK cells lacking T-BET and EOMES also exhibited defective responses to cytokine stimulation. Single-cell RNA-Seq revealed a specific T-box transcriptional program in human NK cells, which was rapidly lost following T-BET and EOMES deletion. Further, T-BET- and EOMES-deleted CD56bright NK cells acquired an innate lymphoid cell precursor-like (ILCP-like) profile with increased expression of the ILC-3-associated TFs RORC and AHR, revealing a role for T-box TFs in maintaining mature NK cell phenotypes and an unexpected role of suppressing alternative ILC lineages. Our study reveals the critical importance of sustained EOMES and T-BET expression to orchestrate mature NK cell function and identity.

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Dispersive currents explain patterns of population connectivity in an ecologically and economically important fish

How to identify the drivers of population connectivity remains a fundamental question in ecology and evolution. Answering this question can be challenging in aquatic environments where dynamic lake and ocean currents coupled with high levels of dispersal and gene flow can decrease the utility of modern population genetic tools. To address this challenge, we used RAD-Seq to genotype 959 yellow perch (Perca flavescens), a species with an ~40-day pelagic larval duration (PLD), collected from 20 sites circumscribing Lake Michigan. We also developed a novel, integrative approach that couples detailed biophysical models with eco-genetic agent-based models to generate "predictive" values of genetic differentiation. By comparing predictive and empirical values of genetic differentiation, we estimated the relative contributions for known drivers of population connectivity (e.g., currents, behavior, PLD). For the main basin populations (i.e., the largest contiguous portion of the lake), we found that high gene flow led to low overall levels of genetic differentiation among populations (F _ST = 0.003). By far the best predictors of genetic differentiation were connectivity matrices that were derived from periods of time when there were strong and highly dispersive currents. Thus, these highly dispersive currents are driving the patterns of population connectivity in the main basin. We also found that populations from the northern and southern main basin are slightly divergent from one another, while those from Green Bay and the main basin are highly divergent (F _ST = 0.11). By integrating biophysical and eco-genetic models with genome-wide data, we illustrate that the drivers of population connectivity can be identified in high gene flow systems.

Assessing whole-exome sequencing data from undiagnosed Brazilian patients to improve the diagnostic yield of inborn errors of immunity

OBJECTIVES

Inborn error of immunity (IEI) comprises a broad group of inherited immunological disorders that usually display an overlap in many clinical manifestations challenging their diagnosis. The identification of disease-causing variants from whole-exome sequencing (WES) data comprises the gold-standard approach to ascertain IEI diagnosis. The efforts to increase the availability of clinically relevant genomic data for these disorders constitute an important improvement in the study of rare genetic disorders. This work aims to make available WES data of Brazilian patients' suspicion of IEI without a genetic diagnosis. We foresee a broad use of this dataset by the scientific community in order to provide a more accurate diagnosis of IEI disorders.

DATA DESCRIPTION

Twenty singleton unrelated patients treated at four different hospitals in the state of Rio de Janeiro, Brazil were enrolled in our study. Half of the patients were male with mean ages of 9 ± 3, while females were 12 ± 10 years old. The WES was performed in the Illumina NextSeq platform with at least 90% of sequenced bases with a minimum of 30 reads depth. Each sample had an average of 20,274 variants, comprising 116 classified as rare pathogenic or likely pathogenic according to American College of Medical Genetics and Genomics and the Association (ACMG) guidelines. The genotype-phenotype association was impaired by the lack of detailed clinical and laboratory information, besides the unavailability of molecular and functional studies which, comprise the limitations of this study. Overall, the access to clinical exome sequencing data is limited, challenging exploratory analyses and the understanding of genetic mechanisms underlying disorders. Therefore, by making these data available, we aim to increase the number of WES data from Brazilian samples despite contributing to the study of monogenic IEI-disorders.

Transcriptional and Histone Acetylation Changes Associated with CRE Elements Expose Key Factors Governing the Regulatory Circuit in the Early Stage of Huntington's Disease Models

Huntington's disease (HD) is a disorder caused by an abnormal expansion of trinucleotide CAG repeats within the huntingtin (Htt) gene. Under normal conditions, the CREB Binding Protein interacts with CREB elements and acetylates Lysine 27 of Histone 3 to direct the expression of several genes. However, mutant Htt causes depletion of CBP, which in turn induces altered histone acetylation patterns and transcriptional deregulation. Here, we have studied a differential expression analysis and H3K27ac variation in 4- and 6-week-old R6/2 mice as a model of juvenile HD. The analysis of differential gene expression and acetylation levels were integrated into Gene Regulatory Networks revealing key regulators involved in the altered transcription cascade. Our results show changes in acetylation and gene expression levels that are related to impaired neuronal development, and key regulators clearly defined in 6-week-old mice are proposed to drive the downstream regulatory cascade in HD. Here, we describe the first approach to determine the relationship among epigenetic changes in the early stages of HD. We determined the existence of changes in pre-symptomatic stages of HD as a starting point for early onset indicators of the progression of this disease.

A pilot study to disentangle the infant gut microbiota composition and identification of bacteria correlates with high fat mass

Background: At birth, the human intestine is colonized by a complex community of microorganisms known as gut microbiota. These complex microbial communities that inhabit the gut microbiota are thought to play a key role in maintaining host physiological homeostasis. For this reason, correct colonization of the gastrointestinal tract in the early stages of life could be fundamental for human health. Furthermore, alterations of the infant microbiota are correlated with the development of human inflammatory diseases and disorders. In this context, the possible relationships between intestinal microbiota and body composition during infancy are of great interest. Methods: In this study, we have performed a pilot study based on 16S rRNA gene profiling and metagenomic approaches on repeatedly measured data on time involving a cohort of 41 Italian newborns, which is aimed to investigate the possible correlation between body fat mass percentage (FM%) and the infant gut microbiota composition. Results and conclusion: The taxonomical analysis of the stool microbiota of each infant included in the cohort allowed the identification of a specific correlation between intestinal bacteria, such as Bifidobacterium and Veillonella, and the increase in FM%. Moreover, the analysis of the infant microbiome's metabolic capabilities suggested that the intestinal microbiome functionally impacts the human host and its possible influence on host physiology.

Multi-omic longitudinal study reveals immune correlates of clinical course among hospitalized COVID-19 patients

The IMPACC cohort, composed of >1,000 hospitalized COVID-19 participants, contains five illness trajectory groups (TGs) during acute infection (first 28 days), ranging from milder (TG1-3) to more severe disease course (TG4) and death (TG5). Here, we report deep immunophenotyping, profiling of >15,000 longitudinal blood and nasal samples from 540 participants of the IMPACC cohort, using 14 distinct assays. These unbiased analyses identify cellular and molecular signatures present within 72 h of hospital admission that distinguish moderate from severe and fatal COVID-19 disease. Importantly, cellular and molecular states also distinguish participants with more severe disease that recover or stabilize within 28 days from those that progress to fatal outcomes (TG4 vs. TG5). Furthermore, our longitudinal design reveals that these biologic states display distinct temporal patterns associated with clinical outcomes. Characterizing host immune responses in relation to heterogeneity in disease course may inform clinical prognosis and opportunities for intervention.

Characteristics of gut microbiota in patients with metabolic associated fatty liver disease

Metabolic associated fatty liver disease (MAFLD) is rising in incidence and is an increasingly common cause of cirrhosis and hepatocellular carcinoma (HCC). Alterations in the gut microbiota have been shown to correlate with the development and progression of MAFLD. However, little is known regarding differences in the gut microbiomes of MAFLD patients and healthy cohorts, and subgroups at the abnormal activity of hepatic enzymes in China. In this study, we enrolled 81 MAFLD patients and 25 healthy volunteers. The fecal microbiota was assessed using 16S rRNA gene sequencing and metagenomic sequencing. The results suggested that Ruminococcus obeum and Alistipes were most enriched in healthy individuals when compared with MAFLD patients. Microbe-set Enrichment Analysis (MSEA) results showed Dorea, Lactobacillus and Megasphaera are enriched in MAFLD group. We also found that Alistipes has negatively related to serum glucose (GLU), gamma-glutamyl transferase (GGT), and alanine aminotransferase (ALT). Moreover, the abundance of Dorea was found to be significantly overrepresented in the MAFLD patients and the degree of enrichment increased with the increasing abnormal liver enzyme. An increase in Dorea, combined with decreases in Alistipes appears to be characteristic of MAFLD patients. Further study of microbiota may provide a novel insight into the pathogenesis of MAFLD as well as a novel treatment strategy.

The compact genome of the sponge Oopsacas minuta (Hexactinellida) is lacking key metazoan core genes

BACKGROUND

Explaining the emergence of the hallmarks of bilaterians is a central focus of evolutionary developmental biology-evodevo-and evolutionary genomics. For this purpose, we must both expand and also refine our knowledge of non-bilaterian genomes, especially by studying early branching animals, in particular those in the metazoan phylum Porifera.

RESULTS

We present a comprehensive analysis of the first whole genome of a glass sponge, Oopsacas minuta, a member of the Hexactinellida. Studying this class of sponge is evolutionary relevant because it differs from the three other Porifera classes in terms of development, tissue organization, ecology, and physiology. Although O. minuta does not exhibit drastic body simplifications, its genome is among the smallest of animal genomes sequenced so far, and surprisingly lacks several metazoan core genes (including Wnt and several key transcription factors). Our study also provides the complete genome of a symbiotic Archaea dominating the associated microbial community: a new Thaumarchaeota species.

CONCLUSIONS

The genome of the glass sponge O. minuta differs from all other available sponge genomes by its compactness and smaller number of encoded proteins. The unexpected loss of numerous genes previously considered ancestral and pivotal for metazoan morphogenetic processes most likely reflects the peculiar syncytial tissue organization in this group. Our work further documents the importance of convergence during animal evolution, with multiple convergent evolution of septate-like junctions, electrical-signaling and multiciliated cells in metazoans.

Nitrate Reductase NarGHJI Modulates Virulence via Regulation of agr Expression in Methicillin-Resistant Staphylococcus aureus Strain USA300 LAC

Staphylococcus aureus is a pathogenic bacterium with a widespread distribution that can cause diverse severe diseases. The membrane-bound nitrate reductase NarGHJI serves respiratory function. However, little is known about its contribution to virulence. In this study, we demonstrated that narGHJI disruption results in the downregulation of virulence genes (e.g., RNAIII, agrBDCA, hla, psmα, and psmβ) and reduces the hemolytic activity of the methicillin-resistant S. aureus (MRSA) strain USA300 LAC. Moreover, we provided evidence that NarGHJI participates in regulating host inflammatory response. A mouse model of subcutaneous abscess and Galleria mellonella survival assay demonstrated that the ΔnarG mutant was significantly less virulent than the wild type. Interestingly, NarGHJI contributes to virulence in an agr-dependent manner, and the role of NarGHJI differs between different S. aureus strains. Our study highlights the novel role of NarGHJI in regulating virulence, thereby providing a new theoretical reference for the prevention and control of S. aureus infection. IMPORTANCE Staphylococcus aureus is a notorious pathogen that poses a great threat to human health. The emergence of drug-resistant strains has significantly increased the difficulty of preventing and treating S. aureus infection and enhanced the pathogenic ability of the bacterium. This indicates the importance of identifying novel pathogenic factors and revealing the regulatory mechanisms through which they regulate virulence. The nitrate reductase NarGHJI is mainly involved in bacterial respiration and denitrification, which can enhance bacterial survival. We demonstrated that narGHJI disruption results in the downregulation of the agr system and agr-dependent virulence genes, suggesting that NarGHJI participates in the regulation of S. aureus virulence in an agr-dependent manner. Moreover, the regulatory approach is strain specific. This study provides a new theoretical reference for the prevention and control of S. aureus infection and reveals new targets for the development of therapeutic drugs.

Combined direct/indirect detection allows identification of DNA termini in diverse sequencing datasets and supports a multiple-initiation-site model for HIV plus-strand synthesis

Replication of genetic material involves the creation of characteristic termini. Determining these termini is important to refine our understanding of the mechanisms involved in maintaining the genomes of cellular organisms and viruses. Here we describe a computational approach combining direct and indirect readouts to detect termini from next-generation short-read sequencing. While a direct inference of termini can come from mapping the most prominent start positions of captured DNA fragments, this approach is insufficient in cases where the DNA termini are not captured, whether for biological or technical reasons. Thus, a complementary (indirect) approach to terminus detection can be applied, taking advantage of the imbalance in coverage between forward and reverse sequence reads near termini. A resulting metric ("strand bias") can be used to detect termini even where termini are naturally blocked from capture or ends are not captured during library preparation (e.g., in tagmentation-based protocols). Applying this analysis to datasets where known DNA termini are present, such as from linear double-stranded viral genomes, yielded distinct strand bias signals corresponding to these termini. To evaluate the potential to analyze a more complex situation, we applied the analysis to examine DNA termini present early after HIV infection in a cell culture model. We observed both the known termini expected based on standard models of HIV reverse transcription (the U5-right-end and U3-left-end termini) as well as a signal corresponding to a previously described additional initiation site for plus-strand synthesis (cPPT [central polypurine tract]). Interestingly, we also detected putative terminus signals at additional sites. The strongest of these are a set that share several characteristics with the previously characterized plus-strand initiation sites (the cPPT and 3' PPT [polypurine tract] sites): (i) an observed spike in directly captured cDNA ends, an indirect terminus signal evident in localized strand bias, (iii) a preference for location on the plus-strand, (iv) an upstream purine-rich motif, and (v) a decrease in terminus signal at late time points after infection. These characteristics are consistent in duplicate samples in two different genotypes (wild type and integrase-lacking HIV). The observation of distinct internal termini associated with multiple purine-rich regions raises a possibility that multiple internal initiations of plus-strand synthesis might contribute to HIV replication.

Impaired histone inheritance promotes tumor progression

Faithful inheritance of parental histones is essential to maintain epigenetic information and cellular identity during cell division. Parental histones are evenly deposited onto the replicating DNA of sister chromatids in a process dependent on the MCM2 subunit of DNA helicase. However, the impact of aberrant parental histone partition on human disease such as cancer is largely unknown. In this study, we construct a model of impaired histone inheritance by introducing MCM2-2A mutation (defective in parental histone binding) in MCF-7 breast cancer cells. The resulting impaired histone inheritance reprograms the histone modification landscapes of progeny cells, especially the repressive histone mark H3K27me3. Lower H3K27me3 levels derepress the expression of genes associated with development, cell proliferation, and epithelial to mesenchymal transition. These epigenetic changes confer fitness advantages to some newly emerged subclones and consequently promote tumor growth and metastasis after orthotopic implantation. In summary, our results indicate that impaired inheritance of parental histones can drive tumor progression.

ddRAD sequencing based genotyping of six indigenous dairy cattle breeds of India to infer existing genetic diversity and population structure

The present investigation aimed to identify genome wide SNPs and to carry out diversity and population structure study using ddRAD-seq based genotyping of 58 individuals of six indigenous milch cattle breeds (Bos indicus) such as Sahiwal, Gir, Rathi, Tharparkar, Red Sindhi and Kankrej of India. A high percentage of reads (94.53%) were mapped to the Bos taurus (ARS-UCD1.2) reference genome assembly. Following filtration criteria, a total of 84,027 high quality SNPs were identified across the genome of 6 cattle breeds with the highest number of SNPs observed in Gir (34,743), followed by Red Sindhi (13,092), Kankrej (12,812), Sahiwal (8956), Tharparkar (7356) and Rathi (7068). Most of these SNPs were distributed in the intronic regions (53.87%) followed by intergenic regions (34.94%) while only 1.23% were located in the exonic regions. Together with analysis of nucleotide diversity (π = 0.373), Tajima's D (D value ranging from - 0.295 to 0.214), observed heterozygosity (HO ranging from 0.464 to 0.551), inbreeding coefficient (FIS ranging from - 0.253 to 0.0513) suggested for the presence of sufficient within breed diversity in the 6 major milch breeds of India. The phylogenetic based structuring, principal component and admixture analysis revealed genetic distinctness as well as purity of almost all of the 6 cattle breeds. Overall, our strategy has successfully identified thousands of high-quality genome wide SNPs that will further enrich the Bos indicus representation basic information about genetic diversity and structure of 6 major Indian milch cattle breeds which should have implications for better management and conservation of valuable indicine cattle diversity.

A multiple super-enhancer region establishes inter-TAD interactions and controls Hoxa function in cranial neural crest

Enhancer-promoter interactions preferentially occur within boundary-insulated topologically associating domains (TADs), limiting inter-TAD interactions. Enhancer clusters in linear proximity, termed super-enhancers (SEs), ensure high target gene expression levels. Little is known about SE topological regulatory impact during craniofacial development. Here, we identify 2232 genome-wide putative SEs in mouse cranial neural crest cells (CNCCs), 147 of which target genes establishing CNCC positional identity during face formation. In second pharyngeal arch (PA2) CNCCs, a multiple SE-containing region, partitioned into Hoxa Inter-TAD Regulatory Element 1 and 2 (HIRE1 and HIRE2), establishes long-range inter-TAD interactions selectively with Hoxa2, that is required for external and middle ear structures. HIRE2 deletion in a Hoxa2 haploinsufficient background results in microtia. HIRE1 deletion phenocopies the full homeotic Hoxa2 knockout phenotype and induces PA3 and PA4 CNCC abnormalities correlating with Hoxa2 and Hoxa3 transcriptional downregulation. Thus, SEs can overcome TAD insulation and regulate anterior Hoxa gene collinear expression in a CNCC subpopulation-specific manner during craniofacial development.

The Role of Repetitive Sequences in Repatterning of Major Ribosomal DNA Clusters in Lepidoptera

Genes for major ribosomal RNAs (rDNA) are present in multiple copies mainly organized in tandem arrays. The number and position of rDNA loci can change dynamically and their repatterning is presumably driven by other repetitive sequences. We explored a peculiar rDNA organization in several representatives of Lepidoptera with either extremely large or numerous rDNA clusters. We combined molecular cytogenetics with analyses of second- and third-generation sequencing data to show that rDNA spreads as a transcription unit and reveal association between rDNA and various repeats. Furthermore, we performed comparative long read analyses among the species with derived rDNA distribution and moths with a single rDNA locus, which is considered ancestral. Our results suggest that satellite arrays, rather than mobile elements, facilitate homology-mediated spread of rDNA via either integration of extrachromosomal rDNA circles or ectopic recombination. The latter arguably better explains preferential spread of rDNA into terminal regions of lepidopteran chromosomes as efficiency of ectopic recombination depends on the proximity of homologous sequences to telomeres.

Quality control of protein synthesis in the early elongation stage

In the early stage of bacterial translation, peptidyl-tRNAs frequently dissociate from the ribosome (pep-tRNA drop-off) and are recycled by peptidyl-tRNA hydrolase. Here, we establish a highly sensitive method for profiling of pep-tRNAs using mass spectrometry, and successfully detect a large number of nascent peptides from pep-tRNAs accumulated in Escherichia coli pthts strain. Based on molecular mass analysis, we found about 20% of the peptides bear single amino-acid substitutions of the N-terminal sequences of E. coli ORFs. Detailed analysis of individual pep-tRNAs and reporter assay revealed that most of the substitutions take place at the C-terminal drop-off site and that the miscoded pep-tRNAs rarely participate in the next round of elongation but dissociate from the ribosome. These findings suggest that pep-tRNA drop-off is an active mechanism by which the ribosome rejects miscoded pep-tRNAs in the early elongation, thereby contributing to quality control of protein synthesis after peptide bond formation.

Altered gut microbiome profile in patients with knee osteoarthritis

Introduction

Osteoarthritis (OA) is a kind of chronic, degenerative disorder with unknown causes. In this study, we aimed to improve our understanding of the gut microbiota profile in patients with knee OA.

Methods

16S rDNA gene sequencing was performed to detect the gut microbiota in fecal samples collected from the patients with OA (n = 32) and normal control (NC, n = 57). Then the metagenomic sequencing was used to identify the genes or functions linked with gut microbial changes at the species level in the fecal samples from patients with OA and NC groups.

Results

The Proteobacteria was identified as dominant bacteria in OA group. We identified 81 genera resulted significantly different in abundance between OA and NC. The abundance of Agathobacter, Ruminococcus, Roseburia, Subdoligranulum, and Lactobacillus showed significant decrease in the OA compared to the NC. The abundance of genera Prevotella_7, Clostridium, Flavonifractor and Klebsiella were increasing in the OA group, and the families Lactobacillaceae, Christensenellaceae, Clostridiaceae_1 and Acidaminococcaceae were increasing in the NC. The metagenomic sequencing showed that the abundance of Bacteroides stercoris, Bacteroides vulgatus and Bacteroides uniformis at the species level were significantly decreasing in the OA, and the abundance of Escherichia coli, Klebsiella pneumoniae, Shigella flexneri and Streptococcus salivarius were significantly increased in OA.

Discussion

The results of our study interpret a comprehensive profile of the gut microbiota in patients with knee OA and offer the evidence that the cartilage-gut-microbiome axis could play a crucial role in underlying the mechanisms and pathogenesis of OA.

Differences in gut metagenomes between dairy workers and community controls: a cross-sectional study

Background

As a nexus of routine antibiotic use and zoonotic pathogen presence, the livestock farming environment is a potential hotspot for the emergence of zoonotic diseases and antibiotic resistant bacteria. Livestock can further facilitate disease transmission by serving as intermediary hosts for pathogens as they undergo evolution prior to a spillover event. In light of this, we are interested in characterizing the microbiome and resistome of dairy workers, whose exposure to the livestock farming environment places them at risk for facilitating community transmission of antibiotic resistant genes and emerging zoonotic diseases.

Results

Using shotgun sequencing, we investigated differences in the taxonomy, diversity and gene presence of the human gut microbiome of 10 dairy farm workers and 6 community controls, supplementing these samples with additional publicly available gut metagenomes. We observed greater abundance of tetracycline resistance genes and prevalence of cephamycin resistance genes in dairy workers' metagenomes, and lower average gene diversity. We also found evidence of commensal organism association with plasmid-mediated tetracycline resistance genes in both dairy workers and community controls (including Faecalibacterium prausnitzii, Ligilactobacillus animalis, and Simiaoa sunii). However, we did not find significant differences in the prevalence of resistance genes or virulence factors overall, nor differences in the taxonomic composition of dairy worker and community control metagenomes.

Conclusions

This study presents the first metagenomics analysis of United States dairy workers, providing insights into potential risks of exposure to antibiotics and pathogens in animal farming environments. Previous metagenomic studies of livestock workers in China and Europe have reported increased abundance and carriage of antibiotic resistance genes in livestock workers. While our investigation found no strong evidence for differences in the abundance or carriage of antibiotic resistance genes and virulence factors between dairy worker and community control gut metagenomes, we did observe patterns in the abundance of tetracycline resistance genes and the prevalence of cephamycin resistance genes that is consistent with previous work.

MicroRNAs from Holarrhena pubescens stems: Identification by small RNA Sequencing and their Potential Contribution to Human Gene Targets

Holarrhena pubescens is an effective medicinal plant from the Apocynaceae family, widely distributed over the Indian subcontinent and extensively used by Ayurveda and ethno-medicine systems without apparent side effects. We postulated that miRNAs, endogenous non-coding small RNAs that regulate gene expression at the post-transcriptional level, may, after ingestion into the human body, contribute to the medicinal properties of plants of this species by inducing regulated human gene expression to modulate. However, knowledge is scarce about miRNA in Holarrhena. In addition, to test the hypothesis on the potential pharmacological properties of miRNA, we performed a high-throughput sequencing analysis using the Next Generation Sequencing Illumina platform; 42,755,236 raw reads have been generated from H. pubescens stems from a library of small RNA isolated, identifying 687 known and 50 new miRNAs led. The novel H. pubescens miRNAs were predicted to regulate specific human genes, and subsequent annotations of gene functions suggested a possible role in various biological processes and signaling pathways, such as Wnt, MAPK, PI3K-Akt, and AMPK signaling pathways and endocytosis. The association of these putative targets with many diseases, including cancer, congenital malformations, nervous system disorders, and cystic fibrosis, has been demonstrated. The top hub proteins STAT3, MDM2, GSK3B, NANOG, IGF1, PRKCA, SNAP25, SRSF1, HTT, and SNCA show their interaction with human diseases, including cancer and cystic fibrosis. To our knowledge, this is the first report of uncovering H. pubescens miRNAs based on high-throughput sequencing and bioinformatics analysis. This study has provided new insight into a potential cross-species control of human gene expression. The potential for miRNA transfer should be evaluated as one possible mechanism of action to account for the beneficial properties of this valuable species.

Substantial rDNA copy number reductions alter timing of development and produce variable tissue-specific phenotypes in C. elegans

The genes that encode ribosomal RNAs are present in several hundred copies in most eukaryotes. These vast arrays of repetitive ribosomal DNA (rDNA) have been implicated not just in ribosome biogenesis, but also aging, cancer, genome stability, and global gene expression. rDNA copy number is highly variable among and within species; this variability is thought to associate with traits relevant to human health and disease. Here we investigate the phenotypic consequences of multicellular life at the lower bounds of rDNA copy number. We use the model Caenorhabditis elegans, which has previously been found to complete embryogenesis using only maternally provided ribosomes. We find that individuals with rDNA copy number reduced to ∼5% of wild type are capable of further development with variable penetrance. Such individuals are sterile and exhibit severe morphological defects, particularly in post-embryonically dividing tissues such as germline and vulva. Developmental completion and fertility are supported by an rDNA copy number ∼10% of wild type, with substantially delayed development. Worms with rDNA copy number reduced to ∼33% of wild type display a subtle developmental timing defect that was absent in worms with higher copy numbers. Our results support the hypothesis that rDNA requirements vary across tissues and indicate that the minimum rDNA copy number for fertile adulthood is substantially less than the lowest naturally observed total copy number. The phenotype of individuals with severely reduced rDNA copy number is highly variable in penetrance and presentation, highlighting the need for continued investigation into the biological consequences of rDNA copy number variation.

Characterization of a novel picornavirus prevalent in experimental rabbits (Oryctolagus cuniculus)

Here, using viral metagenomic method a novel picornavirus (named UJS-2019picorna, GenBank accession number OP821762) was discovered in fecal and blood samples of experimental rabbits (Oryctolagus cuniculus). The complete genome size of UJS-2019picorna is 7832 bp excluding the poly(A)-tail, with GC content of 44.00% and a nucleotide composition of 28.0% A, 28.0% U, 21.5% G, and 22.5% C. The viral genome has a typical picornavirus organization pattern from the 5'-3' direction: VPg-5' UTR-(L)-P1, (VP4-VP2-VP3-VP1)-P2, (2 A-2B-2C)-P3, (3 A-3B-3C-3D)-3' UTR-poly(A). The P1 region of UJS-2019picorna is related to Erbovirus with amino acid identity of 37.31%, while the P2 and P3 regions are the closest to Bopivirus with amino acid identity of 35.66%-39.53%. According to the Picornaviridae Study Group guidelines, UJS-2019picorna should be presumed to be a new genus belonging to the Picornaviridae family. Epidemiologic study revealed that this novel picornavirus was prevalent in a cohort of experimental rabbits, with prevalence rate of 23.68% (9/38) in feces and 18.4% (7/38) in blood samples. Further work is required to elucidate whether this virus is pathogenic to rabbits and whether it has influence on studies using rabbits as experimental animal.

RAS and PP2A activities converge on epigenetic gene regulation

RAS-mediated human cell transformation requires inhibition of the tumor suppressor protein phosphatase 2A (PP2A). However, the phosphoprotein targets and cellular processes in which RAS and PP2A activities converge in human cancers have not been systematically analyzed. Here, we discover that phosphosites co-regulated by RAS and PP2A are enriched on proteins involved in epigenetic gene regulation. As examples, RAS and PP2A co-regulate the same phosphorylation sites on HDAC1/2, KDM1A, MTA1/2, RNF168, and TP53BP1. We validate RAS- and PP2A-elicited regulation of HDAC1/2 chromatin recruitment, of RNF168-TP53BP1 interaction, and of gene expression. Consistent with their known synergistic effects in cancer, RAS activation and PP2A inhibition resulted in epigenetic reporter derepression and activation of oncogenic transcription. Transcriptional derepression by PP2A inhibition was associated with an increase in euchromatin and a decrease in global DNA methylation. Collectively, the results indicate that epigenetic protein complexes constitute a significant point of convergence for RAS hyperactivity and PP2A inhibition in cancer. Furthermore, the work provides an important resource for future studies focusing on phosphoregulation of epigenetic gene regulation in cancer and in other RAS/PP2A-regulated cellular processes.

Complementary gene regulation by NRF1 and NRF2 protects against hepatic cholesterol overload

Hepatic cholesterol overload promotes steatohepatitis. Insufficient understanding of liver stress defense impedes therapy development. Here, we elucidate the role of stress defense transcription factors, nuclear factor erythroid 2 related factor-1 (NRF1) and -2 (NRF2), in counteracting cholesterol-linked liver stress. Using a diet that increases liver cholesterol storage, expression profiles and phenotypes of liver from mice with hepatocyte deficiency of NRF1, NRF2, or both are compared with controls, and chromatin immunoprecipitation sequencing is undertaken to identify target genes. Results show NRF1 and NRF2 co-regulate genes that eliminate cholesterol and mitigate inflammation and oxidative damage. Combined deficiency, but not deficiency of either alone, results in severe steatohepatitis, hepatic cholesterol overload and crystallization, altered bile acid metabolism, and decreased biliary cholesterol. Moreover, therapeutic effects of NRF2-activating drug bardoxolone require NRF1 and are supplemented by NRF1 overexpression. Thus, we discover complementary gene programming by NRF1 and NRF2 that counteract cholesterol-associated fatty liver disease progression.

Differences in Bacterial Small RNAs in Stool Samples from Hypercholesterolemic and Normocholesterolemic Subjects

Cholesterol metabolism is important at the physiological level as well as in several diseases, with small RNA being an element to consider in terms of its epigenetic control. Thus, the aim of this study was to identify differences between bacterial small RNAs present at the gut level in hypercholesterolemic and normocholesterolemic individuals. Twenty stool samples were collected from hypercholesterolemic and normocholesterolemic subjects. RNA extraction and small RNA sequencing were performed, followed by bioinformatics analyses with BrumiR, Bowtie 2, BLASTn, DESeq2, and IntaRNA, after the filtering of the reads with fastp. In addition, the prediction of secondary structures was obtained with RNAfold WebServer. Most of the small RNAs were of bacterial origin and presented a greater number of readings in normocholesterolemic participants. The upregulation of small RNA ID 2909606 associated with Coprococcus eutactus (family Lachnospiraceae) was presented in hypercholesterolemic subjects. In addition, a positive correlation was established between small RNA ID 2149569 from the species Blautia wexlerae and hypercholesterolemic subjects. Other bacterial and archaeal small RNAs that interacted with the LDL receptor (LDLR) were identified. For these sequences, the prediction of secondary structures was also obtained. There were significant differences in bacterial small RNAs associated with cholesterol metabolism in hypercholesterolemic and normocholesterolemic participants.

Characterization of a novel papillomavirus identified from a whale (Delphinapterus leucas) pharyngeal metagenomic library

Here, using viral metagenomic method, a novel whale papillomavirus (temporarily named wPV, GenBank accession number OP856597) was discovered in a whale (Delphinapterus leucas) pharyngeal metagenomic library. The complete genome size of wPV is 7179 bp, with GC content of 54.4% and a nucleotide composition of 23.4% A, 22.3% T, 28.4% G, and 25.9% C. The viral genome has a typical papillomavirus organization pattern, and five ORFs were predicted, including two late genes encoding L1 and L2, and three early genes encoding E1, E2, and E6. Pairwise sequence comparison and phylogenetic analysis based on the L1 gene sequence indicated that wPV may be a novel species within genus Dyodeltapapillomavirus. In addition, the E2 region of wPV was predicted to have a potential recombination event. The discovery of this novel papillomavirus increases our understanding of the viral ecology of marine mammals, providing insights into possible future infectious diseases.

Bacteria require phase separation for fitness in the mammalian gut

Therapeutic manipulation of the gut microbiota holds great potential for human health. The mechanisms bacteria use to colonize the gut therefore present valuable targets for clinical intervention. We now report that bacteria use phase separation to enhance fitness in the mammalian gut. We establish that the intrinsically disordered region (IDR) of the broadly and highly conserved transcription termination factor Rho is necessary and sufficient for phase separation in vivo and in vitro in the human commensal Bacteroides thetaiotaomicron. Phase separation increases transcription termination by Rho in an IDR-dependent manner. Moreover, the IDR is critical for gene regulation in the gut. Our findings expose phase separation as vital for host-commensal bacteria interactions and relevant for novel clinical applications.

Conservation and divergence of canonical and non-canonical imprinting in murids

BACKGROUND

Genomic imprinting affects gene expression in a parent-of-origin manner and has a profound impact on complex traits including growth and behavior. While the rat is widely used to model human pathophysiology, few imprinted genes have been identified in this murid. To systematically identify imprinted genes and genomic imprints in the rat, we use low input methods for genome-wide analyses of gene expression and DNA methylation to profile embryonic and extraembryonic tissues at allele-specific resolution.

RESULTS

We identify 14 and 26 imprinted genes in these tissues, respectively, with 10 of these genes imprinted in both tissues. Comparative analyses with mouse reveal that orthologous imprinted gene expression and associated canonical DNA methylation imprints are conserved in the embryo proper of the Muridae family. However, only 3 paternally expressed imprinted genes are conserved in the extraembryonic tissue of murids, all of which are associated with non-canonical H3K27me3 imprints. The discovery of 8 novel non-canonical imprinted genes unique to the rat is consistent with more rapid evolution of extraembryonic imprinting. Meta-analysis of novel imprinted genes reveals multiple mechanisms by which species-specific imprinted expression may be established, including H3K27me3 deposition in the oocyte, the appearance of ZFP57 binding motifs, and the insertion of endogenous retroviral promoters.

CONCLUSIONS

In summary, we provide an expanded list of imprinted loci in the rat, reveal the extent of conservation of imprinted gene expression, and identify potential mechanisms responsible for the evolution of species-specific imprinting.