InterPro and InterProScan: tools for protein sequence classification and comparison

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.

Chromosome-level genome of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus

Recent conservation efforts to protect rare and endangered aquatic species have intensified. Nevertheless, the ornate spiny lobster (Panulirus ornatus), which is prevalent in the Indo-Pacific waters, has been largely ignored. In the absence of a detailed genomic reference, the conservation and population genetics of this crustacean are poorly understood. Here, We assembled a comprehensive chromosome-level genome for P. ornatus. This genome-among the most detailed for lobsters-spans 2.65 Gb with a contig N50 of 51.05 Mb, and 99.11% of the sequences with incorporated to 73 chromosomes. The ornate spiny lobster genome comprises 65.67% repeat sequences and 22,752 protein-coding genes with 99.20% of the genes functionally annotated. The assembly of the P. ornatus genome provides valuable insights into comparative crustacean genomics and endangered species conservation, and lays the groundwork for future research on the speciation, ecology, and evolution of the ornate spiny lobster.

Chromosome-scale genome assemblies of sexually dimorphic male and female Acrossocheilus fasciatus

Acrossocheilus fasciatus is a stream-dwelling fish species of the Barbinae subfamily. It is valued for its colorfully striped appearance and delicious meat. This species is also characterized by apparent sexual dimorphism and toxic ovum. Biology and aquaculture researches of A. fasciatus are hindered by the lack of a high-quality reference genome. Here, we report chromosome-level genome assemblies of the male and female A. fasciatus. The HiFi-only genome assemblies for both female and male individuals were 899.13 Mb (N50 length of 32.58 Mb) and 885.68 Mb (N50 length of 33.06 Mb), respectively. Notably, a substantial proportion of the assembled sequences, accounting for 96.15% and 98.35% for female and male genomes, respectively, were successfully anchored onto 25 chromosomes utilizing Hi-C data. We annotated the female assembly as a reference genome and identified a total of 400.62 Mb (44.56%) repetitive sequences, 27,392 protein-coding genes, and 35,869 ncRNAs. The high-quality male and female reference genomes will provide genomic resources for developing sex-specific molecular markers, inform single-sex breeding, and elucidate genetic mechanisms of sexual dimorphism.

Single-cell long-read sequencing-based mapping reveals specialized splicing patterns in developing and adult mouse and human brain

RNA isoforms influence cell identity and function. However, a comprehensive brain isoform map was lacking. We analyze single-cell RNA isoforms across brain regions, cell subtypes, developmental time points and species. For 72% of genes, full-length isoform expression varies along one or more axes. Splicing, transcription start and polyadenylation sites vary strongly between cell types, influence protein architecture and associate with disease-linked variation. Additionally, neurotransmitter transport and synapse turnover genes harbor cell-type variability across anatomical regions. Regulation of cell-type-specific splicing is pronounced in the postnatal day 21-to-postnatal day 28 adolescent transition. Developmental isoform regulation is stronger than regional regulation for the same cell type. Cell-type-specific isoform regulation in mice is mostly maintained in the human hippocampus, allowing extrapolation to the human brain. Conversely, the human brain harbors additional cell-type specificity, suggesting gain-of-function isoforms. Together, this detailed single-cell atlas of full-length isoform regulation across development, anatomical regions and species reveals an unappreciated degree of isoform variability across multiple axes.

A new chromosome-scale genome of wild Brassica oleracea provides insights into the domestication of Brassica crops

The cultivated diploid Brassica oleracea is an important vegetable crop, but the genetic basis of its domestication remains largely unclear in the absence of high-quality reference genomes of wild B. oleracea. Here, we report the first chromosome-level assembly of the wild Brassica oleracea L. W03 genome (total genome size, 630.7 Mb; scaffold N50, 64.6 Mb). Using the newly assembled W03 genome, we constructed a gene-based B. oleracea pangenome and identified 29 744 core genes, 23 306 dispensable genes, and 1896 private genes. We re-sequenced 53 accessions, representing six potential wild B. oleracea progenitor species. The results of the population genomic analysis showed that the wild B. oleracea populations had the highest level of diversity and represents the most closely related population to modern-day horticultural B. oleracea. In addition, the WUSCHEL gene was found to play a decisive role in domestication and to be involved in cauliflower and broccoli curd formation. We also illustrate the loss of disease-resistance genes during selection for domestication. Our results provide new insights into the domestication of B. oleracea and will facilitate the future genetic improvement of Brassica crops.

Graph-Based Pangenome of Actinidia chinensis Reveals Structural Variations Mediating Fruit Degreening

Fruit ripening is associated with the degreening process (loss of chlorophyll) that occurs in most fruit species. Kiwifruit is one of the special species whose fruits may maintain green flesh by accumulating a large amount of chlorophyll even after ripening. However, little is known about the genetic variations related to the fruit degreening process. Here, a graph-based kiwifruit pangenome by analyzing 14 chromosome-scale haplotype-resolved genome assemblies from seven representative cultivars or lines in Actinidia chinensis is built. A total of 49,770 non-redundant gene families are identified, with core genes constituting 46.6%, and dispensable genes constituting 53.4%. A total of 84,591 non-redundant structural variations (SVs) are identified. The pangenome graph integrating both reference genome sequences and variant information facilitates the identification of SVs related to fruit color. The SV in the promoter of the AcBCM gene determines its high expression in the late developmental stage of fruits, which causes chlorophyll accumulation in the green-flesh fruits by post-translationally regulating AcSGR2, a key enzyme of chlorophyll catabolism. Taken together, a high-quality pangenome is constructed, unraveled numerous genetic variations, and identified a novel SV mediating fruit coloration and fruit quality, providing valuable information for further investigating genome evolution and domestication, QTL genes function, and genomics-assisted breeding.

The first chromosome-level genome of the stag beetle Dorcus hopei Saunders, 1854 (Coleoptera: Lucanidae)

Stag beetles (Coleoptera: Lucanidae) represent a significant saproxylic assemblage in forest ecosystems and are noted for their enlarged mandibles and male polymorphism. Despite their relevance as ideal models for the study of exaggerated mandibles that aid in attracting mates, the regulatory mechanisms associated with these traits remain understudied, and restricted by the lack of high-quality reference genomes for stag beetles. To address this limitation, we successfully assembled the first chromosome-level genome of a representative species Dorcus hopei. The genome was 496.58 Mb in length, with a scaffold N50 size of 54.61 Mb, BUSCO values of 99.8%, and 96.8% of scaffolds anchored to nine pairs of chromosomes. We identified 285.27 Mb (57.45%) of repeat sequences and annotated 11,231 protein-coding genes. This genome will be a valuable resource for further understanding the evolution and ecology of stag beetles, and provides a basis for studying the mechanisms of exaggerated mandibles through comparative analysis.

Temperature-dependent action of pepper mildew resistance locus O 1 in inducing pathogen immunity and thermotolerance

Plant diseases tend to be more serious under conditions of high-temperature/high-humidity (HTHH) than under moderate conditions, and hence disease resistance under HTHH is an important determinant for plant survival. However, how plants cope with diseases under HTHH remains poorly understood. In this study, we used the pathosystem consisting of pepper (Capsicum annuum) and Ralstonia solanacearum (bacterial wilt) as a model to examine the functions of the protein mildew resistance locus O 1 (CaMLO1) and U-box domain-containing protein 21 (CaPUB21) under conditions of 80% humidity and either 28 °C or 37 °C. Expression profiling, loss- and gain-of-function assays involving virus-induced gene-silencing and overexpression in pepper plants, and protein-protein interaction assays were conducted, and the results showed that CaMLO1 acted negatively in pepper immunity against R. solanacearum at 28 °C but positively at 37 °C. In contrast, CaPUB21 acted positively in immunity at 28 °C but negatively at 37 °C. Importantly, CaPUB21 interacted with CaMLO1 under all of the tested conditions, but only the interaction in response to R. solanacearum at 37 °C or to exposure to 37 °C alone led to CaMLO1 degradation, thereby turning off defence responses against R. solanacearum at 37 °C and under high-temperature stress to conserve resources. Thus, we show that CaMLO1 and CaPUB21 interact with each other and function distinctly in pepper immunity against R. solanacearum in an environment-dependent manner.

Pharmacoinformatics approach for type 2 diabetes mellitus therapeutics using phytocompounds from Costus genus: an in-silico investigation

