Complete Mitochondrial Genome and Its Phylogenetic Position in Red Algae Fushitsunagia catenata from South Korea

The mitogenome is an important tool in taxonomic and evolutionary studies. Only a few complete mitogenomes have been reported for red algae. Herein, we reported the complete mitochondrial genome sequence of Fushitsunagia catenata (Harvey) Filloramo, G.V. and Saunders, G.W. 2016, a monospecific genus. The genome was 25,889 bp in circumference and had a strongly biased AT of 70.4%. It consisted of 2 rRNAs, 23 tRNAs, and 24 protein-coding genes (PCGs). nad5 (1986 bp) was the largest and atp9 (231 bp) was the smallest PCG. All PCGs used ATG as an initiation codon and TAA as a termination codon, except TAG, which was the termination codon used in the sdh3, rps3, and rps11 genes. The general structure and gene content of the present findings were almost identical to those of other red algae genomes, particularly those of the Rhodymeniales order. The maximum likelihood analysis showed that F. catenata was closely related to Rhodymenia pseudopalmata. The mitochondrial genome data presented in this study will enhance our understanding of evolution in Rhodophyta species.

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

EcoGenoRisk: Developing a computational ecological risk assessment tool for synthetic biology

The expanding field of synthetic biology (synbio) supports new opportunities in the design of targeted bioproducts or modified microorganisms. However, this rapid development of synbio products raises concerns surrounding the potential risks of modified microorganisms contaminating unintended environments. These potential invasion risks require new bioinformatic tools to inform the design phase. EcoGenoRisk is a newly constructed computational risk assessment tool for invasiveness that aims to predict where synbio microorganisms may establish a population by screening for habitats of genetically similar microorganisms. The first module of the tool identifies genetically similar microorganisms and potential ecological relationships such as competition, mutualism, and inhibition. In total, 520 archaeal and 32,828 bacterial complete assembly genomes were analyzed to test the specificity and accuracy of the tool as well as to characterize the enzymatic profiles of different taxonomic lineages. Additionally, ecological relationships were analyzed to determine which would result in the greatest potential overlap between shared functional profiles. Notably, competition displayed the significantly highest overlap of shared functions between compared genomes. Overall, EcoGenoRisk is a flexible software pipeline that assists environmental risk assessors to query large databases of known microorganisms and prioritize follow-up bench scale studies.

DandD: Efficient measurement of sequence growth and similarity

Genome assembly databases are growing rapidly. The redundancy of sequence content between a new assembly and previous ones is neither conceptually nor algorithmically easy to measure. We introduce pertinent methods and DandD, a tool addressing how much new sequence is gained when a sequence collection grows. DandD can describe how much structural variation is discovered in each new human genome assembly and when discoveries will level off in the future. DandD uses a measure called δ ("delta"), developed initially for data compression and chiefly dependent on k-mer counts. DandD rapidly estimates δ using genomic sketches. We propose δ as an alternative to k-mer-specific cardinalities when computing the Jaccard coefficient, thereby avoiding the pitfalls of a poor choice of k. We demonstrate the utility of DandD's functions for estimating δ, characterizing the rate of pangenome growth, and computing all-pairs similarities using k-independent Jaccard.

Comparative analysis of the vaginal bacteriome and virome in healthy women living in high-altitude and sea-level areas

The vaginal microbiota plays an important role in the health of the female reproductive tract and is closely associated with various pregnancy outcomes and sexually transmitted diseases. Plenty of internal and external factors have strong influence on the changes in a woman's vaginal microbiome. However, the effect of a high-altitude on female vaginal microbiota has not been described. In this study, we characterized the vaginal bacteriome and virome of 13 and 34 healthy women living in high-altitude and sea-level areas, using whole-metagenome shotgun sequencing of their vaginal mucus samples. The results revealed that the vaginal bacteriomes of high-altitude individuals are featured by a significant increase of species diversity, depletion of Lactobacillus crispatus, and more abundant of some anaerobic bacteria, such as Chlamydia trachomatis, Mageeibacillus indolicus, Dialister micraerophilus, and Sneathia amnii). In addition, the vagina samples of sea-level subjects harbor more Lactobacillus strains, whereas the anaerobic bacteroidetes strains mostly appeared in high-altitude subjects. Identified and assembled 191 virus operational taxonomic units (vOTUs), there were significant differences in the abundance of 107 vOTUs between the two groups. Together, the results of this study raised the understanding of bacteriome and virome in the vagina of women at different elevations, and demonstrated that the vaginal microbiome is related to the high-altitude geographic adaptation.

Draft genome of Brasenia schreberi, a worldwide distributed and endangered aquatic plant

OBJECTIVES

Brasenia is a monotypic genus in the family of Cabombaceae. The only species, B. schreberi, is a macrophyte distributed worldwide. Because it requires good water quality, it is endangered in China and other countries due to the deterioration of aquatic habitats. The young leaves and stems of B. schreberi are covered by thick mucilage, which has high medical value. As an allelopathic aquatic plant, it can also be used in the management of aquatic weeds. Here, we present its assembled and annotated genome to help shed light on medial and allelopathic substrates and facilitate their conservation.

DATA DESCRIPTION

Genomic DNA and RNA extracted from B. schreberi leaf tissues were used for whole genome and RNA sequencing using a Nanopore and/or MGI sequencer. The assembly was 1,055,148,839 bp in length, with 92 contigs and an N50 of 22,379,495 bp. The repetitive elements in the assembly were 555,442,205 bp. A completeness assessment of the assembly with BUSCO and compleasm indicated 88.4 and 90.9% completeness in the Eudicots database and 95.4 and 96.6% completeness in the Embryphyta database. Gene annotation revealed 67,747 genes that coded for 73,344 proteins.

Genomic Insights into High-Altitude Adaptation: A Comparative Analysis of Roscoea alpina and R. purpurea in the Himalayas

Environmental stress at high altitudes drives the development of distinct adaptive mechanisms in plants. However, studies exploring the genetic adaptive mechanisms of high-altitude plant species are scarce. In the present study, we explored the high-altitude adaptive mechanisms of plants in the Himalayas through whole-genome resequencing. We studied two widespread members of the Himalayan endemic alpine genus Roscoea (Zingiberaceae): R. alpina (a selfing species) and R. purpurea (an outcrossing species). These species are distributed widely in the Himalayas with distinct non-overlapping altitude distributions; R. alpina is distributed at higher elevations, and R. purpurea occurs at lower elevations. Compared to R. purpurea, R. alpina exhibited higher levels of linkage disequilibrium, Tajima's D, and inbreeding coefficient, as well as lower recombination rates and genetic diversity. Approximately 96.3% of the genes in the reference genome underwent significant genetic divergence (FST ≥ 0.25). We reported 58 completely divergent genes (FST = 1), of which only 17 genes were annotated with specific functions. The functions of these genes were primarily related to adapting to the specific characteristics of high-altitude environments. Our findings provide novel insights into how evolutionary innovations promote the adaptation of mountain alpine species to high altitudes and harsh habitats.

Chromosome-level genome assembly and demographic history of Euryodendron excelsum in monotypic genus endemic to China

