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

Long-term co-circulation of multiple influenza A viruses in pigs, Guangxi, China

Influenza A viruses (IAVs) pose a persistent potential threat to human health because of the spillover from avian and swine infections. Extensive surveillance was performed in 12 cities of Guangxi, China, during 2018 and 2023. A total of 2540 samples (including 2353 nasal swabs and 187 lung tissues) were collected from 18 pig farms with outbreaks of respiratory disease. From these, 192 IAV-positive samples and 19 genomic sequences were obtained. We found that the H1 and H3 swine influenza A viruses (swIAVs) of multiple lineages and genotypes have continued to co-circulate during that time in this region. Genomic analysis revealed the Eurasian avian-like H1N1 swIAVs (G4) still remained predominant in pig populations. Strikingly, the novel multiple H3N2 genotypes were found to have been generated through the repeated introduction of the early H3N2 North American triple reassortant viruses (TR H3N2 lineage) that emerged in USA and Canada in 1998 and 2005, respectively. Notably, when the matrix gene segment derived from the H9N2 avian influenza virus was introduced into endemic swIAVs, this produced a novel quadruple reassortant H1N2 swIAV that could pose a potential risk for zoonotic infection.

Visualizing metagenomic and metatranscriptomic data: A comprehensive review

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

In silico identification and functional study of long non-coding RNA involved in acute hepatopancreatic necrosis disease caused by Vibrio parahaemolyticus infection in white shrimp, Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) caused by toxin-producing Vibrio parahaemolyticus (VpAHPND) has severely affected shrimp production. Long non-coding RNA (lncRNA), a regulatory non-coding RNA, which can play important function in shrimp disease responses. This study aimed to identify and investigate the role of lncRNA involved in VpAHPND infection in Pacific white shrimp, Litopenaeus vannamei. From a total of 368,736 de novo assembled transcripts, 67,559 were identified as putative lncRNAs, and only 72 putative lncRNAs showed differential expression between VpAHPND-infected and normal shrimp. The six candidate lncRNAs were validated for their expression profiles during VpAHPND infection and tissue distribution using RT-qPCR. The role of lnc2088 in response to VpAHPND infection was investigated through RNA interference. The result indicated that the suppression of lnc2088 expression led to an increase in shrimp mortality after VpAHPND infection. To explore the set of genes involved in lnc2088 knockdown, RNA sequencing was performed. A total of 275 differentially expressed transcripts were identified in the hepatopancreas of lnc2088 knockdown shrimp. The expression profiles of five candidate metabolic and immune-related genes were validated in lnc2088 knockdown and VpAHPND-infected shrimp. The result showed that the expression of ChiNAG was significantly increased, while that of NCBP1, WIPF2, and NFKB1 was significantly downregulated in ds2088-injected shrimp. Additionally, the expression of NFKB1, NCBP1 and WIPF2 was significantly increased, whereas that of ChiNAG and CUL5 were significantly decreased after infection with VpAHPND. Our work identified putative lncRNA profiles in L. vannamei in response to VpAHPND infection and investigated the role of lncRNA in shrimp immunity.

Investigation of parenteral nutrition-induced hepatotoxicity using human liver spheroid co-cultures

Parenteral nutrition (PN) is typically administered to individuals with gastrointestinal dysfunction, a contraindication for enteral feeding, and a need for nutritional therapy. When PN is the only energy source in patients, it is defined as total parenteral nutrition (TPN). TPN is a life-saving approach for different patient populations, both in infants and adults. However, despite numerous benefits, TPN can cause adverse effects, including metabolic disorders and liver injury. TPN-associated liver injury, known as intestinal failure-associated liver disease (IFALD), represents a significant problem affecting up to 90% of individuals receiving TPN. IFALD pathogenesis is complex, depending on the TPN components as well as on the patient's medical conditions. Despite numerous animal studies and clinical observations, the molecular mechanisms driving IFALD remain largely unknown. The present study was set up to elucidate the mechanisms underlying IFALD. For this purpose, human liver spheroid co-cultures were treated with a TPN mixture, followed by RNA sequencing analysis. Subsequently, following exposure to TPN and its single nutritional components, several key events of liver injury, including mitochondrial dysfunction, endoplasmic reticulum stress, oxidative stress, apoptosis, and lipid accumulation (steatosis), were studied using various techniques. It was found that prolonged exposure to TPN substantially changes the transcriptome profile of liver spheroids and affects multiple metabolic and signaling pathways contributing to liver injury. Moreover, TPN and its main components, especially lipid emulsion, induce changes in all key events measured and trigger steatosis.

Characterization of runs of homozygosity islands in American mink using whole-genome sequencing data

The genome-wide analysis of runs of homozygosity (ROH) islands can be an effective strategy for identifying shared variants within a population and uncovering important genomic regions related to complex traits. The current study performed ROH analysis to characterize the genome-wide patterns of homozygosity, identify ROH islands and annotated genes within these candidate regions using whole-genome sequencing data from 100 American mink (Neogale vison). After sequence processing, variants were called using GATK and Samtools pipelines. Subsequent to quality control, 8,373,854 bi-allelic variants identified by both pipelines remained for further analysis. A total of 34,652 ROH segments were identified in all individuals, among which shorter segments (0.3-1 Mb) were abundant throughout the genome, approximately accounting for 84.39% of all ROH. Within these segments, we identified 63 ROH islands housing 156 annotated genes. The genes located in ROH islands were associated with fur quality (EDNRA, FGF2, FOXA2 and SLC24A4), body size/weight (MYLK4, PRIM2, FABP2, EYS and PHF3), immune capacity (IL2, IL21, PTP4A1, SEMA4C, JAK2, CCNA2 and TNIP3) and reproduction (ADAD1, KHDRBS2, INSL6, PGRMC2 and HSPA4L). Furthermore, Gene Ontology and KEGG pathway enrichment analyses revealed 56 and 9 significant terms (FDR-corrected p-value < 0.05), respectively, among which cGMP-PKG signalling pathway, regulation of actin cytoskeleton, and calcium signalling pathway were highlighted due to their functional roles in growth and fur characteristics. This is the first study to present ROH islands in American mink. The candidate genes from ROH islands and functional enrichment analysis suggest possible signatures of selection in response to the mink breeding targets, such as increased body length, reproductive performance and fur quality. These findings contribute to our understanding of genetic characteristics, and provide complementary information to assist with implementation of breeding strategies for genetic improvement in American mink.

In vivo detection of endogenous toxic phenolic compounds of intestine

Phenol and p-cresol are two common toxic small molecules related to various diseases. Existing reports confirmed that high L-tyrosine in the daily diet can increase the concentration of phenolic compounds in blood and urine. L-tyrosine is a common component of protein-rich foods. Some anaerobic bacteria in the gut can convert non-toxic l-tyrosine into these two toxic phenolic compounds, phenol and p-cresol. Existing methods have been constructed for measuring the concentration of phenolic compound in feces. However, there is still a lack of direct visual evidence to measure the phenolic compounds in the intestine. In this study, we aimed to construct a whole-cell biosensor for phenolic compounds detection based on the dmpR, the regulator from the phenol metabolism cluster. The commensal bacterium Citrobacter amalonaticus PS01 was selected and used as the chassis. Compared with the biosensor based on ECN1917, the biosensor PS01[dmpR] could better implant into the mouse gut through gavage and showed a higher sensitive to phenolic compound. And the concentration of phenolic compounds in the intestines could be observed with the help of in vivo imaging system using PS01[dmpR]. This paper demonstrated endogenous phenol synthesis in the gut and the strategy of using commensal bacteria to construct whole-cell biosensors for detecting small molecule compounds in the intestines.

Obesity-enriched gut microbe degrades myo-inositol and promotes lipid absorption

Numerous studies have reported critical roles for the gut microbiota in obesity. However, the specific microbes that causally contribute to obesity and the underlying mechanisms remain undetermined. Here, we conducted shotgun metagenomic sequencing in a Chinese cohort of 631 obese subjects and 374 normal-weight controls and identified a Megamonas-dominated, enterotype-like cluster enriched in obese subjects. Among this cohort, the presence of Megamonas and polygenic risk exhibited an additive impact on obesity. Megamonas rupellensis possessed genes for myo-inositol degradation, as demonstrated in vitro and in vivo, and the addition of myo-inositol effectively inhibited fatty acid absorption in intestinal organoids. Furthermore, mice colonized with M. rupellensis or E. coli heterologously expressing the myo-inositol-degrading iolG gene exhibited enhanced intestinal lipid absorption, thereby leading to obesity. Altogether, our findings uncover roles for M. rupellensis as a myo-inositol degrader that enhances lipid absorption and obesity, suggesting potential strategies for future obesity management.