Type 2 Diabetes Mellitus (T2DM), as a significant health concern globally, particularly in India, underscoring the vital need for effective therapeutics. Current drug therapies for T2DM may have limitations, leading researchers to explore natural products as alternatives. In this study. We have investigated the anti-diabetic compounds from the Costus genus, known as the insulin plant, which is abundant in southern India. The bioinformatics tools and software used for in-silico analysis to identify potential therapeutic compounds and hub genes associated with T2DM in the Indian population that could cut short the in-vitro and in-vivo experimental approaches in near future. The systematic review and combinatorial in-silico analysis revealed IGF2BP2, INS and TCF as the key targets that are associated with T2DM. The compounds stigmasterol, cycloartenol, and diosgenone were explored to be potent among all the 38 phytocompounds from genus Costus with binding energies -8.48, -10.07, and -10.31 kcal/mol against IGF2BP2, INS and TCF. The molecular dynamics (MD) simulation studies of these complexes demonstrated stable and consistent dynamic behavior, particularly in the INS-cycloartenol, IGF2BP2-stigmasterol and TCF7L2-diosgenone complexes. The identified compounds and associated targets represent potential candidates for T2DM therapeutics in the Indian population. The pharmacoinformatics approach presented in the study could streamline the drug discovery process by prioritizing compounds for further experimental validation.Communicated by Ramaswamy H. Sarma.

Degradation of Listeria monocytogenes biofilm by phages belonging to the genus Pecentumvirus

Listeria monocytogenes is a pathogenic foodborne bacterium that is a significant cause of mortality associated with foodborne illness and causes many food recalls attributed to a bacteriological cause. Their ability to form biofilms contributes to the persistence of Listeria spp. in food processing environments. When growing as biofilms, L. monocytogenes are more resistant to sanitizers used in the food industry, such as benzalkonium chloride (BAC), as well as to physical stresses like desiccation and starvation. Lytic phages of Listeria are antagonistic to a broad range of Listeria spp. and may, therefore, have utility in reducing the occurrence of Listeria-associated food recalls by preventing food contamination. We screened nine closely related Listeria phages, including the commercially available Listex P100, for host range and ability to degrade microtiter plate biofilms of L. monocytogenes ATCC 19111 (serovar 1/2a). One phage, CKA15, was selected and shown to rapidly adsorb to its host under conditions relevant to applying the phage in dairy processing environments. Under simulated dairy processing conditions (SDPC), CKA15 caused a 2-log reduction in Lm19111 biofilm bacteria. This work supports the biosanitation potential of phage CKA15 and provides a basis for further investigation of phage-bacteria interactions in biofilms grown under SDPC.

IMPORTANCE

Listeria monocytogenes is a pathogenic bacterium that is especially dangerous for children, the elderly, pregnant women, and immune-compromised people. Because of this, the food industry takes its presence in their plants seriously. Food recalls due to L. monocytogenes are common with a high associated economic cost. In food-processing plants, Listeria spp. typically reside in biofilms, which are structures produced by bacteria that shield them from environmental stressors and are often attached to surfaces. The significance of our work is that we show a bacteriophage-a virus-infecting bacteria-can reduce Listeria counts by two orders of magnitude when the bacterial biofilms were grown under simulated dairy processing conditions. This work provides insights into how phages may be tested and used to develop biosanitizers that are effective but are not harmful to the environment or human health.

A chromosome level genome assembly of Pseudoroegneria Libanotica reveals a key Kcs gene involves in the cuticular wax elongation for drought resistance

BACKGROUND

The genus Pseudoroegneria (Nevski) Löve (Triticeae, Poaceae), whose genome symbol was designed as "St", accounts for more than 60% of perennial Triticeae species. The diploid species Psudoroegneria libanotica (2n = 14) contains the most ancient St genome, exhibited strong drought resistance, and was morphologically covered by cuticular wax on the aerial part. Therefore, the St-genome sequencing data could provide fundamental information for studies of genome evolution and reveal its mechanisms of cuticular wax and drought resistance.

RESULTS

In this study, we reported the chromosome-level genome assembly for the St genome of Pse. libanotica, with a total size of 2.99 Gb. 46,369 protein-coding genes annotated and 71.62% was repeat sequences. Comparative analyses revealed that the genus Pseudoroegneria diverged during the middle and late Miocene. During this period, unique genes, gene family expansion, and contraction in Pse. libanotica were enriched in biotic and abiotic stresses, such as fatty acid biosynthesis which may greatly contribute to its drought adaption. Furthermore, we investigated genes associated with the cuticular wax formation and water deficit and found a new Kcs gene evm.TU.CTG175.54. It plays a critical role in the very long chain fatty acid (VLCFA) elongation from C18 to C26 in Pse. libanotica. The function needs more evidence to be verified.

CONCLUSIONS

We sequenced and assembled the St genome in Triticeae and discovered a new KCS gene that plays a role in wax extension to cope with drought. Our study lays a foundation for the genome diversification of Triticeae species and deciphers cuticular wax formation genes involved in drought resistance.

Whole genome sequencing identifies novel mutations in malaria parasites resistant to artesunate (ATN) and to ATN + mefloquine combination

Introduction

The global evolution of resistance to Artemisinin-based Combination Therapies (ACTs) by malaria parasites, will severely undermine our ability to control this devastating disease.

Methods

Here, we have used whole genome sequencing to characterize the genetic variation in the experimentally evolved Plasmodium chabaudi parasite clone AS-ATNMF1, which is resistant to artesunate + mefloquine.

Results and discussion

Five novel single nucleotide polymorphisms (SNPs) were identified, one of which was a previously undescribed E738K mutation in a 26S proteasome subunit that was selected for under artesunate pressure (in AS-ATN) and retained in AS-ATNMF1. The wild type and mutated three-dimensional (3D) structure models and molecular dynamics simulations of the P. falciparum 26S proteasome subunit Rpn2 suggested that the E738K mutation could change the toroidal proteasome/cyclosome domain organization and change the recognition of ubiquitinated proteins. The mutation in the 26S proteasome subunit may therefore contribute to altering oxidation-dependent ubiquitination of the MDR-1 and/or K13 proteins and/or other targets, resulting in changes in protein turnover. In light of the alarming increase in resistance to artemisin derivatives and ACT partner drugs in natural parasite populations, our results shed new light on the biology of resistance and provide information on novel molecular markers of resistance that may be tested (and potentially validated) in the field.

Genetic dissection of domestication traits in interspecific chickpea populations

Chickpea (Cicer arietinum) is a pulse crop that provides an integral source of nutrition for human consumption. The close wild relatives Cicer reticulatum and Cicer echinospermum harbor untapped genetic diversity that can be exploited by chickpea breeders to improve domestic varieties. Knowledge of genomic loci that control important chickpea domestication traits will expedite the development of improved chickpea varieties derived from interspecific crosses. Therefore, we set out to identify genomic loci underlying key chickpea domestication traits by both association and quantitative trait locus (QTL) mapping using interspecific F2 populations. Diverse phenotypes were recorded for various agronomic traits. A total of 11 high-confidence markers were detected on chromosomes 1, 3, and 7 by both association and QTL mapping; these were associated with growth habit, flowering time, and seed traits. Furthermore, we identified candidate genes linked to these markers, which advanced our understanding of the genetic basis of domestication traits and validated known genes such as the FLOWERING LOCUS gene cluster that regulates flowering time. Collectively, this study has elucidated the genetic basis of chickpea domestication traits, which can facilitate the development of superior chickpea varieties.

Dipterocarpoidae genomics reveal their demography and adaptations to Asian rainforests

Dipterocarpoideae species form the emergent layer of Asian rainforests. They are the indicator species for Asian rainforest distribution, but they are severely threatened. Here, to understand their adaptation and population decline, we assemble high-quality genomes of seven Dipterocarpoideae species including two autotetraploid species. We estimate the divergence time between Dipterocarpoideae and Malvaceae and within Dipterocarpoideae to be 108.2 (97.8‒118.2) and 88.4 (77.7‒102.9) million years ago, and we identify a whole genome duplication event preceding dipterocarp lineage diversification. We find several genes that showed a signature of selection, likely associated with the adaptation to Asian rainforests. By resequencing of two endangered species, we detect an expansion of effective population size after the last glacial period and a recent sharp decline coinciding with the history of local human activities. Our findings contribute to understanding the diversification and adaptation of dipterocarps and highlight anthropogenic disturbances as a major factor in their endangered status.

Chromosomal level genome assemblies of two Malus crabapple cultivars Flame and Royalty