Euryodendron excelsum is in a monotypic genus Euryodendron, endemic to China. It has intermediate morphisms in the Pentaphylacaceae or Theaceae families, which make it distinct. Due to anthropogenic disturbance, E. excelsum is currently found in very restricted and fragmented areas with extremely small populations. Although much research and effort has been applied towards its conservation, its long-term survival mechanisms and evolutionary history remain elusive, especially from a genomic aspect. Therefore, using a combination of long/short whole genome sequencing, RNA sequencing reads, and Hi-C data, we assembled and annotated a high-quality genome for E. excelsum. The genome assembly of E. excelsum comprised 1,059,895,887 bp with 99.66% anchored into 23 pseudo-chromosomes and a 99.0% BUSCO completeness. Comparative genomic analysis revealed the expansion of terpenoid and flavonoid secondary metabolite genes, and displayed a tandem and/or proximal duplication framework of these genes. E. excelsum also displayed genes associated with growth, development, and defence adaptation from whole genome duplication. Demographic analysis indicated that its fluctuations in population size and its recent population decline were related to cold climate changes. The E. excelsum genome assembly provides a highly valuable resource for evolutionary and ecological research in the future, aiding its conservation, management, and restoration.

A chromosome-level genome assembly of anesthetic drug-producing Anisodus acutangulus provides insights into its evolution and the biosynthesis of tropane alkaloids

Tropane alkaloids (TAs), which are anticholinergic agents, are an essential class of natural compounds, and there is a growing demand for TAs with anesthetic, analgesic, and spasmolytic effects. Anisodus acutangulus (Solanaceae) is a TA-producing plant that was used as an anesthetic in ancient China. In this study, we assembled a high-quality, chromosome-scale genome of A. acutangulus with a contig N50 of 7.4 Mb. A recent whole-genome duplication occurred in A. acutangulus after its divergence from other Solanaceae species, which resulted in the duplication of ADC1 and UGT genes involved in TA biosynthesis. The catalytic activities of H6H enzymes were determined for three Solanaceae plants. On the basis of evolution and co-expressed genes, AaWRKY11 was selected for further analyses, which revealed that its encoded transcription factor promotes TA biosynthesis by activating AaH6H1 expression. These findings provide useful insights into genome evolution related to TA biosynthesis and have potential implications for genetic manipulation of TA-producing plants.

Comprehensive analysis of complete chloroplast genome sequence of Plantago asiatica L. (Plantaginaceae)

Plantago asiatica L. is a representative individual species of Plantaginaceae, whose high reputation is owed to its edible and medicinal values. However, the phylogeny and genes of the P. asiatica chloroplast have not yet been well described. Here we report the findings of a comprehensive analysis of the P. asiatica chloroplast genome. The P. asiatica chloroplast genome is 164,992 bp, circular, and has a GC content of 37.98%. The circular genome contains 141 genes, including 8 rRNAs, 38 tRNAs, and 95 protein-coding genes. Seventy-two simple sequence repeats are detected. Comparative chloroplast genome analysis of six related species suggests that a higher similarity exists in the coding region than the non-coding region, and differences in the degree of preservation is smaller between P. asiatica and Plantago depressa than among others. Our phylogenetic analysis illustrates P. asiatica has a relatively close relationship with P. depressa, which was also divided into different clades with Plantago ovata and Plantago lagopus in the genus Plantago. This analysis of the P. asiatica chloroplast genome contributes to an improved deeply understanding of the evolutionary relationships among Plantaginaceae.

Draft genome of Castanopsis chinensis, a dominant species safeguarding biodiversity in subtropical broadleaved evergreen forests

OBJECTIVES

Castanopsis is the third largest genus in the Fagaceae family and is essentially tropical or subtropical in origin. The species in this genus are mainly canopy-dominant trees, and the key components of evergreen broadleaved forests play a crucial role in the maintenance of local biodiversity. Castanopsis chinensis, distributed from South China to Vietnam, is a representative species. It currently suffers from a high disturbance of human activity and climate change. Here, we present its assembled genome to facilitate its preliminary conservation and breeding on the genome level.

DATA DESCRIPTION

The C. chinensis genome was assembled and annotated by Nanopore and MGI whole-genome sequencing and RNA-seq reads using leaf tissues. The assembly was 888,699,661 bp in length, consisting of 133 contigs and a contig N50 of 23,395,510 bp. A completeness assessment of the assembly with Benchmarking Universal Single-Copy Orthologs (BUSCO) indicated a score of 98.3%. Repetitive elements comprised 471,006,885 bp, accounting for 55.9% of the assembled sequences. A total of 51,406 genes that coded for 54,310 proteins were predicted. Multiple databases were used to functionally annotate the protein sequences.

Dominance in self-compatibility between subgenomes of allopolyploid Arabidopsis kamchatica shown by transgenic restoration of self-incompatibility

The evolutionary transition to self-compatibility facilitates polyploid speciation. In Arabidopsis relatives, the self-incompatibility system is characterized by epigenetic dominance modifiers, among which small RNAs suppress the expression of a recessive SCR/SP11 haplogroup. Although the contribution of dominance to polyploid self-compatibility is speculated, little functional evidence has been reported. Here we employ transgenic techniques to the allotetraploid plant A. kamchatica. We find that when the dominant SCR-B is repaired by removing a transposable element insertion, self-incompatibility is restored. This suggests that SCR was responsible for the evolution of self-compatibility. By contrast, the reconstruction of recessive SCR-D cannot restore self-incompatibility. These data indicate that the insertion in SCR-B conferred dominant self-compatibility to A. kamchatica. Dominant self-compatibility supports the prediction that dominant mutations increasing selfing rate can pass through Haldane's sieve against recessive mutations. The dominance regulation between subgenomes inherited from progenitors contrasts with previous studies on novel epigenetic mutations at polyploidization termed genome shock.

Genome assembly of Erythrophleum Fordii, a special "ironwood" tree in China

OBJECTIVES

Erythrophleum is a genus in the Fabaceae family. The genus contains only about 10 species, and it is best known for its hardwood and medical properties worldwide. Erythrophleum fordii Oliv. is the only species of this genus distributed in China. It has superior wood and can be used in folk medicine, which leads to its overexploitation in the wild. For its effective conservation and elucidation of the distinctive genetic traits of wood formation and medical components, we present its first genome assembly.

DATA DESCRIPTION

This work generated ~ 160.8 Gb raw Nanopore whole genome sequencing (WGS) long reads, ~ 126.0 Gb raw MGI WGS short reads and ~ 29.0 Gb raw RNA-seq reads using E. fordii leaf tissues. The de novo assembly contained 864,825,911 bp in the E. fordii genome, with 59 contigs and a contig N50 of 30,830,834 bp. Benchmarking Universal Single-Copy Orthologs (BUSCO) revealed 98.7% completeness of the assembly. The assembly contained 471,006,885 bp (54.4%) repetitive sequences and 28,761 genes that coded for 33,803 proteins. The protein sequences were functionally annotated against multiple databases, facilitating comparative genomic analysis.

De novo assembly and comparative genome analysis for polyhydroxyalkanoates-producing Bacillus sp. BNPI-92 strain

BACKGROUND

Certain Bacillus species play a vital role in polyhydroxyalkanoate (PHA) production. However, most of these isolates did not properly identify to species level when scientifically had been reported.

RESULTS

