Search

Abstract The molecular etiology of idiopathic pulmonary fibrosis (IPF) has been extensively investigated to identify new therapeutic targets. Although anti-inflammatory treatments are not effective for patients with IPF, damaged alveolar epithelial cells play a critical role in lung fibrogenesis. Here, we establish an organoid-based lung fibrosis model using mouse and human lung tissues to assess the direct communication between damaged alveolar type II (AT2)-lineage cells and lung fibroblasts by excluding immune cells. Using this in vitro model and mouse genetics, we demonstrate that bleomycin causes DNA damage and activates p53 signaling in AT2-lineage cells, leading to AT2-to-AT1 transition-like state with a senescence-associated secretory phenotype (SASP). Among SASP-related factors, TGF-β plays an exclusive role in promoting lung fibroblast-to-myofibroblast differentiation. Moreover, the autocrine TGF-β-positive feedback loop in AT2-lineage cells is a critical cellular system in non-inflammatory lung fibrogenesis. These findings provide insights into the mechanism of IPF and potential therapeutic targets.

The diversity of neural stem cells is a hallmark of the cerebral cortex development in gyrencephalic mammals, such as Primates and Carnivora. Among them, ferrets are a good model for mechanistic studies. However, information on their neural progenitor cells (NPC), termed radial glia (RG), is limited. Here, we surveyed the temporal series of single-cell transcriptomes of progenitors regarding ferret corticogenesis and found a conserved diversity and temporal trajectory between human and ferret NPC, despite the large timescale difference. We found truncated RG (tRG) in ferret cortical development, a progenitor subtype previously described in humans. The combination of in silico and in vivo analyses identified that tRG differentiate into both ependymal and astrogenic cells. Via transcriptomic comparison, we predict that this is also the case in humans. Our findings suggest that tRG plays a role in the formation of adult ventricles, thereby providing the architectural bases for brain expansion.

The hagfishes (Myxiniformes) arose from agnathan (jawless vertebrate) lineages and they are one of only two extant cyclostome taxa, together with lampreys (Petromyzontiformes). Even though whole genome sequencing has been achieved for diverse vertebrate taxa, genome-wide sequence information has been highly limited for cyclostomes. Here we sequenced the genome of the inshore hagfish Eptatretus burgeri using DNA extracted from the testis, with a short-read sequencing platform, aiming to reconstruct a high-coverage protein-coding gene catalogue. The obtained genome assembly, scaffolded with mate-pair reads and paired RNA-seq reads, exhibited an N50 scaffold length of 293 Kbp, which allowed the genome-wide prediction of coding genes. This computation resulted in the gene models whose completeness was estimated at the complete coverage of more than 83 % and the partial coverage of more than 93 % by referring to evolutionarily conserved single-copy orthologs. The high contiguity of the assembly and completeness of the gene models promise a high utility in various comparative analyses including phylogenomics and phylome exploration.

With the recent development of ultraviolet light emitting diodes (UV-LEDs), it is now possible to apply UV disinfection, previously used only for medium-to-large sized sewage treatment plants, to Johkasou wastewater, a household wastewater treatment system. However, it is known that coliforms that have been inactivated by UV irradiation regain their activity after of sunlight irradiation, etc. Therefore, in small Johkasou, where discharge channels are often open, photorecovery may occur immediately after UV disinfection. However, no studies have been conducted to investigate photorecovery by sunlight after UV disinfection of wastewater from small Johkasou. In this study, the relationship between the photorecovery of coliforms and the dose of sunlight after UV-LED irradiation of effluent in Johkasou was investigated. As a result, we identified the amount of sunlight that will maximize photoreactivation effect. We also clarified the suppressing effect of an excessive amount of sunlight on the photoreactivation phenomenon. Additionally, since it was confirmed that the photoreactivation effect was suppressed by a certain period of darkness following UV-LED irradiation, it was considered that the photoreactivation effect can also be suppressed by allowing wastewater to reside in darkness for the currently estimated residence time in a disinfection tank.

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.

Measurement(s) sequence_assembly • genome • DNA Technology Type(s) sequence assembly process • DNA sequencing Sample Characteristic - Organism Macaca fascicularis • Callithrix jacchus Machine-accessible metadata file describing the reported data: https://doi.org/10.6084/m9.figshare.14665314