Malus hybrid 'Flame' and Malus hybrid 'Royalty' are representative ornamental crabapples, rich in flavonoids and serving as the preferred materials for studying the coloration mechanism. We generated two sets of high-quality chromosome-level and haplotype-resolved genome of 'Flame' with sizes of 688.2 Mb and 675.7 Mb, and those of 'Royalty' with sizes of 674.1 Mb and 663.6 Mb, all anchored to 17 chromosomes and with a high BUSCO completeness score nearly 99.0%. A total of 47,833 and 47,307 protein-coding genes were annotated in the two haplotype genomes of 'Flame', and the numbers of 'Royalty' were 46,305 and 46,920 individually. The assembled high-quality genomes offer new resources for studying the origin and adaptive evolution of crabapples and the molecular basis of the accumulation of flavonoids and anthocyanins, facilitating molecular breeding of Malus plants.

Computational investigations into structure and function impact of novel mutations identified in targeted exons from ovarian cancer cell lines

The lack of sensitive and specific biomarkers for ovarian cancer leads to late stage diagnosis of the disease in a majority of the cases. Mutation accumulation is the basis for cancer progression, thus identifying mutations is an important step in the disease diagnosis. In the present study, a comprehensive analysis of fifteen Next Generation Sequencing samples from thirteen ovarian cancer cell lines was carried out for the identification of new mutations. The study revealed eight clinically significant novel mutations in six ovarian cancer oncogenes, viz. SMARCA4, ARID1A, PPP2R1A, CTNNB1, DICER1 and PIK3CA. In-depth computational analysis revealed that the mutations affected the structure of the proteins in terms of stability, solvent accessible surface area and molecular dynamics. Moreover, the mutations were present in functionally significant domains of the proteins, thereby adversely affecting the protein functionality. PPI network for SMARCA4, CTNNB1, DICER1, PIK3CA, PPP2R1A and ARID1A showed that these genes were involved in certain significant pathways affecting various hallmarks of cancer. For further validation, in vitro studies were performed that revealed hypermutability of the CTNNB1 gene. Through this study we have identified some key mutations and have analysed their structural and functional impact. The study establishes some key mutations, which can be potentially explored as biomarker and drug target.Communicated by Ramaswamy H. Sarma.

Glutaredoxin regulation of primary root growth is associated with early drought stress tolerance in pearl millet

Seedling root traits impact plant establishment under challenging environments. Pearl millet is one of the most heat and drought tolerant cereal crops that provides a vital food source across the sub-Saharan Sahel region. Pearl millet's early root system features a single fast-growing primary root which we hypothesize is an adaptation to the Sahelian climate. Using crop modeling, we demonstrate that early drought stress is an important constraint in agrosystems in the Sahel where pearl millet was domesticated. Furthermore, we show that increased pearl millet primary root growth is correlated with increased early water stress tolerance in field conditions. Genetics including genome-wide association study and quantitative trait loci (QTL) approaches identify genomic regions controlling this key root trait. Combining gene expression data, re-sequencing and re-annotation of one of these genomic regions identified a glutaredoxin-encoding gene PgGRXC9 as the candidate stress resilience root growth regulator. Functional characterization of its closest Arabidopsis homolog AtROXY19 revealed a novel role for this glutaredoxin (GRX) gene clade in regulating cell elongation. In summary, our study suggests a conserved function for GRX genes in conferring root cell elongation and enhancing resilience of pearl millet to its Sahelian environment.

Conservation genomics analysis reveals recent population decline and possible causes in bumblebee Bombus opulentus

Bumblebees are a genus of pollinators (Bombus) that play important roles in natural ecosystem and agricultural production. Several bumblebee species have been recorded as under population decline, and the proportion of species experiencing population decline within subgenus Thoracobombus is higher than average. Bombus opulentus is 1 species in Thoracobombus, but little is known about its recent population dynamics. Here, we employed conservation genomics methods to investigate the population dynamics of B. opulentus during the recent past and identify the likely environmental factors that may cause population decline. Firstly, we placed the scaffold-level of B. opulentus reference genome sequence onto chromosome-level using Hi-C technique. Then, based on this reference genome and whole-genome resequencing data for 51 B. opulentus samples, we reconstructed the population structure and effective population size (Ne ) trajectories of B. opulentus and identified genes that were under positive selection. Our results revealed that the collected B. opulentus samples could be divided into 2 populations, and 1 of them experienced a recent population decline; the declining population also exhibited lower genetic diversity and higher inbreeding levels. Genes related to high-temperature tolerance, immune response, and detoxication showed signals of positive selection in the declining population, suggesting that climate warming and pathogen/pesticide exposures may contribute to the decline of this B. opulentus population. Taken together, our study provided insights into the demography of B. opulentus populations and highlighted that populations of the same bumblebee species could have contrasting Ne trajectories and population decline could be caused by a combination of various stressors.

Genome and haplotype provide insights into the population differentiation and breeding improvement of Gossypium barbadense

INTRODUCTION

Sea-island cotton (Gossypium barbadense, Gb) is one of the major sources of high-grade natural fiber. Besides the common annual Gb cotton, perennial Gb cotton is also cultivated, but studies on perennial Gb cotton are rare.

OBJECTIVES

We aimed to make a systematic analysis of perennial sea-island cotton and lay a foundation for its utilization in breeding, and try to identify the representative structural variations (SVs) in sea-island cotton, and to reveal the population differentiation and adaptive improvement of sea-island cotton.

METHODS

Through genome assembly of one perennial Gb cotton accession (named Gb_M210936) and comparative genome analysis, variations during Gb cotton domestication were identified by comparing Gb_M210936 with annual Gb accession 3-79 and with wild allotetraploid cotton G. darwinii. Six perennial Gb accessions combining with the resequenced 1,129 cotton accessions were used to conduct population and genetic analysis. Large haplotype blocks (haploblocks), generated from interspecific introgressions and intraspecific inversions, were identified and were used to analyze their effects on population differentiation and agronomic traits of sea-island cotton.

RESULTS

One reference genome of perennial sea-island cotton was assembled. Representative SVs in sea-island cotton were identified, and 31 SVs were found to be associated with agronomic traits. Perennial Gb cotton had a closer kinship with the wild-to-landrace continuum Gb cotton from south America where Gb cotton is originally domesticated. Haploblocks were associated with agronomic traits improvement of sea-island cotton, promoted sea-island cotton differentiation into three subgroups, were suffered from breeding selection, and may drive Gb cotton to be adapted to central Asian.

CONCLUSION

Our study made up the lack of perennial Gb cotton genome, and clarified that exotic introgressions improved the traits of sea-island cotton, promoted the population differentiation, and drove sea-island cotton adaptive to central Asia, which will provide new insights for the genetic breeding improvement of sea-island cottons.

Chromosome-scale genome of the human blood fluke Schistosoma mekongi and its implications for public health

BACKGROUND

Schistosoma mekongi is a human blood fluke causing schistosomiasis that threatens approximately 1.5 million humans in the world. Nonetheless, the limited available S. mekongi genomic resources have hindered understanding of its biology and parasite-host interactions for disease management and pathogen control. The aim of our study was to integrate multiple technologies to construct a high-quality chromosome-level assembly of the S. mekongi genome.

METHODS

The reference genome for S. mekongi was generated through integrating Illumina, PacBio sequencing, 10 × Genomics linked-read sequencing, and high-throughput chromosome conformation capture (Hi-C) methods. In this study, we conducted de novo assembly, alignment, and gene prediction to assemble and annotate the genome. Comparative genomics allowed us to compare genomes across different species, shedding light on conserved regions and evolutionary relationships. Additionally, our transcriptomic analysis focused on genes associated with parasite-snail interactions in S. mekongi infection. We employed gene ontology (GO) enrichment analysis for functional annotation of these genes.

RESULTS

In the present study, the S. mekongi genome was both assembled into 8 pseudochromosomes with a length of 404 Mb, with contig N50 and scaffold N50 lengths of 1168 kb and 46,759 kb, respectively. We detected that 43% of the genome consists of repeat sequences and predicted 9103 protein-coding genes. We also focused on proteases, particularly leishmanolysin-like metalloproteases (M8), which are crucial in the invasion of hosts by 12 flatworm species. Through phylogenetic analysis, it was discovered that the M8 gene exhibits lineage-specific amplification among the genus Schistosoma. Lineage-specific expansion of M8 was observed in blood flukes. Additionally, the results of the RNA-seq revealed that a mass of genes related to metabolic and biosynthetic processes were up-regulated, which might be beneficial for cercaria production.

CONCLUSIONS

This study delivers a high-quality, chromosome-scale reference genome of S. mekongi, enhancing our understanding of the divergence and evolution of Schistosoma. The molecular research conducted here also plays a pivotal role in drug discovery and vaccine development. Furthermore, our work greatly advances the understanding of host-parasite interactions, providing crucial insights for schistosomiasis intervention strategies.

