SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data

Acceleration of genome rearrangement in clitellate annelids

Comparisons of multiple metazoan genomes have revealed the existence of ancestral linkage groups (ALGs), genomic scaffolds sharing sets of orthologous genes that have been inherited from ancestral animals for hundreds of millions of years (Simakov et al. 2022; Schultz et al. 2023) These ALGs have persisted across major animal taxa including Cnidaria, Deuterostomia, Ecdysozoa and Spiralia. Notwithstanding this general trend of chromosome-scale conservation, ALGs have been obliterated by extensive genome rearrangements in certain groups, most notably including Clitellata (oligochaetes and leeches), a group of easily overlooked invertebrates that is of tremendous ecological, agricultural and economic importance (Charles 2019; Barrett 2016). To further investigate these rearrangements, we have undertaken a comparison of 12 clitellate genomes (including four newly sequenced species) and 11 outgroup representatives. We show that these rearrangements began at the base of the Clitellata (rather than progressing gradually throughout polychaete annelids), that the inter-chromosomal rearrangements continue in several clitellate lineages and that these events have substantially shaped the evolution of the otherwise highly conserved Hox cluster.

Growth phase-dependent reorganization of cryptophyte photosystem I antennae

Photosynthetic cryptophytes are eukaryotic algae that utilize membrane-embedded chlorophyll a/c binding proteins (CACs) and lumen-localized phycobiliproteins (PBPs) as their light-harvesting antennae. Cryptophytes go through logarithmic and stationary growth phases, and may adjust their light-harvesting capability according to their particular growth state. How cryptophytes change the type/arrangement of the photosynthetic antenna proteins to regulate their light-harvesting remains unknown. Here we solve four structures of cryptophyte photosystem I (PSI) bound with CACs that show the rearrangement of CACs at different growth phases. We identify a cryptophyte-unique protein, PsaQ, which harbors two chlorophyll molecules. PsaQ specifically binds to the lumenal region of PSI during logarithmic growth phase and may assist the association of PBPs with photosystems and energy transfer from PBPs to photosystems.

A chromosome-level genome assembly of an avivorous bat species (Nyctalus aviator)

Currently, three carnivorous bat species, namely Ia io, Nyctalus lasiopterus, and Nyctalus aviator, are known to actively prey on seasonal migratory birds (hereinafter referred to as "avivorous bats"). However, the absence of reference genomes impedes a thorough comprehension of the molecular adaptations of avivorous bat species. Herein, we present the high-quality chromosome-scale reference genome of N. aviator based on PacBio subreads, DNBSEQ short-reads and Hi-C sequencing data. The genome assembly size of N. aviator is 1.77 Gb, with a scaffold N50 of 102 Mb, of which 99.8% assembly was anchored into 21 pseudo-chromosomes. After masking 635.1 Mb repetitive sequences, a total of 19,412 protein-coding genes were identified, of which 99.3% were functionally annotated. The genome assembly and gene prediction reached 96.1% and 96.1% completeness of Benchmarking Universal Single-Copy Orthologs (BUSCO), respectively. This chromosome-level reference genome of N. aviator fills a gap in the existing information on the genomes of carnivorous bats, especially avivorous ones, and will be valuable for mechanism of adaptations to dietary niche expansion in bat species.

Identification and map-based cloning of an EMS-induced mutation in wheat gene TaSP1 related to spike architecture

KEY MESSAGE

A candidate gene TaSP1 related to spike shape was cloned, and the gene-specific marker was developed to efficiently track the superior haplotype in common wheat. Spike shape, an important factor that affects wheat grain yield, is mainly defined by spike length (SPL), spikelet number (SPN), and compactness. Zhoumai32 mutant 1160 (ZM1160), a mutant obtained from ethyl methane sulfonate (EMS) treatment of hexaploid wheat variety Zhoumai32, was used to identify and clone the candidate gene that conditioned the spike shape. Genetic analysis of an F2 population derived from a cross of ZM1160 and Bainong207 suggested that the compact spike shape in ZM1160 was controlled by a single recessive gene, and therefore, the mutated gene was designated as Tasp1. With polymorphic markers identified through bulked segregant analysis (BSA), the gene was mapped to a 2.65-cM interval flanked by markers YZU0852 and MIS46239 on chromosome 7D, corresponding to a 0.42-Mb physical interval of Chinese spring (CS) reference sequences (RefSeq v1.0). To fine map TaSP1, 15 and seven recombinants were, respectively, screened from 1599 and 1903 F3 plants derived from the heterozygous F2 plants. Finally, TaSP1 was delimited to a 21.9 Kb (4,870,562 to 4,892,493 bp) Xmis48123-Xmis48104 interval. Only one high-confidence gene TraesCS7D02G010200 was annotated in this region, which encodes an unknown protein with a putative vWA domain. Quantitative reverse transcription PCR (qRT-PCR) analysis showed that TraesCS7D02G010200 was mainly expressed in the spike. Haplotype analysis of 655 wheat cultivars using the candidate gene-specific marker Xg010200p2 identified a superior haplotype TaSP1b with longer spike and more spikelet number. TaSP1 is beneficial to the improvement in wheat spike shape.

Metagenomic analysis of individual mosquito viromes reveals the geographical patterns and drivers of viral diversity

Mosquito transmitted viruses are responsible for an increasing burden of human disease. Despite this, little is known about the diversity and ecology of viruses within individual mosquito hosts. Here, using a meta-transcriptomic approach, we determined the viromes of 2,438 individual mosquitoes (81 species), spanning ~4,000 km along latitudes and longitudes in China. From these data we identified 393 viral species associated with mosquitoes, including 7 (putative) species of arthropod-borne viruses (that is, arboviruses). We identified potential mosquito species and geographic hotspots of viral diversity and arbovirus occurrence, and demonstrated that the composition of individual mosquito viromes was strongly associated with host phylogeny. Our data revealed a large number of viruses shared among mosquito species or genera, enhancing our understanding of the host specificity of insect-associated viruses. We also detected multiple virus species that were widespread throughout the country, perhaps reflecting long-distance mosquito dispersal. Together, these results greatly expand the known mosquito virome, linked viral diversity at the scale of individual insects to that at a country-wide scale, and offered unique insights into the biogeography and diversity of viruses in insect vectors.

The rise of baobab trees in Madagascar

The baobab trees (genus Adansonia) have attracted tremendous attention because of their striking shape and distinctive relationships with fauna1. These spectacular trees have also influenced human culture, inspiring innumerable arts, folklore and traditions. Here we sequenced genomes of all eight extant baobab species and argue that Madagascar should be considered the centre of origin for the extant lineages, a key issue in their evolutionary history2,3. Integrated genomic and ecological analyses revealed the reticulate evolution of baobabs, which eventually led to the species diversity seen today. Past population dynamics of Malagasy baobabs may have been influenced by both interspecific competition and the geological history of the island, especially changes in local sea levels. We propose that further attention should be paid to the conservation status of Malagasy baobabs, especially of Adansonia suarezensis and Adansonia grandidieri, and that intensive monitoring of populations of Adansonia za is required, given its propensity for negatively impacting the critically endangered Adansonia perrieri.

Overexpression of the Chrysanthemum lavandulifolium ROS1 gene promotes flowering in Arabidopsis thaliana by reducing the methylation level of CONSTANS

DNA demethylation is involved in the regulation of flowering in plants, yet the underlying molecular mechanisms remain largely unexplored. The RELEASE OF SILENCING 1 (ROS1) gene, encoding a DNA demethyltransferase, plays key roles in many developmental processes. In this study, the ROS1 gene was isolated from Chrysanthemum lavandulifolium, where it was strongly expressed in the leaves, buds and flowers. Overexpression of the ClROS1 gene caused an early flowering phenotype in Arabidopsis thaliana. RNA-seq analysis of the transgenic plants revealed that differentially expressed genes (DEGs) were significantly enriched in the circadian rhythm pathway and that the positive regulator of flowering, CONSTANS (CO), was up-regulated. Additionally, whole-genome bisulphite sequencing (WGBS), PCR following methylation-dependent digestion with the enzyme McrBC, and bisulfite sequencing PCR (BSP) confirmed that the methylation level of the AtCO promoter was reduced, specifically in CG context. Overall, our results demonstrated that ClROS1 accelerates flowering by reducing the methylation level of the AtCO promoter. These findings clarify the epigenetic mechanism by which ClROS1-mediated DNA demethylation regulates flowering.

The high-quality genome of Grona styracifolia uncovers the genomic mechanism of high levels of schaftoside, a promising drug candidate for treatment of COVID-19

Recent study has evidenced that traditional Chinese medicinal (TCM) plant-derived schaftoside shows promise as a potential drug candidate for COVID-19 treatment. However, the biosynthetic pathway of schaftoside in TCM plants remains unknown. In this study, the genome of the TCM herb Grona styracifolia (Osbeck) H.Ohashi & K.Ohashi (GSO), which is rich in schaftoside, was sequenced, and a high-quality assembly of GSO genome was obtained. Our findings revealed that GSO did not undergo recent whole genome duplication (WGD) but shared an ancestral papilionoid polyploidy event, leading to the gene expansion of chalcone synthase (CHS) and isoflavone 2'-hydroxylase (HIDH). Furthermore, GSO-specific tandem gene duplication resulted in the gene expansion of C-glucosyltransferase (CGT). Integrative analysis of the metabolome and transcriptome identified 13 CGTs and eight HIDHs involved in the biosynthetic pathway of schaftoside. Functional studies indicated that CGTs and HIDHs identified here are bona fide responsible for the biosynthesis of schaftoside in GSO, as confirmed through hairy root transgenic system and in vitro enzyme activity assay. Taken together, the ancestral papilionoid polyploidy event expanding CHSs and HIDHs, along with the GSO-specific tandem duplication of CGT, contributes, partially if not completely, to the robust biosynthesis of schaftoside in GSO. These findings provide insights into the genomic mechanisms underlying the abundant biosynthesis of schaftoside in GSO, highlighting the potential of GSO as a source of bioactive compounds for pharmaceutical development.