The chromosome-level genome assembly of the red swamp crayfish Procambarus clarkii

Red swamp crayfish, Procambarus clarkii, is the most cultured freshwater crayfish species. It attracts significant research attention due to its considerable economic importance. However, the limited availability of genome information has impeded further genetic studies and breeding programs. By utilizing Illumina, PacBio, and Hi-C sequencing technologies, we present a more comprehensive and continuous chromosome-level assembly for P. clarkii than the published one. The final genome size is 4.03 Gb, consisting of 2,358 scaffolds with a N50 of 42.87 Mb. Notably, 3.68 Gb, corresponding to 91.42% of the genome, was anchored to 94 chromosomes. The assembly comprises 70.64% repetitive sequences, including 5.21% tandem repeats and 65.40% transposable elements. Additionally, a total of 4,456 non-coding RNAs and 28,852 protein-coding genes were predicted in the P. clarkii genome, with 96.26% of the genes were annotated. This high-quality genome assembly not only represents a significant improvement for the genome of P. clarkii and provides insights into the unique genome evolution, but also offers valuable information for developing freshwater aquaculture and accelerating genetic breeding.

Integrated analyses of the intestinal microbiome and transcriptome in Ningxiang piglets

Ningxiang (NX) pig has been recognized as one of the most famous Chinese indigenous breeds due to its characteristics in stress resistance. However, intestinal microbial feature and gene profiling in NX piglets have not been studied. Here, we compared the intestinal microbiome and transcriptome between NX and Duroc × Landrace × Large white (DLY) piglets and found the high enrichment of several colonic Bacteroides, Prevotella and Clostridium species in NX piglets. Further functional analyses revealed their predominant function in methane, glycolysis and gluconeogenesis metabolism. Our mRNA-sequencing data unraveled the distinct colonic gene expression between these two breeds. In particular, we showed that the improved intestinal function in NX piglets may be determined by enhanced intestinal barrier gene expression and varied immune gene expression through modulating the composition of the gut microbes. Together, our study revealed the intestinal characteristics of NX piglets, providing their potential application in improving breeding strategies and developing dietary interventions.

Ether-Diol Ambiguity: An Inconspicuous Issue in the Structure Elucidation of Oxygenated Natural Products

Tertiary and allylic hydroxyl groups readily eliminate water during positive ion mode mass spectrometry and may show similar NMR spectra to their corresponding ethers. In a routine structure elucidation workflow, these factors can cause researchers to incorrectly assign diol moieties as ethers or vice versa, leading to inaccurate chemical structures. After facing this problem during our work on oxygenated sesquiterpenoids from a Fusarium sp. fungal strain, we became aware of this challenging issue. We examined the literature for oxygenated natural products bearing these functional groups, and with the aid of density functional calculations of NMR chemical shifts, we now report the structures of 15 natural products that should be revised. We further establish that derivatizing sub-micromolar amounts of alcohols to their sulfates can be used to distinguish these from their corresponding ethers using liquid chromatography negative ion mode mass spectrometry. Finally, we isolated lignoren/cyclonerodiol from the Fusarium sp. culture extract and supported its revised identity as cyclonerodiol using this sulfation approach. Our results suggest that ether-diol ambiguity could be a prevalent issue affecting the structure elucidation of oxygenated natural products and highlight the importance of using complementary techniques, such as sulfation with LC-(-)-ESI-MS or density functional calculations of NMR chemical shifts.

Genomic variation, environmental adaptation, and feralization in ramie, an ancient fiber crop

Feralization is an important evolutionary process, but the mechanisms behind it remain poorly understood. Here, we use the ancient fiber crop ramie (Boehmeria nivea (L.) Gaudich.) as a model to investigate genomic changes associated with both domestication and feralization. We first produced a chromosome-scale de novo genome assembly of feral ramie and investigated structural variations between feral and domesticated ramie genomes. Next, we gathered 915 accessions from 23 countries, comprising cultivars, major landraces, feral populations, and the wild progenitor. Based on whole-genome resequencing of these accessions, we constructed the most comprehensive ramie genomic variation map to date. Phylogenetic, demographic, and admixture signal detection analyses indicated that feral ramie is of exoferal or exo-endo origin, i.e., descended from hybridization between domesticated ramie and the wild progenitor or ancient landraces. Feral ramie has higher genetic diversity than wild or domesticated ramie, and genomic regions affected by natural selection during feralization differ from those under selection during domestication. Ecological analyses showed that feral and domesticated ramie have similar ecological niches that differ substantially from the niche of the wild progenitor, and three environmental variables are associated with habitat-specific adaptation in feral ramie. These findings advance our understanding of feralization, providing a scientific basis for the excavation of new crop germplasm resources and offering novel insights into the evolution of feralization in nature.

Defective endomembrane dynamics in Rab27a deficiency impairs nucleic acid sensing and cytokine secretion in immune cells

Endosomal Toll-like receptors (eTLRs) are essential for the sensing of non-self through RNA and DNA detection. Here, using spatiotemporal analysis of vesicular dynamics, super-resolution microscopy studies, and functional assays, we show that endomembrane defects associated with the deficiency of the small GTPase Rab27a cause delayed eTLR ligand recognition, defective early signaling, and impaired cytokine secretion. Rab27a-deficient neutrophils show retention of eTLRs in amphisomes and impaired ligand internalization. Extracellular signal-regulated kinase (ERK) signaling and β2-integrin upregulation, early responses to TLR7 and TLR9 ligands, are defective in Rab27a deficiency. CpG-stimulated Rab27a-deficient neutrophils present increased tumor necrosis factor alpha (TNF-α) secretion and decreased secretion of a selected group of mediators, including interleukin (IL)-10. In vivo, CpG-challenged Rab27a-null mice show decreased production of type I interferons (IFNs) and IFN-γ, and the IFN-α secretion defect is confirmed in Rab27a-null plasmacytoid dendritic cells. Our findings have significant implications for immunodeficiency, inflammation, and CpG adjuvant vaccination.

Prognostic genome and transcriptome signatures in colorectal cancers

Colorectal cancer is caused by a sequence of somatic genomic alterations affecting driver genes in core cancer pathways1. Here, to understand the functional and prognostic impact of cancer-causing somatic mutations, we analysed the whole genomes and transcriptomes of 1,063 primary colorectal cancers in a population-based cohort with long-term follow-up. From the 96 mutated driver genes, 9 were not previously implicated in colorectal cancer and 24 had not been linked to any cancer. Two distinct patterns of pathway co-mutations were observed, timing analyses identified nine early and three late driver gene mutations, and several signatures of colorectal-cancer-specific mutational processes were identified. Mutations in WNT, EGFR and TGFβ pathway genes, the mitochondrial CYB gene and 3 regulatory elements along with 21 copy-number variations and the COSMIC SBS44 signature correlated with survival. Gene expression classification yielded five prognostic subtypes with distinct molecular features, in part explained by underlying genomic alterations. Microsatellite-instable tumours divided into two classes with different levels of hypoxia and infiltration of immune and stromal cells. To our knowledge, this study constitutes the largest integrated genome and transcriptome analysis of colorectal cancer, and interlinks mutations, gene expression and patient outcomes. The identification of prognostic mutations and expression subtypes can guide future efforts to individualize colorectal cancer therapy.

Two high-quality Prototheca zopfii genomes provide new insights into their evolution as obligate algal heterotrophs and their pathogenicity

The majority of the nearly 10,000 described species of green algae are photoautotrophs; however, some species have lost their ability to photosynthesize and become obligate heterotrophs that rely on parasitism for survival. Two high-quality genomes of the heterotrophic algae Prototheca zopfii Pz20 and Pz23 were obtained using short- and long-read genomic as well as transcriptomic data. The genome sizes were 31.2 Mb and 31.3 Mb, respectively, and contig N50 values of 1.99 Mb and 1.26 Mb. Although P. zopfii maintained its plastid genome, the transition to heterotrophy led to a reduction in both plastid and nuclear genome size, including the loss of photosynthesis-related genes from both the nuclear and plastid genomes and the elimination of genes encoding for carotenoid oxygenase and pheophorbide an oxygenase. The loss of genes, including basic leucine-zipper (bZIP) transcription factors, flavin adenine dinucleotide-linked oxidase, and helicase, could have played a role in the transmission of autotrophy to heterotrophs and in the processes of abiotic stress resistance and pathogenicity. A total of 66 (1.37%) and 73 (1.49%) genes were identified as potential horizontal gene transfer events in the two P. zopfii genomes, respectively. Genes for malate synthase and isocitrate lyase, which are horizontally transferred from bacteria, may play a pivotal role in carbon and nitrogen metabolism as well as the pathogenicity of Prototheca and non-photosynthetic organisms. The two high-quality P. zopfii genomes provide new insights into their evolution as obligate heterotrophs and pathogenicity.