From NGS analysis, 5719 genes were predicted in the de novo genome assembly. Based on genome annotation using RAST server, 5,527,513 bp sequences were predicted with 5679 bp number of protein-coding sequence. Its genome sequence contains 35.1% and 156 GC content and contigs, respectively. In RAST server analysis, subsystem (43%) and non-subsystem coverage (57%) were generated. Ortho Venn comparative genome analysis indicated that Bacillus sp. BNPI-92 shared 2930 gene cluster (core gene) with B. cereus ATCC 14579 T (AE016877), B. paranthracis Mn5T (MACE01000012), B. thuringiensis ATCC 10792 T (ACNF01000156), and B. antrics Amen T (AE016879) strains. For our strain, the maximum gene cluster (190) was shared with B. cereus ATCC 14579 T (AE016877). For Ortho Venn pair wise analysis, the maximum overlapping gene clusters thresholds have been detected between Bacillus s p.BNPI-92 and Ba. cereus ATCC 14579 T (5414). Average nucleotide identity (ANI) such as OriginalANI and OrthoANI, in silicon digital DND-DNA hybridization (isDDH), Type (Strain) Genome Server (TYGS), and Genome-Genome Distance Calculator (GGDC) were more essentially related Bacillus sp. BNPI-92 with B. cereus ATCC 14579 T strain. Therefore, based on the combination of RAST annotation, OrthoVenn server, ANI and isDDH result Bacillus sp.BNPI-92 strain was strongly confirmed to be a B. cereus type strain. It was designated as B. cereus BNPI-92 strain. In B. cereus BNPI-92 strain whole genome sequence, PHA biosynthesis encoding genes such as phaP, phaQ, phaR (PHA synthesis repressor phaR gene sequence), phaB/phbB, and phaC were predicted on the same operon. These gene clusters were designated as phaPQRBC. However, phaA was located on other operons.

CONCLUSIONS

This newly obtained isolate was found to be new a strain based on comparative genomic analysis and it was also observed as a potential candidate for PHA biosynthesis.

Draft genome sequence of methicillin-resistant Staphylococcus aureus strain D1418m22 isolated from human wound pus

This study reveals the complete genome sequence of methicillin-resistant Staphylococcus aureus (MRSA) strain d1418m22, sourced from a Karnal, Haryana, human skin wound. Classified as community-associated MRSA, it features a 2.78-MB genome harboring Staphylococcus-specific genetic elements, encompassing 2,625 protein-coding genes and 18 antimicrobial resistance genes.

De novo genome assembly and annotation of the medicinal plant Tinospora cordifolia (Willd.) Miers ex Hook. f. & Thom's

Indian natural climbing shrub Tinospora cordifolia, often known as "Guduchi" and "Amrita," is a highly esteemed medicinal plant in the Indian system of medicine (ISM). It is a member of the Menispermaceae family which consists of a rich source of protein, micronutrients, and rich source of bioactive components which are used in treating various systemic diseases. The current study was designed to know the biological characterization of the plant genome and biosynthesis of plant metabolites essential for its medicinal applications. Tinospora cordifolia's complete genome was sequenced using Illumina HiSeq2500 sequencing technology. The draft genome was assembled through a de novo method. An integrative genome annotation approach was used to perform functional gene prediction. The pathway analysis was carried out using the KEGG database. The total genome size obtained after genome assembly was 894 Mb with an N50 of 9148 bp. The integrative annotation approach resulted in 35,111 protein-coding genes. In addition, genes responsible for the synthesis of syringin, a secondary metabolite found in plants, were identified. In comparison to the standard drug (dopamine, rasagiline, and selegiline), syringin's molecular docking exhibited a greater binding affinity from the range of - 4.3 to - 6.6 kcal/mol for all the targets of Parkinson's disease and for Alzheimer's targets; it has shown the maximum potency from the range of - 6.5 to - 7.4 kcal/mol with respect to the standard drug (donepezil, galantamine, and rivastigmine). This study provides the genomic information of Tinospora cordifolia which is helpful in understanding genomic insights and metabolic pathways connected to the corresponding plant genome and predicts the possible useful effect for the molecular characterization of therapeutic drugs.

Construction and characterization of microsatellite markers for the Schistosoma japonicum isolate from a hilly area of China based on whole genome sequencing

Schistosoma japonicum had once caused the greatest disease burden in China and has still been transmitted in some hilly areas, for example, in Shitai of Anhui province, where rodents are projected to be the main reservoir. This may lead to a critical need of molecular tools with high efficiency in monitoring the dynamic of the rodent-associated S. japonicum, as an appropriate amount of schistosome input can re-establish its life cycle in a place with snails and then result in the re-emergence of schistosomiasis. Therefore, the goal of this study was to develop high polymorphic microsatellites from the whole genome of rodent-associated S. japonicum strain to monitor its transmission dynamic. We sampled the hilly schistosome isolate from Shitai of Anhui in China and sequenced the parasite with the next-generation sequencing technology. The whole genome was assembled with four different approaches. We then developed 71 microsatellite markers at a genome-wide scale throughout two best assembled genomes. Based on their chromosome mapping and the expected length of targeted sequences, we selected 24 markers for the development of multiplex reactions. Two multiplexes composed of 10 loci were finally developed, and their potential was revealed by their successful application on and capturing the genetic diversity of three schistosome populations. The selected 10 markers, each with clear chromosome location and characteristics, will be greatly useful in tracing the dispersal pathways or/and dynamics of the rodent-associated S. japonicum or others in the hilly area of China or elsewhere.

Selection of mutant Listeria phages under food-relevant conditions can enhance application potential

Bacteriophages are viruses that infect and kill bacteria. Currently, phage products are available for the control of the pathogen Listeria monocytogenes in food products in the United States. In this study, we explore whether experimental evolution can be used to generate phages with improved abilities to function under specific food-relevant conditions. Ultra-pasteurized oat and whole milk were chosen as test matrices as they represent different food groups, yet have similar physical traits and macronutrient composition. We showed that (i) wild-type phage LP-125 infection kinetics are different in the two matrices and (ii) LP-125 has a significantly higher burst size in oat milk. From this, we attempted to evolve LP-125 to have improved infection kinetics in whole milk. Ancestral LP-125 was passaged through 10 rounds of amplification in milk conditions. Plaque-purified DNA samples from milk-selected phages were isolated and sequenced, and mutations present in the isolated phages were identified. We found two nonsynonymous substitutions in LP125_108 and LP125_112 genes, which encode putative baseplate-associated glycerophosphoryl diester phosphodiesterase and baseplate protein, respectively. Protein structural modeling showed that the substituted amino acids in the mutant phages are predicted to localize to surface-exposed helices on the corresponding structures, which might affect the surface charge of proteins and their interaction with the bacterial cell. The phage containing the LP125_112 mutation adsorbed significantly faster than the ancestral phage in both oat and whole milk. Follow-up experiments suggest that fat content may be a key factor for the expression of the phenotype of this mutation. IMPORTANCE Bacteriophages are one of the tools available to control the foodborne pathogen, Listeria monocytogenes. Phage products must work under a broad range of food conditions to be an effective control for L. monocytogenes. Here, we show that the experimental evolution of phages can be used to generate new phages with phenotypes useful under specific conditions. We used this approach to select for a mutant phage that more efficiently binds to L. monocytogenes that is grown in whole milk and oat milk. We show that the fat content of these milks is necessary for the expression of this phenotype. Our findings show that experimental evolution can be used to select for improved phages with better performance under specific conditions. This approach has the potential to support the development of condition-specific phage-based biocontrols in the food industry.

Exploring objective feature sets in constructing the evolution relationship of animal genome sequences

BACKGROUND

Exploring evolution regularities of genome sequences and constructing more objective species evolution relationships at the genomic level are high-profile topics. Based on the evolution mechanism of genome sequences proposed in our previous research, we found that only the 8-mers containing CG or TA dinucleotides correlate directly with the evolution of genome sequences, and the relative frequency rather than the actual frequency of these 8-mers is more suitable to characterize the evolution of genome sequences.

RESULT

Therefore, two types of feature sets were obtained, they are the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers. The evolution relationships of mammals and reptiles were constructed by the relative frequency set of CG1 + CG2 8-mers, and two types of evolution relationships of insects were constructed by the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers respectively. Through comparison and analysis, we found that evolution relationships are consistent with the known conclusions. According to the evolution mechanism, we considered that the evolution relationship constructed by CG1 + CG2 8-mers reflects the evolution state of genome sequences in current time, and the evolution relationship constructed by TA1 + TA2 8-mers reflects the evolution state in the early stage.

