Impact of translational regulation on diel expression revealed by time-series ribosome profiling in Arabidopsis.

Plants have developed the ability to adjust to the day/night cycle through the expression of diel genes, which allow them to effectively respond to environmental changes and optimise their growth and development. Diel oscillations also have substantial implications in many physiological processes, including photosynthesis, floral development, and environmental stress responses. The expression of diel genes is regulated by a combination of the circadian clock and responses to environmental cues, such as light and temperature. A great deal of information is available on the transcriptional regulation of diel gene expression. However, the extent to which translational regulation is involved in controlling diel changes in expression is not yet clear. To investigate the impact of translational regulation on diel expression, we conducted Ribo-seq and RNA-seq analyses on a time-series sample of Arabidopsis shoots cultivated under a 12 h light/dark cycle. Our results showed that translational regulation is involved in about 71% of the genes exhibiting diel changes in mRNA abundance or translational activity, including clock genes, many of which are subject to both translational and transcriptional control. They also revealed that the diel expression of glycosylation and ion-transporter-related genes is mainly established through translational regulation. The expression of several diel genes likely subject to translational regulation through upstream open-reading frames was also determined.

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.

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.

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.

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.

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.

Transcriptome Analysis of Rice Root Tips Reveals Auxin, Gibberellin and Ethylene Signaling Underlying Nutritropism.

Nutritropism is a positive tropism toward nutrients in plant roots. An NH4+ gradient is a nutritropic stimulus in rice (Oryza sativa L.). When rice roots are exposed to an NH4+ gradient generated around nutrient sources, root tips bend toward and coil around the sources. The molecular mechanisms are largely unknown. Here, we analyzed the transcriptomes of the inside and outside of bending root tips exhibiting nutritropism to reveal nutritropic signal transduction. Tissues facing the nutrient sources (inside) and away (outside) were separately collected by laser microdissection. Principal component analysis revealed distinct transcriptome patterns between the two tissues. Annotations of 153 differentially expressed genes implied that auxin, gibberellin and ethylene signaling were activated differentially between the sides of the root tips under nutritropism. Exogenous application of transport and/or biosynthesis inhibitors of these phytohormones largely inhibited the nutritropism. Thus, signaling and de novo biosynthesis of the three phytohormones are necessary for nutritropism. Expression patterns of IAA genes implied that auxins accumulated more in the inside tissues, meaning that ammonium stimulus is transduced to auxin signaling in nutritropism similar to gravity stimulus in gravitropism. SAUR and expansin genes, which are known to control cell wall modification and to promote cell elongation in shoot gravitropism, were highly expressed in the inside tissues rather than the outside tissues, and our transcriptome data are unexplainable for differential elongation in root nutritropism.

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.

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.

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.

SDF-1/CXCR4 signal is involved in the induction of Primordial Germ Cell migration in a model marine fish, Japanese anchovy (Engraulis japonicus).

Primordial germ cells (PGCs) are pivotal for gonadal development and reproductive success. Though artificial induction of sterility by targeting PGCs are gaining popularity due to its advantages in fish surrogacy and biodiversity management, it is often skill and time intensive. In this study, we have focused on understanding the role of PGCs and the chemotactic SDF-1/CXCR4 signaling on gonad development of Japanese anchovy (JA, Engraulis japonicus), an upcoming marine model organism with eco-commercial values, with an aim to develop a novel, easy, and versatile gonad sterilization method. Our data showed that PGC migration related genes, i.e., sdf-1a, sdf-1b, cxcr4a, cxcr4b and vasa, are phylogenetically closer relatives of respective herring (Clupea harengus) and zebrafish (Danio rerio) homolog. Subsequently, PGC marking and live tracing experiments confirmed that PGC migration in JA initiates from 16 hours post fertilization (hpf) followed by PGC settlement in the gonadal ridge at 44 hpf. We found that overexpression of zebrafish sdf-1a mRNA in the germ cell suppresses cxcr4a and increases cxcr4b transcription at 8 hpf, dose dependently disrupts PGC migration at 24-48 hpf, induces PGC death and upregulates sdf-1b at 5 days post hatching. 48 h of immersion treatment with CXCR4 antagonist (AMD3100, Abcam) also accelerated PGC mismigration and pushed the PGC away from gonadal ridge in a dose responsive manner, and further when grown to adulthood caused germ cell less gonad formation in some individuals. Cumulatively, our data, for the first time, suggests that JA PGC migration is largely regulated by SDF1/CXCR4 signaling, and modulation of this signaling has strong potential for sterile, germ cell less gonad preparation at a mass scale. However, further in-depth analysis is pertinent to apply this methodology in marine fish species to successfully catapult Japanese anchovy into a true marine fish model.

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.

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.

Modeling the SDF-1/CXCR4 protein using advanced artificial intelligence and antagonist screening for Japanese anchovy.