IMPORTANCE

The genus Prototheca, characterized by its heterotrophic nature and pathogenicity, serves as an exemplary model for investigating pathobiology. The limited understanding of the protothecosis infectious disease is attributed to the lack of genomic resources. Using HiFi long-read sequencing, both nuclear and plastid genomes were generated for two strains of P. zopfii. The findings revealed a concurrent reduction in both plastid and nuclear genome size, accompanied by the loss of genes associated with photosynthesis, carotenoid oxygenase, basic leucine-zipper (bZIP) transcription factors, and others. The analysis of horizontal gene transfer revealed the presence of 1.37% and 1.49% bacterial genes, including malate synthase and isocitrate lyase, which play crucial roles in carbon and nitrogen metabolism, as well as pathogenicity and obligate heterotrophy. The two high-quality P. zopfii genomes represent valuable resources for investigating their adaptation and evolution as obligate heterotrophs, as well as for developing future prevention and treatment strategies against protothecosis.

Telomere-to-telomere genome assembly of sorghum

"Cuohu Bazi" (CHBZ) is an ancient sorghum variety collected from the fields of China, known for its agronomic traits like dwarf stature, early maturation. In this study, we present the first telomere-to-telomere (T2T) and gap-free genome assembly of CHBZ using PacBio HiFi reads, Oxford Nanopore Technologies, and Hi-C data. The assembled genome comprises 724.85 Mb, effectively resolving all 3,913 gaps that were present in the previous sorghum BTx623 reference genome. Notably, the T2T assembly captures 10 centromeres and all 20 telomeres, providing strong support for their integrity. This assembly is of high quality in terms of contiguity (contig N50: 71.1 Mb), completeness (BUSCO score: 99.01%, k-mer completeness: 98.88%), and correctness (QV: 61.60). Repetitive sequences accounted for 70.41% of the genome and a total of 32,855 protein-coding genes have been annotated. Furthermore, 161 CHBZ-specific presence/absence variants genes have been identified when comparing to BTx623 genome. This study provides valuable insights for future research on sorghum genetics, genomics, and evolutionary history.

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.

Establishment of Novel Simple Sequence Repeat (SSR) Markers from Chimonanthus praecox Transcriptome Data and Their Application in the Identification of Varieties

Chimonanthus praecox, a member of the Calycanthaceae family, is a unique, traditional, and famous flowering economic tree species in China. Despite the existence of several varieties, only a few cultivars have been formally named. Currently, expression sequence tag-simple sequence repeat (EST-SSR) markers are extensively used to identify different species and varieties; a large number of microsatellites can be identified from transcriptome databases. A total of 162,638 unigenes were assembled using RNA-seq; 82,778 unigenes were annotated using the Nr, Nt, Swiss-Prot, Pfam, GO, KOG, and KEGG databases. In total, 13,556 SSR loci were detected from 11,691 unigenes, with trinucleotide repeat motifs being the most abundant among the six repeat motifs. To develop the markers, 64,440 pairs of SSR primers with polymorphism potential were designed, and 75 pairs of primers were randomly selected for amplification. Among these markers, seven pairs produced amplified fragments of the expected size with high polymorphism. Using these markers, 12 C. praecox varieties were clustered into two monophyletic clades. Microsatellites in the transcriptome of C. praecox exhibit rich types, strong specificity, and great polymorphism potential. These EST-SSR markers serve as molecular technical methods for identifying different varieties of C. praecox and facilitate the exploration of a large number of candidate genes associated with important traits.

NEXMIF Combined with KIDINS220 Gene Mutation Caused Neurodevelopmental Disorder and Epilepsy: One Case Report

of Medical History: A male infant, 8 months old, was admitted to hospital with cough and fever. The clinical symptoms were found to be mental retardation, obesity, dystonia, movement limitation, and visual retardation. Early development was normal, but after 6 months, the child developed upright head instability, difficulty grasping, and seizures. Symptoms and Signs: The child presents with mental retardation, obesity, increased muscle tone, motor dysfunction, visual impairment, and seizures.

DIAGNOSIS

A whole exon test was performed to detect a neurite extension and migration factor (NEXMIF) gene mutation (NM_001008537.2: c.1042C > T (p. Arg348*)), which is known to be associated with intellectual disability and neurological symptoms. In addition, the test revealed a mutation in the Kinase D interacting substrate of 220 kDa (KIDINS220) gene (NM_020738.2: c.3242_3243insC (p. Leu1082AIafs*5)) with a heterozygous mutation in the father and wild type in the mother.

TREATMENT

The patient was treated with anti-infection, aerosol inhalation, calcium supplement, and anti-epileptic drugs (levetiracetam), and the disease was controlled. Home and hospital rehabilitation is also underway.

CLINICAL OUTCOME

The condition of the child improved after treatment and no seizures occurred again. The patient needs continuous rehabilitation treatment and follow-up observation.

CONCLUSION

For male children with unexplained neurodevelopmental disorders and comorbidities such as obesity, dystonia, and seizures, mutations in related genes such as NEXMIF should be considered. Clinical practice should improve genetic testing as early as possible to provide a basis for genetic counseling.

SMRT sequencing of a full-length transcriptome reveals cold induced alternative splicing in Vitis amurensis root

Alternative splicing enhances diversity at the transcriptional and protein levels that widely involved in plant response to biotic and abiotic stresses. V. amurensis is an extremely cold-tolerant wild grape variety, however, studies on alternative splicing (AS) in amur grape at low temperatures are currently poorly understood. In this study, we analyzed full-length transcriptome and RNA seq data at 0, 2, and 24 h after cold stress in V. amurensis roots. Following quality control and correction, 221,170 high-quality full-length non-concatemer (FLNC) reads were identified. A total of 16,181 loci and 30,733 isoforms were identified. These included 22,868 novel isoforms from annotated genes and 2815 isoforms from 2389 novel genes. Among the distinguished novel isoforms, 673 Long non-coding RNAs (LncRNAs) and 18,164 novel isoforms open reading frame (ORF) region were found. A total of 2958 genes produced 8797 AS events, of which 189 genes were involved in the low-temperature response. Twelve transcription factors show AS during cold treatment and VaMYB108 was selected for initial exploration. Two transcripts, Chr05.63.1 (VaMYB108short) and Chr05.63.2 (VaMYB108normal) of VaMYB108, display up-regulated expression after cold treatment in amur grape roots and are both localized in the nucleus. Only VaMYB108normal exhibits transcriptional activation activity. Overexpression of either VaMYB108short or VaMYB108normal in grape roots leads to increased expression of the other transcript and both increased chilling resistance of amur grape roots. The results improve and supplement the genome annotations and provide insights for further investigation into AS mechanisms during cold stress in V. amurensis.

A chromosome-level genome assembly provides insights into the local adaptation of Tamarix austromongolica in the Yellow River Basin, China

Tamarix austromongolica is endemic to the Yellow River Basin and has adapted to diverse ecological settings in the region, including the arid areas of northwestern China and the saline soil regions of the Yellow River Delta. However, the genetic basis of its local adaptation remains unclear. We report a chromosome-level assembly of the T. austromongolica genome based on PacBio high-fidelity sequencing and Hi-C technology. The 12 pseudochromosomes cover 98.44% of the 1.32 Gb assembly, with a contig N50 of 52.57 Mb and a BUSCO score of 98.2%. The genome comprises 913.6 Mb (68.83%) of repetitive sequences and 22,374 protein-coding genes. Genome evolution analyses suggest that genes under positive selection and significantly expanded gene families have facilitated T. austromongolica's adaptability to diverse environmental factors and high resistance to diseases. Using genotyping-by-sequencing, we conducted population structure and selection analyses of 114 samples from 15 sites. Two genetic groups were identified, and 114 and 289 candidate genes were assigned to the populations of the northwestern and eastern parts of the Yellow River, respectively. Furthermore, we discovered numerous candidate genes associated with high-altitude adaptability and salt tolerance. This research provides valuable genomic resources for the evolutionary study and genetic breeding of tamarisk.