CONCLUSION

Our study provides objective feature sets in constructing evolution relationships at the genomic level.

A chromosome-level reference genome of an aromatic medicinal plant Adenosma buchneroides

Adenosma buchneroides Bonati, belonging to the genus Adenosum (Plantaginaceae), is an aromatic medicinal plant and utilized in traditional Chinese medicine. It has been widely used as plant-based repellents to prevent vector-borne diseases. However, the lack of a reference genome limits the study of conservation management and molecular biology of A. buchneroides. Here, we generated a chromosome-level de novo genome assembly of A. buchneroides which is a high-quality chromosome-scale assembly of aromatic medicinal plant in Plantaginaceae. The genome has a total length of 442.84 Mb with scaffold N50 of 27.98 Mb and 95.55% of the genome assigned to 14 chromosomes. BUSCO assessment yielded a completeness score of 97.2%. Furthermore, we predicted 24,367 protein-coding genes, and 95.79% of them was functionally annotated. The chromosome-scale genome of A. buchneroides will be a significant resource for understanding the genetic basis and evolution of active components biosynthesis, which will facilitate further study and exploit of A. buchneroides.

Genomic and ecological factors shaping specialism and generalism across an entire subphylum

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Paradigms proposed to explain this variation either invoke trade-offs between performance efficiency and breadth or underlying intrinsic or extrinsic factors. We assembled genomic (1,154 yeast strains from 1,049 species), metabolic (quantitative measures of growth of 843 species in 24 conditions), and ecological (environmental ontology of 1,088 species) data from nearly all known species of the ancient fungal subphylum Saccharomycotina to examine niche breadth evolution. We found large interspecific differences in carbon breadth stem from intrinsic differences in genes encoding specific metabolic pathways but no evidence of trade-offs and a limited role of extrinsic ecological factors. These comprehensive data argue that intrinsic factors driving microbial niche breadth variation.

Genome Survey and Chromosome-Level Draft Genome Assembly of Glycine max var. Dongfudou 3: Insights into Genome Characteristics and Protein Deficiencies

Dongfudou 3 is a highly sought-after soybean variety due to its lack of beany flavor. To support molecular breeding efforts, we conducted a genomic survey using next-generation sequencing. We determined the genome size, complexity, and characteristics of Dongfudou 3. Furthermore, we constructed a chromosome-level draft genome and speculated on the molecular basis of protein deficiency in GmLOX1, GmLOX2, and GmLOX3. These findings set the stage for high-quality genome analysis using third-generation sequencing. The estimated genome size is approximately 1.07 Gb, with repetitive sequences accounting for 72.50%. The genome is homozygous and devoid of microbial contamination. The draft genome consists of 916.00 Mb anchored onto 20 chromosomes, with annotations of 46,446 genes and 77,391 transcripts, achieving Benchmarking Single-Copy Orthologue (BUSCO) completeness of 99.5% for genome completeness and 99.1% for annotation. Deletions and substitutions were identified in the three GmLox genes, and they also lack corresponding active proteins. Our proposed approach, involving k-mer analysis after filtering out organellar DNA sequences, is applicable to genome surveys of all plant species, allowing for accurate assessments of size and complexity. Moreover, the process of constructing chromosome-level draft genomes using closely related reference genomes offers cost-effective access to valuable information, maximizing data utilization.

Abundance and Diversification of Repetitive Elements in Decapoda Genomes

Repetitive elements are a major component of DNA sequences due to their ability to propagate through the genome. Characterization of Metazoan repetitive profiles is improving; however, current pipelines fail to identify a significant proportion of divergent repeats in non-model organisms. The Decapoda order, for which repeat content analyses are largely lacking, is characterized by extremely variable genome sizes that suggest an important presence of repetitive elements. Here, we developed a new standardized pipeline to annotate repetitive elements in non-model organisms, which we applied to twenty Decapoda and six other Crustacea genomes. Using this new tool, we identified 10% more repetitive elements than standard pipelines. Repetitive elements were more abundant in Decapoda species than in other Crustacea, with a very large number of highly repeated satellite DNA families. Moreover, we demonstrated a high correlation between assembly size and transposable elements and different repeat dynamics between Dendrobranchiata and Reptantia. The patterns of repetitive elements largely reflect the phylogenetic relationships of Decapoda and the distinct evolutionary trajectories within Crustacea. In summary, our results highlight the impact of repetitive elements on genome evolution in Decapoda and the value of our novel annotation pipeline, which will provide a baseline for future comparative analyses.

Seedability: optimizing alignment parameters for sensitive sequence comparison

Motivation

Most sequence alignment techniques make use of exact k-mer hits, called seeds, as anchors to optimize alignment speed. A large number of bioinformatics tools employing seed-based alignment techniques, such as Minimap2 , use a single value of k per sequencing technology, without a strong guarantee that this is the best possible value. Given the ubiquity of sequence alignment, identifying values of k that lead to more sensitive alignments is thus an important task. To aid this, we present Seedability , a seed-based alignment framework designed for estimating an optimal seed k-mer length (as well as a minimal number of shared seeds) based on a given alignment identity threshold. In particular, we were motivated to make Minimap2 more sensitive in the pairwise alignment of short sequences.

Results

The experimental results herein show improved alignments of short and divergent sequences when using the parameter values determined by Seedability in comparison to the default values of Minimap2 . We also show several cases of pairs of real divergent sequences, where the default parameter values of Minimap2 yield no output alignments, but the values output by Seedability produce plausible alignments.

Availability and implementation

https://github.com/lorrainea/Seedability (distributed under GPL v3.0).

Complete Mitogenome Sequencing, Annotation, and Phylogeny of Grateloupia turuturu, a Red Alga with Intronic cox1 Gene

The mitochondrial genome (mitogenome) is essential for identifying species and tracing genetic variation, gene patterns, and evolutionary studies. Here, the mitogenome of Grateloupia turuturu was sequenced on the Illumina sequencing platform. This circular mitogenome (28,265 bp) contains 49 genes, including three rRNAs, twenty transfer RNAs (tRNAs), and twenty-six protein-coding genes (PCGs). Nucleotide composition indicates biased AT (68.8%) content. A Group II intronic sequence was identified between two exons of the cox1 gene, and this sequence comprises an open reading frame (ORF) that encodes a hypothetical protein. The gene content, annotation, and genetic makeup are identical to those of Halymeniaceae members. The complete mitogenome sequences of the Grateloupia and Polyopes species were used in a phylogenetic analysis, which revealed that these two genera are monophyletic and that G. turuturu and G. elliptica are closely related. This newly constructed mitogenome will help us better understand the general trends in the development of cox1 introns in Halymeniaceae, as well as the evolution of red algal mitogenomes within the Rhodophyta and among diverse algal species.

Seeding with minimized subsequence

MOTIVATION

Modern methods for computation-intensive tasks in sequence analysis (e.g. read mapping, sequence alignment, genome assembly, etc.) often first transform each sequence into a list of short, regular-length seeds so that compact data structures and efficient algorithms can be employed to handle the ever-growing large-scale data. Seeding methods using kmers (substrings of length k) have gained tremendous success in processing sequencing data with low mutation/error rates. However, they are much less effective for sequencing data with high error rates as kmers cannot tolerate errors.

RESULTS