PUMC-MB1 is a novel group 3 medulloblastoma preclinical model, sensitive to PI3K/mTOR dual inhibitor

PURPOSE

Medulloblastoma (MB), a common and heterogeneous posterior fossa tumor in pediatric patients, presents diverse prognostic outcomes. To advance our understanding of MB's intricate biology, the development of novel patient tumor-derived culture MB models with necessary data is still an essential requirement.

METHODS

We continuously passaged PUMC-MB1 in vitro in order to establish a continuous cell line. We examined the in vitro growth using Cell Counting Kit-8 (CCK-8) and in vivo growth with subcutaneous and intracranial xenograft models. The xenografts were investigated histopathologically with Hematoxylin and Eosin (HE) staining and immunohistochemistry (IHC). Concurrently, we explored its molecular features using Whole Genome Sequencing (WGS), targeted sequencing, and RNA sequecing. Guided by bioinformatics analysis, we validated PUMC-MB1's drug sensitivity in vitro and in vivo.

RESULTS

PUMC-MB1, derived from a high-risk MB patient, displayed a population doubling time (PDT) of 48.18 h and achieved 100% tumor growth in SCID mice within 20 days. HE and Immunohistochemical examination of the original tumor and xenografts confirmed the classification of PUMC-MB1 as a classic MB. Genomic analysis via WGS revealed concurrent MYC and OTX2 amplifications. The RNA-seq data classified it within the Group 3 MB subgroup, while according to the WHO classification, it fell under the Non-WNT/Non-SHH MB. Comparative analysis with D283 and D341med identified 4065 differentially expressed genes, with notable enrichment in the PI3K-AKT pathway. Cisplatin, 4-hydroperoxy cyclophosphamide/cyclophosphamide, vincristine, and dactolisib (a selective PI3K/mTOR dual inhibitor) significantly inhibited PUMC-MB1 proliferation in vitro and in vivo.

CONCLUSIONS

PUMC-MB1, a novel Group 3 (Non-WNT/Non-SHH) MB cell line, is comprehensively characterized for its growth, pathology, and molecular characteristics. Notably, dactolisib demonstrated potent anti-proliferative effects with minimal toxicity, promising a potential therapeutic avenue. PUMC-MB1 could serve as a valuable tool for unraveling MB mechanisms and innovative treatment strategies.

Paraflavitalea pollutisoli sp. nov., Pollutibacter soli gen. nov. sp. nov., Polluticoccus soli gen. nov. sp. nov., and Terrimonas pollutisoli sp. nov., four new members of the family Chitinophagaceae from polluted soil

A taxonomic investigation was conducted on four bacterial strains isolated from soil contaminated with polycyclic aromatic hydrocarbons and heavy metals. Phylogenetic analysis revealed that these strains belonged to the family Chitinophagaceae. Examination of the 16S rRNA genes indicated that their sequence identities were below 97.6 % compared to any known and validly nominated bacterial species. The genomes of the four strains ranged from 4.12 to 8.76 Mb, with overall G + C molar contents varying from 41.28 % to 50.39 %. Predominant cellular fatty acids included iso-C15:0, iso-C15:1 G, and iso-C17:0 3-OH. The average nucleotide identity ranged from 66.90 % to 74.63 %, and digital DNA-DNA hybridization was 12.5-12.8 %. Based on the genomic and phenotypic features of the new strains, four novel species and two new genera were proposed within the family Chitinophagaceae. The ecological distributions were investigated by data-mining of NCBI databases, and results showed that additional strains or species of the newly proposed taxa were widely distributed in various environments, including polluted soil and waters. Functional analysis demonstrated that strains H1-2-19XT, JS81T, and JY13-12T exhibited resistance to arsenite (III) and chromate (VI). The proposed names for the four novel species are Paraflavitalea pollutisoli (type strain H1-2-19XT = JCM 36460T = CGMCC 1.61321T), Terrimonas pollutisoli (type strain H1YJ31T = JCM 36215T = CGMCC 1.61343T), Pollutibacter soli (type strain JS81T = JCM 36462T = CGMCC 1.61338T), and Polluticoccus soli (type strain JY13-12T = JCM 36463T = CGMCC 1.61341T).

RNA-Seq Data Analysis: A Practical Guide for Model and Non-Model Organisms

RNA sequencing (RNA-seq) has emerged as a powerful tool for assessing genome-wide gene expression, revolutionizing various fields of biology. However, analyzing large RNA-seq datasets can be challenging, especially for students or researchers lacking bioinformatics experience. To address these challenges, we present a comprehensive guide to provide step-by-step workflows for analyzing RNA-seq data, from raw reads to functional enrichment analysis, starting with considerations for experimental design. This is designed to aid students and researchers working with any organism, irrespective of whether an assembled genome is available. Within this guide, we employ various recognized bioinformatics tools to navigate the landscape of RNA-seq analysis and discuss the advantages and disadvantages of different tools for the same task. Our protocol focuses on clarity, reproducibility, and practicality to enable users to navigate the complexities of RNA-seq data analysis easily and gain valuable biological insights from the datasets. Additionally, all scripts and a sample dataset are available in a GitHub repository to facilitate the implementation of the analysis pipeline. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Analysis of data from a model plant with an available reference genome Basic Protocol 2: Gene ontology enrichment analysis Basic Protocol 3: De novo assembly of data from non-model plants.

Single-cell transcriptomics: background, technologies, applications, and challenges

Single-cell sequencing was developed as a high-throughput tool to elucidate unusual and transient cell states that are barely visible in the bulk. This technology reveals the evolutionary status of cells and differences between populations, helps to identify unique cell subtypes and states, reveals regulatory relationships between genes, targets and molecular mechanisms in disease processes, tumor heterogeneity, the state of the immune environment, etc. However, the high cost and technical limitations of single-cell sequencing initially prevented its widespread application, but with advances in research, numerous new single-cell sequencing techniques have been discovered, lowering the cost barrier. Many single-cell sequencing platforms and bioinformatics methods have recently become commercially available, allowing researchers to make fascinating observations. They are now increasingly being used in various industries. Several protocols have been discovered in this context and each technique has unique characteristics, capabilities and challenges. This review presents the latest advancements in single-cell transcriptomics technologies. This includes single-cell transcriptomics approaches, workflows and statistical approaches to data processing, as well as the potential advances, applications, opportunities and challenges of single-cell transcriptomics technology. You will also get an overview of the entry points for spatial transcriptomics and multi-omics.

Chromosome-level genome assembly of the yellow-cheek carp Elopichthys bambusa

Yellow-cheek carp (Elopichthys bambusa) is a typical large and ferocious carnivorous fish endemic to East Asia, with high growth rate, nutritional value and economic value. In this study, a chromosome-level genome of yellow-cheek carp was generated by combining PacBio reads, Illumina reads and Hi-C data. The genome size is 827.63 Mb with a scaffold N50 size of 33.65 Mb, and 99.51% (823.61 Mb) of the assembled sequences were anchored to 24 pseudo-chromosomes. The genome is predicted to contain 24,153 protein-coding genes, with 95.54% having functional annotations. Repeat elements account for approximately 55.17% of the genomic landscape. The completeness of yellow-cheek carp genome assembly is highlighted by a BUSCO score of 98.4%. This genome will help us understand the genetic diversity of yellow-cheek carp and facilitate its conservation planning.

GWAS determined genetic loci associated with callus induction in oil palm tissue culture

KEY MESSAGE

GWAS identified six loci at 25 kb downstream of WAK2, a crucial gene for cell wall and callus formation, enabling development of a SNP marker for enhanced callus induction potential. Efficient callus induction is vital for successful oil palm tissue culture, yet identifying genomic loci and markers for early detection of genotypes with high potential of callus induction remains unclear. In this study, immature male inflorescences from 198 oil palm accessions (dura, tenera and pisifera) were used as explants for tissue culture. Callus induction rates were collected at one-, two- and three-months after inoculation (C1, C2 and C3) as phenotypes. Resequencing generated 11,475,258 high quality single nucleotide polymorphisms (SNPs) as genotypes. GWAS was then performed, and correlation analysis revealed a positive association of C1 with both C2 (R = 0.81) and C3 (R = 0.50), indicating that C1 could be used as the major phenotype for callus induction rate. Therefore, only significant SNPs (P ≤ 0.05) in C1 were identified to develop markers for screening individuals with high potential of callus induction. Among 21 significant SNPs in C1, LD block analysis revealed six SNPs on chromosome 12 (Chr12) potentially linked to callus formation. Subsequently, 13 SNP markers were identified from these loci and electrophoresis results showed that marker C-12 at locus Chr12_12704856 can be used effectively to distinguish the GG allele, which showed the highest probability (69%) of callus induction. Furthermore, a rapid SNP variant detection method without electrophoresis was established via qPCR-based melting curve analysis. Our findings facilitated marker-assisted selection for specific palms with high potential of callus induction using immature male inflorescence as explant, aiding ortet palm selection in oil palm tissue culture.