Gut microbial interactions based on network construction and bacterial pairwise cultivation

Association networks are widely applied for the prediction of bacterial interactions in studies of human gut microbiomes. However, the experimental validation of the predicted interactions is challenging due to the complexity of gut microbiomes and the limited number of cultivated bacteria. In this study, we addressed this challenge by integrating in vitro time series network (TSN) associations and co-cultivation of TSN taxon pairs. Fecal samples were collected and used for cultivation and enrichment of gut microbiome on YCFA agar plates for 13 days. Enriched cells were harvested for DNA extraction and metagenomic sequencing. A total of 198 metagenome-assembled genomes (MAGs) were recovered. Temporal dynamics of bacteria growing on the YCFA agar were used to infer microbial association networks. To experimentally validate the interactions of taxon pairs in networks, we selected 24 and 19 bacterial strains from this study and from the previously established human gut microbial biobank, respectively, for pairwise co-cultures. The co-culture experiments revealed that most of the interactions between taxa in networks were identified as neutralism (51.67%), followed by commensalism (21.67%), amensalism (18.33%), competition (5%) and exploitation (3.33%). Genome-centric analysis further revealed that the commensal gut bacteria (helpers and beneficiaries) might interact with each other via the exchanges of amino acids with high biosynthetic costs, short-chain fatty acids, and/or vitamins. We also validated 12 beneficiaries by adding 16 additives into the basic YCFA medium and found that the growth of 66.7% of these strains was significantly promoted. This approach provides new insights into the gut microbiome complexity and microbial interactions in association networks. Our work highlights that the positive relationships in gut microbial communities tend to be overestimated, and that amino acids, short-chain fatty acids, and vitamins are contributed to the positive relationships.

Induction of let-7d-5p miRNA modulates aortic smooth muscle inflammatory signaling and phenotypic switching

BACKGROUND AND AIMS

Activation of vascular smooth muscle cell inflammation is recognised as an important early driver of vascular disease. We have previously identified the let-7 miRNA family as important regulators of inflammation in in vitro and in vivo models of atherosclerosis. Here we investigated a dual statin/let-7d-5p miRNA combination therapy approach to target human aortic SMC (HAoSMC) activation and inflammation.

METHODS

In vitro studies using primary HAoSMCs were performed to investigate the effects of let-7d-5p miRNA overexpression and inhibition. HAoSMCs were treated with combinations of the inflammatory cytokine tumor necrosis factor-α (TNF-α), and atorvastatin or lovastatin. HAoSMC Bulk RNA-seq transcriptomics of HAoSMCs revealed downstream regulatory networks modulated by let-7d-5p miRNA overexpression and statins. Proteome profiler cytokine array, Western blotting and quantitative PCR analyses were performed on HAoSMCs to validate key findings.

RESULTS

Let-7d-5p overexpression significantly attenuated TNF-α-induced upregulation of IL-6, ICAM1, VCAM1, CCL2, CD68, MYOCD gene expression in HAoSMCs (p<0.05). Statins (atorvastatin, lovastatin) significantly attenuated inflammatory gene expression and upregulated Let-7d levels in HAoSMCs (p<0.05). Bulk RNA-seq analysis of a dual Let-7d-5p overexpression/statin therapy in HAoSMCs revealed that let-7d-5p activation and statins converge on key inflammatory pathways (IL-6, IL-1β, TNF-α, IFN-γ). Let-7d-5p overexpression led to reduced expression of the ox-LDL receptor OLR1, and this was associated with lower ox-LDL uptake in HAoSMCs. In silico analysis of smooth muscle cell phenotypic switching shows that overexpression of let-7d-5p in HAoSMCs maintains a contractile phenotype.

CONCLUSIONS

Targeting the Let-7 network alongside statins can modulate HAoSMC activation and attenuate key inflammatory pathway signals.

Dietary α-linolenic acid supplementation enhances resistance to Salmonella Typhimurium challenge in chickens by altering the intestinal mucosal barrier integrity and cecal microbes

Salmonella is an important foodborne pathogen. Given the ban on the use of antibiotics during the egg-laying period in China, finding safe and effective alternatives to antibiotics to reduce Salmonella enterica subsp. enterica serovar Typhimurium (S. Typhimurium) infections in chickens is essential for the prevention and control of this pathogen and the protection of human health. Numerous studies have shown that unsaturated fatty acids have a positive effect on intestinal inflammation and resistance to infection by intestinal pathogens. Here we investigated the protective effect of α-linolenic acid (ALA) against S. Typhimurium infection in chickens and further explored its mechanism of action. We added different proportions of ALA to the feed and observed the effect of ALA on S. Typhimurium colonization using metagenomic sequencing technology and physiological index measurements. The role of gut flora on S. Typhimurium colonization was subsequently verified by fecal microbiota transplantation (FMT). We found that ALA protects chickens from S. Typhimurium infection by reducing intestinal inflammation through remodeling the gut microbiota, up-regulating the expression of ileocecal barrier-related genes, and maintaining the integrity of the intestinal epithelium. Our data suggest that supplementation of feed with ALA may be an effective strategy to alleviate S. Typhimurium infection in chickens.

Isolation and Genomic Characterization of Kadipiro Virus from Mosquitoes in Yunnan, China

Background: Kadipiro virus (KDV) is a species of the new 12 segmented RNA virus grouped under the genus Seadornavirus within the Reoviridae family. It has previously been isolated or detected from mosquito, Odonata, and bat feces in Indonesia, China, and Denmark, respectively. Here, we describe the isolation and characterization of a viral strain from mosquitoes in Yunnan Province, China. Methods: Mosquitoes were collected overnight using light traps in Shizong county, on July 17, 2023. Virus was isolated from the mosquito homogenate and grown using baby hamster kidney and Aedes albopictus (C6/36) cells. Preliminary identification of the virus was performed by agarose gel electrophoresis (AGE). The full-genome sequences of the strain were determined by full-length amplification of cDNAs and sequenced using next-generation sequencing. Results: We isolated a viral strain (SZ_M48) from mosquitoes (Culex tritaeniorhynchus Giles) that caused cytopathogenic effects in C6/36 cells. AGE analysis indicated a genome consisting of 12 segments of double-stranded RNA that demonstrated a "6-5-1" pattern, similar to the migrating bands of KDV. Phylogenetic analysis based on the full-genome sequence revealed that SZ_M48 is more clustered with KDV isolates from Hubei and Shangdong in China than with Indonesian and Danish strains. The identity between SZ_M48 and SDKL1625 (Shandong, China) is slightly lower than that of QTM27331 (Hubei, China), and the identity with JKT-7075 (Indonesia) and 21164-6/M.dau/DK (Denmark) is the lowest. Conclusion: The full-genome sequence of the new KDV strain described in this study may be useful for surveillance of the evolutionary characteristics of KDVs. Moreover, these findings extend the knowledge about the genomic diversity, potential vectors, and the distribution of KDVs in China.

Complete genome sequence of Alloalcanivorax xenomutans HF10, an alkane-degrading bacterium under hypersaline conditions

We report the complete genome sequence of Alloalcanivorax xenomutans HF10, an alkane-degrading strain isolated from the sediments of ocean in Xiamen, China, with a high salt tolerance potential of more than 10%. Its genome is composed of a 4.76-Mb chromosome.

High-throughput single-molecule long-read RNA sequencing analysis of tissue-specific genes and isoforms in lettuce (Lactuca sativa L.)

Lettuce is one of the most widely cultivated and consumed dicotyledonous vegetables globally. Despite the availability of its reference genome sequence, lettuce gene annotation remains incomplete, impeding comprehensive research and the broad application of genomic resources. Long-read RNA isoform sequencing (Iso-Seq) offers substantial advantages for analyzing RNA alternative splicing and aiding gene annotation, yet it faces throughput limitations. We present the HIT-ISOseq method tailored for bulk sample analysis, significantly enhancing RNA sequencing throughput on the PacBio platform by concatenating cDNA. Here we show, HIT-ISOseq generates 3-4 cDNA molecules per CCS read in lettuce, yielding 15.7 million long reads per PacBio Sequel II SMRT Cell 8 M. We validate its effectiveness in analyzing six lettuce tissue samples, including roots, stems, and leaves, revealing tissue-specific gene expression patterns and RNA isoforms. Leveraging diverse tissue long-read RNA sequencing, we refine the transcript annotation of the lettuce reference genome, expanding its GO and KEGG annotation repertoire. Collectively, this study serves as a foundational reference for genome annotation and the analysis of multi-sample isoform expression, utilizing high-throughput long-read transcriptome sequencing.