RNA virus populations are not clonal; rather, they comprise a mutant swarm in which sequences are slightly different from the master sequence. Genetic diversity within a population (intrapopulation genetic diversity) is critical for RNA viruses to survive under environmental stresses. Disinfection has become an important practice in the control of pathogenic viruses; however, the impact of intrapopulation genetic diversity on the sensitivity of disinfection, defined as –log10 (postdisinfected infectious titer/predisinfected titer), has not been elucidated. In this study, we serially passaged populations of rhesus rotavirus. We demonstrated that populations with reduced chlorine sensitivity emerged at random and independently of chlorine exposure. Sequencing analysis revealed that compared with sensitive populations, less-sensitive ones had higher non-synonymous genetic diversity of the outer capsid protein gene, suggesting that changes in the amino acid sequences of the outer capsid protein were the main factors influencing chlorine sensitivity. No common mutations were found among less-sensitive populations, indicating that rather than specific mutations, the diversity of the outer capsid protein itself was associated with the disinfection sensitivity and that the disinfection sensitivity changed stochastically. Simulation results suggest that the disinfection sensitivity of a genetically diverse population is destabilized if cooperative viral clusters including multiple sequences are formed. These results advocate that any prevention measures leading to low intrapopulation genetic diversity are important to prevent the spread and evolution of pathogenic RNA viruses in society.

Abstract Background Vertebrates are characterized by possession of hypobranchial muscles (HBMs). Cyclostomes, or modern jawless vertebrates, possess a rudimentary and superficial HBM lateral to the pharynx, whereas the HBM in jawed vertebrates is internalized and anteroposteriorly specified. Precursor cells of the HBM, marked by expression of Lbx1, originate from somites and undergo extensive migration before becoming innervated by the hypoglossal nerve. How the complex form of HBM arose in evolution is relevant to the establishment of the vertebrate body plan, but despite having long been assumed to be similar to that of limb muscles, modification of developmental mechanisms of HBM remains enigmatic. Results Here we characterize the expression of Lbx genes in lamprey and hagfish (cyclostomes) and catshark (gnathostome; jawed vertebrates). We show that the expression patterns of the single cyclostome Lbx homologue, Lbx-A, do not resemble the somitic expression of mammalian Lbx1. Disruption of Lbx-A revealed that LjLbx-A is required for the formation of both HBM and body wall muscles, likely due to the insufficient extension of precursor cells rather than to hindered muscle differentiation. Both homologues of Lbx in the catshark were expressed in the somitic muscle primordia, unlike in amniotes. During catshark embryogenesis, Lbx2 is expressed in the caudal HBM as well as in the abdominal rectus muscle, similar to lamprey Lbx-A, whereas Lbx1 marks the rostral HBM and pectoral fin muscle. Conclusions We conclude that the vertebrate HBM primarily emerged as a specialized somatic muscle to cover the pharynx, and the anterior internalized HBM of the gnathostomes is likely a novelty added rostral to the cyclostome-like HBM, for which duplication and functionalization of Lbx genes would have been a prerequisite.

Lake Tazawa in Akita Prefecture, Japan, is an acidic lake caused by the inflow of low-pH water from Tamagawa Hot Springs. Because of this anthropogenic acidification, many indigenous fish species and aquatic ecosystem have perished. Although several counter-measures such as artificial neutralization have been implemented against the acidic waters of Tamagawa River flowing into Lake Tazawa, the pH level of the lake remained low. Therefore, a bioassay evaluation of the low-pH lake water is necessary for ecosystem restoration. In this study, short-term chronic toxicity tests were applied to observe biological responses of three trophic levels of aquatic organisms using water samples of Lake Tazawa. As a result, Lake Tazawa water has toxic effect on all tested aquatic organisms. Especially, the fish embryo/larvae was directly affected by low-pH water, while algae and crustaceans were affected not only by low pH but other factors such as fluorine ion and metals. In order to improve the water quality management of Lake Tazawa for ecosystem restoration, further efforts are needed to reduce multiple toxic substances in addition to the pH neutralization of lake water.