We propose SubseqHash, a strategy that uses subsequences, rather than substrings, as seeds. Formally, SubseqHash maps a string of length n to its smallest subsequence of length k, k < n, according to a given order overall length-k strings. Finding the smallest subsequence of a string by enumeration is impractical as the number of subsequences grows exponentially. To overcome this barrier, we propose a novel algorithmic framework that consists of a specifically designed order (termed ABC order) and an algorithm that computes the minimized subsequence under an ABC order in polynomial time. We first show that the ABC order exhibits the desired property and the probability of hash collision using the ABC order is close to the Jaccard index. We then show that SubseqHash overwhelmingly outperforms the substring-based seeding methods in producing high-quality seed-matches for three critical applications: read mapping, sequence alignment, and overlap detection. SubseqHash presents a major algorithmic breakthrough for tackling the high error rates and we expect it to be widely adapted for long-reads analysis.

AVAILABILITY AND IMPLEMENTATION

SubseqHash is freely available at https://github.com/Shao-Group/subseqhash.

The genome of Acorus deciphers insights into early monocot evolution

Acorales is the sister lineage to all the other extant monocot plants. Genomic resource enhancement of this genus can help to reveal early monocot genomic architecture and evolution. Here, we assemble the genome of Acorus gramineus and reveal that it has ~45% fewer genes than the majority of monocots, although they have similar genome size. Phylogenetic analyses based on both chloroplast and nuclear genes consistently support that A. gramineus is the sister to the remaining monocots. In addition, we assemble a 2.2 Mb mitochondrial genome and observe many genes exhibit higher mutation rates than that of most angiosperms, which could be the reason leading to the controversies of nuclear genes- and mitochondrial genes-based phylogenetic trees existing in the literature. Further, Acorales did not experience tau (τ) whole-genome duplication, unlike majority of monocot clades, and no large-scale gene expansion is observed. Moreover, we identify gene contractions and expansions likely linking to plant architecture, stress resistance, light harvesting, and essential oil metabolism. These findings shed light on the evolution of early monocots and genomic footprints of wetland plant adaptations.

Dataset from genome sequencing, assembly and mining of microsatellite markers in barred-button quail (Turnix suscitator)

Turnix suscitator (barred-button quail) is a member of the primitive genus Turnix in the highly diverse order of shore birds Charadriiformes. Absence of genome scale data of T. suscitator has limited our understanding about its systematics, taxonomic and evolutionary history as well has hindered the characterization of genome wide microsatellite markers of the same. Hence we generated whole genome short read sequences of T. suscitator, created a high quality assembly and mined genome-wide microsatellite markers from the same. A total of 34142524 reads were sequenced with an estimated genome size of 817 mb. SPAdes assembly consisted of 320761 total contigs and an estimated N50 value of 907 base pairs. Krait identified a total of 77028 microsatellite motifs covering 0.64% of the total sequences in the SPAdes assembly. Further the whole genome sequence and genome wide microsatellites dataset of T. suscitator will facilitate future genomic/evolutionary studies of Turnix species.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates change at the macroevolutionary level is a major question in evolutionary biology. With >66,000 described species, rove beetles (Staphylinidae) comprise the largest metazoan family. Their exceptional radiation has been coupled to pervasive biosynthetic innovation whereby numerous lineages bear defensive glands with diverse chemistries. Here, we combine comparative genomic and single-cell transcriptomic data from across the largest rove beetle clade, Aleocharinae. We retrace the functional evolution of two novel secretory cell types that together comprise the tergal gland-a putative catalyst behind Aleocharinae's megadiversity. We identify key genomic contingencies that were critical to the assembly of each cell type and their organ-level partnership in manufacturing the beetle's defensive secretion. This process hinged on evolving a mechanism for regulated production of noxious benzoquinones that appears convergent with plant toxin release systems, and synthesis of an effective benzoquinone solvent that weaponized the total secretion. We show that this cooperative biosynthetic system arose at the Jurassic-Cretaceous boundary, and that following its establishment, both cell types underwent ∼150 million years of stasis, their chemistry and core molecular architecture maintained almost clade-wide as Aleocharinae radiated globally into tens of thousands of lineages. Despite this deep conservation, we show that the two cell types have acted as substrates for the emergence of adaptive, biochemical novelties-most dramatically in symbiotic lineages that have infiltrated social insect colonies and produce host behavior-manipulating secretions. Our findings uncover genomic and cell type evolutionary processes underlying the origin, functional conservation and evolvability of a chemical innovation in beetles.

Photoreceptor genes in a trechine beetle, Trechiama kuznetsovi, living in the upper hypogean zone

To address how organisms adapt to a new environment, subterranean organisms whose ancestors colonized subterranean habitats from surface habitats have been studied. Photoreception abilities have been shown to have degenerated in organisms living in caves and calcrete aquifers. Meanwhile, the organisms living in a shallow subterranean environment, which are inferred to reflect an intermediate stage in an evolutionary pathway to colonization of a deeper subterranean environment, have not been studied well. In the present study, we examined the photoreception ability in a trechine beetle, Trechiama kuznetsovi, which inhabits the upper hypogean zone and has a vestigial compound eye. By de novo assembly of genome and transcript sequences, we were able to identify photoreceptor genes and phototransduction genes. Specifically, we focused on opsin genes, where one long wavelength opsin gene and one ultraviolet opsin gene were identified. The encoded amino acid sequences had neither a premature stop codon nor a frameshift mutation, and appeared to be subject to purifying selection. Subsequently, we examined the internal structure of the compound eye and nerve tissue in the adult head, and found potential photoreceptor cells in the compound eye and nerve bundle connected to the brain. The present findings suggest that T. kuznetsovi has retained the ability of photoreception. This species represents a transitional stage of vision, in which the compound eye regresses, but it may retain the ability of photoreception using the vestigial eye.

Haplotype resolved chromosome level genome assembly of Citrus australis reveals disease resistance and other citrus specific genes

Recent advances in genome sequencing and assembly techniques have made it possible to achieve chromosome level reference genomes for citrus. Relatively few genomes have been anchored at the chromosome level and/or are haplotype phased, with the available genomes of varying accuracy and completeness. We now report a phased high-quality chromosome level genome assembly for an Australian native citrus species; Citrus australis (round lime) using highly accurate PacBio HiFi long reads, complemented with Hi-C scaffolding. Hifiasm with Hi-C integrated assembly resulted in a 331 Mb genome of C. australis with two haplotypes of nine pseudochromosomes with an N50 of 36.3 Mb and 98.8% genome assembly completeness (BUSCO). Repeat analysis showed that more than 50% of the genome contained interspersed repeats. Among them, LTR elements were the predominant type (21.0%), of which LTR Gypsy (9.8%) and LTR copia (7.7%) elements were the most abundant repeats. A total of 29 464 genes and 32 009 transcripts were identified in the genome. Of these, 28 222 CDS (25 753 genes) had BLAST hits and 21 401 CDS (75.8%) were annotated with at least one GO term. Citrus specific genes for antimicrobial peptides, defense, volatile compounds and acidity regulation were identified. The synteny analysis showed conserved regions between the two haplotypes with some structural variations in Chromosomes 2, 4, 7 and 8. This chromosome scale, and haplotype resolved C. australis genome will facilitate the study of important genes for citrus breeding and will also allow the enhanced definition of the evolutionary relationships between wild and domesticated citrus species.

Molecular evaluation of the metabolism of estrogenic di(2-ethylhexyl) phthalate in Mycolicibacterium sp

BACKGROUND

Di(2-ethylhexyl) phthalate (DEHP) is a widely detected plasticizer and a priority pollutant of utmost concern for its adverse impact on humans, wildlife and the environment. To eliminate such toxic burden, biological processes are the most promising ways to combat rampant environmental insults under eco-friendly conditions. The present study investigated the biochemical and molecular assessment of the catabolic potential of Mycolicibacterium sp. strain MBM in the assimilation of estrogenic DEHP.

RESULTS