A chromosome-level genome assembly provides insights into Cornus wilsoniana evolution, oil biosynthesis, and floral bud development

Cornus wilsoniana W. is a woody oil plant with high oil content and strong hypolipidemic effects, making it a valuable species for medicinal, landscaping, and ecological purposes in China. To advance genetic research on this species, we employed PacBio together with Hi-C data to create a draft genome assembly for C. wilsoniana. Based on an 11-chromosome anchored chromosome-level assembly, the estimated genome size was determined to be 843.51 Mb. The N50 contig size and N50 scaffold size were calculated to be 4.49 and 78.00 Mb, respectively. Furthermore, 30 474 protein-coding genes were annotated. Comparative genomics analysis revealed that C. wilsoniana diverged from its closest species ~12.46 million years ago (Mya). Furthermore, the divergence between Cornaceae and Nyssaceae occurred >62.22 Mya. We also found evidence of whole-genome duplication events and whole-genome triplication γ, occurring at ~44.90 and 115.86 Mya. We further inferred the origins of chromosomes, which sheds light on the complex evolutionary history of the karyotype of C. wilsoniana. Through transcriptional and metabolic analysis, we identified two FAD2 homologous genes that may play a crucial role in controlling the oleic to linoleic acid ratio. We further investigated the correlation between metabolites and genes and identified 33 MADS-TF homologous genes that may affect flower morphology in C. wilsoniana. Overall, this study lays the groundwork for future research aimed at identifying the genetic basis of crucial traits in C. wilsoniana.

Chromosome level genome assembly of colored calla lily (Zantedeschia elliottiana)

The colored calla lily is an ornamental floral plant native to southern Africa, belonging to the Zantedeschia genus of the Araceae family. We generated a high-quality chromosome-level genome of the colored calla lily, with a size of 1,154 Mb and a contig N50 of 42 Mb. We anchored 98.5% of the contigs (1,137 Mb) into 16 pseudo-chromosomes, and identified 60.18% of the sequences (694 Mb) as repetitive sequences. Functional annotations were assigned to 95.1% of the predicted protein-coding genes (36,165). Additionally, we annotated 469 miRNAs, 1,652 tRNAs, 10,033 rRNAs, and 1,677 snRNAs. Furthermore, Gypsy-type LTR retrotransposons insertions in the genome are the primary factor causing significant genome size variation in Araceae species. This high-quality genome assembly provides valuable resources for understanding genome size differences within the Araceae family and advancing genomic research on colored calla lily.

Multiple roles of wheat ferritin genes during stress treatment and TaFER5D-1 as a positive regulator in response to drought and salt tolerance

Ferritin not only regulates the plant's iron content but also plays a significant role in the plant's development and resistance to oxidative damage. However, the role of the FER family in wheat has not been systematically elucidated. In this study, 39 FERs identified from wheat and its ancestral species were clustered into two subgroups, and gene members from the same group contain relatively conservative protein models. The structural analyses indicated that the gene members from the same group contained relatively conserved protein models. The cis-acting elements and expression patterns analysis suggested that TaFERs might play an important role combating to abiotic and biotic stresses. In the transcriptional analysis, the TaFER5D-1 gene was found to be significantly up-regulated under drought and salt stresses and was, therefore, selected to further explore the biological functions Moreover, the GFP expression assay revealed the subcellular localization of TaFER5D-1 proteins in the chloroplast, nucleus, membrane and cytoplasm. Over-expression of TaFER5D-1 in transgenic Arabidopsis lines conferred greater tolerance to drought and salt stress. According to the qRT-PCR data, TaFER5D-1 gene over-expression increased the expression of genes related to root development (Atsweet-17 and AtRSL4), iron storage (AtVIT1 and AtYSL1), and stress response (AtGolS1 and AtCOR47). So it is speculated that TaFER5D-1 could improve stress tolerance by promoting root growth, iron storage, and stress-response ability. Thus, the current study provides insight into the role of TaFER genes in wheat.

Genome assembly of Polygala tenuifolia provides insights into its karyotype evolution and triterpenoid saponin biosynthesis

Polygala tenuifolia is a perennial medicinal plant that has been widely used in traditional Chinese medicine for treating mental diseases. However, the lack of genomic resources limits the insight into its evolutionary and biological characterization. In the present work, we reported the P. tenuifolia genome, the first genome assembly of the Polygalaceae family. We sequenced and assembled this genome by a combination of Illumnina, PacBio HiFi, and Hi-C mapping. The assembly includes 19 pseudochromosomes covering ~92.68% of the assembled genome (~769.62 Mb). There are 36 463 protein-coding genes annotated in this genome. Detailed comparative genome analysis revealed that P. tenuifolia experienced two rounds of whole genome duplication that occurred ~39-44 and ~18-20 million years ago, respectively. Accordingly, we systematically reconstructed ancestral chromosomes of P. tenuifolia and inferred its chromosome evolution trajectories from the common ancestor of core eudicots to the present species. Based on the transcriptomics data, enzyme genes and transcription factors involved in the synthesis of triterpenoid saponin in P. tenuifolia were identified. Further analysis demonstrated that whole-genome duplications and tandem duplications play critical roles in the expansion of P450 and UGT gene families, which contributed to the synthesis of triterpenoid saponins. The genome and transcriptome data will not only provide valuable resources for comparative and functional genomic researches on Polygalaceae, but also shed light on the synthesis of triterpenoid saponin.

Lampshade web spider Ectatosticta davidi chromosome-level genome assembly provides evidence for its phylogenetic position

The spider of Ectatosticta davidi, belonging to the lamp-shade web spider family, Hypochilidae, which is closely related to Hypochilidae and Filistatidae and recovered as sister of the rest Araneomorphs spiders. Here we show the final assembled genome of E. davidi with 2.16 Gb in 15 chromosomes. Then we confirm the evolutionary position of Hypochilidae. Moreover, we find that the GMC gene family exhibit high conservation throughout the evolution of true spiders. We also find that the MaSp genes of E. davidi may represent an early stage of MaSp and MiSp genes in other true spiders, while CrSp shares a common origin with AgSp and PySp but differ from MaSp. Altogether, this study contributes to addressing the limited availability of genomic sequences from Hypochilidae spiders, and provides a valuable resource for investigating the genomic evolution of spiders.

Characterization of crAss-like phage isolates highlights Crassvirales genetic heterogeneity and worldwide distribution

Crassvirales (crAss-like phages) are an abundant group of human gut-specific bacteriophages discovered in silico. The use of crAss-like phages as human fecal indicators is proposed but the isolation of only seven cultured strains of crAss-like phages to date has greatly hindered their study. Here, we report the isolation and genetic characterization of 25 new crAss-like phages (termed crAssBcn) infecting Bacteroides intestinalis, belonging to the order Crassvirales, genus Kehishuvirus and, based on their genomic variability, classified into six species. CrAssBcn phage genomes are similar to ΦCrAss001 but show genomic and aminoacidic differences when compared to other crAss-like phages of the same family. CrAssBcn phages are detected in fecal metagenomes around the world at a higher frequency than ΦCrAss001. This study increases the known crAss-like phage isolates and their abundance and heterogeneity open the question of what member of the Crassvirales group should be selected as human fecal marker.

Chromosome-level genome assembly of the sea cucumber Apostichopus japonicus

Sea cucumber is a morphologically diverse and ecologically important clade of echinoderms. The sea cucumber Apostichopus japonicus is the most economically valuable species of sea cucumber. The initial assembly of the A. japonicus genome was released in 2017. However, this genome assembly is fragmented and lacks relative position information of genes on chromosomes. In this study, we produced a high-quality chromosome-level genome of A. japonicus using Pacbio HiFi long-reads and Hi-C sequencing data. The assembled A. japonicus genome spanned 671.60 Mb with a contig N50 size of 17.20 Mb and scaffold N50 size of 29.65 Mb. A total of 99.9% of the assembly was anchored to 23 chromosomes. In total, 19,828 genes were annotated, and 97.2% of BUSCO genes were fully represented. This high-quality genome of A. japonicus will not only aid in the development of sustainable aquaculture practices, but also lay a foundation for a deeper understanding of their genetic makeup, evolutionary history, and ecological adaptation.

Chromosome-level genome and multi-omics analyses provide insights into the geo-herbalism properties of Alpinia oxyphylla

Introduction

Alpinia oxyphylla Miquel (A. oxyphylla), one of the "Four Famous South Medicines" in China, is an essential understory cash crop that is planted widely in the Hainan, Guangdong, Guangxi, and Fujian provinces. Particularly, A. oxyphylla from Hainan province is highly valued as the best national product for geo-herbalism and is an important indicator of traditional Chinese medicine efficacy. However, the molecular mechanism underlying the formation of its quality remains unspecified.