Comparative Genomics of Lotus japonicus Reveals Insights into Proanthocyanidin Accumulation and Abiotic Stress Response

Lotus japonicus, is an important perennial model legume, has been widely used for studying biological processes such as symbiotic nitrogen fixation, proanthocyanidin (PA) biosynthesis, and abiotic stress response. High-quality L. japonicus genomes have been reported recently; however, the genetic basis of genes associated with specific characters including proanthocyanidin distribution in most tissues and tolerance to stress has not been systematically explored yet. Here, based on our previous high-quality L. japonicus genome assembly and annotation, we compared the L. japonicus MG-20 genome with those of other legume species. We revealed the expansive and specific gene families enriched in secondary metabolite biosynthesis and the detection of external stimuli. We suggested that increased copy numbers and transcription of PA-related genes contribute to PA accumulation in the stem, petiole, flower, pod, and seed coat of L. japonicus. Meanwhile, According to shared and unique transcription factors responding to five abiotic stresses, we revealed that MYB and AP2/ERF play more crucial roles in abiotic stresses. Our study provides new insights into the key agricultural traits of L. japonicus including PA biosynthesis and response to abiotic stress. This may provide valuable gene resources for legume forage abiotic stress resistance and nutrient improvement.

A naturally isolated symbiotic bacterium suppresses flavivirus transmission by Aedes mosquitoes

The commensal microbiota of the mosquito gut plays a complex role in determining the vector competence for arboviruses. In this study, we identified a bacterium from the gut of field Aedes albopictus mosquitoes named Rosenbergiella sp. YN46 (Rosenbergiella_YN46) that rendered mosquitoes refractory to infection with dengue and Zika viruses. Inoculation of 1.6 × 103 colony forming units (CFUs) of Rosenbergiella_YN46 into A. albopictus mosquitoes effectively prevents viral infection. Mechanistically, this bacterium secretes glucose dehydrogenase (RyGDH), which acidifies the gut lumen of fed mosquitoes, causing irreversible conformational changes in the flavivirus envelope protein that prevent viral entry into cells. In semifield conditions, Rosenbergiella_YN46 exhibits effective transstadial transmission in field mosquitoes, which blocks transmission of dengue virus by newly emerged adult mosquitoes. The prevalence of Rosenbergiella_YN46 is greater in mosquitoes from low-dengue areas (52.9 to ~91.7%) than in those from dengue-endemic regions (0 to ~6.7%). Rosenbergiella_YN46 may offer an effective and safe lead for flavivirus biocontrol.

Dissecting Holistic Metabolic Acclimatization of Mucor circinelloides WJ11 Defective in Carotenoid Biosynthesis

Mucor circinelloides WJ11 is a lipid-producing strain with industrial potential. A holistic approach using gene manipulation and bioprocessing development has improved lipid production and the strain's economic viability. However, the systematic regulation of lipid accumulation and carotenoid biosynthesis in M. circinelloides remains unknown. To dissect the metabolic mechanism underlying lipid and carotenoid biosynthesis, transcriptome analysis and reporter metabolites identification were implemented between the wild-type (WJ11) and ΔcarRP WJ11 strains of M. circinelloides. As a result, transcriptome analysis revealed 10,287 expressed genes, with 657 differentially expressed genes (DEGs) primarily involved in amino acid, carbohydrate, and energy metabolism. Integration with a genome-scale metabolic model (GSMM) identified reporter metabolites in the ΔcarRP WJ11 strain, highlighting metabolic pathways crucial for amino acid, energy, and nitrogen metabolism. Notably, the downregulation of genes associated with carotenoid biosynthesis and acetyl-CoA generation suggests a coordinated relationship between the carotenoid and fatty acid biosynthesis pathways. Despite disruptions in the carotenoid pathway, lipid production remains stagnant due to reduced acetyl-CoA availability, emphasizing the intricate metabolic interplay. These findings provide insights into the coordinated relationship between carotenoid and fatty acid biosynthesis in M. circinelloides that are valuable in applied research to design optimized strains for producing desired bioproducts through emerging technology.

Integrated transcriptomic and metabolomic data reveal the cold stress responses molecular mechanisms of two coconut varieties

Among tropical fruit trees, coconut holds significant edible and economic importance. The natural growth of coconuts faces a challenge in the form of low temperatures, which is a crucial factor among adverse environmental stresses impacting their geographical distribution. Hence, it is essential to enhance our comprehension of the molecular mechanisms through which cold stress influences various coconut varieties. We employed analyses of leaf growth morphology and physiological traits to examine how coconuts respond to low temperatures over 2-hour, 8-hour, 2-day, and 7-day intervals. Additionally, we performed transcriptome and metabolome analyses to identify the molecular and physiological shifts in two coconut varieties displaying distinct sensitivities to the cold stress. As the length of cold stress extended, there was a prominent escalation within the soluble protein (SP), proline (Pro) concentrations, the activity of peroxidase (POD) and superoxide dismutase (SOD) in the leaves. Contrariwise, the activity of glutathione peroxidase (GSH) underwent a substantial reduction during this period. The widespread analysis of metabolome and transcriptome disclosed a nexus of genes and metabolites intricately cold stress were chiefly involved in pathways centered around amino acid, flavonoid, carbohydrate and lipid metabolism. We perceived several stress-responsive metabolites, such as flavonoids, carbohydrates, lipids, and amino acids, which unveiled considerably, lower in the genotype subtle to cold stress. Furthermore, we uncovered pivotal genes in the amino acid biosynthesis, antioxidant system and flavonoid biosynthesis pathway that presented down-regulation in coconut varieties sensitive to cold stress. This study broadly enriches our contemporary perception of the molecular machinery that contributes to altering levels of cold stress tolerance amid coconut genotypes. It also unlocks several unique prospects for exploration in the areas of breeding or engineering, aiming to identifying tolerant and/or sensitive coconut varieties encompassing multi-omics layers in response to cold stress conditions.

Unveiling Gene Expression Dynamics during Early Embryogenesis in Cynoglossus semilaevis: A Transcriptomic Perspective

Commencing with sperm-egg fusion, the early stages of metazoan development include the cleavage and formation of blastula and gastrula. These early embryonic events play a crucial role in ontogeny and are accompanied by a dramatic remodeling of the gene network, particularly encompassing the maternal-to-zygotic transition. Nonetheless, the gene expression dynamics governing early embryogenesis remain unclear in most metazoan lineages. We conducted transcriptomic profiling on two types of gametes (oocytes and sperms) and early embryos (ranging from the four-cell to the gastrula stage) of an economically valuable flatfish-the Chinese tongue sole Cynoglossus semilaevis (Pleuronectiformes: Cynoglossidae). Comparative transcriptome analysis revealed that large-scale zygotic genome activation (ZGA) occurs in the blastula stage, aligning with previous findings in zebrafish. Through the comparison of the most abundant transcripts identified in each sample and the functional analysis of co-expression modules, we unveiled distinct functional enrichments across different gametes/developmental stages: actin- and immune-related functions in sperms; mitosis, transcription inhibition, and mitochondrial function in oocytes and in pre-ZGA embryos (four- to 1000-cell stage); and organ development in post-ZGA embryos (blastula and gastrula). These results provide insights into the intricate transcriptional regulation of early embryonic development in Cynoglossidae fish and expand our knowledge of developmental constraints in vertebrates.

The complete Chloroplast genome of Stachys geobombycis and comparative analysis with related Stachys species

Herb genomics, at the forefront of traditional Chinese medicine research, combines genomics with traditional practices, facilitating the scientific validation of ancient remedies. This integration enhances public understanding of traditional Chinese medicine's efficacy and broadens its scope in modern healthcare. Stachys species encompass annual or perennial herbs or small shrubs, exhibiting simple petiolate or sessile leaves. Despite their wide-ranging applications across various fields, molecular data have been lacking, hindering the precise identification and taxonomic elucidation of Stachys species. To address this gap, we assembled the complete chloroplast (CP) genome of Stachys geobombycis and conducted reannotation and comparative analysis of seven additional species within the Stachys genus. The findings demonstrate that the CP genomes of these species exhibit quadripartite structures, with lengths ranging from 14,523 to 150,599 bp. Overall, the genome structure remains relatively conserved, hosting 131 annotated genes, including 87 protein coding genes, 36 tRNA genes, and 8 rRNA genes. Additionally, 78 to 98 SSRs and long repeat sequences were detected , and notably, 6 highly variable regions were identified as potential molecular markers in the CP genome through sequence alignment. Phylogenetic analysis based on Bayesian inference and maximum likelihood methods strongly supported the phylogenetic position of the genus Stachys as a member of Stachydeae tribe. Overall, this comprehensive bioinformatics study of Stachys CP genomes lays the groundwork for phylogenetic classification, plant identification, genetic engineering, evolutionary studies, and breeding research concerning medicinal plants within the Stachys genus.

Gonadal transcriptome sequencing reveals sexual dimorphism in expression profiling of sex-related genes in Asian arowana (Scleropages formosus)