Wastewater reclamation and reuse have been practically applied to water-stressed regions, but waterborne pathogens remaining in insufficiently treated wastewater are of concern. Sanitation Safety Planning adopts the hazard analysis and critical control point (HACCP) approach to manage human health risks upon exposure to reclaimed wastewater. HACCP requires a predetermined reference value (critical limit: CL) at critical control points (CCPs), in which specific parameters are monitored and recorded in real time. A disinfection reactor of a wastewater treatment plant (WWTP) is regarded as a CCP, and one of the CCP parameters is the disinfection intensity (e.g., initial disinfectant concentration and contact time), which is proportional to the log reduction value (LRV) of waterborne pathogens. However, the achievable LRVs are not always stable because the disinfection intensity is affected by water quality parameters, which vary among WWTPs. In this study, we established models for projecting virus LRVs using ozone, in which water quality and operational parameters were used as explanatory variables. For the model construction, we used five machine learning algorithms and found that automatic relevance determination with interaction terms resulted in better prediction performances for norovirus and rotavirus LRVs. Poliovirus and coxsackievirus LRVs were predicted well by a Bayesian ridge with interaction terms and lasso with quadratic terms, respectively. The established models were relatively robust to predict LRV using new datasets that were out of the range of the training data used here, but it is important to collect LRV datasets further to make the models more predictable and flexible for newly obtained datasets. The modeling framework proposed here can help WWTP operators and risk assessors determine the appropriate CL to protect human health in wastewater reclamation and reuse.

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Abstract The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in animal genomes. Currently, our knowledge of its binding property is confined mainly to mammals. In this study, we identified CTCF homologs in extant jawless fishes and performed ChIP-seq for the CTCF protein in the Arctic lamprey. Our phylogenetic analysis suggests that the lamprey lineage experienced gene duplication that gave rise to its unique paralog, designated CTCF2, which is independent from the previously recognized duplication between CTCF and CTCFL. The ChIP-seq analysis detected comparable numbers of CTCF binding sites between lamprey, chicken, and human, and revealed that the lamprey CTCF protein binds to the two-part motif, consisting of core and upstream motifs previously reported for mammals. These findings suggest that this mode of CTCF binding was established in the last common ancestor of extant vertebrates (more than 500 million years ago). We analyzed CTCF binding inside Hox clusters, which revealed a reinforcement of CTCF binding in the region spanning Hox1-4 genes that is unique to lamprey. Our study provides not only biological insights into the antiquity of CTCF-based epigenomic regulation known in mammals but also a technical basis for comparative epigenomic studies encompassing the whole taxon Vertebrata.

With advances in sequencing technology, a vast amount of genomic sequence information has become available. However, annotating biological functions particularly of non-protein-coding regions in genome sequences without experiments is still a challenging task. Recently deep learning-based methods were shown to have the ability to predict gene regulatory regions from genome sequences, promising to aid the interpretation of genomic sequence data. Here, we report an improvement of the prediction accuracy for gene regulatory regions by using the design of convolution layers that efficiently process genomic sequence information, and developed a software, DeepGMAP, to train and compare different deep learning-based models (https://github.com/koonimaru/DeepGMAP). First, we demonstrate that our convolution layers, termed forward- and reverse-sequence scan (FRSS) layers, integrate both forward and reverse strand information, and enhance the power to predict gene regulatory regions. Second, we assessed previous studies and identified problems associated with data structures that caused overfitting. Finally, we introduce visualization methods to examine what the program learned. Together, our FRSS layers improve the prediction accuracy for gene regulatory regions.

Bivalves serve as an important aquaculture product, as they are the source of essential fatty acids, such as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), in our diet. However, their cultivation in the wild can be affected by fouling organisms that, in turn, affect their EPA and DHA content. The effects of fouling organisms on the EPA and DHA contents of cultivated bivalves have not been well documented. We examined the effects of fouling organisms on the EPA and DHA contents and condition index of cultured oysters, Crassostrea gigas, in an aquaculture system. We sampled two-year-old oysters from five sites in Shizugawa Bay, Japan, in August 2014. Most of the fouling organisms were sponges, macroalgae, and Mytilus galloprovincialis. A significant negative relationship existed between the DHA content in C. gigas and the presence of sponges and macroalgae. A lower C. gigas EPA content corresponded to a higher M. galloprovincialis fouling mass and a lower C. gigas condition index. This can be explained by dietary competition between C. gigas and M. galloprovincialis for diatoms, which were the main producer of EPA in our study sites. Our findings indicate that fouling organisms likely reduce the EPA and DHA content in cultivated oysters. Therefore, our results suggest that the current efforts to remove fouling organisms from oyster clusters is an effective strategy to enhance the content of EPA and DHA in oysters.

