Identification of BAY61-3606 Derivatives With Improved Activity in Splicing Modulation That Induces Inclusion of Cassette Exons Similar to the Splicing Factor 3B Subunit 1 Mutation.

Splicing modulation by a small compound offers therapeutic potential for diseases caused by splicing abnormality. However, only a few classes of compounds that can modulate splicing have been identified. We previously identified BAY61-3606, a multiple kinase inhibitor, as a compound that relaxes the splicing fidelity at the 3' splice site recognition. We have also reported the synthesis of derivatives of BAY61-3606. In this study, we tested those compounds for their splicing modulation capabilities and identified two contrasting compounds. These compounds were further investigated for their effects on the whole transcriptome, and analysis of changes in transcription and splicing revealed that the highly active derivative in the splicing reporter assay also showed significantly higher activity in modulating the splicing of endogenously expressed genes. Particularly, cassette exon inclusion was highly upregulated by this compound, and clustering analysis revealed that these effects resembled those in splicing factor 3b subunit 1 (SF3B1) K700E mutant cells but contrasted with those of the splicing inhibitor H3B-8800. Additionally, a group of serine/arginine-rich (SR) protein genes was identified as representatively affected, likely via modulation of poison exon inclusion. This finding could guide further analysis of the mode of action of these compounds on splicing, which could be valuable for developing drugs for diseases associated with splicing abnormalities.

Multi-omics analysis using antibody-based in situ biotinylation technique suggests the mechanism of Cajal body formation.

Membrane-less subcellular compartments play important roles in various cellular functions. Although techniques exist to identify components of cellular bodies, a comprehensive method for analyzing both static and dynamic states has not been established. Here, we apply an antibody-based in situ biotinylation proximity-labeling technique to identify components of static and dynamic nuclear bodies. Using this approach, we comprehensively identify DNA, RNA, and protein components of Cajal bodies (CBs) and then clarify their interactome. By inhibiting transcription, we capture dynamic changes in CBs. Our analysis reveals that nascent small nuclear RNAs (snRNAs) transcribed in CBs contribute to CB formation by assembling RNA-binding proteins, including frontotemporal dementia-related proteins, RNA-binding motif proteins, and heterogeneous nuclear ribonucleoproteins.

Mettl1-dependent m7G tRNA modification is essential for maintaining spermatogenesis and fertility in Drosophila melanogaster.

Modification of guanosine to N7-methylguanosine (m7G) in the variable loop region of tRNA is catalyzed by the METTL1/WDR4 heterodimer and stabilizes target tRNA. Here, we reveal essential functions of Mettl1 in Drosophila fertility. Knockout of Mettl1 (Mettl1-KO) causes no major effect on the development of non-gonadal tissues, but abolishes the production of elongated spermatids and mature sperm, which is fully rescued by expression of a Mettl1-transgene, but not a catalytic-dead Mettl1 transgene. This demonstrates that Mettl1-dependent m7G is required for spermatogenesis. Mettl1-KO results in a loss of m7G modification on a subset of tRNAs and decreased tRNA abundance. Ribosome profiling shows that Mettl1-KO led to ribosomes stalling at codons decoded by tRNAs that were reduced in abundance. Mettl1-KO also significantly reduces the translation efficiency of genes involved in elongated spermatid formation and sperm stability. Germ cell-specific expression of Mettl1 rescues disrupted m7G tRNA modification and tRNA abundance in Mettl1-KO testes but not in non-gonadal tissues. Ribosome stalling is much less detectable in non-gonadal tissues than in Mettl1-KO testes. These findings reveal a developmental role for m7G tRNA modification and indicate that m7G modification-dependent tRNA abundance differs among tissues.

Mitochondrial phylogeny and distribution of cytoplasmic male sterility-associated genes in Beta vulgaris.