A detailed biochemical study revealed an initial hydrolytic pathway of degradation for DEHP followed by the assimilation of hydrolyzed phthalic acid and 2-ethylhexanol to TCA cycle intermediates. Besides the inducible nature of DEHP-catabolic enzymes, strain MBM can efficiently utilize various low- and high-molecular-weight phthalate diesters and can grow under moderately halotolerant conditions. Whole genome sequence analysis exhibited a genome size of 6.2 Mb with a GC content of 66.51% containing 6,878 coding sequences, including multiple genes, annotated as relevant to the catabolism of phthalic acid esters (PAEs). Substantiating the annotated genes through transcriptome assessment followed by RT-qPCR analysis, the possible roles of upregulated genes/gene clusters in the metabolism of DEHP were revealed, reinforcing the biochemical pathway of degradation at the molecular level.

CONCLUSIONS

A detailed co-relation of biochemical, genomic, transcriptomic and RT-qPCR analyses highlights the PAE-degrading catabolic machineries in strain MBM. Further, due to functional attributes in the salinity range of both freshwater and seawater, strain MBM may find use as a suitable candidate in the bioremediation of PAEs.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

Draft genome and transcriptome of Nepenthes mirabilis, a carnivorous plant in China

OBJECTIVES

Nepenthes belongs to the monotypic family Nepenthaceae, one of the largest carnivorous plant families. Nepenthes species show impressive adaptive radiation and suffer from being overexploited in nature. Nepenthes mirabilis is the most widely distributed species and the only Nepenthes species that is naturally distributed within China. Herein, we reported the genome and transcriptome assemblies of N. mirabilis. The assemblies will be useful resources for comparative genomics, to understand the adaptation and conservation of carnivorous species.

DATA DESCRIPTION

This work produced ~ 139.5 Gb N. mirabilis whole genome sequencing reads using leaf tissues, and ~ 21.7 Gb and ~ 27.9 Gb of raw RNA-seq reads for its leaves and flowers, respectively. Transcriptome assembly obtained 339,802 transcripts, in which 79,758 open reading frames (ORFs) were identified. Function analysis indicated that these ORFs were mainly associated with proteolysis and DNA integration. The assembled genome was 691,409,685 bp with 159,555 contigs/scaffolds and an N50 of 10,307 bp. The BUSCO assessment of the assembled genome and transcriptome indicated 91.1% and 93.7% completeness, respectively. A total of 42,961 genes were predicted in the genome identified, coding for 45,461 proteins. The predicted genes were annotated using multiple databases, facilitating future functional analyses of them. This is the first genome report on the Nepenthaceae family.

Exploring the sorghum race level diversity utilizing 272 sorghum accessions genomic resources

Due to evolutionary divergence, sorghum race populations exhibit significant genetic and morphological variation. A k-mer-based sorghum race sequence comparison identified the conserved k-mers of all 272 accessions from sorghum and the race-specific genetic signatures identified the gene variability in 10,321 genes (PAVs). To understand sorghum race structure, diversity and domestication, a deep learning-based variant calling approach was employed in a set of genotypic data derived from a diverse panel of 272 sorghum accessions. The data resulted in 1.7 million high-quality genome-wide SNPs and identified selective signature (both positive and negative) regions through a genome-wide scan with different (iHS and XP-EHH) statistical methods. We discovered 2,370 genes associated with selection signatures including 179 selective sweep regions distributed over 10 chromosomes. Co-localization of these regions undergoing selective pressure with previously reported QTLs and genes revealed that the signatures of selection could be related to the domestication of important agronomic traits such as biomass and plant height. The developed k-mer signatures will be useful in the future to identify the sorghum race and for trait and SNP markers for assisting in plant breeding programs.

Stoichiometric balance ratio of cellobiose and gentiobiose induces cellulase production in Talaromyces cellulolyticus

BACKGROUND

The exact mechanism by which fungal strains sense insoluble cellulose is unknown, but research points to the importance of transglycosylation products generated by fungi during cellulose breakdown. Here, we used multi-omics approach to identify the transglycosylation metabolites and determine their function in cellulase induction in a model strain, Talaromyces cellulolyticus MTCC25456.

RESULTS

Talaromyces sp. is a novel hypercellulolytic fungal strain. Based on genome scrutiny and biochemical analysis, we predicted the presence of cellulases on the surface of its spores. We performed metabolome analysis to show that these membrane-bound cellulases act on polysaccharides to form a mixture of disaccharides and their transglycosylated derivatives. Inevitably, a high correlation existed between metabolite data and the KEGG enrichment analysis of differentially expressed genes in the carbohydrate metabolic pathway. Analysis of the contribution of the transglycosylation product mixtures to cellulase induction revealed a 57% increase in total cellulase. Further research into the metabolites, using in vitro induction tests and response surface methodology, revealed that Talaromyces sp. produces cell wall-breaking enzymes in response to cellobiose and gentiobiose as a stimulant. Precisely, a 2.5:1 stoichiometric ratio of cellobiose to gentiobiose led to a 2.4-fold increase in cellulase synthesis. The application of the optimized inducers in cre knockout strain significantly increased the enzyme output.

CONCLUSION

This is the first study on the objective evaluation and enhancement of cellulase production using optimized inducers. Inducer identification and genetic engineering boosted the cellulase production in the cellulolytic fungus Talaromyces sp.

Introgressive hybridization in the west Pacific pen shells (genus Atrina): Restricted interspecies gene flow within the genome

A compelling interest in marine biology is to elucidate how species boundaries between sympatric free-spawning marine invertebrates such as bivalve molluscs are maintained in the face of potential hybridization. Hybrid zones provide the natural resources for us to study the underlying genetic mechanisms of reproductive isolation between hybridizing species. Against this backdrop, we examined the occurrence of introgressive hybridization (introgression) between two bivalves distributed in the western Pacific margin, Atrina japonica and Atrina lischkeana, based on single-nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing. Using 1066 ancestry-informative SNP sites, we also investigated the extent of introgression within the genome to search for SNP sites with reduced interspecies gene flow. A series of our individual-level clustering analyses including the principal component analysis, Bayesian model-based clustering, and triangle plotting based on ancestry-heterozygosity relationships for an admixed population sample from the Seto Inland Sea (Japan) consistently suggested the presence of specimens with varying degrees of genomic admixture, thereby implying that the two species are not completely isolated. The Bayesian genomic cline analysis identified 10 SNP sites with reduced introgression, each of which was located within a genic region or an intergenic region physically close to a functional gene. No, or very few, heterozygotes were observed at these sites in the hybrid zone, suggesting that selection acts against heterozygotes. Accordingly, we raised the possibility that the SNP sites are within genomic regions that are incompatible between the two species. Our finding of restricted interspecies gene flow at certain genomic regions gives new insight into the maintenance of species boundaries in hybridizing broadcast-spawning molluscs.

An improved reference genome for Trifolium subterraneum L. provides insight into molecular diversity and intra-specific phylogeny