Asia arowana (Scleropages formosus) is an ornamental fish with high economic value, while its sex determination mechanism is still poorly understood. By far, no morphological evidence or molecular marker has been developed for effective distinguishment of genders, which poses a critical challenge to our captive breeding efforts. In this study, we sequenced gonadal transcriptomes of adult Asian arowanas and revealed differential expression profiling of sex-related genes. Based on the comparative transcriptomics analysis of testes (n = 3) and ovaries (n = 3), we identified a total of 8,872 differentially expressed genes (DEGs) and 18,490 differentially expressed transposable elements (TEs) between male and female individuals. Interestingly, the expression of TEs usually has been more significantly testis-biased than related coding genes. As expected, several genes related to females (such as foxl2 and cyp19a1a) are significantly transcribed in the ovary, and some genes related to male gonad development (such as dmrt1, gsdf and amh) are highly expressed in the testis. This sexual dimorphism is valuable for ascertaining the differential expression patterns of sex-related genes and enriching the genetic resources of this economically important species. These valuable genetic materials thereby provide instructive references for gender identification and one-to-one breeding practices so as to expand fish numbers for a rapid elevation of economic value.

Draft genome of Microsporidia sp. MB-a malaria-blocking microsporidian symbiont of the Anopheles arabiensis

We report the draft whole-genome assembly of Microsporidia sp. MB, a symbiotic malaria-transmission-blocking microsporidian isolated from Anopheles arabiensis in Kenya. The whole-genome sequence of Microsporidia sp. MB has a length of 5,908,979 bp, 2,335 contigs, and an average GC content of 31.12%.

Deep-sea in situ and laboratory multi-omics provide insights into the sulfur assimilation of a deep-sea Chloroflexota bacterium

Chloroflexota bacteria are abundant and globally distributed in various deep-sea ecosystems. It has been reported based on metagenomics data that two deep-sea Chloroflexota lineages (the SAR202 group and Dehalococcoidia class) have the potential to drive sulfur cycling. However, the absence of cultured Chloroflexota representatives is a significant bottleneck toward understanding their contribution to the deep-sea sulfur cycling. In this study, we find that Phototrophicus methaneseepsis ZRK33 isolated from deep-sea sediment has a heterotrophic lifestyle and can assimilate sulfate and thiosulfate. Using combined physiological, genomic, proteomic, and in situ transcriptomic methods, we find that strain ZRK33 can perform assimilatory sulfate reduction in both laboratory and deep-sea conditions. Metabolism of sulfate or thiosulfate by strain ZRK33 significantly promotes the transport and degradation of various macromolecules and thereby stimulates the energy production. In addition, metagenomic results show that genes associated with assimilatory and dissimilatory sulfate reduction are ubiquitously distributed in the metagenome-assembled genomes of Chloroflexota members derived from deep-sea sediments. Metatranscriptomic results also show that the expression levels of related genes are upregulated, strongly suggesting that Chloroflexota bacteria may play undocumented roles in deep-sea sulfur cycling.

IMPORTANCE

The cycling of sulfur is one of Earth's major biogeochemical processes and is closely related to the energy metabolism of microorganisms living in the deep-sea cold seep and hydrothermal vents. To date, some of the members of Chloroflexota are proposed to play a previously unrecognized role in sulfur cycling. However, the sulfur metabolic characteristics of deep-sea Chloroflexota bacteria have never been reported, and remain to be verified in cultured deep-sea representatives. Here, we show that the deep-sea Chloroflexota bacterium ZRK33 can perform sulfate assimilation in both laboratory and deep-sea conditions, which expands our knowledge of the sulfur metabolic potential of deep-sea Chloroflexota bacteria. We also show that the genes associated with assimilatory and dissimilatory sulfate reduction ubiquitously distribute in the deep-sea Chloroflexota members, providing hints to the roles of Chloroflexota bacteria in deep-sea sulfur biogeochemical cycling.

Substantial viral diversity in bats and rodents from East Africa: insights into evolution, recombination, and cocirculation

BACKGROUND

Zoonotic viruses cause substantial public health and socioeconomic problems worldwide. Understanding how viruses evolve and spread within and among wildlife species is a critical step when aiming for proactive identification of viral threats to prevent future pandemics. Despite the many proposed factors influencing viral diversity, the genomic diversity and structure of viral communities in East Africa are largely unknown.

RESULTS

Using 38.3 Tb of metatranscriptomic data obtained via ultradeep sequencing, we screened vertebrate-associated viromes from 844 bats and 250 rodents from Kenya and Uganda collected from the wild. The 251 vertebrate-associated viral genomes of bats (212) and rodents (39) revealed the vast diversity, host-related variability, and high geographic specificity of viruses in East Africa. Among the surveyed viral families, Coronaviridae and Circoviridae showed low host specificity, high conservation of replication-associated proteins, high divergence among viral entry proteins, and frequent recombination. Despite major dispersal limitations, recurrent mutations, cocirculation, and occasional gene flow contribute to the high local diversity of viral genomes.

CONCLUSIONS

The present study not only shows the landscape of bat and rodent viromes in this zoonotic hotspot but also reveals genomic signatures driven by the evolution and dispersal of the viral community, laying solid groundwork for future proactive surveillance of emerging zoonotic pathogens in wildlife. Video Abstract.

Differential expression and co-localization of transcriptional factors during callus transition to differentiation for shoot organogenesis in the water fern Ceratopteris richardii

BACKGROUND AND AIMS

In flowering plants, regeneration can be achieved by a variety of approaches, and different sets of transcriptional factors are involved in these processes. However, regeneration in taxa other than flowering plants remains a mystery. Ceratopteris richardii is a representative fern capable of both direct and indirect organogenesis, and we aimed to investigate the genetics underlying the transition from callus proliferation to differentiation.

METHODS

Morphological and histological analyses were used to determine the type of regeneration involved. RNA sequencing and differential gene expression were used to investigate how the callus switches from proliferation to differentiation. Phylogenetic analysis and RNA in situ hybridization were used to understand whether transcriptional factors are involved in this transition.

KEY RESULTS

The callus formed on nascent leaves and subsequently developed the shoot pro-meristem and shoot meristem, thus completing indirect de novo shoot organogenesis in C. richardii. Genes were differentially expressed during the callus transition from proliferation to differentiation, indicating a role for photosynthesis, stimulus response and transmembrane signalling in this transition and the involvement of almost all cell layers that make up the callus. Transcriptional factors were either downregulated or upregulated, which were generally in many-to-many orthology with genes known to be involved in callus development in flowering plants, suggesting that the genetics of fern callus development are both conserved and divergent. Among them, an STM-like, a PLT-like and an ethylene- and salt-inducible ERF gene3-like gene were expressed simultaneously in the vasculature but not in the other parts of the callus, indicating that the vasculature played a role in the callus transition from proliferation to differentiation.

CONCLUSIONS

Indirect de novo shoot organogenesis could occur in ferns, and the callus transition from proliferation to differentiation required physiological changes, differential expression of transcriptional factors and involvement of the vasculature.

Intergenic sequences harboring potential enhancer elements contribute to Axenfeld-Rieger syndrome by regulating PITX2

Recent studies have uncovered that noncoding sequence variants may relate to Axenfeld-Rieger syndrome (ARS), a rare developmental anomaly with genetic heterogeneity. However, how these genomic regions are functionally and structurally associated with ARS is still unclear. In this study, we performed genome-wide linkage analysis and whole-genome sequencing in a Chinese family with ARS and identified a heterozygous deletion of about 570 kb (termed LOH-1) in the intergenic sequence between paired-like homeodomain transcription factor 2 (PITX2) and family with sequence similarity 241 member A. Knockout of LOH-1 homologous sequences caused ARS phenotypes in mice. RNA-Seq and real-time quantitative PCR revealed a significant reduction in Pitx2 gene expression in LOH-1-/- mice, while forkhead box C1 expression remained unchanged. ChIP-Seq and bioinformatics analysis identified a potential enhancer region (LOH-E1) within LOH-1. Deletion of LOH-E1 led to a substantial downregulation of the PITX2 gene. Mechanistically, we found a sequence (hg38 chr4:111,399,594-111,399,691) that is on LOH-E1 could regulate PITX2 by binding to RAD21, a critical component of the cohesin complex. Knockdown of RAD21 resulted in reduced PITX2 expression. Collectively, our findings indicate that a potential enhancer sequence that is within LOH-1 may regulate PITX2 expression remotely through cohesin-mediated loop domains, leading to ARS when absent.

Unraveling the molecular landscape of Ataxia Telangiectasia: Insights into Neuroinflammation, immune dysfunction, and potential therapeutic target

BACKGROUND

Ataxia Telangiectasia (AT) is a genetic disorder characterized by compromised DNA repair, cerebellar degeneration, and immune dysfunction. Understanding the molecular mechanisms driving AT pathology is crucial for developing targeted therapies.

METHODS

In this study, we conducted a comprehensive analysis to elucidate the molecular mechanisms underlying AT pathology. Using publicly available RNA-seq datasets comparing control and AT samples, we employed in silico transcriptomics to identify potential genes and pathways. We performed differential gene expression analysis with DESeq2 to reveal dysregulated genes associated with AT. Additionally, we constructed a Protein-Protein Interaction (PPI) network to explore the interactions between proteins implicated in AT.

RESULTS

The network analysis identified hub genes, including TYROBP and PCP2, crucial in immune regulation and cerebellar function, respectively. Furthermore, pathway enrichment analysis unveiled dysregulated pathways linked to AT pathology, providing insights into disease progression.