Cytoplasmic male sterility (CMS) is a mitochondrial-encoded trait that confers reproductive defects in males but not in females or any vegetative function. Why CMS is so often found in plants should be investigated from the viewpoint of mitochondrial phylogeny. Beta vulgaris, including the wild subspecies maritima and cultivated subspecies vulgaris (e.g., sugar beet), is known to be mitochondrially polymorphic, from which multiple CMS mitochondria have been found, but their evolutionary relationship has been obscure. We first refined the B. vulgaris reference mitochondrial genome to conduct a more accurate phylogenetic study. We identified mitochondrial single-nucleotide polymorphic sites from 600 B. vulgaris accessions. Principal component analysis, hierarchical clustering analysis, and creation of a phylogenetic tree consistently suggested that B. vulgaris mitochondria can be classified into several groups whose geographical distribution tends to be biased toward either the Atlantic or Mediterranean coasts. We examined the distribution of CMS-associated mitochondrial genes from Owen, E- and G-type CMS mitochondria. About one-third of cultivated beets had Owen-type CMS, which reflects the prevalence of using Owen-type CMS in hybrid breeding. Occurrence frequencies for each of the three CMS genes in wild beet were less than 4%. CMS genes were tightly associated with specific mitochondrial groups that are phylogenetically distinct, suggesting their independent origin. However, homologous sequences of the Owen type CMS gene occurred in several different mitochondrial groups, for which an intricate explanation is necessary. Whereas the origin of cultivated beet had been presumed to be Greece, we found an absence of Owen-type mitochondria in Greek accessions.

Exploring Concurrent Approach for Respiratory Epidemiological Surveillance and Symptom Screening (CARES): a new strategy for preventing respiratory infection outbreaks in long-term care facilities.

OBJECTIVE: The Concurrent Approach for Respiratory Epidemiological Surveillance and Symptom Screening (CARES) combines symptom screening to prevent external respiratory infections by managing staff and visitor health and surveillance to monitor the infection rates within the facility and take timely decisions on activity programs based on outbreak prevalence. This study examines the efficacy of the CARES strategy in preventing respiratory infection outbreaks in long-term care facilities. DESIGN: Prospective cohort study utilizing historical controls including an intervention and control group. SETTING: This study was conducted in two long-term care facility wards (total: 110 beds) in Japan. PATIENTS: We enrolled patients aged 1-72 years from the target ward with severe intellectual and physical disabilities from October 1, 2018, to March 31, 2019, and from October 1, 2017, to March 31, 2018, for the intervention (n = 104) and control (n = 98) groups, respectively. The study included all admitted patients and excluded those hospitalized or discharged during the study period. INTERVENTION: The total number of days that activity programs, new admissions, and visitations were cancelled in the two groups was compared before and after the introduction of CARES. RESULTS: CARES reduced the duration of new admission cancellations and visitation cancellations by 16 and 23 days, respectively (α = 0.1, P value < 0.001). Additionally, the maximum duration of activity program cancellations was reduced by 2 days. Furthermore, five cases of presenteeism were prevented. CONCLUSIONS: CARES improves patients' quality of life by continuation of activity programs, new patient admission, and ongoing visitations.

Induction of male-like mandibles in XX individuals of a stag beetle by gene knockdown of a feminizer gene transformer.

Males and females share most of the genome, but many animals show different phenotypes between the sexes, known as sexual dimorphism. Many insect species show extreme sexual dimorphism, including beetles with "weapon traits" represented by extremely developed horns and mandibles. Existing studies of sex-specific development of beetle weapon traits suggest that sex-specific gene expression plays an important role. On the other hand, contributions of the Y-chromosome, which may potentially carry genes necessary for male development, to weapon trait expression have not been examined. In holometabolous insects, including beetles, the feminizing gene transformer (tra) is roughly conserved in its feminizing function. Only females express a functional isoform of Tra, which causes female differentiation. Knocking down tra in females leads to male tissue differentiation, enabling us to analyze male phenotypes in individuals lacking a Y-chromosome (XX-males). In this study, we investigate whether the Y-chromosome is necessary for stag beetles to express male-specific weapon traits by comparing tra-knockdown-induced XX-males with natural XY males. We show that XX-males could express weapons (enlarged mandibles) as in XY-males. These results suggest that the Y-chromosome does not have a major role in weapon trait expression in this species.

Genome-wide analysis of histone modifications can contribute to the identification of candidate cis-regulatory regions in the threespine stickleback fish.