Genome-wide association studies of body size traits in Tibetan sheep

BACKGROUND

Elucidating the genetic variation underlying phenotypic diversity will facilitate improving production performance in livestock species. The Tibetan sheep breed in China holds significant historical importance, serving as a fundamental pillar of Qinghai's animal husbandry sector. The Plateau-type Tibetan sheep, comprising 90% of the province's population, are characterized by their tall stature and serve as the primary breed among Tibetan sheep. In contrast, Zhashijia sheep exhibit larger size and superior meat quality. These two species provide an excellent model for elucidating the genetic basis of body size variation. Therefore, this study aims to conduct a comprehensive genome-wide association study on these two Tibetan sheep breeds to identify single nucleotide polymorphism loci and regulatory genes that influence body size traits in Tibetan sheep.

RESULT

In this study, the phenotypic traits of body weight, body length, body height, chest circumference, chest depth, chest width, waist angle width, and pipe circumference were evaluated in two Tibetan sheep breeds: Plateau-type sheep and Zhashijia Tibetan sheep. Whole genome sequencing generated 48,215,130 high-quality SNPs for genome-wide association study. Four methods were applied and identified 623 SNPs significantly associated with body size traits. The significantly associated single nucleotide polymorphisms identified in this study are located near or within 111 candidate genes. These genes exhibit enrichment in the cAMP and Rap1 signaling pathways, significantly affecting animal growth, and body size. Specifically, the following genes were associated: ASAP1, CDK6, FRYL, NAV2, PTPRM, GPC6, PTPRG, KANK1, NTRK2 and ADCY8.

CONCLUSION

By genome-wide association study, we identified 16 SNPs and 10 candidate genes associated with body size traits in Tibetan sheep, which hold potential for application in genomic selection breeding programs in sheep. Identifying these candidate genes will establish a solid foundation for applying molecular marker-assisted selection in sheep breeding and improve our understanding of body size control in farmed animals.

Transoral Robotic Surgery and Human Papillomavirus Infection: Impact on Oropharyngeal Cancer Prognosis

Background/Objective: The incidence of oropharyngeal cancer (OPC) remains significant, with a rising prevalence of HPV-positive (HPV+) cases, underscoring the growing importance of appropriate treatment approaches for this condition. While HPV+ OPC typically exhibits a more favorable prognosis than HPV-negative (HPV-) OPC, certain HPV+ OPC patients still face adverse outcomes. This study aimed to assess the effectiveness of TORS versus traditional surgery in treating OPC patients and investigate the prognostic implications of specific variants in the HPV genome. Methods: The clinical information, including pathological features, treatments, and outcomes (death), of 135 OPC patients treated with traditional surgery from 2008 to 2018 (the non-TORS group) and 130 OPC patients treated with TORS from 2017 to 2021 (the TORS group) were obtained from Sun Yat-sen University Cancer Center (SYSUCC). A comparative analysis of 3-year overall survival (OS) was performed between these two groups. Furthermore, we conducted next-generation sequencing for the HPV16 genome of the 68 HPV+ OPC cases to characterize single-nucleotide variations (SNVs) in the HPV16 genome and evaluate its association with HPV+ OPC patient survival. Results: The comparative analysis of 3-year OS between the two groups (TORS vs. non-TORS) revealed a significant prognostic improvement in the TORS group for OPC patients with a T1-T2 classification (89.3% vs. 72.0%; p = 1.1 × 10-2), stages I-II (92.1% vs. 82.2%; p = 4.6 × 10-2), and stages III-IV (82.8% vs. 62.2%; p = 5.7 × 10-2) and for HPV- patients (85.5% vs. 33.3%; p < 1.0 × 10-6). Furthermore, three SNVs (SNV1339A>G, SNV1950A>C, and SNV4298A>G) in the HPV16 genome were identified as being associated with worse survival. These SNVs could alter protein interactions and weaken the binding affinity for MHC-II, promoting viral amplification and immune evasion. Conclusions: TORS exhibited a superior prognosis to traditional surgery in OPC patients. Additionally, identifying specific SNVs within the HPV16 genome provided potential prognostic markers for HPV+ OPC. These significant findings hold clinical relevance for treatment decision-making and prognostic assessment in patients with OPC.

Water wisteria genome reveals environmental adaptation and heterophylly regulation in amphibious plants

Heterophylly is a phenomenon whereby an individual plant dramatically changes leaf shape in response to the surroundings. Hygrophila difformis (Acanthaceae; water wisteria), has recently emerged as a model plant to study heterophylly because of its striking leaf shape variation in response to various environmental factors. When submerged, H. difformis often develops complex leaves, but on land it develops simple leaves. Leaf complexity is also influenced by other factors, such as light density, humidity, and temperature. Here, we sequenced and assembled the H. difformis chromosome-level genome (scaffold N50: 60.43 Mb, genome size: 871.92 Mb), which revealed 36 099 predicted protein-coding genes distributed over 15 pseudochromosomes. H. difformis diverged from its relatives during the Oligocene climate-change period and expanded gene families related to its amphibious habit. Genes related to environmental stimuli, leaf development, and other pathways were differentially expressed in submerged and terrestrial conditions, possibly modulating morphological and physiological acclimation to changing environments. We also found that auxin plays a role in H. difformis heterophylly. Finally, we discovered candidate genes that respond to different environmental conditions and elucidated the role of LATE MERISTEM IDENTITY 1 (LMI1) in heterophylly. We established H. difformis as a model for studying interconnections between environmental adaptation and morphogenesis.

Whole-Genome Metagenomic Analysis of the Oral Microbiota in Patients with Obstructive Sleep Apnea Comorbid with Major Depressive Disorder

Background

Obstructive sleep apnea (OSA) patients commonly experience high rates of depression. This study aims to examine the oral microbiota characteristics of OSA and those with comorbid major depressive disorder (OSA+MDD) patients.

Methods

Participants were enrolled from Aug 2022 to Apr 2023. Polysomnography, psychiatrist interviews, and scales were used to diagnose OSA and MDD. Oral samples were collected from participants by rubbing swabs on buccal mucosa, palate, and gums. Oral microbiota was analyzed via whole-genome metagenomics and bioinformatic analysis followed sequencing. Venous blood was drawn to detect plasma inflammatory factor levels.

Results

The study enrolled 33 OSA patients, 28 OSA+MDD patients, and 28 healthy controls. Significant differences were found in 8 phyla, 229 genera, and 700 species of oral microbiota among the three groups. Prevotellaceae abundance in the OSA and OSA+MDD groups was significantly lower than that in healthy controls. Linear discriminant analysis effect size (LEfSe) analysis showed that Streptococcaceae and Actinobacteria were the characteristic oral microbiota of the OSA and OSA+MDD groups, respectively. KEGG analysis indicates 30 pathways were changed in the OSA and OSA+MDD groups compared with healthy controls, and 23 pathways were changed in the OSA group compared with the OSA+MDD group. Levels of IL-6 in the OSA+MDD group were significantly higher than in the healthy group, correlating positively with the abundance of Schaalia, Campylobacter, Fusobacterium, Alloprevotella, and Candidatus Nanosynbacter in the oral, as well as with Hamilton Anxiety Rating Scale and Hamilton Depression Rating Scale scores.

Conclusion

Significant differences in oral microbiota populations and gene function were observed among the three groups. OSA patients were characterized by a decreased abundance of Prevotellaceae and an increased abundance of Streptococcaceae. OSA+MDD patients had an increased abundance of Actinobacteria. IL-6 might regulate the relationship between depression and the oral microbiota in OSA+MDD patients.

Mar1, a high mobility group box protein, regulates n-alkane adsorption and cell morphology of the dimorphic yeast Yarrowia lipolytica