CONCLUSION

Our integrated approach offers a holistic understanding of the complex molecular landscape of AT and identifies potential targets for therapeutic intervention. By combining transcriptomic analysis with network-based methods, we provide valuable insights into the underlying mechanisms of AT pathogenesis.

Genomics comparisons of three chromosome-level mudskipper genome assemblies reveal molecular clues for water-to-land evolution and adaptation

INTRODUCTION

Mudskippers are a large group of amphibious fishes that have developed many morphological and physiological capacities to live on land. Genomics comparisons of chromosome-level genome assemblies of three representative mudskippers, Boleophthalmus pectinirostris (BP), Periophthalmus magnuspinnatus (PM) and P. modestus (PMO), may be able to provide novel insights into the water-to-land evolution and adaptation.

METHODS

Two chromosome-level genome assemblies for BP and PM were respectively sequenced by an integration of PacBio, Nanopore and Hi-C sequencing. A series of standard assembly and annotation pipelines were subsequently performed for both mudskippers. We also re-annotated the PMO genome, downloaded from NCBI, to obtain a redundancy-reduced annotation. Three-way comparative analyses of the three mudskipper genomes in a large scale were carried out to discover detailed genomic differences, such as different gene sizes, and potential chromosomal fission and fusion events. Comparisons of several representative gene families among the three amphibious mudskippers and some other teleosts were also performed to find some molecular clues for terrestrial adaptation.

RESULTS

We obtained two high-quality haplotype genome assemblies with 23 and 25 chromosomes for BP and PM respectively. We also found two specific chromosome fission events in PM. Ancestor chromosome analysis has discovered a common fusion event in mudskipper ancestor. This fusion was then retained in all the three mudskipper species. A loss of some SCPP (secretory calcium-binding phosphoprotein) genes were identified in the three mudskipper genomes, which could lead to reduction of scales for a part-time terrestrial residence. The loss of aanat1a gene, encoding an important enzyme (arylalkylamine N-acetyltransferase 1a, AANAT1a) for dopamine metabolism and melatonin biosynthesis, was confirmed in PM but not in PMO (as previously reported existence in BP), suggesting a better air vision of PM than both PMO and BP. Such a tiny variation within the genus Periophthalmus exemplifies to prove a step-by-step evolution for the mudskippers' water-to-land adaptation.

CONCLUSION

These high-quality mudskipper genome assemblies will become valuable genetic resources for in-depth discovery of genomic evolution for the terrestrial adaptation of amphibious fishes.

Genomics of post-bottleneck recovery in the northern elephant seal

Populations and species are threatened by human pressure, but their fate is variable. Some depleted populations, such as that of the northern elephant seal (Mirounga angustirostris), recover rapidly even when the surviving population was small. The northern elephant seal was hunted extensively and taken by collectors between the early 1800s and 1892, suffering an extreme population bottleneck as a consequence. Recovery was rapid and now there are over 200,000 individuals. We sequenced 260 modern and 8 historical northern elephant seal nuclear genomes to assess the impact of the population bottleneck on individual northern elephant seals and to better understand their recovery. Here we show that inbreeding, an increase in the frequency of alleles compromised by lost function, and allele frequency distortion, reduced the fitness of breeding males and females, as well as the performance of adult females on foraging migrations. We provide a detailed investigation of the impact of a severe bottleneck on fitness at the genomic level and report on the role of specific gene systems.

Effects of Cold Plasma Pretreatment on the Synthesis of Polysaccharide from Pleurotus ostreatus

Fungal polysaccharides have attracted wide attention because of their medical pharmaceutical and health care value. So far, many efforts have been made in strain improvement to produce polysaccharides on a large scale at low cost. Here, a novel cold plasma-induced strain improvement technology was employed to pretreat Pleurotus ostreatus CGMCC 5.374 by radio-frequency (RF) low-vacuum cold plasma (LVCP) for the purpose of obtaining a high-yield polysaccharide strain. The optimum pretreatment conditions including discharge power, treatment time, and working pressure were determined by single factor and orthogonal experiment in succession. Furthermore, transcriptome analysis was conducted to study the effects of RF-LVCP on cell metabolism and proliferation. Results showed that under the optimal condition of discharge power of 130 W, treatment time of 25 s and working pressure of 140 Pa, polysaccharide content in mycelium was increased by 3.16% after 6 days in comparison to the original strain. Transcriptome analysis showed that RF-LVCP is helpful for specific gene transcription profiles, Gene Ontology (GO) and KEGG pathways, of which the differentially expressed genes (DEGs) were mainly involve with the up-regulation of polysaccharide transport, physiology, synthesis and metabolism, as well as the down-regulation of polysaccharide hydrolysis and macromolecular degradation.

Comparative genomics of a novel Erwinia species associated with the Highland midge (Culicoides impunctatus)

Erwinia (Enterobacterales: Erwiniaceae) are a group of cosmopolitan bacteria best known as the causative agents of various plant diseases. However, other species in this genus have been found to play important roles as insect endosymbionts supplementing the diet of their hosts. Here, I describe Candidatus Erwinia impunctatus (Erwimp) associated with the Highland midge Culicoides impunctatus (Diptera: Ceratopogonidae), an abundant biting pest in the Scottish Highlands. The genome of this new Erwinia species was assembled using hybrid long and short read techniques, and a comparative analysis was undertaken with other members of the genus to understand its potential ecological niche and impact. Genome composition analysis revealed that Erwimp is similar to other endophytic and ectophytic species in the genus and is unlikely to be restricted to its insect host. Evidence for an additional plant host includes the presence of a carotenoid synthesis operon implicated as a virulence factor in plant-associated members in the sister genus Pantoea. Unique features of Erwimp include several copies of intimin-like proteins which, along with signs of genome pseudogenization and a loss of certain metabolic pathways, suggests an element of host restriction seen elsewhere in the genus. Furthermore, a screening of individuals over two field seasons revealed the absence of the bacteria in Culicoides impunctatus during the second year indicating this microbe-insect interaction is likely to be transient. These data suggest that Culicoides impunctatus may have an important role to play beyond a biting nuisance, as an insect vector transmitting Erwimp alongside any conferred impacts to surrounding biota.

Transcriptomics aids in uncovering the metabolic shifts and molecular machinery of Schizochytrium limacinum during biotransformation of hydrophobic substrates to docosahexaenoic acid

BACKGROUND

Biotransformation of waste oil into value-added nutraceuticals provides a sustainable strategy. Thraustochytrids are heterotrophic marine protists and promising producers of omega (ω) fatty acids. Although the metabolic routes for the assimilation of hydrophilic carbon substrates such as glucose are known for these microbes, the mechanisms employed for the conversion of hydrophobic substrates are not well established. Here, thraustochytrid Schizochytrium limacinum SR21 was investigated for its ability to convert oils (commercial oils with varying fatty acid composition and waste cooking oil) into ω-3 fatty acid; docosahexaenoic acid (DHA).

RESULTS

Within 72 h SR21 consumed ~ 90% of the oils resulting in enhanced biomass (7.5 g L- 1) which was 2-fold higher as compared to glucose. Statistical analysis highlights C16 fatty acids as important precursors of DHA biosynthesis. Transcriptomic data indicated the upregulation of multiple lipases, predicted to possess signal peptides for secretory, membrane-anchored and cytoplasmic localization. Additionally, transcripts encoding for mitochondrial and peroxisomal β-oxidation along with acyl-carnitine transporters were abundant for oil substrates that allowed complete degradation of fatty acids to acetyl CoA. Further, low levels of oxidative biomarkers (H2O2, malondialdehyde) and antioxidants were determined for hydrophobic substrates, suggesting that SR21 efficiently mitigates the metabolic load and diverts the acetyl CoA towards energy generation and DHA accumulation.

CONCLUSIONS

The findings of this study contribute to uncovering the route of assimilation of oil substrates by SR21. The thraustochytrid employs an intricate crosstalk among the extracellular and intracellular molecular machinery favoring energy generation. The conversion of hydrophobic substrates to DHA can be further improved using synthetic biology tools, thereby providing a unique platform for the sustainable recycling of waste oil substrates.

Transcriptomic and biochemical analysis of metabolic remodeling in Bacillus subtilis MSC4 under Benzo[a]pyrene stress

Polyaromatic benzo[a]pyrene (B[a]P) is a toxic carcinogenic environmental pollutant, and the use of microorganisms to remediate B[a]P contamination is considered to be one of the most effective strategies. However, there is still a gap in studying the metabolic remodeling of microorganisms under B[a]P stress. In this study, our systematically investigated the effects of B[a]P on the metabolism of Bacillus subtilis MSC4 based on transcriptomic, molecular and biochemical analyses. The results showed that in response to B[a]P stress, MSC4 formed more biofilm matrix and endospores, the structure of the endospores also was changed, which led to a reduction in their resistance and made them more difficult to germinate. In addition to an increase in glycolysis activity, the activities of tricarboxylic acid cycle, pentose phosphate pathway and the electron transport chain were decreased. B[a]P stress forced MSC4 to strengthen arginine synthesis, urea cycle, and urea decomposition, meanwhile, synthesize more ribonucleotides. The activity of DNA replication, transcription activities and the expression of multiple ribosomal protein genes were reduced. Moreover, all of the reported enzymes involved in B[a]P degradation showed decreased transcript abundance, and the degradation of B[a]P caused significant up-regulation of the gene expression of the acid inducible enzyme OxdC and the synthesis of acetoin. In addition, the cytotoxicity of B[a]P to bacteria was directly displayed in four aspects: increased intracellular level of reactive oxygen species (ROS), elevated cell membrane permeability, up-regulation of the cell envelope stress-sensing two-component system LiaRS, and downregulation of siderophores biosynthesis. Finally, B[a]P also caused morphological changes in the cells, with some cells exhibiting significant deformation and concavity. These findings provide effective research directions for targeted improvement the cellular activity of B[a]P-degrading strains, and is beneficial for further application of microorganisms to remediate B[a]P -contaminated soils.