Methods

To this end, we employed a multi-omics approach to investigate the authentic quality formation of A. oxyphylla.

Results

In this study, we present a high-quality chromosome-level genome assembly of A. oxyphylla, with contig N50 of 76.96 Mb and a size of approximately 2.08Gb. A total of 38,178 genes were annotated, and the long terminal repeats were found to have a high frequency of 61.70%. Phylogenetic analysis demonstrated a recent whole-genome duplication event (WGD), which occurred before A. oxyphylla's divergence from W. villosa (~14 Mya) and is shared by other species from the Zingiberaceae family (Ks, ~0.3; 4DTv, ~0.125). Further, 17 regions from four provinces were comprehensively assessed for their metabolite content, and the quality of these four regions varied significantly. Finally, genomic, metabolic, and transcriptomic analyses undertaken on these regions revealed that the content of nootkatone in Hainan was significantly different from that in other provinces.

Discussion

Overall, our findings provide novel insights into germplasm conservation, geo-herbalism evaluation, and functional genomic research for the medicinal plant A. oxyphylla.

Gap-free genome assembly of anadromous Coilia nasus

The Chinese tapertail anchovy, Coilia nasus, is a socioeconomically important anadromous fish that migrates from near ocean waters to freshwater to spawn every spring. The analysis of genomic architecture and information of C. nasus were hindered by the previously released versions of reference genomes with gaps. Here, we report the assembly of a chromosome-level gap-free genome of C. nasus by incorporating high-coverage and accurate long-read sequence data with multiple assembly strategies. All 24 chromosomes were assembled without gaps, representing the highest completeness and assembly quality. We assembled the genome with a size of 851.67 Mb and used BUSCO to estimate the completeness of the assembly as 92.5%. Using a combination of de novo prediction, protein homology and RNA-seq annotation, 21,900 genes were functionally annotated, representing 99.68% of the total predicted protein-coding genes. The availability of gap-free reference genomes for C. nasus will provide the opportunity for understanding genome structure and function, and will also lay a solid foundation for further management and conservation of this important species.

Functional annotation of a hugely expanded nanos repertoire in Lytechinus variegatus, the green sea urchin

Nanos genes encode essential RNA-binding proteins involved in germline determination and germline stem cell maintenance. When examining diverse classes of echinoderms, typically three, sometimes four, nanos genes are present. In this analysis, we identify and annotate nine nanos orthologs in the green sea urchin, Lytechinus variegatus (Lv). All nine genes are transcribed and grouped into three distinct classes. Class one includes the germline Nanos, with one member: Nanos2. Class two includes Nanos3-like genes, with significant sequence similarity to Nanos3 in the purple sea urchin, Strongylocentrotus purpuratus (Sp), but with wildly variable expression patterns. The third class includes several previously undescribed nanos zinc-finger genes that may be the result of duplications of Nanos2. All nine nanos transcripts occupy unique genomic loci and are expressed with unique temporal profiles during development. Importantly, here we describe and characterize the unique genomic location, conservation, and phylogeny of the Lv ortholog of the well-studied Sp Nanos2. However, in addition to the conserved germline functioning Nanos2, the green sea urchin appears to be an outlier in the echinoderm phyla with eight additional nanos genes. We hypothesize that this expansion of nanos gene members may be the result of a previously uncharacterized L1-class transposon encoded on the opposite strand of a nanos2 pseudogene present on chromosome 12 in this species. The expansion of nanos genes described here represents intriguing insights into germline specification and nanos evolution in this species of sea urchin.

Population genomic analysis provides evidence of the past success and future potential of South China tiger captive conservation

BACKGROUND

Among six extant tiger subspecies, the South China tiger (Panthera tigris amoyensis) once was widely distributed but is now the rarest one and extinct in the wild. All living South China tigers are descendants of only two male and four female wild-caught tigers and they survive solely in zoos after 60 years of effective conservation efforts. Inbreeding depression and hybridization with other tiger subspecies were believed to have occurred within the small, captive South China tiger population. It is therefore urgently needed to examine the genomic landscape of existing genetic variation among the South China tigers.

RESULTS

In this study, we assembled a high-quality chromosome-level genome using long-read sequences and re-sequenced 29 high-depth genomes of the South China tigers. By combining and comparing our data with the other 40 genomes of six tiger subspecies, we identified two significantly differentiated genomic lineages among the South China tigers, which harbored some rare genetic variants introgressed from other tiger subspecies and thus maintained a moderate genetic diversity. We noticed that the South China tiger had higher F_ROH values for longer runs of homozygosity (ROH > 1 Mb), an indication of recent inbreeding/founder events. We also observed that the South China tiger had the least frequent homozygous genotypes of both high- and moderate-impact deleterious mutations, and lower mutation loads than both Amur and Sumatran tigers. Altogether, our analyses indicated an effective genetic purging of deleterious mutations in homozygous states from the South China tiger, following its population contraction with a controlled increase in inbreeding based on its pedigree records.

CONCLUSIONS

The identification of two unique founder/genomic lineages coupled with active genetic purging of deleterious mutations in homozygous states and the genomic resources generated in our study pave the way for a genomics-informed conservation, following the real-time monitoring and rational exchange of reproductive South China tigers among zoos.

Single-cell long-read mRNA isoform regulation is pervasive across mammalian brain regions, cell types, and development

RNA isoforms influence cell identity and function. Until recently, technological limitations prevented a genome-wide appraisal of isoform influence on cell identity in various parts of the brain. Using enhanced long-read single-cell isoform sequencing, we comprehensively analyze RNA isoforms in multiple mouse brain regions, cell subtypes, and developmental timepoints from postnatal day 14 (P14) to adult (P56). For 75% of genes, full-length isoform expression varies along one or more axes of phenotypic origin, underscoring the pervasiveness of isoform regulation across multiple scales. As expected, splicing varies strongly between cell types. However, certain gene classes including neurotransmitter release and reuptake as well as synapse turnover, harbor significant variability in the same cell type across anatomical regions, suggesting differences in network activity may influence cell-type identity. Glial brain-region specificity in isoform expression includes strong poly(A)-site regulation, whereas neurons have stronger TSS regulation. Furthermore, developmental patterns of cell-type specific splicing are especially pronounced in the murine adolescent transition from P21 to P28. The same cell type traced across development shows more isoform variability than across adult anatomical regions, indicating a coordinated modulation of functional programs dictating neural development. As most cell-type specific exons in P56 mouse hippocampus behave similarly in newly generated data from human hippocampi, these principles may be extrapolated to human brain. However, human brains have evolved additional cell-type specificity in splicing, suggesting gain-of-function isoforms. Taken together, we present a detailed single-cell atlas of full-length brain isoform regulation across development and anatomical regions, providing a previously unappreciated degree of isoform variability across multiple scales of the brain.

Genome assembly of the Brassicaceae diploid Orychophragmus violaceus reveals complex whole-genome duplication and evolution of dihydroxy fatty acid metabolism

Orychophragmus violaceus is a Brassicaceae species widely cultivated in China, particularly as a winter cover crop in northern China because of its low-temperature tolerance and low water demand. Recently, O. violaceus has also been cultivated as a potential industrial oilseed crop because of its abundant 24-carbon dihydroxy fatty acids (diOH-FAs), which contribute to superior high-temperature lubricant properties. In this study, we performed de novo assembly of the O. violaceus genome. Whole-genome synteny analysis of the genomes of its relatives demonstrated that O. violaceus is a diploid that has undergone an extra whole-genome duplication (WGD) after the Brassicaceae-specific α-WGD event, with a basic chromosome number of x = 12. Formation of diOH-FAs is hypothesized to have occurred after the WGD event. Based on the genome and the transcriptome data from multiple stages of seed development, we predicted that OvDGAT1-1 and OvDGAT1-2 are candidate genes for the regulation of diOH-FA storage in O. violaceus seeds. These results may greatly facilitate the development of heat-tolerant and eco-friendly plant-based lubricants using O. violaceus seed oil and improve our understanding of the genomic evolution of Brassicaceae.

The genome of a vestimentiferan tubeworm (Ridgeia piscesae) provides insights into its adaptation to a deep-sea environment

BACKGROUND