BACKGROUND: Cis-regulatory mutations often underlie phenotypic evolution. However, because identifying the locations of promoters and enhancers in non-coding regions is challenging, we have fewer examples of identified causative cis-regulatory mutations that underlie naturally occurring phenotypic variations than of causative amino acid-altering mutations. Because cis-regulatory elements have epigenetic marks of specific histone modifications, we can detect cis-regulatory elements by mapping and analyzing them. Here, we investigated histone modifications and chromatin accessibility with cleavage under targets and tagmentation (CUT&Tag) and assay for transposase-accessible chromatin-sequencing (ATAC-seq). RESULTS: Using the threespine stickleback (Gasterosteus aculeatus) as a model, we confirmed that the genes for which nearby regions showed active marks, such as H3K4me1, H3K4me3, and high chromatin accessibility, were highly expressed. In contrast, the expression levels of genes for which nearby regions showed repressive marks, such as H3K27me3, were reduced, suggesting that our chromatin analysis protocols overall worked well. Genomic regions with peaks of histone modifications showed higher nucleotide diversity within and between populations. By comparing gene expression in the gills of the marine and stream ecotypes, we identified several insertions and deletions (indels) with transposable element fragments in the candidate cis-regulatory regions. CONCLUSIONS: Thus, mapping and analyzing histone modifications can help identify cis-regulatory elements and accelerate the identification of causative mutations in the non-coding regions underlying naturally occurring phenotypic variations.

Effects of subchronic and mild social defeat stress on the intestinal microbiota and fecal bile acid composition in mice.

The gut microbiota plays a crucial role in both the pathogenesis and alleviation of host depression by modulating the brain-gut axis. We have developed a murine model of human depression called the subchronic and mild social defeat stress (sCSDS) model, which impacts not only behavior but also the host gut microbiota and gut metabolites, including bile acids. In this study, we utilized liquid chromatography/mass spectrometry (LC/MS) to explore the effects of sCSDS on the mouse fecal bile acid profile. sCSDS mice exhibited significantly elevated levels of deoxycholic acid (DCA) and lithocholic acid (LCA) in fecal extracts, leading to a notable increase in total bile acids and 7α-dehydroxylated secondary bile acids. Consequently, a noteworthy negative correlation was identified between the abundances of DCA and LCA and the social interaction score, an indicator of susceptibility in stressed mice. Furthermore, analysis of the colonic microbiome unveiled a negative correlation between the abundance of CDCA and Turicibacter. Additionally, DCA and LCA exhibited positive correlations with Oscillospiraceae and Lachnospiraceae but negative correlations with the Eubacterium coprostanoligenes group. These findings suggest that sCSDS impacts the bidirectional interaction between the gut microbiota and bile acids and is associated with reduced social interaction, a behavioral indicator of susceptibility in stressed mice.

Metagenome-assembled genome sequence of Spiroplasma phoeniceum, assembled from the hindgut of Locusta migratoria, a migration grasshopper species.

Spiroplasma phoeniceum is a plant pathogen and a mesophilic microaerophile. Here, we report the metagenome-assembled genome (MAG) sequence of S. phoeniceum binned from hindgut contents of the wild-type male Locusta migratoria, a grasshopper species. The MAG sequence comprises 1,059,205 bp in 91 contigs with a 26.3% of GC content.

Odd-Paired is Involved in Morphological Divergence of Snail-Feeding Beetles.

Body shape and size diversity and their evolutionary rates correlate with species richness at the macroevolutionary scale. However, the molecular genetic mechanisms underlying the morphological diversification across related species are poorly understood. In beetles, which account for one-fourth of the known species, adaptation to different trophic niches through morphological diversification appears to have contributed to species radiation. Here, we explored the key genes for the morphological divergence of the slender to stout body shape related to divergent feeding methods on large to small snails within the genus Carabus. We show that the zinc-finger transcription factor encoded by odd-paired (opa) controls morphological variation in the snail-feeding ground beetle Carabus blaptoides. Specifically, opa was identified as the gene underlying the slender to stout morphological difference between subspecies through genetic mapping and functional analysis via gene knockdown. Further analyses revealed that changes in opa cis-regulatory sequences likely contributed to the differences in body shape and size between C. blaptoides subspecies. Among opa cis-regulatory sequences, single nucleotide polymorphisms on the transcription factor binding sites may be associated with the morphological differences between C. blaptoides subspecies. opa was highly conserved in a wide range of taxa, especially in beetles. Therefore, opa may play an important role in adaptive morphological divergence in beetles.