Feasibility of functional precision medicine for guiding treatment of relapsed or refractory pediatric cancers

Children with rare, relapsed or refractory cancers often face limited treatment options, and few predictive biomarkers are available that can enable personalized treatment recommendations. The implementation of functional precision medicine (FPM), which combines genomic profiling with drug sensitivity testing (DST) of patient-derived tumor cells, has potential to identify treatment options when standard-of-care is exhausted. The goal of this prospective observational study was to generate FPM data for pediatric patients with relapsed or refractory cancer. The primary objective was to determine the feasibility of returning FPM-based treatment recommendations in real time to the FPM tumor board (FPMTB) within a clinically actionable timeframe (<4 weeks). The secondary objective was to assess clinical outcomes from patients enrolled in the study. Twenty-five patients with relapsed or refractory solid and hematological cancers were enrolled; 21 patients underwent DST and 20 also completed genomic profiling. Median turnaround times for DST and genomics were within 10 days and 27 days, respectively. Treatment recommendations were made for 19 patients (76%), of whom 14 received therapeutic interventions. Six patients received subsequent FPM-guided treatments. Among these patients, five (83%) experienced a greater than 1.3-fold improvement in progression-free survival associated with their FPM-guided therapy relative to their previous therapy, and demonstrated a significant increase in progression-free survival and objective response rate compared to those of eight non-guided patients. The findings from our proof-of-principle study illustrate the potential for FPM to positively impact clinical care for pediatric and adolescent patients with relapsed or refractory cancers and warrant further validation in large prospective studies. ClinicalTrials.gov registration: NCT03860376 .

Whole transcriptome scanning and validation of negatively related genes in UC-MSCs

Background

Human umbilical cord mesenchymal stem cells (UC-MSCs) are one of the most extensively researched stem cell types due to their potential for multi-lineage differentiation, secretion of regenerative factors, modulations of immunological activities, and the release of regenerative substances and influence immunological processes. Since UC-MSCs must be cultivated on a large scale for clinical use, selecting the appropriate storing passage, such as the usage-based passage of UC-MSCs, is critical for long-term autologous or allogeneic usage. Long-term cultivation of stem cells, on the other hand, causes them to lose their pluripotent differentiation capacity. As a result, distinguishing between high and low passages of UC-MSCs and identifying the particular variations associated with stem cells and their modes of action is essential for regenerative medicine. Therefore, we investigated the biological features and transcriptional changes of UC-MSCs over passages.

Methods

UC-MSCs were isolated from the tissues of the human umbilical cord, and UC-MSCs from five passages (P1, P3, P5, P10 and P15) with three repetitions were compared and identified based on morphology, cell markers, differentiation capacity, and aging-related characteristics. It was previously assumed that the phenotype of cells before the P10 passage was stable, defined as early passage, and that culture could be continued until the 15th passage, defined as late passage. Next, the five passages of UC-MSCs were sequenced using high-throughput complete transcriptome sequencing. Fuzzy C-Means Clustering (FCM) and Weighted Gene Co-expression Network Analysis (WGCNA) were used to find hub genes, and gene silencing was performed to investigate the impact of missing genes on the stemness of UC-MSC cells.

Results

UC-MSCs of different passages displayed similar surface markers, including CD73, CD105, CD90, CD34, CD45 and HLA-DR. However, the proliferation time of late-phase UC-MSCs was longer than that of early-phase UC-MSCs, and the expression of the senescence-associated (SA)-β-gal staining marker was higher. At the same time, pluripotency markers (NANOG, OCT4, SOX2 and KIF4A) were down-regulated, and the multi-differentiation potential was reduced. Meanwhile, KIFC1 and UBE2C were down-regulated in late-phase UC-MSCs, which were involved in the maintenance of stemness.

Conclusions

KIFC1 and UBE2C were highly expressed in early-UC-MSCs and showed a downward gradient trend with cell expansion in vitro. They regulated UC-MSC proliferation, colony sphere formation, multiple differentiation, stemness maintenance, and other biological manifestations. Therefore, they are anticipated to be new biomarkers for UC-MSCs quality identification in regenerative medicine applications.

A high-quality chromosome-level genome assembly of the Chinese medaka Oryzias sinensis

Oryzias sinensis, also known as Chinese medaka or Chinese ricefish, is a commonly used animal model for aquatic environmental assessment in the wild as well as gene function validation or toxicology research in the lab. Here, a high-quality chromosome-level genome assembly of O. sinensis was generated using single-tube long fragment read (stLFR) reads, Nanopore long-reads, and Hi-C sequencing data. The genome is 796.58 Mb, and a total of 712.17 Mb of the assembled sequences were anchored to 23 pseudo-chromosomes. A final set of 22,461 genes were annotated, with 98.67% being functionally annotated. The Benchmarking Universal Single-Copy Orthologs (BUSCO) benchmark of genome assembly and gene annotation reached 95.1% (93.3% single-copy) and 94.6% (91.7% single-copy), respectively. Furthermore, we also use ATAC-seq to uncover chromosome transposase-accessibility as well as related genome area function enrichment for Oryzias sinensis. This study offers a new improved foundation for future genomics research in Chinese medaka.

Transcriptome-wide analysis of the differences between MCF7 cells cultured in DMEM or αMEM

MCF7 cells have been used as an experimental model for breast cancer for decades. Typically, a culture medium is designed to supply cells with the nutrients essential for their continuous proliferation. Each medium has a specific nutritional composition. Therefore, cells cultured in different media may exhibit differences in their metabolism. However, only a few studies have investigated the effects of media on cells. In this study, we compared the effects of Dulbecco's modified Eagle medium (DMEM) and minimum essential medium alpha modification (αMEM) on MCF7 cells. The two media differentially affected the morphology, cell cycle, and proliferation of MCF7 cells, but had no effect on cell death. Replacement of DMEM with αMEM led to a decrease in ATP production and an increase in reactive oxygen species production, but did not affect the cell viability. RNA-sequencing and bioinformatic analyses revealed 721 significantly upregulated and 1247 downregulated genes in cells cultured in αMEM for 48 h compared with that in cells cultured in DMEM. The enriched gene ontology terms were related to mitosis and cell proliferation. Kyoto encyclopedia of genes and genomes analysis revealed cell cycle and DNA replication as the top two significant pathways. MCF7 cells were hypoxic when cultured in αMEM. These results show that the culture medium considerably affects cultured cells. Thus, the stability of the culture system in a study is very important to obtain reliable results.

Comparative Metabolomics and Transcriptome Studies of Two Forms of Rhododendron chrysanthum Pall. under UV-B Stress

Rhododendron chrysanthum Pall. (R. chrysanthum), a plant with UV-B resistance mechanisms that can adapt to alpine environments, has gained attention as an important plant resource with the ability to cope with UV-B stress. In this experiment, R. chrysanthums derived from the same origin were migrated to different culture environments (artificial climate chamber and intelligent artificial incubator) to obtain two forms of R. chrysanthum. After UV-B irradiation, 404 metabolites and 93,034 unigenes were detected. Twenty-six of these different metabolites were classified as UV-B-responsive metabolites. Glyceric acid is used as a potential UV-B stress biomarker. The domesticated Rhododendron chrysanthum Pall. had high amino acid and SOD contents. The study shows that the domesticated Rhododendron chrysanthum Pall. has significant UV-B resistance. The transcriptomics results show that the trends of DEGs after UV-B radiation were similar for both forms of R. chrysanthum: cellular process and metabolic process accounted for a higher proportion in biological processes, cellular anatomical entity accounted for the highest proportion in the cellular component, and catalytic activity and binding accounted for the highest proportion in the molecular function category. Through comparative study, the forms of metabolites resistant to UV-B stress in plants can be reflected, and UV-B radiation absorption complexes can be screened for application in future specific practices. Moreover, by comparing the differences in response to UV-B stress between the two forms of R. chrysanthum, references can be provided for cultivating domesticated plants with UV-B stress resistance characteristics. Research on the complex mechanism of plant adaptation to UV-B will be aided by these results.

Overexpression of SLIM1 transcription factor accelerates vegetative development in Arabidopsis thaliana

The transcription factor Sulfur Limitation 1 (SLIM1) belongs to the plant-specific Ethylene Insenstive3-Like transcription factor family and is known to coordinate gene expression in response to sulfur deficiency. However, the roles of SLIM1 in nutrient-sufficient conditions have not been characterized. Employing constitutive SLIM1 overexpression (35S::SLIM1) and CRISPR/Cas9 mutant plants (slim1-cr), we identified several distinct phenotypes in nutrient-sufficient conditions in Arabidopsis thaliana. Overexpression of SLIM1 results in plants with approximately twofold greater rosette area throughout vegetative development. 35S::SLIM1 plants also bolt earlier and exhibit earlier downregulation of photosynthesis-associated genes and earlier upregulation of senescence-associated genes than Col-0 and slim1-cr plants. This suggests that overexpression of SLIM1 accelerates development in A. thaliana. Genome-wide differential gene expression analysis relative to Col-0 at three time points with slim1-cr and two 35S::SLIM1 lines allowed us to identify 1,731 genes regulated directly or indirectly by SLIM1 in vivo.