SDF-1/CXCR4 chemokine signaling are indispensable for cell migration, especially the Primordial Germ Cell (PGC) migration towards the gonadal ridge during early development. We earlier found that this signaling is largely conserved in the Japanese anchovy (Engraulis japonicus, EJ), and a mere treatment of CXCR4 antagonist, AMD3100, leads to germ cell depletion and thereafter gonad sterilization. However, the effect of AMD3100 was limited. So, in this research, we scouted for CXCR4 antagonist with higher potency by employing advanced artificial intelligence deep learning-based computer simulations. Three potential candidates, AMD3465, WZ811, and LY2510924, were selected and in vivo validation was conducted using Japanese anchovy embryos. We found that seven transmembrane motif of EJ CXCR4a and EJ CXCR4b were extremely similar with human homolog while the CXCR4 chemokine receptor N terminal (PF12109, essential for SDF-1 binding) was missing in EJ CXCR4b. 3D protein analysis and cavity search predicted the cavity in EJ CXCR4a to be five times larger (6,307 Å³) than that in EJ CXCR4b (1,241 Å³). Docking analysis demonstrated lower binding energy of AMD3100 and AMD3465 to EJ CXCR4a (Vina score -9.6) and EJ CXCR4b (Vina score -8.8), respectively. Furthermore, we observed significant PGC mismigration in microinjected AMD3465 treated groups at 10, 100 and 1 × 105 nM concentration in 48 h post fertilized embryos. The other three antagonists showed various degrees of PGC dispersion, but no significant effect compared to their solvent control at tested concentrations was observed. Cumulatively, our results suggests that AMD3645 might be a better candidate for abnormal PGC migration in Japanese anchovy and warrants further investigation.

Alpha-glycosyl isoquercitrin alleviates subchronic social defeat stress-induced depression symptoms by modulating the microbiota-gut-brain axis in mice.

AIMS: Increasing evidence suggests a link between gut microbial dysbiosis and the pathogenesis of depression. Alpha-glycosyl isoquercitrin (AGIQ), consisting of isoquercitrin and its glycosylated quercetin, has beneficial effects on the gut microbiome and brain function. Here, we detected the potential antidepressant impact of a four-week administration of AGIQ and its underlying mechanisms using a mouse model of depression. MAIN METHODS: Male C57BL/6 mice were orally administered AGIQ (0.05 % or 0.5 % in drinking water) for 28 days; subchronic social defeat stress was performed in the last 10 days. Behavior tests were conducted to assess anxiety and depressive-like behaviors. Additionally, evaluations encompassed 5-hydroxytryptamine (5-HT) levels, the gut microbiota composition, lipopolysaccharide (LPS) concentrations, short-chain fatty acids levels, and intestinal barrier integrity changes. KEY FINDINGS: AGIQ significantly alleviated depression-like behaviors and increased hippocampal 5-HT levels. Further, AGIQ mitigated stress-induced gut microbial abnormalities and reduced the levels of LPS in the serum, which affected the relative gene expression levels of 5-HT biosynthesis enzymes in vitro. Furthermore, AGIQ reversed the reduced butyrate levels in cecal contents and improved the impaired intestinal barrier by increasing the expression of colonic zonula occluden-1 (ZO-1) and occludin, thereby decreasing LPS leakage. SIGNIFICANCE: Our results suggest that AGIQ could improve stress-induced depression by regulating the gut microbiome, which inhibits LPS production and maintains the gut barrier. This is the first report on the potential effect of AGIQ on depression via the gut microbiota-brain axis, shedding new light on treatment options.

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.

QTL analysis of femaleness in monoecious spinach and fine mapping of a major QTL using an updated version of chromosome-scale pseudomolecules.

Although spinach is predominantly dioecious, monoecious plants with varying proportions of female and male flowers are also present. Recently, monoecious inbred lines with highly female and male conditions have been preferentially used as parents for F1-hybrids, rather than dioecious lines. Accordingly, identifying the loci for monoecism is an important issue for spinach breeding. We here used long-read sequencing and Hi-C technology to construct SOL_r2.0_pseudomolecule, a set of six pseudomolecules of spinach chromosomes (total length: 879.2 Mb; BUSCO complete 97.0%) that are longer and more genetically complete than our previous version of pseudomolecules (688.0 Mb; 81.5%). Three QTLs, qFem2.1, qFem3.1, and qFem6.1, responsible for monoecism were mapped to SOL_r2.0_pseudomolecule. qFem3.1 had the highest LOD score and corresponded to the M locus, which was previously identified as a determinant of monoecious expression, by genetic analysis of progeny from female and monoecious plants. The other QTLs were shown to modulate the ratio of female to male flowers in monoecious plants harboring a dominant allele of the M gene. Our findings will enable breeders to efficiently produce highly female- and male-monoecious parental lines for F1-hybrids by pyramiding the three QTLs. Through fine-mapping, we narrowed the candidate region for the M locus to a 19.5 kb interval containing three protein-coding genes and one long non-coding RNA gene. Among them, only RADIALIS-like-2a showed a higher expression in the reproductive organs, suggesting that it might play a role in reproductive organogenesis. However, there is no evidence that it is involved in the regulation of stamen and pistil initiation, which are directly related to the floral sex differentiation system in spinach. Given that auxin is involved in reproductive organ formation in many plant species, genes related to auxin transport/response, in addition to floral organ formation, were identified as candidates for regulators of floral sex-differentiation from qFem2.1 and qFem6.1.

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.

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.

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.