Characterization of the far-red light absorbing light-harvesting chlorophyll a/b binding complex, a derivative of the distinctive Lhca gene family in green algae.

Prasiola crispa, an aerial green alga, exhibits remarkable adaptability to the extreme conditions of Antarctica by forming layered colonies capable of utilizing far-red light for photosynthesis. Despite a recent report on the structure of P. crispa's unique light-harvesting chlorophyll (Chl)-binding protein complex (Pc-frLHC), which facilitates far-red light absorption and uphill excitation energy transfer to photosystem II, the specific genes encoding the subunits of Pc-frLHC have not yet been identified. Here, we report a draft genome sequence of P. crispa strain 4113, originally isolated from soil samples on Ongul Island, Antarctica. We obtained a 92 Mbp sequence distributed in 1,045 scaffolds comprising 10,244 genes, reflecting 87.1% of the core eukaryotic gene set. Notably, 26 genes associated with the light-harvesting Chl a/b binding complex (LHC) were identified, including four Pc-frLHC genes, with similarity to a noncanonical Lhca gene with four transmembrane helices, such as Ot_Lhca6 in Ostreococcus tauri and Cr_LHCA2 in Chlamydomonas reinhardtii. A comparative analysis revealed that Pc-frLHC shares homology with certain Lhca genes found in Coccomyxa and Trebouxia species. This similarity indicates that Pc-frLHC has evolved from an ancestral Lhca gene with four transmembrane helices and branched out within the Trebouxiaceae family. Furthermore, RNA-seq analysis conducted during the initiation of Pc-frLHC gene induction under red light illumination indicated that Pc-frLHC genes were induced independently from other genes associated with photosystems or LHCs. Instead, the genes of transcription factors, helicases, chaperones, heat shock proteins, and components of blue light receptors were identified to coexpress with Pc-frLHC. Those kinds of information could provide insights into the expression mechanisms of Pc-frLHC and its evolutional development.

W chromosome sequences of two bombycid moths provide an insight into the origin of Fem.

Fem is a W-linked gene that encodes a piRNA precursor, and its product, Fem piRNA, is a master factor of female determination in Bombyx mori. Fem has low similarity to any known sequences, and the origin of Fem remains unclear. So far, two hypotheses have been proposed for the origin of Fem: The first hypothesis is that Fem is an allele of Masc, which assumes that the W chromosome was originally a homologous chromosome of the Z chromosome. The second hypothesis is that Fem arose by the transposition of Masc to the W chromosome. To explore the origin of Fem, we determined the W chromosome sequences of B. mori and, as a comparison, a closely relative bombycid species of Trilocha varians with a Fem-independent sex determination system. To our surprise, although the sequences of W and Z chromosomes show no homology to each other, a few pairs of homologues are shared by W and Z chromosomes, indicating the W chromosome of both species originated from Z chromosome. In addition, the W chromosome of T. varians lacks Fem, while the W chromosome of B. mori has over 100 copies of Fem. The high-quality assembly of the W chromosome of B. mori arose the third hypothesis about the origin of Fem: Fem is a chimeric sequence of multiple transposons. More than half of one transcriptional unit of Fem shows a significant homology to RTE-BovB. Moreover, the Fem piRNA-producing region could correspond to the boundary of the two transposons, gypsy and satellite DNA.

Chimera RNA transcribed from integrated HPV18 genome with adjacent host genomic region promotes oncogenic gene expression through condensate formation.