Vestimentifera (Polychaeta, Siboglinidae) is a taxon of deep-sea worm-like animals living in deep-sea hydrothermal vents, cold seeps, and organic falls. The morphology and lifespan of Ridgeia piscesae, which is the only vestimentiferan tubeworm species found in the hydrothermal vents on the Juan de Fuca Ridge, vary greatly according to endemic environment. Recent analyses have revealed the genomic basis of adaptation in three vent- and seep-dwelling vestimentiferan tubeworms. However, the evolutionary history and mechanism of adaptation in R. piscesae, a unique species in the family Siboglinidae, remain to be investigated.

RESULT

We assembled a draft genome of R. piscesae collected at the Cathedral vent of the Juan de Fuca Ridge. Comparative genomic analysis showed that vent-dwelling tubeworms with a higher growth rate had smaller genome sizes than seep-dwelling tubeworms that grew much slower. A strong positive correlation between repeat content and genome size but not intron size and the number of protein-coding genes was identified in these deep-sea tubeworm species. Evolutionary analysis revealed that Ridgeia pachyptila and R. piscesae, the two tubeworm species that are endemic to hydrothermal vents of the eastern Pacific, started to diverge between 28.5 and 35 million years ago. Four genes involved in cell proliferation were found to be subject to positive selection in the genome of R. piscesae.

CONCLUSION

Ridgeia pachyptila and R. piscesae started to diverge after the formation of the Gorda/Juan de Fuca/Explorer ridge systems and the East Pacific Rise. The high growth rates of vent-dwelling tubeworms might be derived from their small genome sizes. Cell proliferation is important for regulating the growth rate in R. piscesae.

Comparative genomics of five Valsa species gives insights on their pathogenicity evolution

Valsa is a genus of ascomycetes within the Valsaceae family. This family includes many wood destructive pathogens such as the well known Valsa mali and Valsa pyri which cause canker diseases in fruit trees and threaten the global fruit production. Lack of genomic information of this family is impeding our understandings about their evolution and genetic basis of their pathogenicity divergence. Here, we report genome assemblies of Valsa malicola, Valsa persoonii, and Valsa sordida which represent close relatives of Valsa mali and Valsa pyri with different host preferences. Comparative genomics analysis revealed that segmental rearrangements, inversions, and translocations frequently occurred among Valsa spp. genomes. Gene families that exhibited gene copy expansions tended to be associated with secondary metabolism, transmembrane transport, and pyrophosphatase activities. Orthologous genes in regions lost synteny exhibited significantly higher rate of synonymous substitution (KS) than those in regions retained synteny. Moreover, among these genes, membrane transporter families associated with antidrug (MFS, DHA) activities and nutrient transportation (SP and APCs) activities were significantly over-represented. Lineage specific synonymous substitution (KS) and nonsynonymous substitution (KA) analysis based on the phylogeny constructed from 11 fungal species identified a set of genes with selection signatures in Valsa clade and these genes were significantly enriched in functions associated with fatty acid beta-oxidation, DNA helicase activity, and ATPase activity. Furthermore, unique genes that possessed or retained by each of the five Valsa species are more likely part of the secondary metabolic (SM) gene clusters. SM gene clusters conserved across five Valsa species showed various degrees of diversification in both identity and completeness. All 11 syntenically conserved SM clusters showed differential expression during the infection of apple branch with Valsa mali suggesting involvements of secondary metabolism in the pathogenicity of Valsa species.

The Chromosome-Level Genome of Hestina assimilis (Lepidoptera: Nymphalidae) Reveals the Evolution of Saprophagy-Related Genes in Brush-Footed Butterflies

Most butterflies feed on nectar, while some saprophagous butterflies forage on various non-nectar foods. To date, little is known about the genomic and molecular shifts associated with the evolution of the saprophagous feeding strategy. Here, we assembled the high-quality chromosome-level genome of Hestina assimilis to explore its saprophagous molecular and genetic mechanisms. This chromosome-level genome of H. assimilis is 412.82 Mb, with a scaffold N50 of 15.70 Mb. In total, 98.11% of contigs were anchored to 30 chromosomes. Compared with H. assimilis and other Nymphalidae butterflies, the genes of metabolism and detoxification experienced expansions. We annotated 80 cytochrome P450 (CYP) genes in the H. assimilis genome, among which genes belonging to the CYP4 subfamily were significantly expanded (p < 0.01). These P450 genes were unevenly distributed and mainly concentrated on chromosomes 6-9. We identified 33 olfactory receptor (OR), 20 odorant-binding protein (OBP), and six gustatory receptor (GR) genes in the H. assimilis genome, which were fewer than in the nectarivorous Danaus plexippus. A decreased number of OBP, OR, and GR genes implied that H. assimilis should resort less to olfaction and gustation than their nectarivorous counterparts, which need highly specialized olfactory and gustatory functions. Moreover, we found one site under positive selection occurred in residue 996 (phenylalanine) of GR genes exclusive to H. assimilis, which is conservative in most lineages. Our study provides support for the adaptive evolution of feeding habits in butterflies.

Computational studies reveal co-occurrence of two mutations in IL7R gene of high-grade serous carcinoma patients

Major cause of mortality in ovarian cancer can be attributed to a lack of specific and sensitive biomarkers for diagnosis and prognosis of the disease. Uncovering the mutations in genes involved in crucial oncogenic pathways is a key step in discovery and development of novel biomarkers. Whole exome sequencing (WES) is a powerful method for the detection of cancer driver mutations. The present work focuses on identifying functionally damaging mutations in patients with high-grade serous ovarian carcinoma (HGSC) through computational analysis of WES. In this study, WES data of HGSC patients was retrieved from the genomic literature available in sequence read archive, the variants were identified and comprehensive structural and functional analysis was performed. Interestingly, I66T and V138I mutations were found to be co-occurring in the IL7R gene in four out of five HGSC patient samples investigated in this study. The V138I mutation was located in the fibronectin type-3 domain and computationally assessed to be causing disruptive effects on the structure and dynamics of IL7R protein. This mutation was found to be co-occurring with the neutral I66T mutation in the same domain which compensated the disruptive effects of V138I variant. These comprehensive studies point to a hitherto unexplored significant role of the IL7R gene in ovarian carcinoma. It is envisaged that the work will lay the foundation for the development of a novel biomarker with potential application in molecular profiling and in estimation of the disease prognosis.Communicated by Ramaswamy H. Sarma.

Gapless genome assembly of East Asian finless porpoise

In recent years, conservation efforts have increased for rare and endangered aquatic wildlife, especially cetaceans. However, the East Asian finless porpoise (Neophocaena asiaeorientalis sunameri), which has a wide distribution in China, has received far less attention and protection. As an endangered small cetacean, the lack of a chromosomal-level reference for the East Asian finless porpoise limits our understanding of its population genetics and conservation biology. To address this issue, we combined PacBio HiFi long reads and Hi-C sequencing data to generate a gapless genome of the East Asian finless porpoise that is approximately 2.5 Gb in size over its 21 autosomes and two sex chromosomes (X and Y). A total of 22,814 protein-coding genes were predicted where ~97.31% were functionally annotated. This high-quality genome assembly of East Asian finless porpoise will not only provide new resources for the comparative genomics of cetaceans and conservation biology of threatened species, but also lay a foundation for more speciation, ecology, and evolutionary studies. Measurement(s) Neophocaena asiaeorientalis sunameri • Gapless genome assembly • sequence annotation Technology Type(s) MGISEQ. 2000 • PacBio HiFi Sequencing • Hi-C Sample Characteristic - Organism Neophocaena asiaeorientalis sunameri Sample Characteristic - Environment seawater Sample Characteristic - Location Yellow Sea near Lianyungang City, Jiangsu Province, China.

Computationally Guided Enzymatic Studies on Schizochytrium-Sourced Malonyl-CoA:ACP Transacylase

The malonyl-CoA:ACP transacylase (MAT) domain is responsible for the selection and incorporation of malonyl building blocks in the biosynthesis of polyunsaturated fatty acids (PUFAs) in eukaryotic microalgae (Schizochytrium) and marine bacteria (Moritella marina, Photobacterium profundum, and Shewanella). Elucidation of the structural basis underlying the substrate specificity and catalytic mechanism of the MAT will help to improve the yield and quality of PUFAs. Here, a methodology guided by molecular dynamics simulations was carried out to identify and mutate specificity-conferring residues within the MAT domain of Schizochytrium. Combining mutagenesis, cell-free protein synthesis, and in vitro biochemical assay, we dissected nearby interactions and molecular mechanisms relevant for binding and catalysis and found that the reorientation of the Ser154 C^β-O^γ bond establishes distinctive proton-transfer chains (His153-Ser154 and Asn235-His153-Ser154) for catalysis. Gln66 can be replaced by tyrosine to shorten the distance between His153 (N^ε2) and Ser154 (O^γ), which facilitates a faster proton-transfer rate, allowing better use of acyl substrates than the wild type. Furthermore, we screened a mutant that displayed an 18.4% increase in PUFA accumulation. These findings provide important insights into the study of MAT through protein engineering and will benefit dissecting the molecular mechanisms of other PUFA-related catalytic domains.