Subterranean clover (Trifolium subterraneum L., Ts) is a geocarpic, self-fertile annual forage legume with a compact diploid genome (n = x = 8, 544 Mb/1C). Its resilience and climate adaptivity have made it an economically important species in Mediterranean and temperate zones. Using the cultivar Daliak, we generated higher resolution sequence data, created a new genome assembly TSUd_3.0, and conducted molecular diversity analysis for copy number variant (CNV) and single-nucleotide polymorphism (SNP) among 36 cultivars. TSUd_3.0 substantively improves prior genome assemblies with new Hi-C and long-read sequence data, covering 531 Mb, containing 41,979 annotated genes and generating a 94.4% BUSCO score. Comparative genomic analysis among select members of the tribe Trifolieae indicated TSUd 3.0 corrects six assembly-error inversion/duplications and confirmed phylogenetic relationships. Its synteny with T. pratense, T. repens, Medicago truncatula and Lotus japonicus genomes were assessed, with the more distantly related T. repens and M. truncatula showing higher levels of co-linearity with Ts than between Ts and its close relative T. pratense. Resequencing of 36 cultivars discovered 7,789,537 SNPs subsequently used for genomic diversity assessment and sequence-based clustering. Heterozygosity estimates ranged from 1% to 21% within the 36 cultivars and may be influenced by admixture. Phylogenetic analysis supported subspecific genetic structure, although it indicates four or five groups, rather than the three recognized subspecies. Furthermore, there were incidences where cultivars characterized as belonging to a particular subspecies clustered with another subspecies when using genomic data. These outcomes suggest that further investigation of Ts sub-specific classification using molecular and morpho-physiological data is needed to clarify these relationships. This upgraded reference genome, complemented with comprehensive sequence diversity analysis of 36 cultivars, provides a platform for future gene functional analysis of key traits, and genome-based breeding strategies for climate adaptation and agronomic performance. Pangenome analysis, more in-depth intra-specific phylogenomic analysis using the Ts core collection, and functional genetic and genomic studies are needed to further augment knowledge of Trifolium genomes.

Chromosome-level genomes of multicellular algal sisters to land plants illuminate signaling network evolution

The filamentous and unicellular algae of the class Zygnematophyceae are the closest algal relatives of land plants. Inferring the properties of the last common ancestor shared by these algae and land plants allows us to identify decisive traits that enabled the conquest of land by plants. We sequenced four genomes of filamentous Zygnematophyceae (three strains of Zygnema circumcarinatum and one strain of Z. cylindricum) and generated chromosome-scale assemblies for all strains of the emerging model system Z. circumcarinatum. Comparative genomic analyses reveal expanded genes for signaling cascades, environmental response, and intracellular trafficking that we associate with multicellularity. Gene family analyses suggest that Zygnematophyceae share all the major enzymes with land plants for cell wall polysaccharide synthesis, degradation, and modifications; most of the enzymes for cell wall innovations, especially for polysaccharide backbone synthesis, were gained more than 700 million years ago. In Zygnematophyceae, these enzyme families expanded, forming co-expressed modules. Transcriptomic profiling of over 19 growth conditions combined with co-expression network analyses uncover cohorts of genes that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

The chromosome-level genome and key genes associated with mud-dwelling behavior and adaptations of hypoxia and noxious environments in loach (Misgurnus anguillicaudatus)

BACKGROUND

The loach (Misgurnus anguillicaudatus), the most widely distributed species of the family Cobitidae, displays a mud-dwelling behavior and intestinal air-breathing, inhabiting the muddy bottom of extensive freshwater habitats. However, lack of high-quality reference genome seriously limits the interpretation of the genetic basis of specialized adaptations of the loach to the adverse environments including but not limited to the extreme water temperature, hypoxic and noxious mud environment.

RESULTS

This study generated a 1.10-Gb high-quality, chromosome-anchored genome assembly, with a contig N50 of 3.83 Mb. Multiple comparative genomic analyses found that proto-oncogene c-Fos (fos), a regulator of bone development, is positively selected in loach. Knockout of fos (ID: Mis0086400.1) led to severe osteopetrosis and movement difficulties, combined with the comparison results of bone mineral density, supporting the hypothesis that fos is associated with loach mud-dwelling behavior. Based on genomic and transcriptomic analysis, we identified two key elements involved in the intestinal air-breathing of loach: a novel gene (ID: mis0158000.1) and heat shock protein beta-1 (hspb1). The flavin-containing monooxygenase 5 (fmo5) genes, central to xenobiotic metabolism, undergone expansion in loach and were identified as differentially expressed genes in a drug stress trial. A fmo5-/- (ID: Mis0185930.1) loach displayed liver and intestine injury, indicating the importance of this gene to the adaptation of the loach to the noxious mud.

CONCLUSIONS

Our work provides valuable insights into the genetic basis of biological adaptation to adverse environments.

Independent recruitment of FRUITFULL-like transcription factors in the convergent origins of vernalization-responsive grass flowering

Flowering in response to low temperatures (vernalization) has evolved multiple times independently across angiosperms as an adaptation to match reproductive development with the short growing season of temperate habitats. Despite the context of a generally conserved flowering time network, evidence suggests that the genes underlying vernalization responsiveness are distinct across major plant clades. Whether different or similar mechanisms underlie vernalization-induced flowering at narrower (e.g., family-level) phylogenetic scales is not well understood. To test the hypothesis that vernalization responsiveness has evolved convergently in temperate species of the grass family (Poaceae), we carried out flowering time experiments with and without vernalization in several representative species from different subfamilies. We then determined the likelihood that vernalization responsiveness evolved through parallel mechanisms by quantifying the response of Pooideae vernalization pathway FRUITFULL (FUL)-like genes to extended periods of cold. Our results demonstrate that vernalization-induced flowering has evolved multiple times independently in at least five grass subfamilies, and that different combinations of FUL-like genes have been recruited to this pathway on several occasions.

Characterization and genome-wide sequence analysis of an ectomycorrhizal fungus Pisolithus albus, a potential source for reclamation of degraded lands

Pisolithus albus is a ubiquitous ectomycorrhizal fungus that establishes symbiosis with a wide range of woody plants around the globe. The symbiotic association of this fungus plays a crucial role in the nutrient cycling of their host plants and enables them to thrive in adverse environmental conditions. Based on its ecological importance and lack of genomic studies, whole-genome sequencing was carried out to analyze P. albus sequences through an Illumina HiSeq X system. The functional annotations were performed against various databases to explore genomic patterns and traits possibly attributing to its specialization. Comparative genomics of P. albus with phylogenetically related Pisolithus microcarpus and Pisolithus tinctorius (only available genomes of Pisolithus at NCBI till now) led to the identification of their unique and shared basic functional and stress adaptation capabilities. The de novo assembled genome of 56.15 Mb with 91.8% BUSCO completeness is predicted to encode 23,035 genes. The study is aimed to generate solid genomic data resources for P. albus, forming the theoretical basis for future transcriptomic, proteomic and metabolomic studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03483-5.

SPRDA: a link prediction approach based on the structural perturbation to infer disease-associated Piwi-interacting RNAs

piRNA and PIWI proteins have been confirmed for disease diagnosis and treatment as novel biomarkers due to its abnormal expression in various cancers. However, the current research is not strong enough to further clarify the functions of piRNA in cancer and its underlying mechanism. Therefore, how to provide large-scale and serious piRNA candidates for biological research has grown up to be a pressing issue. In this study, a novel computational model based on the structural perturbation method is proposed to predict potential disease-associated piRNAs, called SPRDA. Notably, SPRDA belongs to positive-unlabeled learning, which is unaffected by negative examples in contrast to previous approaches. In the 5-fold cross-validation, SPRDA shows high performance on the benchmark dataset piRDisease, with an AUC of 0.9529. Furthermore, the predictive performance of SPRDA for 10 diseases shows the robustness of the proposed method. Overall, the proposed approach can provide unique insights into the pathogenesis of the disease and will advance the field of oncology diagnosis and treatment.

Genome Survey Sequencing of the Mole Cricket Gryllotalpa orientalis