The dimorphic yeast Yarrowia lipolytica possesses an excellent ability to utilize n-alkane as a sole carbon and energy source. Although there are detailed studies on the enzymes that catalyze the reactions in the metabolic processes of n-alkane in Y. lipolytica, the molecular mechanism underlying the incorporation of n-alkane into the cells remains to be elucidated. Because Y. lipolytica adsorbs n-alkane, we postulated that Y. lipolytica incorporates n-alkane through direct interaction with it. We isolated and characterized mutants defective in adsorption to n-hexadecane. One of the mutants harbored a nonsense mutation in MAR1 (Morphology and n-alkane Adsorption Regulator 1) encoding a protein containing a high mobility group box. The deletion mutant of MAR1 exhibited defects in adsorption to n-hexadecane and filamentous growth on solid media, whereas the strain that overexpressed MAR1 exhibited hyperfilamentous growth. Fluorescence microscopic observations suggested that Mar1 localizes in the nucleus. RNA-sequencing analysis revealed the alteration of the transcript levels of several genes, including those encoding transcription factors and cell surface proteins, by the deletion of MAR1. These findings suggest that MAR1 is involved in the transcriptional regulation of the genes required for n-alkane adsorption and cell morphology transition.IMPORTANCEYarrowia lipolytica, a dimorphic yeast capable of assimilating n-alkane as a carbon and energy source, has been extensively studied as a promising host for bioconversion of n-alkane into useful chemicals and bioremediation of soil and water contaminated by petroleum. While the metabolic pathway of n-alkane in this yeast and the enzymes involved in this pathway have been well characterized, the molecular mechanism to incorporate n-alkane into the cells is yet to be fully understood. Due to the ability of Y. lipolytica to adsorb n-alkane, it has been hypothesized that Y. lipolytica incorporates n-alkane through direct interaction with it. In this study, we identified a gene, MAR1, which plays a crucial role in the transcriptional regulation of the genes necessary for the adsorption to n-alkane and the transition of the cell morphology in Y. lipolytica. Our findings provide valuable insights that could lead to advanced applications of Y. lipolytica in n-alkane bioconversion and bioremediation.

Prioritizing Maize Metabolic Gene Regulators through Multi-Omic Network Integration

Elucidating gene regulatory networks is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. Here, we integrated 46 co-expression networks, 283 protein-DNA interaction (PDI) assays, and 16 million SNPs used to identify expression quantitative trait loci (eQTL) to construct TF-target networks. In total, we analyzed ∼4.6M interactions to generate four distinct types of TF-target networks: co-expression, PDI, trans -eQTL, and cis -eQTL combined with PDIs. To functionally annotate TFs based on their target genes, we implemented three different network integration strategies. We evaluated the effectiveness of each strategy through TF loss-of function mutant inspection and random network analyses. The multi-network integration allowed us to identify transcriptional regulators of several biological processes. Using the topological properties of the fully integrated network, we identified potential functionally redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

GRAPHICAL ABSTRACT

Diversity and phage sensitivity to phages of porcine enterotoxigenic Escherichia coli

Enterotoxigenic Escherichia coli (ETEC) is a diverse and poorly characterized E. coli pathotype that causes diarrhea in humans and animals. Phages have been proposed for the veterinary biocontrol of ETEC, but effective solutions require understanding of porcine ETEC diversity that affects phage infection. Here, we sequenced and analyzed the genomes of the PHAGEBio ETEC collection, gathering 79 diverse ETEC strains isolated from European pigs with post-weaning diarrhea (PWD). We identified the virulence factors characterizing the pathotype and several antibiotic resistance genes on plasmids, while phage resistance genes and other virulence factors were mostly chromosome encoded. We experienced that ETEC strains were highly resistant to Enterobacteriaceae phage infection. It was only by enrichment of numerous diverse samples with different media and conditions, using the 41 ETEC strains of our collection as hosts, that we could isolate two lytic phages that could infect a large part of our diverse ETEC collection: vB_EcoP_ETEP21B and vB_EcoS_ETEP102. Based on genome and host range analyses, we discussed the infection strategies of the two phages and identified components of lipopolysaccharides ( LPS) as receptors for the two phages. Our detailed computational structural analysis highlights several loops and pockets in the tail fibers that may allow recognition and binding of ETEC strains, also in the presence of O-antigens. Despite the importance of receptor recognition, the diversity of the ETEC strains remains a significant challenge for isolating ETEC phages and developing sustainable phage-based products to address ETEC-induced PWD.IMPORTANCEEnterotoxigenic Escherichia coli (ETEC)-induced post-weaning diarrhea is a severe disease in piglets that leads to weight loss and potentially death, with high economic and animal welfare costs worldwide. Phage-based approaches have been proposed, but available data are insufficient to ensure efficacy. Genome analysis of an extensive collection of ETEC strains revealed that phage defense mechanisms were mostly chromosome encoded, suggesting a lower chance of spread and selection by phage exposure. The difficulty in isolating lytic phages and the molecular and structural analyses of two ETEC phages point toward a multifactorial resistance of ETEC to phage infection and the importance of extensive phage screenings specifically against clinically relevant strains. The PHAGEBio ETEC collection and these two phages are valuable tools for the scientific community to expand our knowledge on the most studied, but still enigmatic, bacterial species-E. coli.

Multi-omics analysis unveils immunosuppressive microenvironment in the occurrence and development of multiple pulmonary lung cancers

Multiple pulmonary lung cancers (MPLCs) are frequently encountered on computed tomography (CT) scanning of chest, yet their intrinsic characteristics associated with genomic features and radiological or pathological textures that may lead to distinct clinical outcomes remain largely unexplored. A total of 27 pulmonary nodules covering different radiological or pathological textures as well as matched adjacent normal tissues and blood samples were collected from patients diagnosed with MPLCs. Whole-exome sequencing (WES) and whole-transcriptome sequencing were performed. The molecular and immune features of MPLCs associated with distinct radiological or pathological textures were comprehensively investigated. Genomics analysis unveiled the distinct branches of pulmonary nodules originating independently within the same individual. EGFR and KRAS mutations were found to be prevalent in MPLCs, exhibiting mutual exclusivity. The group with KRAS mutations exhibited stronger immune signatures compared to the group with EGFR mutations. Additionally, MPLCs exhibited a pronounced immunosuppressive microenvironment, which was particularly distinct when compared with normal tissues. The expression of the FDSCP gene was specifically observed in MPLCs. When categorizing MPLCs based on radiological or pathological characteristics, a progressive increase in mutation accumulation was observed, accompanied by heightened chromatin-level instability as ground-glass opacity component declined or invasive progression occurred. A close association with the immunosuppressive microenvironment was also observed during the progression of pulmonary nodules. Notably, the upregulation of B cell and regulatory T cell marker genes occurred progressively. Immune cell abundance analysis further demonstrated a marked increase in exhausted cells and regulatory T cells during the progression of pulmonary nodules. These results were further validated by independent datasets including nCounter RNA profiling, single-cell RNA sequencing, and spatial transcriptomic datasets. Our study provided a comprehensive representation of the diverse landscape of MPLCs originating within the same individual and emphasized the significant influence of the immunosuppressive microenvironment in the occurrence and development of pulmonary nodules. These findings hold great potential for enhancing the clinical diagnosis and treatment strategies for MPLCs.

Integration of metagenomics and metaproteomics in the intestinal lavage fluids benefits construction of discriminative model and discovery of biomarkers for HBV liver diseases

Intestinal lavage fluid (IVF) containing the mucosa-associated microbiota instead of fecal samples was used to study the gut microbiota using different omics approaches. Focusing on the 63 IVF samples collected from healthy and hepatitis B virus-liver disease (HBV-LD), a question is prompted whether omics features could be extracted to distinguish these samples. The IVF-related microbiota derived from the omics data was classified into two enterotype sets, whereas the genomics-based enterotypes were poorly overlapped with the proteomics-based one in either distribution of microbiota or of IVFs. There is lack of molecular features in these enterotypes to specifically recognize healthy or HBV-LD. Running machine learning against the omics data sought the appropriate models to discriminate the healthy and HBV-LD IVFs based on selected genes or proteins. Although a single omics dataset is basically workable in such discrimination, integration of the two datasets enhances discrimination efficiency. The protein features with higher frequencies in the models are further compared between healthy and HBV-LD based on their abundance, bringing about three potential protein biomarkers. This study highlights that integration of metaomics data is beneficial for a molecular discriminator of healthy and HBV-LD, and reveals the IVF samples are valuable for microbiome in a small cohort.

Chromosome-level genome assembly and annotation of xerophyte secretohalophyte Reaumuria soongarica