We have examined the role of Fam60a, a gene highly expressed in embryonic stem cells, in mouse development. Fam60a interacts with components of the Sin3a-Hdac transcriptional corepressor complex, and most Fam60a–/– embryos manifest hypoplasia of visceral organs and die in utero. Fam60a is recruited to the promoter regions of a subset of genes, with the expression of these genes being either up- or down-regulated in Fam60a–/– embryos. The DNA methylation level of the Fam60a target gene Adhfe1 is maintained at embryonic day (E) 7.5 but markedly reduced at E9.5 in Fam60a–/– embryos, suggesting that DNA demethylation is enhanced in the mutant. Examination of genome-wide DNA methylation identified several differentially methylated regions, which were preferentially hypomethylated, in Fam60a–/– embryos. Our data suggest that Fam60a is required for proper embryogenesis, at least in part as a result of its regulation of DNA methylation at specific gene promoters.

The identification of molecular networks at the system level in mammals is accelerated by next-generation mammalian genetics without crossing, which requires both the efficient production of whole-body biallelic knockout (KO) mice in a single generation and high-performance phenotype analyses. Here, we show that the triple targeting of a single gene using the CRISPR/Cas9 system achieves almost perfect KO efficiency (96%–100%). In addition, we developed a respiration-based fully automated non-invasive sleep phenotyping system, the Snappy Sleep Stager (SSS), for high-performance (95.3% accuracy) sleep/wake staging. Using the triple-target CRISPR and SSS in tandem, we reliably obtained sleep/wake phenotypes, even in double-KO mice. By using this system to comprehensively analyze all of the N-methyl-D-aspartate (NMDA) receptor family members, we found Nr3a as a short-sleeper gene, which is verified by an independent set of triple-target CRISPR. These results demonstrate the application of mammalian reverse genetics without crossing to organism-level systems biology in sleep research.

In de novo genome sequencing, mate-pair reads are crucial for scaffolding assembled contigs. However, preparation of mate-pair libraries is not a trivial task, even when using one of the latest approaches, the Nextera Mate Pair Sample Prep Kit from Illumina. To reduce cost and enhance library yield and fidelity when using this kit, we have modified the manufacturer's protocol based on (i) variable tagmentation conditions, (ii) intensive DNA shearing to decrease library insert length, and (iii) sequencing on an Illumina HiSeq with >150 cycles. Finally, we provide additional suggestions for further improvement in the application of this kit.

To better understand the process of karyotype evolution in Galloanserae (Galliformes and Anseriformes), we performed comparative chromosome painting with chicken chromosome-specific DNA probes and FISH mapping of the 18S-28S ribosomal RNA (rRNA) genes, telomeric (TTAGGG)n repeats, and cDNA clones of 37 genes for three anserid species, the domestic duck (Anas platyrhynchos), Muscovy duck (Cairina moschata), and Chinese goose (Anser cygnoides). Each chicken probe of chromosomes 1-9 and the Z chromosome painted a single pair of chromosomes in the three species except for the chromosome 4 probe, which painted acrocentric chromosome 4 and a pair of microchromosomes. The 18S-28S rRNA genes were localized to four pairs of microchromosomes in the domestic duck and Muscovy duck, and eight pairs of microchromosomes in the Chinese goose. The (TTAGGG)n repeats were localized to both telomeric ends of all chromosomes in the three species, and were additionally located in the interstitial region of the short arm of chromosome 1 in the domestic duck and in the centromeric region of chromosome 1 in the Muscovy duck. Comparative gene mapping of 37 chicken chromosome 1-4-linked and Z-linked genes revealed high chromosome homologies among three anserid species and also between the chicken and the three anserid species, although there were several chromosome rearrangements: pericentric inversion in chromosome 2, pericentric and paracentric inversions or centromere repositioning in the Z chromosome between the chicken and the anserid species, and pericentric inversions in chromosome 4 and the Z chromosome between the Chinese goose and the two duck species. These results collectively suggest that karyotypes have been highly conserved in Anseriformes and that chromosome rearrangements also occurred less frequently between Galliformes and Anseriformes after they diverged around 100 million years ago.