Multi-omics provides new insights into the domestication and improvement of dark jute (Corchorus olitorius)

Jute (Corchorus sp.) is the most important bast fiber crop worldwide; however, the mechanisms underlying domestication and improvement remain largely unknown. We performed multi-omics analysis by integrating de novo sequencing, resequencing, and transcriptomic and epigenetic sequencing to clarify the domestication and improvement of dark jute Corchorus olitorius. We demonstrated that dark jute underwent early domestication and a relatively moderate genetic bottleneck during improvement breeding. A genome-wide association study of 11 important agronomic traits identified abundant candidate loci. We characterized the selective sweeps in the two breeding stages of jute, prominently, soil salinity differences played an important role in environmental adaptation during domestication, and the strongly selected genes for improvement had an increased frequency of favorable haplotypes. Furthermore, we speculated that an encoding auxin/indole-3-acetic acid protein COS07g_00652 could enhance the flexibility and strength of the stem to improve fiber yield. Our study not only provides valuable genetic resources for future fiber breeding in jute, but also is of great significance for reviewing the genetic basis of early crop breeding.

Genome-Wide Identification and Functional Analysis of the GASA Gene Family Responding to Multiple Stressors in Canavalia rosea

In plants, the Gibberellic Acid-Stimulated Arabidopsis (GASA) gene family is unique and responds to ubiquitous stress and hormones, playing important regulatory roles in the growth and development of plants, as well as in the resistance mechanisms to biotic and abiotic stress. In this study, a total of 23 CrGASAs were characterized in C. rosea using a genome-wide approach, and their phylogenetic relationships, gene structures, conserved motifs, chromosomal locations, gene duplications, and promoter regions were systematically analyzed. Expression profile analysis derived from transcriptome data showed that CrGASAs are expressed at higher levels in the flowers or fruit than in the leaves, vines, and roots. The expression of CrGASAs also showed habitat- and environmental-stress-regulated patterns in C. rosea analyzed by transcriptome and quantitative reverse transcription PCR (qRT-PCR). The heterologous induced expression of some CrGASAs in yeast enhanced the tolerance to H2O2, and some CrGASAs showed elevated heat tolerance and heavy metal (HM) Cd/Cu tolerance. These findings will provide an important foundation to elucidate the biological functions of CrGASA genes, especially their role in the ecological adaptation of specific plant species to tropical islands and reefs in C. rosea.

Genome-wide identification and expression analysis of SBP-box gene family reveal their involvement in hormone response and abiotic stresses in Chrysanthemum nankingense

SQUAMOSA promoter-binding-protein (SBP)-box family proteins are a class of plant-specific transcription factors, and widely regulate the development of floral and leaf morphology in plant growth and involve in environment and hormone signal response. In this study, we isolated and identified 21 non-redundant SBP-box genes in Chrysanthemum nankingense with bioinformatics analysis. Sequence alignments of 21 CnSBP proteins discovered a highly conserved SBP domain including two zinc finger-like structures and a nuclear localization signal region. According to the amino acid sequence alignments, 67 SBP-box genes from Arabidopsis thaliana, rice, Artemisia annua and C. nankingense were clustered into eight groups, and the motif and gene structure analysis also sustained this classification. The gene evolution analysis indicated the CnSBP genes experienced a duplication event about 10 million years ago (Mya), and the CnSBP and AtSPL genes occurred a divergence at 24 Mya. Transcriptome data provided valuable information for tissue-specific expression profiles of the CnSBPs, which highly expressed in floral tissues and differentially expressed in leaf, root and stem organs. Quantitative Real-time Polymerase Chain Reaction data showed expression patterns of the CnSBPs under exogenous hormone and abiotic stress treatments, separately abscisic acid, salicylic acid, gibberellin A3, methyl jasmonate and ethylene spraying as well as salt and drought stresses, indicating that the candidate CnSBP genes showed differentiated spatiotemporal expression patterns in response to hormone and abiotic stresses. Our study provides a systematic genome-wide analysis of the SBP-box gene family in C. nankingense. In general, it provides a fundamental theoretical basis that SBP-box genes may regulate the resistance of stress physiology in chrysanthemum via exogenous hormone pathways.

A metabolic, phylogenomic and environmental atlas of diatom plastid transporters from the model species Phaeodactylum

Diatoms are an important group of algae, contributing nearly 40% of total marine photosynthetic activity. However, the specific molecular agents and transporters underpinning the metabolic efficiency of the diatom plastid remain to be revealed. We performed in silico analyses of 70 predicted plastid transporters identified by genome-wide searches of Phaeodactylum tricornutum. We considered similarity with Arabidopsis thaliana plastid transporters, transcriptional co-regulation with genes encoding core plastid metabolic pathways and with genes encoded in the mitochondrial genomes, inferred evolutionary histories using single-gene phylogeny, and environmental expression trends using Tara Oceans meta-transcriptomics and meta-genomes data. Our data reveal diatoms conserve some of the ion, nucleotide and sugar plastid transporters associated with plants, such as non-specific triose phosphate transporters implicated in the transport of phosphorylated sugars, NTP/NDP and cation exchange transporters. However, our data also highlight the presence of diatom-specific transporter functions, such as carbon and amino acid transporters implicated in intricate plastid-mitochondria crosstalk events. These confirm previous observations that substrate non-specific triose phosphate transporters (TPT) may exist as principal transporters of phosphorylated sugars into and out of the diatom plastid, alongside suggesting probable agents of NTP exchange. Carbon and amino acid transport may be related to intricate metabolic plastid-mitochondria crosstalk. We additionally provide evidence from environmental meta-transcriptomic/meta- genomic data that plastid transporters may underpin diatom sensitivity to ocean warming, and identify a diatom plastid transporter (J43171) whose expression may be positively correlated with temperature.

Exploitation of a Bacterium-Encoded Lytic Transglycosylase by a Human Oral Lytic Phage To Facilitate Infection

Bacteriophages (phages) are an integral part of the human oral microbiome. Their roles in modulating bacterial physiology and shaping microbial communities have been discussed but remain understudied due to limited isolation and characterization of oral phage. Here, we report the isolation of LC001, a lytic phage targeting human oral Schaalia odontolytica (formerly known as Actinomyces odontolyticus) strain XH001. We showed that LC001 attached to and infected surface-grown, but not planktonic, XH001 cells, and it displayed remarkable host specificity at the strain level. Whole-genome sequencing of spontaneous LC001-resistant, surface-grown XH001 mutants revealed that the majority of the mutants carry nonsense or frameshift mutations in XH001 gene APY09_05145 (renamed ltg-1), which encodes a putative lytic transglycosylase (LT). The mutants are defective in LC001 binding, as revealed by direct visualization of the significantly reduced attachment of phage particles to the XH001 spontaneous mutants compared that to the wild type. Meanwhile, targeted deletion of ltg-1 produced a mutant that is defective in LC001 binding and resistant to LC001 infection even as surface-grown cells, while complementation of ltg-1 in the mutant background restored the LC001-sensitive phenotype. Intriguingly, similar expression levels of ltg-1 were observed in surface-grown and planktonic XH001, which displayed LC001-binding and nonbinding phenotypes, respectively. Furthermore, the overexpression of ltg-1 failed to confer an LC001-binding and -sensitive phenotype to planktonic XH001. Thus, our data suggested that rather than directly serving as a phage receptor, ltg-1-encoded LT may increase the accessibility of phage receptor, possibly via its enzymatic activity, by cleaving the peptidoglycan structure for better receptor exposure during peptidoglycan remodeling, a function that can be exploited by LC001 to facilitate infection. IMPORTANCE The evidence for the presence of a diverse and abundant phage population in the host-associated oral microbiome came largely from metagenomic analysis or the observation of virus-like particles within saliva/plaque samples, while the isolation of oral phage and investigation of their interaction with bacterial hosts are limited. Here, we report the isolation of LC001, the first lytic phage targeting oral Schaalia odontolytica. Our study suggested that LC001 may exploit the host bacterium-encoded lytic transglycosylase function to gain access to the receptor, thus facilitating its infection.

Chromosomal-scale genome assembly of the near-extinction big-head schizothorcin (Aspiorhynchus laticeps)