The mole cricket Gryllotalpa orientalis is an evolutionarily, medicinal, and agriculturally significant insect that inhabits underground environments and is distributed globally. This study measured genome size by flow cytometry and k-mer based on low-coverage sequencing, and nuclear repetitive elements were also identified. The haploid genome size estimate is 3.14 Gb by flow cytometry, 3.17 Gb, and 3.77 Gb-based two k-mer methods, respectively, which is well within the range previously reported for other species of the suborder Ensifera. 56% of repetitive elements were found in G. orientalis, similar to 56.83% in Locusta migratoria. However, the great size of repetitive sequences could not be annotated to specific repeat element families. For the repetitive elements that were annotated, Class I-LINE retrotransposon elements were the most common families and more abundant than satellite and Class I-LTR. These results based on the newly developed genome survey could be used in the taxonomic study and whole genome sequencing to improve the understanding of the biology of G. orientalis.

Isolation and Characterization of Bacteria from Natural Hot Spring and Insights into the Thermophilic Cellulase Production

A thermophilic cellulase-producing bacterium, Bacillus velezensis strain MRC 5958, from Bakra natural hot Springs, India was characterized through genome sequencing. It has a genome size of 4,467,129 bp and a GC content of 45.7%. A cellulase purified from its fermentation broth has a molecular weight of about 18 kDa. The optimum temperature and pH for carboxymethyl cellulase activity were at 55 °C and pH ~ 7.0. The enzyme is stable over a wide range of temperatures from 30 °C to 70 °C with maximum activity observed at 48 h of incubation. The strain produces cellulase on alkali-treated sugarcane bagasse, rice straw, rice husk, rice bran, and sawdust. The sugarcane bagasse exhibited the most effective carbon source for cellulase production at (85 U/ml) followed by rice bran (68 U/ml), rice husk (60 U/ml), rice straw (48 U/ml), and sawdust (39 U/ml). Therefore, this strain can be a potential thermostable cellulase-producing candidate for converting the waste biomass into biofuel and other industrial enzymes.

Transformation of 2, 4, 6-trinitrotoluene by Stenotrophomonas strain SG1 under aerobic and anaerobic conditions

TNT, or 2,4,6-trinitrotoluene, is a common explosive that can contaminate soil and groundwater in production sites, military training areas, and disposal locations. The compound is highly toxic; therefore, there is an urgent need to rehabilitate the impacted environments. Harnessing the microbial ability to biodegrade TNT into environmentally harmless compound(s) is one approach to remediating contaminated sites. In our study, we report on the genomic and metabolic ability of Stenotrophomonas strain SG1 to degrade TNT under aerobic and anaerobic conditions. The bacterial strain SG1 was first isolated as a contaminant from a culture of Diaphorobacter sp. strain DS2 over minimal media supplemented with TNT. The draft genome assembly of strain SG1 is ∼4.7 Mb and is distributed among 358 contigs. The homology search against the custom database of enzymes responsible for TNT biodegradation revealed the presence of three N-ethylmaleimide reductases (NemA) with a defined KEGG ortholog and KEGG pathway of TNT degradation. The presence of respiratory nitrate reductases has also been mapped, which supports denitrification under anaerobic conditions. Experimentally, the TNT transformation rate accelerated when carbon sources, such as sodium acetate, sodium citrate, sodium succinate, sucrose, and glucose (final concentration of 5 mM), were added. Citrate promoted the highest growth and TNT transformation ratio (88.35%) in 120 h. With the addition of 5 mM ammonium chloride, TNT completely disappeared in the citrate and sucrose-containing treatments in 120 h. However, higher biomass was obtained in the sucrose and glucose-containing treatments in 120 h. During incubation, the formation of amino dinitrotoluene isomers, dinitrotoluene isomers, trinitrobenzene, azoxy isomers, diaryl hydroxylamines, and corresponding secondary amines was confirmed by GC/MS and UPLC/MS. 2-Amino-4-nitrotoluene, 4-amino-2-nitrotoluene, and 2-amino-6-nitrotoluene were also identified in the culture supernatant by GC/MS. Under anaerobic conditions, TNT completely disappeared in the citrate and citrate plus nitrate treatments. Since the strain shows the ability to remove TNT, this research should be useful in basic research and practical applications for removing TNT from wastewater.

K-Mer-Based Genome Size Estimation in Theory and Practice

Recent advances in sequencing technologies have made genome sequencing of non-model organisms with very large and complex genomes possible. The data can be used to estimate diverse genome characteristics, including genome size, repeat content, and levels of heterozygosity. K-mer analysis is a powerful biocomputational approach with a wide range of applications, including estimation of genome sizes. However, interpretation of the results is not always straightforward. Here, I review k-mer-based genome size estimation, focusing specifically on k-mer theory and peak calling in k-mer frequency histograms. I highlight common pitfalls in data analysis and result interpretation, and provide a comprehensive overview on current methods and programs developed to conduct these analyses.

Epigenetics and island-mainland divergence in an insectivorous small mammal

Geographically isolated populations, specifically island-mainland counterparts, tend to exhibit phenotypic variation in many species. The so-called island syndrome occurs when different environmental pressures lead to insular divergence from mainland populations. This phenomenon can be seen in an island population of Nova Scotia masked shrews (Sorex cinereus), which have developed a specialized feeding habit and digestive enzyme compared to their mainland counterparts. Epigenetic modifications, such as DNA methylation (DNAm), can impact phenotypes by altering gene expression without changing the DNA sequence. Here, we used a de novo masked shrew genome assembly and a mammalian methylation array profiling 37 thousand conserved CpGs to investigate morphological and DNA methylation patterns between island and mainland populations. Island shrews were morphologically and epigenetically different than their mainland counterparts, exhibiting a smaller body size. A gene ontology enrichment analyses of differentially methylated CpGs implicated developmental and digestive system related pathways. Based on our shrew epigenetic clock, island shrews might also be aging faster than their mainland counterparts. This study provides novel insight on phenotypic and epigenetic divergence in island-mainland mammal populations and suggests an underlying role of methylation in island-mainland divergence.

spa Typing of Methicillin-Resistant Staphylococcus aureus Based on Whole-Genome Sequencing: the Impact of the Assembler

Sequencing of the spa gene of methicillin-resistant Staphylococcus aureus (MRSA) is used for assigning spa types to e.g., detect transmission and control outbreaks. Traditionally, spa typing is performed by Sanger sequencing but has in recent years been replaced by whole-genome sequencing (WGS) in some laboratories. Spa typing by WGS involves de novo assembly of millions of short sequencing reads into larger contiguous sequences, from which the spa type is then determined. The choice of assembly program therefore potentially impacts the spa typing result. In this study, WGS of 1,754 MRSA isolates was followed by de novo assembly using the assembly programs SPAdes (with two different sets of parameters) and SKESA. The spa types were assigned and compared to the spa types obtained by Sanger sequencing, regarding the latter as the correct spa types. SPAdes with the two different settings resulted in assembly of the correct spa type for 84.8% and 97.6% of the isolates, respectively, while SKESA assembled the correct spa type in 98.6% of cases. The misassembled spa types were generally two spa repeats shorter than the correct spa type and mainly included spa types with repetition of the same repeats. WGS-based spa typing is thus very accurate compared to Sanger sequencing, when the best assembly program for this purpose is used. IMPORTANCE spa typing of methicillin-resistant Staphylococcus aureus (MRSA) is widely used by clinicians, infection control workers, and researchers both in local outbreak investigations and as an easy way to communicate and compare MRSA types between laboratories and countries. Traditionally, spa types are determined by Sanger sequencing, but in recent years a whole-genome sequencing (WGS)-based approach has become increasingly used. In this study, we compared spa typing by WGS using different methods for assembling the genome from short sequencing reads and compared to Sanger sequencing as the gold standard. We find substantial differences in correct assembly of spa types between the assembly methods. Our findings are therefore important for the quality of WGS based spa typing data being exchanged by clinical microbiology laboratories.

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.