Reaumuria soongarica is a xerophytic shrub belonging to the Tamaricaceae family. The species is widely distributed in the deserts of Central Asia and is characterized by its remarkable adaptability to saline and barren desert environments. Using PacBio long-read sequencing and Hi-C technologies, we assembled a chromosome-level genome of R. soongarica. The genome assembly has a size of 1.28 Gb with a scaffold N50 of 116.15 Mb, and approximately 1.25 Gb sequences were anchored in 11 pseudo-chromosomes. A completeness assessment of the assembled genome revealed a BUSCO score of 97.5% and an LTR Assembly Index of 12.37. R. soongarica genome had approximately 60.07% repeat sequences. In total, 21,791 protein-coding genes were predicted, of which 95.64% were functionally annotated. This high-quality genome will serve as a foundation for studying the genomic evolution and adaptive mechanisms to arid-saline environments in R. soongarica, facilitating the exploration and utilization of its unique genetic resources.

Transcriptome Analysis Reveals Genes and Pathways Associated with Drought Tolerance of Early Stages in Sweet Potato (Ipomoea batatas (L.) Lam.)

The yield of sweet potato [Ipomoea batatas (L.) Lam] can be easily threatened by drought stress. Typically, early stages like the seedling stage and tuber-root expansion stage are more vulnerable to drought stress. In this study, a highly drought-tolerant sweet potato cultivar "WanSu 63" was subjected to drought stress at both the seedling stage (15 days after transplanting, 15 DAT) and the tuber-root expansion stage (45 DAT). Twenty-four cDNA libraries were constructed from leaf segments and root tissues at 15 and 45 DAT for Next-Generation Sequencing. A total of 663, 063, and 218 clean reads were obtained and then aligned to the reference genome with a total mapped ratio greater than 82.73%. A sum of 7119, 8811, 5463, and 930 differentially expressed genes were identified from leaves in 15 days (L15), roots in 15 days (R15), leaves in 45 days (L45), and roots in 45 days (R45), respectively, in drought stress versus control. It was found that genes encoding heat shock proteins, sporamin, LEA protein dehydrin, ABA signaling pathway protein gene NCED1, as well as a group of receptor-like protein kinases genes were enriched in differentially expressed genes. ABA content was significantly higher in drought-treated tissues than in the control. The sweet potato biomass declined sharply to nearly one-quarter after drought stress. In conclusion, this study is the first to identify the differentially expressed drought-responsive genes and signaling pathways in the leaves and roots of sweet potato at the seedling and root expansion stages. The results provide potential resources for drought resistance breeding of sweet potato.

Complete genome sequence of Rhodococcus sp. strain PD04, a pyridine-degrading bacterium under hypersaline condition

We report the complete genome sequence of a pyridine-degrading Rhodococcus sp. strain PD04 under 4% salinity environment, isolated from wastewater of coking plant. The genome is 6.07 Mb with 5,767 annotated gene coding sequences.

Kinship analysis and pedigree reconstruction by RAD sequencing in cattle

Kinship and pedigree, used for estimating inbreeding, heritability, selection, and gene flow, are useful for breeding and animal conservation. However, as the size of crossbred populations increases, inaccurate generation and parentage assignment in livestock farms increase. Restriction-site-associated DNA sequencing is a cost-effective platform for single nucleotide polymorphism (SNP) discovery and genotyping. Here, we performed a kinship analysis and pedigree reconstruction for Angus and Xiangxi yellow cattle. A total of 975 cattle, including 923 offspring with 24 known sires and 28 known dams, were sampled and subjected to SNP discovery and genotyping. The identified SNP panel included 7,305 SNPs capturing the maximum difference between paternal and maternal genome information, allowing us to distinguish F1 from F2 generations with 90% accuracy. In conclusion, we provided a low-cost and efficient SNP panel for kinship analyses and the improvement of local genetic resources, which are valuable for breed improvement, local resource utilization, and conservation.

Complete genome sequence of Pseudomonas aeruginosa L3, a bacterium with Mn(II) oxidative property isolated from heavy metal-contaminated soils

Pseudomonas aeruginosa L3, isolated from heavy metal-contaminated soils, possesses the ability of Mn(II) oxidation. To further enhance the understanding of genes involved in Mn(II) oxidation, the complete genome of this strain was sequenced and annotated, which has a total size of 6.39 Mb with a G + C content of 66.39%.

A chromosome-level genome assembly and evolutionary analysis of Coregonus ussuriensis Berg

Coregonus ussuriensis Berg, distributed widely in cold waters above 45° N latitude, is a savored freshwater whitefish that has been included in the list of endangered animals as a consequence of overfishing. Lack of genomic information seriously hampers evolutionary and genetic research on C. ussuriensis warranting the need to assemble a high-quality reference genome to promote its genetic breeding. We assembled and constructed a reference chromosome-level C. ussuriensis genome (sequence length, 2.51 Gb; contig N50 length, 4.27 Mb) using PacBio sequencing and Hi-C assembly technology, 3,109 contigs were assembled into scaffolds, resulting in a genome assembly with 40 chromosomes and a scaffold N50 length of 62.20 Mb. In addition, 43,320 protein-coding genes were annotated. The peak Ks position in the species comparison reflects the whole-genome replication event of C. ussuriensis. This chromosome-level genome provides reference data for further studies on the molecular breeding of C. ussuriensis.

Genome sequences of two Bacillus anthracis strains utilized as veterinary vaccines in China

In this report, we present the complete genome sequences of two Bacillus anthracis strains utilized as veterinary vaccines in China. The sequencing was conducted using a hybrid assembly methodology that combined Illumina short reads and PacBio long reads. This approach provides a high-quality representative sequence for the strains mentioned above.

Chromosome-level genome assembly of Aquilaria yunnanensis

Aquilaria yunnanensis is an endangered agarwood-producing tree currently listed on the IUCN Red List of Threatened Species. The agarwood it produces has important medicinal and economic value, but its population has sharply declined due to human destruction and habitat reduction. Therefore, obtaining genomic information on A. yunnanensis is beneficial for its protection work. We assembled a chromosome-level reference genome of A. yunnanensis by using BGI short reads, PacBio HiFi long reads, coupled with Hi-C technology. The final genome assembly of A. yunnanensis is 847.04 Mb, with N50 size of 99.68 Mb, in which 805.49 Mb of the bases were anchored on eight pseudo-chromosomes. Two gapless pseudo-chromosomes were detected in the assembly. A total of 27,955 protein-coding genes as well as 74.65% repetitive elements were annotated. These findings may provide valuable resources in conservation, functional genomics, and molecular breeding of A. yunnanensis, as well as the molecular phylogenetics and evolutionary patterns in Aquilaria.

Prolonged Sleep Deprivation Induces a Reprogramming of Circadian Rhythmicity with the Hepatic Metabolic Transcriptomic Profile

Sleep disturbances can disrupt the overall circadian rhythm. However, the impact of sleep deprivation on the circadian rhythm of the liver and its underlying mechanisms still requires further exploration. In this study, we subjected male mice to 5 days of sleep deprivation and performed liver transcriptome sequencing analysis at various time points within a 24-h period. Subsequently, we monitored the autonomic activity and food intake in these male mice for six days post-sleep deprivation. We observed alterations in sleep-wake and feeding rhythms in the first two days following sleep deprivation. Additionally, we also observed a decrease in 24-h serum-glucose levels. Liver transcriptome sequencing has shown that sleep deprivation induces the rhythmic transcription of a large number of genes, or alters the rhythmic properties of genes, which were then significantly enriched in the carbohydrate, lipid, and protein metabolism pathways. Our findings suggest that under conditions of prolonged sleep deprivation, the expression of metabolic-related genes in the liver was reset, leading to changes in the organism's metabolic state to ensure energy supply to sustain prolonged wakefulness.

Transcriptomic changes in the posterior pallium of male zebra finches associated with social niche conformance

Animals plastically adjust their physiological and behavioural phenotypes to conform to their social environment-social niche conformance. The degree of sexual competition is a critical part of the social environment to which animals adjust their phenotypes, but the underlying genetic mechanisms are poorly understood. We conducted a study to investigate how differences in sperm competition risk affect the gene expression profiles of the testes and two brain areas (posterior pallium and optic tectum) in breeding male zebra finches (Taeniopygia castanotis). In this pre-registered study, we investigated a large sample of 59 individual transcriptomes. We compared two experimental groups: males held in single breeding pairs (low sexual competition) versus those held in two pairs (elevated sexual competition) per breeding cage. Using weighted gene co-expression network analysis (WGCNA), we observed significant effects of the social treatment in all three tissues. However, only the treatment effects found in the pallium were confirmed by an additional randomisation test for statistical robustness. Likewise, the differential gene expression analysis revealed treatment effects only in the posterior pallium (ten genes) and optic tectum (six genes). No treatment effects were found in the testis at the single gene level. Thus, our experiments do not provide strong evidence for transcriptomic adjustment specific to manipulated sperm competition risk. However, we did observe transcriptomic adjustments to the manipulated social environment in the posterior pallium. These effects were polygenic rather than based on few individual genes with strong effects. Our findings are discussed in relation to an accompanying paper using the same animals, which reports behavioural results consistent with the results presented here.