Elucidation of the genetic determination of body weight and size in Chinese local chicken breeds by large-scale genomic analyses

BACKGROUND

Body weight and size are important economic traits in chickens. While many growth-related quantitative trait loci (QTLs) and candidate genes have been identified, further research is needed to confirm and characterize these findings. In this study, we investigate genetic and genomic markers associated with chicken body weight and size. This study provides new insights into potential markers for genomic selection and breeding strategies to improve meat production in chickens.

METHODS

We performed whole-genome resequencing of and Wenshang Barred (WB) chickens (n = 596) and three additional breeds with varying body sizes (Recessive White (RW), WB, and Luxi Mini (LM) chickens; (n = 50)). We then used selective sweeps of mutations coupled with genome-wide association study (GWAS) to identify genomic markers associated with body weight and size.

RESULTS

We identified over 9.4 million high-quality single nucleotide polymorphisms (SNPs) among three chicken breeds/lines. Among these breeds, 287 protein-coding genes exhibited positive selection in the RW and WB populations, while 241 protein-coding genes showed positive selection in the LM and WB populations. Genomic heritability estimates were calculated for 26 body weight and size traits, including body weight, chest breadth, chest depth, thoracic horn, body oblique length, keel length, pelvic width, shank length, and shank circumference in the WB breed. The estimates ranged from 0.04 to 0.67. Our analysis also identified a total of 2,522 genome-wide significant SNPs, with 2,474 SNPs clustered around two genomic regions. The first region, located on chromosome 4 (7.41-7.64 Mb), was linked to body weight after ten weeks and body size traits. LCORL, LDB2, and PPARGC1A were identified as candidate genes in this region. The other region, located on chromosome 1 (170.46-171.53 Mb), was associated with body weight from four to eighteen weeks and body size traits. This region contained CAB39L and WDFY2 as candidate genes. Notably, LCORL, LDB2, and PPARGC1A showed highly selective signatures among the three breeds of chicken with varying body sizes.

CONCLUSION

Overall this study provides a comprehensive map of genomic variants associated with body weight and size in chickens. We propose two genomic regions, one on chromosome 1 and the other on chromosome 4, that could helpful for developing genome selection breeding strategies to enhance meat yield in chickens.

Magnesium oxide nanoparticles reduce clubroot by regulating plant defense response and rhizosphere microbial community of tumorous stem mustard (Brassica juncea var. tumida)

Clubroot, caused by Plasmodiophora brassicae, is a major disease that significantly impairs the yield of cruciferous crops and causes significant economic losses across the globe. The prevention of clubroot, especially in tumorous stem mustard (without resistant varieties), are is limited and primarily relies on fungicides. Engineered nanoparticles have opened up new avenues for the management of plant diseases, but there is no report on their application in the prevention of clubroot. The results showed that the control efficacy of 500 mg/L MgO NPs against clubroot was 54.92%. However, when the concentration was increased to 1,500 and 2,500 mg/L, there was no significant change in the control effect. Compared with CK, the average fresh and dry weight of the aerial part of plants treated with MgO NPs increased by 392.83 and 240.81%, respectively. Compared with the F1000 treatment, increases were observed in the content of soil available phosphorus (+16.72%), potassium (+9.82%), exchangeable magnesium (+24.20%), and water-soluble magnesium (+20.64%) in the 1,500 mg/L MgO NPs treatment. The enzyme-linked immune sorbent assay (ELISA) results showed that the application of MgO NPs significantly increased soil peroxidase (POD, +52.69%), alkaline protease (AP, +41.21%), alkaline phosphatase (ALP, +79.26%), urease (+52.69%), and sucrase (+56.88%) activities; And also increased plant L-phenylalanine ammonla-lyase (PAL, +70.49%), polyphenol oxidase (PPO, +36.77%), POD (+38.30%), guaiacol peroxidase (POX, +55.46%) activities and salicylic acid (SA, +59.86%) content. However, soil and plant catalase (CAT, -27.22 and - 19.89%, respectively), and plant super oxidase dismutase (SOD, -36.33%) activities were significantly decreased after the application of MgO NPs. The metagenomic sequencing analysis showed that the MgO NPs treatments significantly improved the α-diversity of the rhizosphere soil microbial community. The relative abundance of beneficial bacteria genera in the rhizosphere soil, including Pseudomonas, Sphingopyxis, Acidovorax, Variovorax, and Bosea, was significantly increased. Soil metabolic functions, such as oxidative phosphorylation (ko00190), carbon fixation pathways in prokaryotes (ko00720), indole alkaloid biosynthesis (ko00901), and biosynthesis of various antibiotics (ko00998) were significantly enriched. These results suggested that MgO NPs might control clubroot by promoting the transformation and utilization of soil nutrients, stimulating plant defense responses, and enriching soil beneficial bacteria.

Characteristics of germline DNA damage response gene mutations in ovarian cancer in Southwest China

DNA damage response (DDR) pathways are responsible for repairing endogenous or exogenous DNA damage to maintain the stability of the cellular genome, including homologous recombination repair (HRR) pathway, mismatch repair (MMR) pathway, etc. In ovarian cancer, current studies are focused on HRR genes, especially BRCA1/2, and the results show regional and population differences. To characterize germline mutations in DDR genes in ovarian cancer in Southwest China, 432 unselected ovarian cancer patients underwent multi-gene panel testing from October 2016 to October 2020. Overall, deleterious germline mutations in DDR genes were detected in 346 patients (80.1%), and in BRCA1/2 were detected in 126 patients (29.2%). The prevalence of deleterious germline mutations in BRCA2 is higher than in other studies (patients are mainly from Eastern China), and so is the mismatch repair genes. We identified three novel BRCA1/2 mutations, two of which probably deleterious (BRCA1 p.K1622* and BRCA2 p.L2987P). Furthermore, we pointed out that deleterious mutations of FNACD2 and RECQL4 are potential ovarian cancer susceptibility genes and may predispose carriers to ovarian cancer. In conclusion, our study highlights the necessity of comprehensive germline mutation detection of DNA damage response genes in ovarian cancer patients, which is conducive to patient management and genetic counseling.

Lose-lose consequences of bacterial community-driven invasions in soil

BACKGROUND

Community-driven invasion, also known as community coalescence, occurs widely in natural ecosystems. Despite that, our knowledge about the process and mechanisms controlling community-driven invasion in soil ecosystems is lacking. Here, we performed a set of coalescence experiments in soil microcosms and assessed impacts up to 60 days after coalescence by quantifying multiple traits (compositional, functional, and metabolic) of the invasive and coalescent communities.

RESULTS

Our results showed that coalescences significantly triggered changes in the resident community's succession trajectory and functionality (carbohydrate metabolism), even when the size of the invasive community is small (~ 5% of the resident density) and 99% of the invaders failed to survive. The invasion impact was mainly due to the high suppression of constant residents (65% on average), leading to a lose-lose situation where both invaders and residents suffered with coalescence. Our results showed that surviving residents could benefit from the coalescence, which supports the theory of "competition-driven niche segregation" at the microbial community level. Furthermore, the result showed that both short- and long-term coalescence effects were predicted by similarity and unevenness indexes of compositional, functional, and metabolic traits of invasive communities. This indicates the power of multi-level traits in monitoring microbial community succession. In contrast, the varied importance of different levels of traits suggests that competitive processes depend on the composition of the invasive community.

CONCLUSIONS

Our results shed light on the process and consequence of community coalescences and highlight that resource competition between invaders and residents plays a critical role in soil microbial community coalescences. These findings provide valuable insights for understanding and predicting soil microbial community succession in frequently disturbed natural and agroecosystems. Video Abstract.

Stimulation of PSTPIP1 to trigger proinflammatory responses in asymptomatic SARS-CoV-2 infections

Background

A hyperinflammatory response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection gravely worsens the clinical progression of coronavirus disease 2019 (COVID-19). Although the undesirable effects of inflammasome activation have been correlated to the severity of COVID-19, the mechanisms of this process in the asymptomatic infection and disease progression have not yet been clearly elucidated.

Methods

We performed strand-specific RNA sequencing in 39 peripheral blood mononuclear cell (PBMC) samples from asymptomatic individuals（n = 10）, symptomatic patients（n = 16） and healthy donors（n = 13）.

Results

Dysregulation of pyrin inflammasomes along with the proline-serine-threonine phosphatase-interacting protein 1 (PSTPIP1) gene was identified in SARS-COV-2 infection. Notably, the PSTPIP1 expression level showed a significant negative correlation with an adjacent long-noncoding RNA (lncRNA) RP11-797A18.6 in the asymptomatic individuals compared with the healthy controls. In addition, a decline in the nuclear factor kappa B subunit 1 (NFKB1) gene expression was observed in asymptomatic infection, followed by a rise in the mild and moderate disease stages, suggesting that altered NFKB1 expression and associated proinflammatory signals may trigger a disease progression.