Most cervical cancers are caused by human papillomavirus (HPV) infection. In HeLa cells, the HPV18 viral genome is integrated at chromosome 8q24.21 and activates transcription of the proto-oncogene c-Myc. However, the mechanism of how the integrated HPV genome and its transcribed RNAs exhibit transcription activation function has not been fully elucidated. In this study, we found that HPV18 transcripts contain an enhancer RNA-like function to activate proximal genes including CCAT1-5L and c-Myc. We showed that the human genome-integrated HPV18 genes are activated by transcription coregulators including BRD4 and Mediator. The transcribed HPV18 RNAs form a liquid-like condensate at chromosome 8q24.21 locus, which in turn accumulates RNA polymerase II. Moreover, we focused on a relatively uncharacterized transcript from the upstream region of CCAT1, named URC. The URC RNA is transcribed as a chimera RNA with HPV18 and is composed of the 3'-untranslated region of the HPV18 transcript. We experimentally showed that the URC contributes to stabilization of HPV18 RNAs by supplying a polyadenylation site for the HPV18 transcript. Our findings suggest that integrated HPV18 at 8q24.21 locus produces HPV18-URC chimera RNA and promotes tumorigenesis through RNA-based condensate formation.

Behaviors of nucleosomes with mutant histone H4s in euchromatic domains of living human cells.

Since Robert Feulgen first stained DNA in the cell, visualizing genome chromatin has been a central issue in cell biology to uncover how chromatin is organized and behaves in the cell. To approach this issue, we have developed single-molecule imaging of nucleosomes, a basic unit of chromatin, to unveil local nucleosome behavior in living cells. In this study, we investigated behaviors of nucleosomes with various histone H4 mutants in living HeLa cells to address the role of H4 tail acetylation, including H4K16Ac and others, which are generally associated with more transcriptionally active chromatin regions. We ectopically expressed wild-type (wt) or mutated H4s (H4K16 point; H4K5,8,12,16 quadruple; and H4 tail deletion) fused with HaloTag in HeLa cells. Cells that expressed wtH4-Halo, H4K16-Halo mutants, and multiple H4-Halo mutants had euchromatin-concentrated distribution. Consistently, the genomic regions of the wtH4-Halo nucleosomes corresponded to Hi-C contact domains (or topologically associating domains, TADs) with active chromatin marks (A-compartment). Utilizing single-nucleosome imaging, we found that none of the H4 deacetylation or acetylation mimicked H4 mutants altered the overall local nucleosome motion. This finding suggests that H4 mutant nucleosomes embedded in the condensed euchromatic domains with excess endogenous H4 nucleosomes cannot cause an observable change in the local motion. Interestingly, H4 with four lysine-to-arginine mutations displayed a substantial freely diffusing fraction in the nucleoplasm, whereas H4 with a truncated N-terminal tail was incorporated in heterochromatic regions as well as euchromatin. Our study indicates the power of single-nucleosome imaging to understand individual histone/nucleosome behavior reflecting chromatin environments in living cells.

Loss of One X and the Y Chromosome Changes the Configuration of the X Inactivation Center in the Genus Tokudaia.

INTRODUCTION: X chromosome inactivation (XCI) is an essential mechanism for dosage compensation between females and males in mammals. In females, XCI is controlled by a complex, conserved locus termed the X inactivation center (Xic), in which the lncRNA Xist is the key regulator. However, little is known about the Xic in species with unusual sex chromosomes. The genus Tokudaia includes three rodent species endemic to Japan. Tokudaia osimensis and Tokudaia tokunoshimensis lost the Y chromosome (XO/XO), while Tokudaia muenninki (TMU) acquired a neo-X region by fusion of the X chromosome and an autosome (XX/XY). We compared the gene location and structure in the Xic among Tokudaia species. METHODS: Gene structure of nine genes in Xic was predicted, and the gene location and genome sequences of Xic were compared between mouse and Tokudaia species. The expression level of the gene was confirmed by transcripts per million calculation using RNA-seq data. RESULTS: Compared to mouse, the Xic gene order and location were conserved in Tokudaia species. However, remarkable structure changes were observed in lncRNA genes, Xist and Tsix, in the XO/XO species. In Xist, important functional repeats, B-, C-, D-, and E-repeats, were partially or completely lost due to deletions in these species. RNA-seq data showed that female-specific expression patterns of Xist and Tsix were confirmed in TMU, however, not in the XO/XO species. Additionally, three deletions and one inversion were confirmed in the intergenic region between Jpx and Ftx in the XO/XO species. CONCLUSION: Our findings indicate that even if the Xist and Tsix lncRNAs are expressed, they are incapable of producing a successful and lasting XCI in the XO/XO species. We hypothesized that the significant structure change in the intergenic region of Jpx-Ftx resulted in the inability to perform the XCI, and, as a result, a lack of Xist expression. Our results collectively suggest that structural changes in the Xic occurred in the ancestral lineage of XO/XO species, likely due to the loss of one X chromosome and the Y chromosome as a consequence of the degradation of the XCI system.