The big-head schizothorcin (Aspiorhynchus laticeps) is an endemic and near-extinction freshwater fish in Xinjiang, China. In this study, a chromosome-scale genome assembly of A. laticeps was generated using PacBio and Hi-C techniques. The PacBio sequencing data resulted in a 1.58 Gb assembly with a contig N50 of 1.27 Mb. Using Hi-C scaffolding approach, 88.38% of the initial assembled sequences were anchored and oriented into a chromosomal-scale assembly. The final assembly consisted of 25 pseudo-chromosomes that yielded 1.37 Gb of sequence, with a scaffold N50 of 44.02 Mb. BUSCO analysis showed a completeness score of 93.7%. The genome contained 48,537 predicted protein-coding genes and 58.31% of the assembly was annotated as repetitive sequences. Whole genome duplication events were further confirmed using 4dTv analysis. The genome assembly of A. laticeps should be valuable and important to understand the genetic adaptation and endangerment process of this species, which could lead to more effective management and conservation of the big-head schizothorcin and related freshwater fish species.

Measurement(s)	Genome Assembly Sequence
Technology Type(s)	Next Generation Sequencing
Sample Characteristic - Organism	Aspiorhynchus laticeps
Sample Characteristic - Location	Xinjiang Autonomous Region

Multiple roles of wheat calmodulin genes during stress treatment and TaCAM2-D as a positive regulator in response to drought and salt tolerance

Calmodulin (CaM) and calmodulin-like (CML) proteins are the most prominent calcium (Ca²⁺) sensing proteins involved in Ca²⁺-signaling processes. However, the function of these calcium sensors in wheat remains unclear. In this study, 15 TaCAMs and 113 TaCMLs were identified from the wheat reference genome. The analysis of cis-acting elements and expression patterns showed that TaCAMs might play an important role in response to abiotic and biotic stresses. TaCAM2-D gene was found to be significantly upregulated under drought and salt stresses, and thus, it was selected to further explore the biological function. Moreover, TaCAM2-D was observed to be localized in the nucleus, membrane and cytoplasm. Overexpression of TaCAM2-D in Arabidopsis conferred greater tolerance to drought and salt. The prediction analysis, the yeast two-hybrid analysis, and bimolecular fluorescence complementation assay indicated that TaCAM2-D interacted with TaMPK8, which is one of the wheat mitogen-activated protein kinases. Thus, the current study provides insights into the understanding of the TaCAM and TaCML genes in wheat.

Association of CARD14 Single-Nucleotide Polymorphisms with Psoriasis

Psoriasis is an immune-mediated chronic and painful disease characterized by red raised patches of inflamed skin that may have desquamation, silvery-white scales, itching and cracks. The susceptibility of developing psoriasis depends on multiple factors, with a complex interplay between genetic and environmental factors. Studies have suggested an association between autosomal dominant CARD14 (caspase recruitment domain-containing protein 14) gain-of-function mutations with the pathophysiology of psoriasis. In this study, non-synonymous single-nucleotide polymorphisms (nsSNPs) of CARD14 gene were assessed to determine their association with psoriasis in Pakistani population. A total of 123 subjects (63 patients with psoriasis and 60 normal controls) were included in this study. DNA was extracted from blood, and PCR analysis was performed followed by Sanger sequencing for 18 CARD14 specific nsSNPs (14 previously reported and the 4 most pathogenic nsSNPs identified using bioinformatics analysis). Among the 18 tested SNPs, only 2 nsSNP, rs2066965 (R547S) and rs34367357 (V585I), were found to be associated with psoriasis. Furthermore, rs2066965 heterozygous genotype was found to be more prevalent in patients with joint pain. Additionally, the 3D structure of CARD14 protein was predicted using alpha-fold2. NMSim web server was used to perform coarse grind simulations of wild-type CARD14 and two mutated structures. R547S increases protein flexibility, whereas V353I is shown to promote CARD14-induced NF-kappa B activation. This study confirms the association between two CARD14 nsSNPs, rs2066965 and rs34367357 with psoriasis in a Pakistani population, and could be helpful in identifying the role of CARD14 gene variants as potential genetic markers in patients with psoriasis.

R gene triplication confers European fodder turnip with improved clubroot resistance

Clubroot is one of the most important diseases for many important cruciferous vegetables and oilseed crops worldwide. Different clubroot resistance (CR) loci have been identified from only limited species in Brassica, making it difficult to compare and utilize these loci. European fodder turnip ECD04 is considered one of the most valuable resources for CR breeding. To explore the genetic and evolutionary basis of CR in ECD04, we sequenced the genome of ECD04 using de novo assembly and identified 978 candidate R genes. Subsequently, the 28 published CR loci were physically mapped to 15 loci in the ECD04 genome, including 62 candidate CR genes. Among them, two CR genes, CRA3.7.1 and CRA8.2.4, were functionally validated. Phylogenetic analysis revealed that CRA3.7.1 and CRA8.2.4 originated from a common ancestor before the whole-genome triplication (WGT) event. In clubroot susceptible Brassica species, CR-gene homologues were affected by transposable element (TE) insertion, resulting in the loss of CR function. It can be concluded that the current functional CR genes in Brassica rapa and non-functional CR genes in other Brassica species were derived from a common ancestral gene before WGT. Finally, a hypothesis for CR gene evolution is proposed for further discussion.

Genome Analyses of Two Blueberry Pathogens: Diaporthe amygdali CAA958 and Diaporthe eres CBS 160.32

The genus Diaporthe includes pathogenic species distributed worldwide and affecting a wide variety of hosts. Diaporthe amygdali and Diaporthe eres have been found to cause cankers, dieback, or twig blights on economically important crops such as soybean, almond, grapevine, and blueberry. Despite their importance as plant pathogens, the strategies of species of Diaporthe to infect host plants are poorly explored. To provide a genomic basis of pathogenicity, the genomes of D. amygdali CAA958 and D. eres CBS 160.32 were sequenced and analyzed. Cellular transporters involved in the transport of toxins, ions, sugars, effectors, and genes implicated in pathogenicity were detected in both genomes. Hydrolases and oxidoreductases were the most prevalent carbohydrate-active enzymes (CAZymes). However, analyses of the secreted proteins revealed that the secretome of D. eres CBS 160.32 is represented by 5.4% of CAZymes, whereas the secreted CAZymes repertoire of D. amygdali CAA958 represents 29.1% of all secretomes. Biosynthetic gene clusters (BGCs) encoding compounds related to phytotoxins and mycotoxins were detected in D. eres and D. amygdali genomes. The core gene clusters of the phytotoxin Fusicoccin A in D. amygdali are reported here through a genome-scale assembly. Comparative analyses of the genomes from 11 Diaporthe species revealed an average of 874 CAZymes, 101 secondary metabolite BGCs, 1640 secreted proteins per species, and genome sizes ranging from 51.5 to 63.6 Mbp. This study offers insights into the overall features and characteristics of Diaporthe genomes. Our findings enrich the knowledge about D. eres and D. amygdali, which will facilitate further research into the pathogenicity mechanisms of these species.

Comparative proteomics reveals biochemical changes in Salvia miltiorrhiza Bunge during sweating processing

ETHNOPHARMACOLOGICAL RELEVANCE

Salvia miltiorrhiza Bunge is a bulk medicinal material used in traditional Chinese medicine, that can cure cardiovascular diseases, neurasthenia, and other conditions. Sweating is a frequently used method of processing S. miltiorrhiza for medical applications. We previously demonstrated changes to the metabolic profile of linoleic acid, glyoxylate, and dicarboxylate after Sweating. However, this alteration has not been explained at the molecular level.

MATERIALS AND METHODS

Fresh roots of Salvia miltiorrhiza Bunge were treated by the Sweating processing, and then the tandem mass tag technique was used to compare the proteome difference between Sweating S. miltiorrhiza and non-Sweating S. miltiorrhiza.

RESULTS

We identified a total of 850 differentially expressed proteins after Sweating treatment in S. miltiorrhiza, including 529 upregulated proteins and 321 downregulated proteins. GO enrichment analysis indicated that these differentially expressed proteins are involved in external encapsulating structure, cell wall, oxidoreductase activity, ligase activity, and others. Further analysis showed that CYP450, the pathogenesis-related protein Bet v 1 family, and the peroxidase domain were the major protein domains. KEGG enrichment identified 18 pathways, of which phenylpropanoid biosynthesis is the most important one related to the metabolite profile and is the principal chemical component of S. miltiorrhiza.

CONCLUSION

This study addressed potential molecular mechanisms in S. miltiorrhiza after Sweating, and the findings provide reasons for the changes in biochemical properties and metabolites changes which might cause pharmacological variation at the proteome level.