Safety of embryo cryopreservation: insights from mid-term placental transcriptional changes

BACKGROUND

In recent years, with benefits from the continuous improvement of clinical technology and the advantage of fertility preservation, the application of embryo cryopreservation has been growing rapidly worldwide. However, amidst this growth, concerns about its safety persist. Numerous studies have highlighted the elevated risk of perinatal complications linked to frozen embryo transfer (FET), such as large for gestational age (LGA) and hypertensive disorders during pregnancy. Thus, it is imperative to explore the potential risk of embryo cryopreservation and its related mechanisms.

METHODS

Given the strict ethical constraints on clinical samples, we employed mouse models in this study. Three experimental groups were established: the naturally conceived (NC) group, the fresh embryo transfer (Fresh-ET) group, and the FET group. Blastocyst formation rates and implantation rates were calculated post-embryo cryopreservation. The impact of FET on fetal growth was evaluated upon fetal and placental weight. Placental RNA-seq was conducted, encompassing comprehensive analyses of various comparisons (Fresh-ET vs. NC, FET vs. NC, and FET vs. Fresh-ET).

RESULTS

Reduced rates of blastocyst formation and implantation were observed post-embryo cryopreservation. Fresh-ET resulted in a significant decrease in fetal weight compared to NC group, whereas FET reversed this decline. RNA-seq analysis indicated that the majority of the expression changes in FET were inherited from Fresh-ET, and alterations solely attributed to embryo cryopreservation were moderate. Unexpectedly, certain genes that showed alterations in Fresh-ET tended to be restored in FET. Further analysis suggested that this regression may underlie the improvement of fetal growth restriction in FET. The expression of imprinted genes was disrupted in both FET and Fresh-ET groups.

CONCLUSION

Based on our experimental data on mouse models, the impact of embryo cryopreservation is less pronounced than other in vitro manipulations in Fresh-ET. However, the impairment of the embryonic developmental potential and the gene alterations in placenta still suggested it to be a risky operation.

Impact of Methyl Jasmonate on Terpenoid Biosynthesis and Functional Analysis of Sesquiterpene Synthesis Genes in Schizonepeta tenuifolia

This study investigates the impact of methyl jasmonate (MeJA) on the volatile oil composition of Schizonepeta tenuifolia and elucidates the function of the StTPS45 gene, a key player in terpenoid biosynthesis. The effect of different concentrations of MeJA (0, 50, 100, 200, and 300 μmol/L) on the growth of S. tenuifolia adventitious bud clusters was analyzed over a 20 d period. Using gas chromatography-mass spectrometry (GC-MS), 17 compounds were identified from the adventitious bud clusters of S. tenuifolia. Significant changes in the levels of major monoterpenes, including increased contents of (+)-limonene and (+)-menthone, were observed, particularly at higher concentrations of MeJA. Analysis of transcriptome data from three groups treated with 0, 100, and 300 μmol/L MeJA revealed significant changes in the gene expression profiles following MeJA treatment. At 100 μmol/L MeJA, most terpene synthase (TPS) genes were overexpressed. Additionally, gene expression and functional predictions suggested that StTPS45 acts as germacrene D synthase. Therefore, StTPS45 was cloned and expressed in Escherichia coli, and enzyme activity assays confirmed its function as a germacrene D synthase. Molecular docking and structural prediction of StTPS45 further suggested specific interactions with farnesyl diphosphate (FPP), aligning with its role in the terpenoid synthesis pathway. These findings provide valuable insights into the modulation of secondary metabolite pathways by jasmonate signaling and underscore the potential of genetic engineering approaches to enhance the production of specific terpenoids in medicinal plants.

A near complete genome assembly of the East Friesian sheep genome

Advancements in sequencing have enabled the assembly of numerous sheep genomes, significantly advancing our understanding of the link between genetic variation and phenotypic traits. However, the genome of East Friesian sheep (Ostfriesisches Milchschaf), a key high-yield milk breed, remains to be fully assembled. Here, we constructed a near-complete and gap-free East Friesian genome assembly using PacBio HiFi, ultra-long ONT and Hi-C sequencing. The resulting genome assembly spans approximately 2.96 Gb, with a contig N50 length of 104.1 Mb and only 164 unplaced sequences. Remarkably, our assembly has captured 41 telomeres and 24 centromeres. The assembled sequence is of high quality on completeness (BUSCO score: 97.1%) and correctness (QV: 69.1). In addition, a total of 24,580 protein-coding genes were predicted, of which 97.2% (23,891) carried at least one conserved functional domain. Collectively, this assembly provides not only a near T2T gap-free genome, but also provides a valuable genetic resource for comparative genome studies of sheep and will serve as an important tool for the sheep research community.

Transcriptome and metabolome analyses reveal regulatory networks associated with nutrition synthesis in sorghum seeds

Cereal seeds are vital for food, feed, and agricultural sustainability because they store and provide essential nutrients to human and animal food and feed systems. Unraveling molecular processes in seed development is crucial for enhancing cereal grain yield and quality. We analyze spatiotemporal transcriptome and metabolome profiles during sorghum seed development in the inbred line 'BTx623'. Morphological and molecular analyses identify the key stages of seed maturation, specifying starch biosynthesis onset at 5 days post-anthesis (dpa) and protein at 10 dpa. Transcriptome profiling from 1 to 25 dpa reveal dynamic gene expression pathways, shifting from cellular growth and embryo development (1-5 dpa) to cell division, fatty acid biosynthesis (5-25 dpa), and seed storage compounds synthesis in the endosperm (5-25 dpa). Network analysis identifies 361 and 207 hub genes linked to starch and protein synthesis in the endosperm, respectively, which will help breeders enhance sorghum grain quality. The availability of this data in the sorghum reference genome line establishes a baseline for future studies as new pangenomes emerge, which will consider copy number and presence-absence variation in functional food traits.

Whole-genome analysis reveals distinct adaptation signatures to diverse environments in Chinese domestic pigs

BACKGROUND

Long-term natural and artificial selection has resulted in many genetic footprints within the genomes of pig breeds across distinct agroecological zones. Nevertheless, the mechanisms by which these signatures contribute to phenotypic diversity and facilitate environmental adaptation remain unclear.

RESULTS

Here, we leveraged whole-genome sequencing data from 82 individuals from 6 domestic pig breeds originating in tropical, high-altitude, and frigid regions. Population genetic analysis suggested that habitat isolation significantly shaped the genetic diversity and contributed to population stratification in local Chinese pig breeds. Analysis of selection signals revealed regions under selection for adaptation in tropical (55.5 Mb), high-altitude (43.6 Mb), and frigid (17.72 Mb) regions. The potential functions of the selective sweep regions were linked to certain complex traits that might play critical roles in different geographic environments, including fat coverage in frigid environments and blood indicators in tropical and high-altitude environments. Candidate genes under selection were significantly enriched in biological pathways involved in environmental adaptation. These pathways included blood circulation, protein degradation, and inflammation for adaptation to tropical environments; heart and lung development, hypoxia response, and DNA damage repair for high-altitude adaptation; and thermogenesis, cold-induced vasodilation (CIVD), and the cell cycle for adaptation to frigid environments. By examining the chromatin state of the selection signatures, we identified the lung and ileum as two candidate functional tissues for environmental adaptation. Finally, we identified a mutation (chr1: G246,175,129A) in the cis-regulatory region of ABCA1 as a plausible promising variant for adaptation to tropical environments.

CONCLUSIONS

In this study, we conducted a genome-wide exploration of the genetic mechanisms underlying the adaptability of local Chinese pig breeds to tropical, high-altitude, and frigid environments. Our findings shed light on the prominent role of cis-regulatory elements in environmental adaptation in pigs and may serve as a valuable biological model of human plateau-related disorders and cardiovascular diseases.