Identification of gene isoforms and their switching events between male and female embryos of the parthenogenetic crustacean Daphnia magna.

The cladoceran crustacean Daphnia exhibits phenotypic plasticity, a phenomenon that leads to diverse phenotypes from one genome. Alternative usage of gene isoforms has been considered a key gene regulation mechanism for controlling different phenotypes. However, to understand the phenotypic plasticity of Daphnia, gene isoforms have not been comprehensively analyzed. Here we identified 25,654 transcripts derived from the 9710 genes expressed during environmental sex determination of Daphnia magna using the long-read RNA-Seq with PacBio Iso-Seq. We found that 14,924 transcripts were previously unidentified and 5713 genes produced two or more isoforms. By a combination of Illumina short-read RNA-Seq, we detected 824 genes that implemented switching of the highest expressed isoform between females and males. Among the 824 genes, we found isoform switching of an ortholog of CREB-regulated transcription coactivator, a major regulator of carbohydrate metabolism in animals, and a correlation of this switching event with the sexually dimorphic expression of carbohydrate metabolic genes. These results suggest that a comprehensive catalog of isoforms may lead to understanding the molecular basis for environmental sex determination of Daphnia. We also infer the applicability of the full-length isoform analyses to the elucidation of phenotypic plasticity in Daphnia.

Genome sequence and cell biological toolbox of the highly regenerative, coenocytic green feather alga Bryopsis.

Green feather algae (Bryopsidales) undergo a unique life cycle in which a single cell repeatedly executes nuclear division without cytokinesis, resulting in the development of a thallus (>100 mm) with characteristic morphology called coenocyte. Bryopsis is a representative coenocytic alga that has exceptionally high regeneration ability: extruded cytoplasm aggregates rapidly in seawater, leading to the formation of protoplasts. However, the genetic basis of the unique cell biology of Bryopsis remains poorly understood. Here, we present a high-quality assembly and annotation of the nuclear genome of Bryopsis sp. (90.7 Mbp, 27 contigs, N50 = 6.7 Mbp, 14 034 protein-coding genes). Comparative genomic analyses indicate that the genes encoding BPL-1/Bryohealin, the aggregation-promoting lectin, are heavily duplicated in Bryopsis, whereas homologous genes are absent in other ulvophyceans, suggesting the basis of regeneration capability of Bryopsis. Bryopsis sp. possesses >30 kinesins but only a single myosin, which differs from other green algae that have multiple types of myosin genes. Consistent with this biased motor toolkit, we observed that the bidirectional motility of chloroplasts in the cytoplasm was dependent on microtubules but not actin in Bryopsis sp. Most genes required for cytokinesis in plants are present in Bryopsis, including those in the SNARE or kinesin superfamily. Nevertheless, a kinesin crucial for cytokinesis initiation in plants (NACK/Kinesin-7II) is hardly expressed in the coenocytic part of the thallus, possibly underlying the lack of cytokinesis in this portion. The present genome sequence lays the foundation for experimental biology in coenocytic macroalgae.

Physical and functional interaction among Irf8 enhancers during dendritic cell differentiation.