Conclusions

Overall, our results indicate that PSTPIP1-dependent pyrin inflammasomes-mediated pyroptosis and NF-κB activation might be potential preventive targets for COVID-19 disease development and progression.

Thermal acclimation uncovers a simple genetic basis of adaptation to high temperature in a cosmopolitan pest

Understanding a population's fitness heterogeneity and genetic basis of thermal adaptation is essential for predicting the responses to global warming. We examined the thermotolerance and genetic adaptation of Plutella xylostella to exposure to hot temperatures. The population fitness parameters of the hot-acclimated DBM strains varied in the thermal environments. Using genome scanning and transcription profiling, we find a number of genes potentially involved in thermal adaptation of DBM. Editing two ABCG transporter genes, PxWhite and PxABCG, confirmed their role in altering cuticle permeability and influencing thermal responses. Our results demonstrate that SNP mutations in genes and changes in gene expression can allow DBM to rapidly adapt to thermal environment. ABCG transporter genes play an important role in thermal adaptation of DBM. This work improves our understanding of genetic adaptation mechanisms of insects to thermal stress and our capacity to predict the effects of rising global temperatures on ectotherms.

Structural and functional characterization of genes PYL-PP2C-SnRK2s in the ABA signalling pathway of Cucurbita pepo

BACKGROUND

The core regulation of the abscisic acid (ABA) signalling pathway comprises the multigenic families PYL, PP2C, and SnRK2. In this work, we conducted a genome-wide study of the components of these families in Cucurbita pepo.

RESULTS

The bioinformatic analysis of the C. pepo genome resulted in the identification of 19 CpPYL, 102 CpPP2C and 10 CpSnRK2 genes. The investigation of gene structure and protein motifs allowed to define 4 PYL, 13 PP2C and 3 SnRK2 subfamilies. RNA-seq analysis was used to determine the expression of these gene families in different plant organs, as well as to detect their differential gene expression during germination, and in response to ABA and cold stress in leaves. The specific tissue expression of some gene members indicated the relevant role of some ABA signalling genes in plant development. Moreover, their differential expression under ABA treatment or cold stress revealed those ABA signalling genes that responded to ABA, and those that were up- or down-regulated in response to cold stress. A reduced number of genes responded to both treatments. Specific PYL-PP2C-SnRK2 genes that had potential roles in germination were also detected, including those regulated early during the imbibition phase, those regulated later during the embryo extension and radicle emergence phase, and those induced or repressed during the whole germination process.

CONCLUSIONS

The outcomes of this research open new research lines for agriculture and for assessing gene function in future studies.

Light and dopamine impact two circadian neurons to promote morning wakefulness

In both mammals and flies, circadian brain neurons orchestrate physiological oscillations and behaviors like wake and sleep; these neurons can be subdivided by morphology and by gene expression patterns. Recent single-cell sequencing studies identified 17 Drosophila circadian neuron groups. One of these include only two lateral neurons (LNs), which are marked by the expression of the neuropeptide ion transport peptide (ITP). Although these two ITP+ LNs have long been grouped with five other circadian evening activity cells, inhibiting the two neurons alone strongly reduces morning activity; this indicates that they are prominent morning neurons. As dopamine signaling promotes activity in Drosophila like in mammals, we considered that dopamine might influence this morning activity function. Moreover, the ITP+ LNs express higher mRNA levels than other LNs of the type 1-like dopamine receptor Dop1R1. Consistent with the importance of Dop1R1, CRISPR/Cas9 mutagenesis of this receptor only in the two ITP+ LNs renders flies significantly less active in the morning, and ex vivo live imaging shows that dopamine increases cAMP levels in these two neurons; cell-specific mutagenesis of Dop1R1 eliminates this cAMP response to dopamine. Notably, the response is more robust in the morning, reflecting higher morning Dop1R1 mRNA levels in the two neurons. As morning levels are not elevated in constant darkness, this suggests light-dependent upregulation of morning Dop1R1 transcript levels. Taken together with enhanced morning cAMP response to dopamine, the data indicate how light stimulates morning wakefulness in flies, which mimics the important effect of light on morning wakefulness in humans.

Novel insights into aerobic 17β-estradiol degradation by enriched microbial communities from mangrove sediments

Various persistent organic pollutants (POPs) including estrogens are often enriched in mangrove regions. This research investigated the estrogens pollution levels in six mangroves located in the Southern China. The estrogen levels were found to be in the range of 5.3-24.9 ng/g dry weight, suggesting that these mangroves had been seriously contaminated. The bacterial communities under estrogen stress were further enriched by supplementing 17β-estradiol (E2) as the sole carbon source. The enriched bacterial communities showed an excellent E2 degradation capacity > 95 %. These communities were able to transform E2 into estrone (E1), 4-hydroxy-estrone, and keto-estrone, etc. 16 S rDNA sequencing and metagenomics analysis revealed that bacterial taxa Oleiagrimonas, Pseudomonas, Terrimonas, and Nitratireductor etc. were the main contributors to estrogen degradation. Moreover, the genes involved in E2 degradation were enriched in the microbial communities, including the genes encoding 17β-hydroxysteroid dehydrogenase, estrone 4-hydroxylase, etc. Finally, the analyses of functional genes and binning genomes demonstrated that E2 was degraded by bacterial communities via dehydrogenation into E1 by 17β-hydroxysteroid dehydrogenase. E1 was then catabolically converted to 3aα-H-4α(3'-propanoate)- 7aβ-methylhexahydro-1,5-indanedione via 4,5-seco pathway. Alternatively, E1 could also be hydroxylated to keto-estrone, followed by B-ring cleavage. This study provides novel insights into the biodegradation of E2 by the bacterial communities in estrogen-contaminated mangroves.

An integrated transcriptome-microbiome host relationship associated with paraben toxicity in the brackish water flea Diaphanosoma celebensis

Parabens, a group of alkyl esters of p-hydroxybenzoic acid, have been found in aquatic systems in particular, leading to concerns about their potential impact on ecosystems. This study investigated the effects of three commonly used parabens, methylparaben (MeP), ethylparaben (EtP), and propylparaben (PrP), on the brackish water flea Diaphanosoma celebensis. The results showed that PrP had the most adverse impact on survival rates, followed by EtP and MeP, while MeP and EtP induced significant adverse effects on reproductive performance. A transcriptome analysis revealed significant differential gene expression patterns in response to paraben exposure, with MeP associated with the most significant effects. MeP and EtP exposure produced greater disruption in the microbiota of D. celebensis than did PrP compared with control groups, and we identified eight key microbiota, including Ruegeria and Roseovarius. Correlation analysis between transcriptome and microbiome data revealed key interactions between specific microbiota and host gene expression. Certain microbial taxa were associated with specific genes (e.g. cuticle related genes) and toxicological pathways, shedding light on the complex molecular response and in vivo toxicity effects of parabens. These findings contribute to a deeper understanding of the molecular mechanisms underlying paraben toxicity and highlight the importance of considering the ecological impact of chemical contaminants in aquatic ecosystems.

In situ community transcriptomics illuminates CO2-fixation potentials and supporting roles of phagotrophy and proton pump in plankton in a subtropical marginal sea

Lineage-wise physiological activities of plankton communities in the ocean are important but challenging to characterize. Here, we conducted whole-assemblage metatranscriptomic profiling at continental shelf and slope sites in the South China Sea to investigate carbon fixation potential in different lineages. RuBisCO expression, the proxy of Calvin carbon fixation (CCF) potential, was mainly contributed by Bacillariophyta, Chlorophyta, Cyanobacteria, and Haptophyta, which was differentially affected by environmental factors among lineages. CCF potential exhibited positive or negative correlations with phagotrophy gene expression, suggesting phagotrophy possibly enhances or complements CCF. Our data also reveal significant non-Calvin carbon fixation (NCF) potential, as indicated by the active expression of genes in all five currently recognized NCF pathways, mainly contributed by Flavobacteriales, Alteromonadales, and Oceanospirillales. Furthermore, in Flavobacteriales, Alteromonadales, Pelagibacterales, and Rhodobacterales, NCF potential was positively correlated with proton-pump rhodopsin (PPR) expression, suggesting that NCF might be energetically supported by PPR. The novel insights into the lineage-differential potential of carbon fixation, widespread mixotrophy, and PPR as an energy source for NCF lay a methodological and informational foundation for further research to understand carbon fixation and the trophic landscape in the ocean.IMPORTANCEMarine plankton plays an important role in global carbon cycling and climate regulation. Phytoplankton and cyanobacteria fix CO2 to produce organic compounds using solar energy and mainly by the Calvin cycle, whereas autotrophic bacteria and archaea may fix CO2 by non-Calvin cycle carbon fixation pathways. How active individual lineages are in carbon fixation and mixotrophy, and what energy source bacteria may employ in non-Calvin carbon fixation, in a natural plankton assemblage are poorly understood and underexplored. Using metatranscriptomics, we studied carbon fixation in marine plankton with lineage resolution in tropical marginal shelf and slope areas. Based on the sequencing results, we characterized the carbon fixation potential of different lineages and assessed Calvin- and non-Calvin- carbon fixation activities and energy sources. Data revealed a high number of unigenes (4.4 million), lineage-dependent differential potentials of Calvin carbon fixation and responses to environmental conditions, major contributors of non-Calvin carbon fixation, and their potential energy source.