The production of type 1 conventional dendritic cells (cDC1s) requires high expression of the transcription factor IRF8. Three enhancers at the Irf8 3' region function in a differentiation stage-specific manner. However, whether and how these enhancers interact physically and functionally remains unclear. Here, we show that the Irf8 3' enhancers directly interact with each other and contact the Irf8 gene body during cDC1 differentiation. The +56 kb enhancer, which functions from multipotent progenitor stages, activates the other 3' enhancers through an IRF8-dependent transcription factor program, that is, in trans. Then, the +32 kb enhancer, which operates in cDC1-committed cells, reversely acts in cis on the other 3' enhancers to maintain the high expression of Irf8. Indeed, mice with compound heterozygous deletion of the +56 and +32 kb enhancers are unable to generate cDC1s. These results illustrate how multiple enhancers cooperate to induce a lineage-determining transcription factor gene during cell differentiation.

Decoding host-microbiome interactions through co-expression network analysis within the non-human primate intestine.

The gut microbiome affects the health status of the host through complex interactions with the host's intestinal wall. These host-microbiome interactions may spatially vary along the physical and chemical environment of the intestine, but these changes remain unknown. This study investigated these intricate relationships through a gene co-expression network analysis based on dual transcriptome profiling of different intestinal sites-cecum, transverse colon, and rectum-of the primate common marmoset. We proposed a gene module extraction algorithm based on the graph theory to find tightly interacting gene modules of the host and the microbiome from a vast co-expression network. The 27 gene modules identified by this method, which include both host and microbiome genes, not only produced results consistent with previous studies regarding the host-microbiome relationships, but also provided new insights into microbiome genes acting as potential mediators in host-microbiome interplays. Specifically, we discovered associations between the host gene FBP1, a cancer marker, and polysaccharide degradation-related genes (pfkA and fucI) coded by Bacteroides vulgatus, as well as relationships between host B cell-specific genes (CD19, CD22, CD79B, and PTPN6) and a tryptophan synthesis gene (trpB) coded by Parabacteroides distasonis. Furthermore, our proposed module extraction algorithm surpassed existing approaches by successfully defining more functionally related gene modules, providing insights for understanding the complex relationship between the host and the microbiome.IMPORTANCEWe unveiled the intricate dynamics of the host-microbiome interactions along the colon by identifying closely interacting gene modules from a vast gene co-expression network, constructed based on simultaneous profiling of both host and microbiome transcriptomes. Our proposed gene module extraction algorithm, designed to interpret inter-species interactions, enabled the identification of functionally related gene modules encompassing both host and microbiome genes, which was challenging with conventional modularity maximization algorithms. Through these identified gene modules, we discerned previously unrecognized bacterial genes that potentially mediate in known relationships between host genes and specific bacterial species. Our findings underscore the spatial variations in host-microbiome interactions along the colon, rather than displaying a uniform pattern throughout the colon.

Synchronized Expansion and Contraction of Olfactory, Vomeronasal, and Taste Receptor Gene Families in Hystricomorph Rodents.

Chemical senses, including olfaction, pheromones, and taste, are crucial for the survival of most animals. There has long been a debate about whether different types of senses might influence each other. For instance, primates with a strong sense of vision are thought to have weakened olfactory abilities, although the oversimplified trade-off theory is now being questioned. It is uncertain whether such interactions between different chemical senses occur during evolution. To address this question, we examined four receptor gene families related to olfaction, pheromones, and taste: olfactory receptor (OR), vomeronasal receptor type 1 and type 2 (V1R and V2R), and bitter taste receptor (T2R) genes in Hystricomorpha, which is morphologically and ecologically the most diverse group of rodents. We also sequenced and assembled the genome of the grasscutter, Thryonomys swinderianus. By examining 16 available genome assemblies alongside the grasscutter genome, we identified orthologous gene groups among hystricomorph rodents for these gene families to separate the gene gain and loss events in each phylogenetic branch of the Hystricomorpha evolutionary tree. Our analysis revealed that the expansion or contraction of the four gene families occurred synchronously, indicating that when one chemical sense develops or deteriorates, the others follow suit. The results also showed that V1R/V2R genes underwent the fastest evolution, followed by OR genes, and T2R genes were the most evolutionarily stable. This variation likely reflects the difference in ligands of V1R/V2Rs, ORs, and T2Rs: species-specific pheromones, environment-based scents, and toxic substances common to many animals, respectively.