Bridger: a new framework for de novo transcriptome assembly using RNA-seq data

The double-banded false coralsnake Erythrolamprus bizona (Dipsadidae, Xenodontinae, Xenodontini) has a metalloproteinase-rich venom with proteolytic activity towards azocasein and α-fibrinogen

The venom of the double-banded false coralsnake, Erythrolamprus bizona, is proteolytic and attenuates neuromuscular contractile activity in vitro. Here, we examined the Duvernoy's venom gland histology and general composition of E. bizona venom using a combination of chromatographic, electrophoretic, enzymatic and proteomic analyses. Histologically, the venom gland consisted of serous epithelium-lined secretory tubules and a supralabial gland that stained positively for mucopolysaccharide. SDS-PAGE showed that the venom had a simple composition, with proteins in the range of 15-60 kDa. This simple composition was confirmed by RP-HPLC that revealed 15 main protein peaks. The venom (1-10 μg) was highly proteolytic towards azocasein, but was devoid of esterase, phospholipase (PLA2), and L-amino acid oxidase activities. The venom also degraded casein and gelatin in zymographic assays, with activity towards gelatin being particularly potent and detected over the range of 18.7 ng-30 μg of venom; gelatinolytic activity was also detected in four of the RP-HPLC peaks. The venom (10 μg) selectively degraded the α-chain of fibrinogen. All proteolytic activity was inhibited by EDTA (metalloproteinase inhibitor) but not by AEBSF (serine proteinase inhibitor). SDS-PAGE followed by in-gel digestion of the main electrophoretic bands coupled with LC-MS/MS analysis revealed the presence of five toxin families: C-type lectin-like proteins (CTL), cysteine-rich secretory proteins (CRiSP), phospholipase B (PLB), snake venom matrix metalloproteinases (svMMP), and snake venom metalloproteinases (SVMP). These findings extend our knowledge of the toxinology of E. bizona and suggest that the local manifestations (pain, edema, erythema, and ecchymosis) seen in human envenomation by this species are probably mediated by venom metalloproteinases.

VirDiG: a de novo transcriptome assembler for coronavirus

Motivation

The discontinuous transcription mechanism of coronaviruses contributes to their adaptation to different host environments and plays a critical role in their lifecycle. Accurate assembly of coronavirus transcripts is vital for understanding the virus's biological traits and developing precise prevention and treatment strategies. However, existing de novo assembly algorithms are primarily designed for alternative splicing events in eukaryotes and are not suitable for assembling coronavirus transcriptome, which consists of both genomic RNA and subgenomic mRNAs. Coronavirus transcriptome reconstruction from short reads remains a challenging problem.

Results

In this work, we present VirDiG, a de novo transcriptome assembler specifically designed for coronaviruses. VirDiG utilizes a discontinuous graph to facilitate accurate transcript assembly by incorporating information from paired-end reads, sequence depth, and start and stop codons. Experimental results from both simulated and real datasets show that VirDiG exhibits significant advantages in reconstructing the transcriptome of coronaviruses when compared to traditional de novo assemblers tailored for classical eukaryotic transcriptome assembly.

Availability and implementation

VirDiG is freely available at https://github.com/Limh616/VirDiG.git.

Chromosome-level genome assembly of the threatened ornamental plant Hibiscus yunnanensis

Hibiscus yunnanensis S.Y. Hu is an endangered species of the genus Hibiscus (Malvaceae), which has high potential economic value. However, the absence of a high-quality reference genome impedes the study of the ecology and molecular biology of H. yunnanensis. Here, we present a high-quality chromosome-level assembly of H. yunnanensis using BGI-DIPSEQ, Nanopore, and Hi-C sequencing. The assembled genome size is 2.2 Gb with a contig N50 of 12.1 Mb and a scaffold N50 of 137.1 Mb. Approximately 99.2% of the assembly is anchored into 17 pseudochromosomes, and a BUSCO analysis indicates a completeness score of 99.6%. Furthermore, we identify 42,085 protein-coding genes, of which 96.4% are functionally annotated. This genome resource provides a foundation for future studies on unique traits, including drought-tolerant, savanna-adapted, and long-flowering traits. Its ability to flower in winter, along with its automatic selfing and lack of delayed inbreeding depression, makes it an excellent model for studying style curvature mechanism and its adaptive significance in the Malvaceae.

The removal of high Se(IV) and Cd(II) concentrations in sulfur autotrophic reactor based on the "hibernation-like microbial survival strategy"

The removal of selenite (Se(IV)) and cadmium (Cd(II)) from low-carbon wastewater presents significant challenges. However, the addition of external organic carbon sources is limited in application due to the high cost and potential for secondary pollution. This study introduced a "hibernation-like microbial survival strategy", enabling efficient removal of Se(IV) and Cd(II) in sulfur autotrophic reactor, with S0 acting as the electron donor. The removal efficiencies of 5-120 mg/L Se(IV) and 50 mg/L Cd(II) were higher than 99 % in phase I-IV, and the nanoparticles formed in sulfur autotrophic reactor were available for recycling. The analysis of X-ray photoelectron spectroscopy confirmed that the removal pathways of Se(IV) and Cd(II) were biological reduction, adsorption, and biosynthesis. The decreased ratio of actual to theoretical sulfate concentrations indicated the weakened sulfur disproportionation trend in sulfur autotrophic reactor. The formation of autotrophic-heterotrophic symbiont was beneficial for promoting electron transfer, material exchange, and information flow. Microorganisms strategically decreased metabolic activity to reduce extra energy consumption under Se(IV) and Cd(II) stress, which was manifested in the decreased extracellular DNA, extracellular polymeric substances, and electron transfer system activity. Furthermore, microorganisms reduced the secretion of nicotinamide adenine dinucleotide, cytochrome c, and cyt-c oxidase on the premise of ensuring the required electron flux. The "hibernation-like microbial survival strategy" was proposed to explain the removal of Se(IV) and Cd(II) in sulfur autotrophic reactor, expanding the potential application of sulfur autotrophy in environmental engineering.

Cov-trans: an efficient algorithm for discontinuous transcript assembly in coronaviruses

BACKGROUND

Discontinuous transcription allows coronaviruses to efficiently replicate and transmit within host cells, enhancing their adaptability and survival. Assembling viral transcripts is crucial for virology research and the development of antiviral strategies. However, traditional transcript assembly methods primarily designed for variable alternative splicing events in eukaryotes are not suitable for the viral transcript assembly problem. The current algorithms designed for assembling viral transcripts often struggle with low accuracy in determining the transcript boundaries. There is an urgent need to develop a highly accurate viral transcript assembly algorithm.

RESULTS

In this work, we propose Cov-trans, a reference-based transcript assembler specifically tailored for the discontinuous transcription of coronaviruses. Cov-trans first identifies canonical transcripts based on discontinuous transcription mechanisms, start and stop codons, as well as reads alignment information. Subsequently, it formulates the assembly of non-canonical transcripts as a path extraction problem, and introduces a mixed integer linear programming to recover these non-canonical transcripts.

CONCLUSION

Experimental results show that Cov-trans outperforms other assemblers in both accuracy and recall, with a notable strength in accurately identifying the boundaries of transcripts. Cov-trans is freely available at https://github.com/computer-Bioinfo/Cov-trans.git .

Comprehensive Analysis of the Influence of Technical and Biological Variations on De Novo Assembly of RNA-Seq Datasets

De novo assembly of transcriptomes from species without reference genome remains a common problem in functional genomics. While methods and algorithms for transcriptome assembly are continually being developed and published, the quality of de novo assemblies using short reads depends on the complexity of the transcriptome and is limited by several types of errors. One problem to overcome is the research gap regarding the best method to use in each study to obtain high-quality de novo assembly. Currently, there are no established protocols for solving the assembly problem considering the transcriptome complexity. In addition, the accuracy of quality metrics used to evaluate assemblies remains unclear. In this study, we investigate and discuss how different variables accounting for the complexity of RNA-Seq data influence assembly results independently of the software used. For this purpose, we simulated transcriptomic short-read sequence datasets from high-quality full-length predicted transcript models with varying degrees of complexity. Subsequently, we conducted de novo assemblies using different assembly programs, and compared and classified the results using both reference-dependent and independent metrics. These metrics were assessed both individually and combined through multivariate analysis. The degree of alternative splicing and the fragment size of the paired-end reads were identified as the variables with the greatest influence on the assembly results. Moreover, read length and fragment size had different influences on the reconstruction of longer and shorter transcripts. These results underscore the importance of understanding the composition of the transcriptome under study, and making experimental design decisions related to the need to work with reads and fragments of different sizes. In addition, the choice of assembly software will positively impact the final assembly outcome. This selection will affect the completeness of represented genes and assembled isoforms, as well as contribute to error reduction.

Predicting Promoters in Multiple Prokaryotes with Prompt

Promoters are important cis-regulatory elements for the regulation of gene expression, and their accurate predictions are crucial for elucidating the biological functions and potential mechanisms of genes. Many previous prokaryotic promoter prediction methods are encouraging in terms of the prediction performance, but most of them focus on the recognition of promoters in only one or a few bacterial species. Moreover, due to ignoring the promoter sequence motifs, the interpretability of predictions with existing methods is limited. In this work, we present a generalized method Prompt (Promoters in multiple prokaryotes) to predict promoters in 16 prokaryotes and improve the interpretability of prediction results. Prompt integrates three methods including RSK (Regression based on Selected k-mer), CL (Contrastive Learning) and MLP (Multilayer Perception), and employs a voting strategy to divide the datasets into high-confidence and low-confidence categories. Results on the promoter prediction tasks in 16 prokaryotes show that the accuracy (Accuracy, Matthews correlation coefficient) of Prompt is greater than 80% in highly credible datasets of 16 prokaryotes, and is greater than 90% in 12 prokaryotes, and Prompt performs the best compared with other existing methods. Moreover, by identifying promoter sequence motifs, Prompt can improve the interpretability of the predictions. Prompt is freely available at https://github.com/duqimeng/PromptPrompt , and will contribute to the research of promoters in prokaryote.

Chemical-defensome and whole-transcriptome expression of the silverside fish Basilichthys microlepidotus in response to chronic pollution in the Maipo River basin, Central Chile

Pollution is a major global concern affecting biodiversity, particularly of freshwater species. Populations have developed mechanisms to deal with pollution, such as the chemical defensome, which is a set of genes involved in maintaining internal stability. Pollution significantly affects the Maipo River basin in Chile. This area is home to the endemic silverside fish Basilichthys microlepidotus, whose populations are affected by pollution to different degrees. We assessed gene expression in the liver and gill of this species, focusing on whole-transcriptome and chemical-defensome levels, to identify both independent and shared mechanisms in response to pollution. The results showed that 14-18 genes were consistently expressed differently among populations in polluted areas. These genes were primarily involved in liver cell mitosis and in responses to organic chemicals and carcinogenic processes. Genes expressed differently in the gill were more abundant in immune system biological processes. All populations consistently downregulated chemical-defensome genes in the liver. In differentially expressed chemical-defensome genes, shared biological processes included virus response, cellular redox homeostasis and transport, organic cyclic compound response and DNA-templated transcription regulation. Studying chemical-defensome genes can help reveal common ways that pollution builds up over time, and examining the whole transcriptome can elucidate the context in which this response develops.

Differential expression analyses and detection of SNP loci associated with environmental variables: Are salinity and temperature factors involved in population differentiation and speciation in Odontesthes?

Environmental factors play a key role in individual adaptation to different local conditions. Because of this, studies about the physiological and genetic responses of individuals exposed to different natural environments offer clues about mechanisms involved in population differentiation, and as a subsequent result, speciation. Marine environments are especially suited to survey this kind of phenomena because they commonly harbor species adapted to different local conditions along a geographic continuum. Silversides belonging to Odontesthes are commonly distributed in tropical and temperate regions of South America and exhibit noticeable phenotypic plasticity, which allows them to adapt to contrasting environments. In this study, the genetic expression of O. argentinensis sampled along the Uruguayan Atlantic coast and estuarine adjacent areas was investigated. In addition, the correlation between individual genotypes and environmental variables was also analysed in O. argentinensis and O. bonariensis. Results obtained suggest a differential expression pattern of low magnitude among individuals from the different areas sampled and a correlation between several SNP loci and environmental variables. The analyses carried out did not show a clear differentiation among individuals sampled along different salinity regimens, but enriched GOTerms seem to be driven by water oxygen content. On the other hand, a total of 46 SNPs analysed in O. argentinensis and O. bonariensis showed a correlation with salinity and temperature. Although none of the correlated SNPs and corresponding genes from our both analyses were directly associated with hypoxia, genes related to the cardiovascular system and muscle cell differentiation were found. All these genes are interesting candidates for future studies since they are closely related to the differentially expressed genes. Although salinity was also mentioned as an important parameter limiting introgression between O. argentinensis and O. bonariensis, it was found that salinity does not drive differential expression in O. argentinensis, but rather oxygen levels. Moreover, salinity does not directly affect the structure and genetic divergence of the populations, they appear to be structured based on their degree of isolation and geographical distance between them. Further studies, like genome-wide analyses, could help to elucidate additional genes adapted to the different environments in these silverside species.

Chromosome-level assembly for the complex genome of land hermit crab Coenobita brevimanus

Land hermit crabs are a group of shell-carrying crabs that have evolved remarkable terrestrial adaptations in behavior, morphology, physiology, and biochemistry. However, the genetic mechanisms underlying these adaptations remain unclear. In addition, usually it is very difficult to get good genome assemblies for crustaceans. In this study, we managed to assemble the first chromosome-level genome for a land hermit crab (Coenobita brevimanus) with careful manual curation. The final assembly spans 4.74 Gb, with the contig N50 of 1.75 Mb and scaffold N50 of 42.95 Mb, encompassing 117 chromosomes that account for 96.54% of the genome. The evaluations including genome BUSCO (95.26%), Merqury qv (35.88) and the mapping ratio of pair-end short reads (99.48%) showed the high-continuity of C. brevimanus genome assembly, making it the genome with the highest quality in crustaceans with genome size bigger than 3 Gb. The availability of this chromosome-scale genome of crustaceans represents a valuable resource for the land hermit crab, which represents an independent water-to-land transition evolutionary event in the animal kingdom.

Type 2 vomeronasal receptor expression in the olfactory organ of African lungfish, Protopterus annectens

The olfactory organ of tetrapods, with few exceptions, comprises the main and accessory organs: olfactory epithelium (OE) and vomeronasal organ (VNO). Unlike tetrapods, teleost fish lack a VNO. However, lungfish, a type of sarcopterygian fish closely related to tetrapods, possesses a lamellar OE similar to the OE of teleosts and a recess epithelium (RecE) resembling the amphibian VNO. The RecE has been hypothesized as a primordial VNO. Olfactory receptors in tetrapods are distinctively expressed in the OE and VNO. For instance, type 2 vomeronasal receptors (V2Rs) in Xenopus are categorized into those exclusively expressed in the OE and those solely expressed in the VNO. It remains unclear whether V2Rs are differentially expressed between the lamellar OE and RecE in lungfish. This study investigated V2R expression in the lamellar OE and RecE of the African lungfish, Protopterus annectens. P. annectens V2Rs were categorized into three groups: those exclusively expressed in the lamellar OE, those exclusively expressed in the RecE, and those expressed in both the lamellar OE and RecE. V2Rs exclusively expressed in the RecE and those expressed in both the lamellar OE and RecE formed a distinct clade in the phylogenetic tree, whereas others were solely expressed in the lamellar OE. These findings suggest that lungfish V2R expression represents an intermediate stage toward complete segregation between V2Rs expressed in the OE and those expressed in the VNO.

All tangled up: Unraveling phylogenetics and reticulate evolution in the vining ferns, Lygodium (Schizaeales)

PREMISE

Reticulate evolution, often accompanied by polyploidy, is prevalent in plants, and particularly in the ferns. Resolving the resulting non-bifurcating histories remains a major challenge for plant phylogenetics. Here, we present a phylogenomic investigation into the complex evolutionary history of the vining ferns, Lygodium (Lygodiaceae, Schizaeales).

METHODS

Using a targeted enrichment approach with the GoFlag 408 flagellate land plant probe set, we generated large nuclear and plastid sequence datasets for nearly all taxa in the genus and constructed the most comprehensive phylogeny of the family to date using concatenated maximum likelihood and coalescence approaches. We integrated this phylogeny with cytological and spore data to explore karyotype evolution and generate hypotheses about the origins of putative polyploids and hybrids.

RESULTS

Our data and analyses support the origins of several putative allopolyploids (e.g., L. cubense, L. heterodoxum) and hybrids (e.g., L. ×fayae) and also highlight the potential prevalence of autopolyploidy in this clade (e.g., L. articulatum, L. flexuosum, and L. longifolium).

CONCLUSIONS

Our robust phylogenetic framework provides valuable insights into dynamic reticulate evolution in this clade and demonstrates the utility of target-capture data for resolving these complex relationships.

Phylogenomics recovers multiple origins of portable case making in caddisflies (Insecta: Trichoptera), nature's underwater architects

Caddisflies (Trichoptera) are among the most diverse groups of freshwater animals with more than 16 000 described species. They play a fundamental role in freshwater ecology and environmental engineering in streams, rivers and lakes. Because of this, they are frequently used as indicator organisms in biomonitoring programmes. Despite their importance, key questions concerning the evolutionary history of caddisflies, such as the timing and origin of larval case making, remain unanswered owing to the lack of a well-resolved phylogeny. Here, we estimated a phylogenetic tree using a combination of transcriptomes and targeted enrichment data for 207 species, representing 48 of 52 extant families and 174 genera. We calibrated and dated the tree with 33 carefully selected fossils. The first caddisflies originated approximately 295 million years ago in the Permian, and major suborders began to diversify in the Triassic. Furthermore, we show that portable case making evolved in three separate lineages, and shifts in diversification occurred in concert with key evolutionary innovations beyond case making.

The chromosome-level genomes of the herbal magnoliids Warburgia ugandensis and Saururus chinensis

Warburgia ugandensis and Saururus chinensis are two of the most important medicinal plants in magnoliids and are widely utilized in traditional Kenya and Chinese medicine, respectively. The absence of higher-quality reference genomes has hindered research on the medicinal compound biosynthesis mechanisms of these plants. We report the chromosome-level genome assemblies of W. ugandensis and S. chinensis, and generated 1.13 Gb and 0.53 Gb genomes from 74 and 27 scaffolds, respectively, using BGI-DIPSEQ, Nanopore, and Hi-C sequencing. The scaffold N50 lengths were 82.97 Mb and 48.53 Mb, and the assemblies were anchored to 14 and 11 chromosomes of W. ugandensis and S. chinensis, respectively. In total, 24,739 and 20,561 genes were annotated, and 98.5% and 98% of the BUSCO genes were fully represented, respectively. The chromosome-level genomes of W. ugandensis and S. chinensis will be valuable resources for understanding the genetics of these medicinal plants, studying the evolution of magnoliids and angiosperms and conserving plant genetic resources.

Comprehensive phylogeny of Pieridae butterflies reveals strong correlation between diversification and temperature

Temperature is thought to be a key factor influencing global species richness patterns. We investigate the link between temperature and diversification in the butterfly family Pieridae by combining next generation DNA sequences and published molecular data with fine-grained distribution data. We sampled nearly 600 pierid butterfly species to infer the most comprehensive molecular phylogeny of the family and curated a distribution dataset of more than 800,000 occurrences. We found strong evidence that species in environments with more stable daily temperatures or cooler maximum temperatures in the warm seasons have higher speciation rates. Furthermore, speciation and extinction rates decreased in tandem with global temperatures through geological time, resulting in a constant net diversification.

Cloning and characterization of luciferase from an Asian firefly Pygoluciola qingyu and its comparison with other beetle luciferases

The bioluminescence system of luminescent beetles has extensive applications in biological imaging, protein labeling and drug screening. To explore wild luciferases with excellent catalytic activity and thermal stability, we cloned the luciferase of Pygoluciola qingyu, one species living in areas of high temperature and with strong bioluminescence, by combining transcriptomic sequencing and reverse transcription polymerase chain reaction (RT-PCR). The total length of luciferase gene is 1638 bp and the luciferase consists 544 amino acids. The recombinant P. qingyu luciferase was produced in vitro and its characteristics were compared with those of eight luciferases from China firefly species and two commercial luciferases. Compared with these luciferases, the P. qingyu luciferase shows the highest luminescence activity at room temperature (about 25-28 ℃) with similar KM value for D-luciferin and ATP to the Photinus pyralis luciferase. The P. qingyu luciferase activity was highest at 35 ℃ and can keep high activity at 30-40 ℃, which suggests the potential of P. qingyu luciferase for in vivo and cell application. Our results provide new insights into P. qingyu luciferase and give a new resource for the application of luciferases.

Venom Composition of Neglected Bothropoid Snakes from the Amazon Rainforest: Ecological and Toxinological Implications

Snake venoms have evolved in several families of Caenophidae, and their toxins have been assumed to be biochemical weapons with a role as a trophic adaptation. However, it remains unclear how venom contributes to the success of venomous species for adaptation to different environments. Here we compared the venoms from Bothrocophias hyoprora, Bothrops taeniatus, Bothrops bilineatus smaragdinus, Bothrops brazili, and Bothrops atrox collected in the Amazon Rainforest, aiming to understand the ecological and toxinological consequences of venom composition. Transcriptomic and proteomic analyses indicated that the venoms presented the same toxin groups characteristic from bothropoids, but with distinct isoforms with variable qualitative and quantitative abundances, contributing to distinct enzymatic and toxic effects. Despite the particularities of each venom, commercial Bothrops antivenom recognized the venom components and neutralized the lethality of all species. No clear features could be observed between venoms from arboreal and terrestrial habitats, nor in the dispersion of the species throughout the Amazon habitats, supporting the notion that venom composition may not shape the ecological or toxinological characteristics of these snake species and that other factors influence their foraging or dispersal in different ecological niches.

A vertebrate-wide catalogue of T1R receptors reveals diversity in taste perception

Taste is a vital chemical sense for feeding behaviour. In mammals, the umami and sweet taste receptors comprise three members of the taste receptor type 1 (T1R/TAS1R) family: T1R1, T1R2 and T1R3. Because their functional homologues exist in teleosts, only three TAS1R genes generated by gene duplication are believed to have been inherited from the common ancestor of bony vertebrates. Here, we report five previously uncharacterized TAS1R members in vertebrates, TAS1R4, TAS1R5, TAS1R6, TAS1R7 and TAS1R8, based on genome-wide survey of diverse taxa. We show that mammalian and teleost fish TAS1R2 and TAS1R3 genes are paralogues. Our phylogenetic analysis suggests that the bony vertebrate ancestor had nine TAS1Rs resulting from multiple gene duplications. Some TAS1Rs were lost independently in descendent lineages resulting in retention of only three TAS1Rs in mammals and teleosts. Combining functional assays and expression analysis of non-teleost fishes we show that the novel T1Rs form heterodimers in taste-receptor cells and recognize a broad range of ligands such as essential amino acids, including branched-chain amino acids, which have not been previously considered as T1R ligands. This study reveals diversity of taste sensations in both modern vertebrates and their ancestors, which might have enabled vertebrates to adapt to diverse habitats on Earth.

A reference transcriptome for walnut anthracnose pathogen, Ophiognomonia leptostyla, guides the discovery of candidate virulence genes

Despite the economic losses due to the walnut anthracnose, Ophiognomonia leptostyla is an orphan fungus with respect to genomic resources. In the present study, the transcriptome of O. leptostyla was assembled for the first time. RNA sequencing was conducted for the fungal mycelia grown in a liquid media, and the inoculated leaf samples of walnut with the fungal conidia sampled at 48, 96 and 144 h post inoculation (hpi). The completeness, correctness, and contiguity of the de novo transcriptome assemblies generated with Trinity, Oases, SOAPdenovo-Trans and Bridger were compared to identify a single superior reference assembly. In most of the assessment criteria including N50, Transrate score, number of ORFs with known description in gene bank, the percentage of reads mapped back to the transcript (RMBT), BUSCO score, Swiss-Prot coverage bin and RESM-EVAL score, the Bridger assembly was the superior and thus used as a reference for profiling the O. leptostyla transcriptome in liquid media vs. during walnut infection. The k-means clustering of transcripts resulted in four distinct transcription patterns across the three sampling time points. Most of the detected CAZy transcripts had elevated transcription at 96 hpi that is hypothetically concurrent with the start of intracellular growth. The in-silico analysis revealed 103 candidate effectors of which six were members of Necrosis and Ethylene Inducing Like Protein (NLP) gene family belonging to three distinct k-means clusters. This study provided a complex and temporal pattern of the CAZys and candidate effectors transcription during six days post O. leptostyla inoculation on walnut leaves, introducing a list of candidate virulence genes for validation in future studies.

Transcriptome analysis of flower colour reveals the correlation between SNP and differential expression genes in Phalaenopsis

BACKGROUND

Phalaenopsis is an important ornamental plant that has great economic value in the world flower market as one of the most popular flower resources.

OBJECTIVE

To investigate the flower colour formation of Phalaenopsis at the transcription level, the genes involved in flower color formation were identified from RNA-seq in this study.

METHODS

In this study, white and purple petals of Phalaenopsis were collected and analyzed to obtained (1) differential expression genes (DEGs) between white and purple flower color and (2) the association between single nucleotide polymorphisms (SNP) mutations and DEGs at the transcriptome level.

RESULTS

The results indicated that a total of 1,175 DEGs were identified, and 718 and 457 of them were up- and down-regulated genes, respectively. Gene Ontology and pathway enrichment showed that the biosynthesis of the secondary metabolites pathway was key to color formation, and the expression of 12 crucial genes (C4H, CCoAOMT, F3'H, UA3'5'GT, PAL, 4CL, CCR, CAD, CALDH, bglx, SGTase, and E1.11.17) that are involved in the regulation of flower color in Phalaenopsis.

CONCLUSION

This study reported the association between the SNP mutations and DEGs for color formation at RNA level, and provides a new insight to further investigate the gene expression and its relationship with genetic variants from RNA-seq data in other species.

StringFix: an annotation-guided transcriptome assembler improves the recovery of amino acid sequences from RNA-Seq reads

BACKGROUND

Reconstruction of amino acid sequences from assembled transcriptome is of interest in personalized medicine, for example, to predict drug-target (or protein-protein) interaction considering individual's genomic variations. Most of the existing transcriptome assemblers, however, seems not well suited for this purpose.

METHODS

In this work, we present StringFix, an annotation guided transcriptome assembly and protein sequence reconstruction software tool that takes genome-aligned reads and the annotations associated to the reference genome as input. The tool 'fixes' the pre-annotated transcript sequence by taking small variations into account, finally to produce possible amino acid sequences that are likely to exist in the test tissue.

RESULTS

The results show that, using outputs from existing reference-based assemblers as the input GTF-guide, StringFix could reconstruct amino acid sequences more precisely with higher sensitivity than direct generation using the recovered transcripts from all the assemblers we tested.

CONCLUSION

By using StringFix with the existing reference-based assemblers, one can recover not only a novel transcripts and isoforms but also the possible amino acid sequence stemming from them.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

Comprehensive phylogenomic time tree of bryophytes reveals deep relationships and uncovers gene incongruences in the last 500 million years of diversification

PREMISE

Bryophytes form a major component of terrestrial plant biomass, structuring ecological communities in all biomes. Our understanding of the evolutionary history of hornworts, liverworts, and mosses has been significantly reshaped by inferences from molecular data, which have highlighted extensive homoplasy in various traits and repeated bursts of diversification. However, the timing of key events in the phylogeny, patterns, and processes of diversification across bryophytes remain unclear.

METHODS

Using the GoFlag probe set, we sequenced 405 exons representing 228 nuclear genes for 531 species from 52 of the 54 orders of bryophytes. We inferred the species phylogeny from gene tree analyses using concatenated and coalescence approaches, assessed gene conflict, and estimated the timing of divergences based on 29 fossil calibrations.

RESULTS

The phylogeny resolves many relationships across the bryophytes, enabling us to resurrect five liverwort orders and recognize three more and propose 10 new orders of mosses. Most orders originated in the Jurassic and diversified in the Cretaceous or later. The phylogenomic data also highlight topological conflict in parts of the tree, suggesting complex processes of diversification that cannot be adequately captured in a single gene-tree topology.

CONCLUSIONS

We sampled hundreds of loci across a broad phylogenetic spectrum spanning at least 450 Ma of evolution; these data resolved many of the critical nodes of the diversification of bryophytes. The data also highlight the need to explore the mechanisms underlying the phylogenetic ambiguity at specific nodes. The phylogenomic data provide an expandable framework toward reconstructing a comprehensive phylogeny of this important group of plants.

Multifaceted synergistic facilitation mechanism of conductive polymers in promoting selenite bioreduction and biological detoxification

Here, polypyrrole (PPY) was first used to the bioreduction of toxic selenite, while the acceleration effect and mechanism were explored. Experiment results suggested that PPY could enhance the selenite bioreduction from 0.42 to 1.04 mg/(L·h). The tests of electrochemical analysis and cytochrome c (cyt-c) content confirmed that PPY promoted the intracellular/intracellular electron transfer of Shewanella oneidensis·MR-1 in selenite bioreduction process. The enhancement of metabolic activity by PPY contributed to biological detoxification, which was manifested in the increased extracellular polymeric substances (EPS), adenosine triphosphate (ATP), electron transfer system activity (ETSA), membrane permeability and enzyme activity. Transcriptome analysis of DEGs, KEGG pathway enrichment and GO functional classification verified that the environmental adaptability of Shewanella oneidensis·MR-1 was enhanced with the addition of PPY. The transmission electron microscopy (TEM) images indicated that PPY promoted the biosynthesis of selenium nanoparticles (SeNPs), which was beneficial to reduce cell damage. Combined with the above results, a multifaceted synergistic facilitation mechanism based on "conductive cross-linking network" was elaborated from electron transfer, microbial metabolism and environmental adaptability. This study shed light the effect of conductive polymers (CPs) on selenite bioreduction and provided new insights into the bioremediation of toxic pollutants.

Evolution and development of fruits of Erycina pusilla and other orchid species

Fruits play a crucial role in seed dispersal. They open along dehiscence zones. Fruit dehiscence zone formation has been intensively studied in Arabidopsis thaliana. However, little is known about the mechanisms and genes involved in the formation of fruit dehiscence zones in species outside the Brassicaceae. The dehiscence zone of A. thaliana contains a lignified layer, while dehiscence zone tissues of the emerging orchid model Erycina pusilla include a lipid layer. Here we present an analysis of evolution and development of fruit dehiscence zones in orchids. We performed ancestral state reconstructions across the five orchid subfamilies to study the evolution of selected fruit traits and explored dehiscence zone developmental genes using RNA-seq and qPCR. We found that erect dehiscent fruits with non-lignified dehiscence zones and a short ripening period are ancestral characters in orchids. Lignified dehiscence zones in orchid fruits evolved multiple times from non-lignified zones. Furthermore, we carried out gene expression analysis of tissues from different developmental stages of E. pusilla fruits. We found that fruit dehiscence genes from the MADS-box gene family and other important regulators in E. pusilla differed in their expression pattern from their homologs in A. thaliana. This suggests that the current A. thaliana fruit dehiscence model requires adjustment for orchids. Additionally, we discovered that homologs of A. thaliana genes involved in the development of carpel, gynoecium and ovules, and genes involved in lipid biosynthesis were expressed in the fruit valves of E. pusilla, implying that these genes may play a novel role in formation of dehiscence zone tissues in orchids. Future functional analysis of developmental regulators, lipid identification and quantification can shed more light on lipid-layer based dehiscence of orchid fruits.

Expression of type 1 vomeronasal receptors in the olfactory organ of the African lungfish, Protopterus dolloi

The vomeronasal organ is an olfactory organ found in amphibians and higher vertebrates. Type 1 vomeronasal receptors, one of the major olfactory receptors in vertebrates, are expressed in the vomeronasal organ in mammals. In amphibians and fish, they are expressed in the olfactory epithelium. The lungfish, which is the species of fish most closely related to amphibians, has a primitive vomeronasal organ: the recess epithelium. Expression of type 1 vomeronasal receptors has been reported in both the olfactory epithelium and the recess epithelium in three species of African lungfish and one species of South American lungfish. However, a previous study suggested that in the African lungfish Protopterus dolloi these receptors are expressed only in the olfactory epithelium. In this study, we identified 21 type 1 vomeronasal receptor genes in P. dolloi and examined the expression sites in the olfactory organ. In P. dolloi, most cells expressing the type 1 vomeronasal receptor were distributed in the olfactory epithelium, but a few were also found in the recess epithelium. This implies that the functions of the olfactory epithelium and the primitive vomeronasal organ are incompletely separated, and that all extant African and South American lungfish share this trait.

Comparing morphological and secretory aspects of cephalic glands among the New World coral snakes brings novel insights on their biological roles

Oral and other cephalic glands have been surveyed by several studies with distinct purposes. Despite the wide diversity and medical relevance of the New World coral snakes, studies focusing on understanding the biological roles of the glands within this group are still scarce. Specifically, the venom glands of some coral snakes were previously investigated but all other cephalic glands remain uncharacterized. In this sense, performing morphological and molecular analysis of these glands may help better understand their biological role. Here, we studied the morphology of the venom, infralabial, rictal, and harderian glands of thirteen species of Micrurus and Micruroides euryxanthus. We also performed a molecular characterization of these glands from selected species of Micrurus using transcriptomic and proteomic approaches. We described substantial morphological variation in the cephalic glands of New World coral snakes and structural evidence for protein-secreting cells in the inferior rictal glands. Our molecular analysis revealed that the venom glands, as expected, are majorly devoted to toxin production, however, the infralabial and inferior rictal glands also expressed some toxin genes at low to medium levels, despite the marked morphological differences. On the other hand, the harderian glands were dominated by the expression of lipocalins, but do not produce toxins. Our integrative analysis, including the prediction of biological processes and pathways, helped decipher some important traits of cephalic glands and better understand their biology.

Genome Sequencing Provides Novel Insights into Mudflat Burrowing Adaptations in Eel Goby Taenioides sp. (Teleost: Amblyopinae)

Amblyopinae is one of the lineage of bony fish that preserves amphibious traits living in tidal mudflat habitats. In contrast to other active amphibious fish, Amblyopinae species adopt a seemly more passive lifestyle by living in deep burrows of mudflat to circumvent the typical negative effects associated with terrestriality. However, little is known about the genetic origin of these mudflat deep-burrowing adaptations in Amblyopinae. Here we sequenced the first genome of Amblyopinae species, Taenioides sp., to elucidate their mudflat deep-burrowing adaptations. Our results revealed an assembled genome size of 774.06 Mb with 23 pseudochromosomes anchored, which predicted 22,399 protein-coding genes. Phylogenetic analyses indicated that Taenioides sp. diverged from the active amphibious fish of mudskipper approximately 28.3 Ma ago. In addition, 185 and 977 putative gene families were identified to be under expansion, contraction and 172 genes were undergone positive selection in Taenioides sp., respectively. Enrichment categories of top candidate genes under significant expansion and selection were mainly associated with hematopoiesis or angiogenesis, DNA repairs and the immune response, possibly suggesting their involvement in the adaptation to the hypoxia and diverse pathogens typically observed in mudflat burrowing environments. Some carbohydrate/lipid metabolism, and insulin signaling genes were also remarkably alterated, illustrating physiological remolding associated with nutrient-limited subterranean environments. Interestingly, several genes related to visual perception (e.g., crystallins) have undergone apparent gene losses, pointing to their role in the small vestigial eyes development in Taenioides sp. Our work provide valuable resources for understanding the molecular mechanisms underlying mudflat deep-burrowing adaptations in Amblyopinae, as well as in other tidal burrowing teleosts.

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

Optimizing an efficient ensemble approach for high-quality de novo transcriptome assembly of Thymus daenensis

Non-erroneous and well-optimized transcriptome assembly is a crucial prerequisite for authentic downstream analyses. Each de novo assembler has its own algorithm-dependent pros and cons to handle the assembly issues and should be specifically tested for each dataset. Here, we examined efficiency of seven state-of-art assemblers on ~ 30 Gb data obtained from mRNA-sequencing of Thymus daenensis. In an ensemble workflow, combining the outputs of different assemblers associated with an additional redundancy-reducing step could generate an optimized outcome in terms of completeness, annotatability, and ORF richness. Based on the normalized scores of 16 benchmarking metrics, EvidentialGene, BinPacker, Trinity, rnaSPAdes, CAP3, IDBA-trans, and Velvet-Oases performed better, respectively. EvidentialGene, as the best assembler, totally produced 316,786 transcripts, of which 235,730 (74%) were predicted to have a unique protein hit (on uniref100), and also half of its transcripts contained an ORF. The total number of unique BLAST hits for EvidentialGene was approximately three times greater than that of the worst assembler (Velvet-Oases). EvidentialGene could even capture 17% and 7% more average BLAST hits than BinPacker and Trinity. Although BinPacker and CAP3 produced longer transcripts, the EvidentialGene showed a higher collinearity between transcript size and ORF length. Compared with the other programs, EvidentialGene yielded a higher number of optimal transcript sets, further full-length transcripts, and lower possible misassemblies. Our finding corroborates that in non-model species, relying on a single assembler may not give an entirely satisfactory result. Therefore, this study proposes an ensemble approach of accompanying EvidentialGene pipelines to acquire a superior assembly for T. daenensis.

The Venom Composition of the Snake Tribe Philodryadini: 'Omic' Techniques Reveal Intergeneric Variability among South American Racers

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas, Chlorosoma and Xenoxybelis. The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas, while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis, respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.

Type 1 vomeronasal receptor expression in juvenile and adult lungfish olfactory organ

Lungfish are the most closely related fish to tetrapods. The olfactory organ of lungfish contains lamellae and abundant recesses at the base of lamellae. Based on the ultrastructural and histochemical characteristics, the lamellar olfactory epithelium (OE), covering the surface of lamellae, and the recess epithelium, contained in the recesses, are thought to correspond to the OE of teleosts and the vomeronasal organ (VNO) of tetrapods. With increasing body size, the recesses increase in number and distribution range in the olfactory organ. In tetrapods, the expression of olfactory receptors is different between the OE and VNO; for instance, the type 1 vomeronasal receptor (V1R) is expressed only in the OE in amphibians and mainly in the VNO in mammals. We recently reported that V1R-expressing cells are contained mainly in the lamellar OE but also rarely in the recess epithelium in the olfactory organ of lungfish of approximately 30 cm body length. However, it is unclear whether the distribution of V1R-expressing cells in the olfactory organ varies during development. In this study, we compared the expression of V1Rs in the olfactory organs between juveniles and adults of the African lungfish Protopterus aethiopicus and South American lungfish, Lepidosiren paradoxa. The density of V1R-expressing cells was higher in the lamellae than in the recesses in all specimens evaluated, and this pattern was more pronounced in juveniles than adults. In addition, the juveniles showed a higher density of V1R-expressing cells in the lamellae compared with the adults. Our results imply that differences in lifestyle between juveniles and adults are related to differences in the density of V1R-expressing cells in the lamellae of lungfish.

A Guide to Sequencing for Long Repetitive Regions

Full-length analysis of genes with highly repetitive sequences is challenging in two respects: assembly algorithm and sequencing accuracy. The de Bruijn graph often used in short-read assembly cannot distinguish adjacent repeat units. On the other hand, the accuracy of long reads is not yet high enough to identify each and every repeat unit. In this chapter, I present an example of a strategy to solve these problems and obtain the full length of long repeats by combining the extraction and assembly of repeat units based on overlap-layout-consensus and scaffolding by long reads.

Transcriptome alterations of radish shoots exposed to cadmium can be interpreted in the context of leaf senescence

Till now few transcriptome studies have described shoot responses of heavy metal (HM)-sensitive plants to excess Cd and still a unifying model of Cd action is lacking. Using RNA-seq technique, the transcriptome responses of radish (Raphanus sativus L.) leaves to Cd stress were investigated in plants raised hydroponically under control and 5.0 mg L^-1 Cd. The element was mainly accumulated in roots and led to declined biomass and photosynthetic pigments, increased H₂O₂ and lipid peroxidation, and the accumulation of sugars, protein thiols, and phytochelatins. Out of 524 differentially expressed genes (DEGs), 244 and 280 upregulated and downregulated ones were assigned to 82 and 115 GO terms, respectively. The upregulated DEGs were involved in osmotic regulation, protein metabolism, chelators, and carbohydrate metabolisms, whereas downregulated DEGs were related to photosynthesis, response to oxidative stress, glucosinolate, and secondary metabolite biosynthesis. Our transcriptome data suggest that Cd triggers ROS production and photosynthesis decline associated with increased proteolysis through ubiquitin-proteasome system (UPS)- and chloroplast-proteases and in this way brings about re-mobilization of N and C stores into amino acids and sugars. Meanwhile, declined glucosinolate metabolism in favor of chelator synthesis and upregulation of dehydrins as inferred from transcriptome analysis confers shoots some tolerance to the HM-derived ionic/osmotic imbalances. Thus, the induction of leaf senescence might be a major long-term response of HM-sensitive plants to Cd toxicity.

De novo transcriptome sequencing and gene co-expression reveal a genomic basis for drought sensitivity and evidence of a rapid local adaptation on Atlas cedar (Cedrus atlantica)

Introduction

Understanding the adaptive capacity to current climate change of drought-sensitive tree species is mandatory, given their limited prospect of migration and adaptation as long-lived, sessile organisms. Knowledge about the molecular and eco-physiological mechanisms that control drought resilience is thus key, since water shortage appears as one of the main abiotic factors threatening forests ecosystems. However, our current background is scarce, especially in conifers, due to their huge and complex genomes.

Methods

Here we investigated the eco-physiological and transcriptomic basis of drought response of the climate change-threatened conifer Cedrus atlantica. We studied C. atlantica seedlings from two locations with contrasting drought conditions to investigate a local adaptation. Seedlings were subjected to experimental drought conditions, and were monitored at immediate (24 hours) and extended (20 days) times. In addition, post-drought recovery was investigated, depicting two contrasting responses in both locations (drought resilient and non-resilient). Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of drought resilience and investigate a rapid local adaptation of C. atlantica.

Results

De novo transcriptome assembly was performed for the first time in this species, providing differences in gene expression between the immediate and extended treatments, as well as among the post-drought recovery phenotypes. Weighted gene co-expression network analysis showed a regulation of stomatal closing and photosynthetic activity during the immediate drought, consistent with an isohydric dynamic. During the extended drought, growth and flavonoid biosynthesis inhibition mechanisms prevailed, probably to increase root-to-shoot ratio and to limit the energy-intensive biosynthesis of secondary metabolites. Drought sensitive individuals failed in metabolism and photosynthesis regulation under drought stress, and in limiting secondary metabolite production. Moreover, genomic differences (SNPs) were found between drought resilient and sensitive seedlings, and between the two studied locations, which were mostly related to transposable elements.

Discussion

This work provides novel insights into the transcriptomic basis of drought response of C. atlantica, a set of candidate genes mechanistically involved in its drought sensitivity and evidence of a rapid local adaptation. Our results may help guide conservation programs for this threatened conifer, contribute to advance drought-resilience research and shed light on trees' adaptive potential to current climate change.

1000 spider silkomes: Linking sequences to silk physical properties

Spider silks are among the toughest known materials and thus provide models for renewable, biodegradable, and sustainable biopolymers. However, the entirety of their diversity still remains elusive, and silks that exceed the performance limits of industrial fibers are constantly being found. We obtained transcriptome assemblies from 1098 species of spiders to comprehensively catalog silk gene sequences and measured the mechanical, thermal, structural, and hydration properties of the dragline silks of 446 species. The combination of these silk protein genotype-phenotype data revealed essential contributions of multicomponent structures with major ampullate spidroin 1 to 3 paralogs in high-performance dragline silks and numerous amino acid motifs contributing to each of the measured properties. We hope that our global sampling, comprehensive testing, integrated analysis, and open data will provide a solid starting point for future biomaterial designs.

Enhancement of hexavalent chromium reduction by Shewanella oneidensis MR-1 in presence of copper nanoparticles via stimulating bacterial extracellular electron transfer and environmental adaptability

The bioreduction of hexavalent chromium (Cr(VI)) depends highly on bacterial activity, while the release of copper nanoparticles (Cu NPs) poses threats to microorganisms in the environment. This work demonstrated that Cr(VI) reduction efficiency of Shewanella oneidensis MR-1 was remarkably enhanced by 83.7% under 20 mg/L Cu NPs exposure. Cu NPs improved the electron migration capacity of Shewanella oneidensis MR-1 by enhancing bioelectrochemical performance and flavin mononucleotide secretion. Moreover, key genes related to extracellular electron transfer pathways, including direct electron transfer through outer-membrane proteins, flavin-mediated electron transfer, and conductive flagella, were generally upregulated under Cu NPs exposure. In addition, environmental adaptability of Shewanella oneidensis MR-1 was enhanced under Cu NPs exposure by improving environmental information processing and energy and reducing power production, promoting Cr(VI) reduction by Shewanella oneidensis MR-1. This work indicated that Cu NPs could enhance Cr(VI) reduction by Shewanella oneidensis MR-1 through regulating extracellular electron transfer and environmental adaptability.

Normalized Workflow to Optimize Hybrid De Novo Transcriptome Assembly for Non-Model Species: A Case Study in Lilium ledebourii (Baker) Boiss

A high-quality transcriptome is required to advance numerous bioinformatics workflows. Nevertheless, the effectuality of tools for de novo assembly and real precision assembled transcriptomes looks somewhat unexplored, particularly for non-model organisms with complicated (very long, heterozygous, polyploid) genomes. To disclose the performance of various transcriptome assembly programs, this study built 11 single assemblies and analyzed their performance on some significant reference-free and reference-based criteria. As well as to reconfirm the outputs of benchmarks, 55 BLAST were performed and compared using 11 constructed transcriptomes. Concisely, normalized benchmarking demonstrated that Velvet–Oases suffer from the worst results, while the EvidentialGene strategy can provide the most comprehensive and accurate transcriptome of Lilium ledebourii (Baker) Boiss. The BLAST results also confirmed the superiority of EvidentialGene, so it could capture even up to 59% more (than Velvet–Oases) unique gene hits. To promote assembly optimization, with the help of normalized benchmarking, PCA and AHC, it is emphasized that each metric can only provide part of the transcriptome status, and one should never settle for just a few evaluation criteria. This study supplies a framework for benchmarking and optimizing the efficiency of assembly approaches to analyze RNA-Seq data and reveals that selecting an inefficient assembly strategy might result in less identification of unique gene hits.

A chromosome-level genome of Brachymystax tsinlingensis provides resources and insights into salmonids evolution

Brachymystax tsinlingensis Li, 1966 is an endangered freshwater fish with economic, ecological, and scientific values. Study of the genome of B. tsinlingensis might be particularly insightful given that this is the only Brachymystax species with genome. We present a high-quality chromosome-level genome assembly and protein-coding gene annotation for B. tsinlingensis with Illumina short reads, Nanopore long reads, Hi-C sequencing reads, and RNA-seq reads from 5 tissues/organs. The final chromosome-level genome size is 2,031,709,341 bp with 40 chromosomes. We found that the salmonids have a unique GC content and codon usage, have a slower evolutionary rate, and possess specific positively selected genes. We also confirmed the salmonids have undergone a whole-genome duplication event and a burst of transposon-mediated repeat expansion, and lost HoxAbβ Hox cluster, highly expressed genes in muscle may partially explain the migratory habits of B. tsinlingensis. The high-quality B. tsinlingensis assembled genome could provide a valuable reference for the study of other salmonids as well as aid the conservation of this endangered species.

Proteotranscriptomics - A facilitator in omics research

Applications in omics research, such as comparative transcriptomics and proteomics, require the knowledge of the species-specific gene sequence and benefit from a comprehensive high-quality annotation of the coding genes to achieve high coverage. While protein-coding genes can in simple cases be detected by scanning the genome for open reading frames, in more complex genomes exonic sequences are separated by introns. Despite advances in sequencing technologies that allow for ever-growing numbers of genomes, the quality of many of the provided genome assemblies do not reach reference quality. These non-contiguous assemblies with gaps and the necessity to predict splice sites limit accurate gene annotation from solely genomic data. In contrast, the transcriptome only contains transcribed gene regions, is devoid of introns and thus provides the optimal basis for the identification of open reading frames. The additional integration of proteomics data to validate predicted protein-coding genes further enriches for accurate gene models. This review outlines the principles of the proteotranscriptomics approach, discusses common challenges and suggests methods for improvement.

TransMeta simultaneously assembles multisample RNA-seq reads

Assembling RNA-seq reads into full-length transcripts is crucial in transcriptomic studies and poses computational challenges. Here we present TransMeta, a simple and robust algorithm that simultaneously assembles RNA-seq reads from multiple samples. TransMeta is designed based on the newly introduced vector-weighted splicing graph model, which enables accurate reconstruction of the consensus transcriptome via incorporating a cosine similarity-based combing strategy and a newly designed label-setting path-searching strategy. Tests on both simulated and real data sets show that TransMeta consistently outperforms PsiCLASS, StringTie2 plus its merge mode, and Scallop plus TACO, the most popular tools, in terms of precision and recall under a wide range of coverage thresholds at the meta-assembly level. Additionally, TransMeta consistently shows superior performance at the individual sample level.

Time-series transcriptomic screening of factors contributing to the cross-tolerance to UV radiation and anhydrobiosis in tardigrades

BACKGROUND

Tardigrades are microscopic animals that are capable of tolerating extreme environments by entering a desiccated state of suspended animation known as anhydrobiosis. While antioxidative stress proteins, antiapoptotic pathways and tardigrade-specific intrinsically disordered proteins have been implicated in the anhydrobiotic machinery, conservation of these mechanisms is not universal within the phylum Tardigrada, suggesting the existence of overlooked components.

RESULTS

Here, we show that a novel Mn-dependent peroxidase is an important factor in tardigrade anhydrobiosis. Through time-series transcriptome analysis of Ramazzottius varieornatus specimens exposed to ultraviolet light and comparison with anhydrobiosis entry, we first identified several novel gene families without similarity to existing sequences that are induced rapidly after stress exposure. Among these, a single gene family with multiple orthologs that is highly conserved within the phylum Tardigrada and enhances oxidative stress tolerance when expressed in human cells was identified. Crystallographic study of this protein suggested Zn or Mn binding at the active site, and we further confirmed that this protein has Mn-dependent peroxidase activity in vitro.

CONCLUSIONS

Our results demonstrated novel mechanisms for coping with oxidative stress that may be a fundamental mechanism of anhydrobiosis in tardigrades. Furthermore, localization of these sets of proteins mainly in the Golgi apparatus suggests an indispensable role of the Golgi stress response in desiccation tolerance.

Chromosome-level genome assembly of the black widow spider Latrodectus elegans illuminates composition and evolution of venom and silk proteins

BACKGROUND

The black widow spider has both extraordinarily neurotoxic venom and three-dimensional cobwebs composed of diverse types of silk. However, a high-quality reference genome for the black widow spider was still unavailable, which hindered deep understanding and application of the valuable biomass.

FINDINGS

We assembled the Latrodectus elegans genome, including a genome size of 1.57 Gb with contig N50 of 4.34 Mb and scaffold N50 of 114.31 Mb. Hi-C scaffolding assigned 98.08% of the genome to 14 pseudo-chromosomes, and with BUSCO, completeness analysis revealed that 98.4% of the core eukaryotic genes were completely present in this genome. Annotation of this genome identified that repetitive sequences account for 506.09 Mb (32.30%) and 20,167 protein-coding genes, and specifically, we identified 55 toxin genes and 26 spidroins and provide preliminary analysis of their composition and evolution.

CONCLUSIONS

We present the first chromosome-level genome assembly of a black widow spider and provide substantial toxin and spidroin gene resources. These high-qualified genomic data add valuable resources from a representative spider group and contribute to deep exploration of spider genome evolution, especially in terms of the important issues on the diversification of venom and web-weaving pattern. The sequence data are also firsthand templates for further application of the spider biomass.

An NGS-Based Phylogeny of Orthotricheae (Orthotrichaceae, Bryophyta) With the Proposal of the New Genus Rehubryum From Zealandia

Phylogenomic data increase the possibilities of resolving the evolutionary and systematic relationships among taxa. This is especially valuable in groups with few and homoplasious morphological characters, in which systematic and taxonomical delimitations have been traditionally difficult. Such is the case of several lineages within Bryophyta, like Orthotrichaceae, the second most diverse family of mosses. Members of tribe Orthotricheae are common in temperate and cold regions, as well as in high tropical mountains. In extratropical areas, they represent one of the main components of epiphytic communities, both in dry and oceanic or hyperoceanic conditions. The epiphytic environment is considered a hostile one for plant development, mainly due to its low capacity of moisture retention. Thus, the diversification of the Orthotrichaceae in this environment could be seen as striking. Over the last two decades, great taxonomic and systematic progresses have led to a rearrangement at the generic level in this tribe, providing a new framework to link environment to patterns of diversification. Here, we use nuclear loci targeted with the GoFlag 408 enrichment probe set to generate a well-sampled phylogeny with well-supported suprageneric taxa and increasing the phylogenetic resolution within the two recognized subtribes. Specifically, we show that several genera with Ulota-like morphology jointly constitute an independent lineage. Within this lineage, the recently described Atlantichella from Macaronesia and Western Europe appears as the sister group of Ulota bellii from Zealandia. This latter species is here segregated in the new genus Rehubryum. Assessment of the ecological and biogeographical affinities of the species within the phylogenetic framework suggests that niche adaptation (including climate and substrate) may be a key evolutionary driver that shaped the high diversification of Orthotricheae.

Partitioning RNAs by length improves transcriptome reconstruction from short-read RNA-seq data

The accuracy of methods for assembling transcripts from short-read RNA sequencing data is limited by the lack of long-range information. Here we introduce Ladder-seq, an approach that separates transcripts according to their lengths before sequencing and uses the additional information to improve the quantification and assembly of transcripts. Using simulated data, we show that a kallisto algorithm extended to process Ladder-seq data quantifies transcripts of complex genes with substantially higher accuracy than conventional kallisto. For reference-based assembly, a tailored scheme based on the StringTie2 algorithm reconstructs a single transcript with 30.8% higher precision than its conventional counterpart and is more than 30% more sensitive for complex genes. For de novo assembly, a similar scheme based on the Trinity algorithm correctly assembles 78% more transcripts than conventional Trinity while improving precision by 78%. In experimental data, Ladder-seq reveals 40% more genes harboring isoform switches compared to conventional RNA sequencing and unveils widespread changes in isoform usage upon m6A depletion by Mettl14 knockout.

FrogCap: A modular sequence capture probe-set for phylogenomics and population genetics for all frogs, assessed across multiple phylogenetic scales

Despite the prevalence of high-throughput sequencing in phylogenetics, many relationships remain difficult to resolve because of conflicting signal among genomic regions. Selection of different types of molecular markers from different genomic regions is required to overcome these challenges. For evolutionary studies in frogs, we introduce the publicly available FrogCap suite of genomic resources, which is a large collection of ~15,000 markers that unifies previous genetic sequencing efforts. FrogCap is designed to be modular, such that subsets of markers and SNPs can be selected based on the desired phylogenetic scale. FrogCap uses a variety of marker types that include exons and introns, ultraconserved elements, and previously sequenced Sanger markers, which span up to 10,000 bp in alignment lengths; in addition, we demonstrate potential for SNP-based analyses. We tested FrogCap using 121 samples distributed across five phylogenetic scales, comparing probes designed using a consensus- or exemplar genome-based approach. Using the consensus design is more resilient to issues with sensitivity, specificity, and missing data than picking an exemplar genome sequence. We also tested the impact of different bait kit sizes (20,020 vs. 40,040) on depth of coverage and found triple the depth for the 20,020 bait  kit. We observed sequence capture success (i.e., missing data, sequenced markers/bases, marker length, and informative sites) across phylogenetic scales. The incorporation of different marker types is effective for deep phylogenetic relationships and shallow population genetics studies. Having demonstrated FrogCap's utility and modularity, we conclude that these new resources are efficacious for high-throughput sequencing projects across variable timescales.

Differences in PLA2 Constitution Distinguish the Venom of Two Endemic Brazilian Mountain Lanceheads, Bothrops cotiara and Bothrops fonsecai

Interspecific differences in snake venom compositions can result from distinct regulatory mechanisms acting in each species. However, comparative analyses focusing on identifying regulatory elements and patterns that led to distinct venom composition are still scarce. Among venomous snakes, Bothrops cotiara and Bothrops fonsecai represent ideal models to complement our understanding of the regulatory mechanisms of venom production. These recently diverged species share a similar specialized diet, habitat, and natural history, but each presents a distinct venom phenotype. Here, we integrated data from the venom gland transcriptome and miRNome and the venom proteome of B. fonsecai and B. cotiara to better understand the regulatory mechanisms that may be acting to produce differing venom compositions. We detected not only the presence of similar toxin isoforms in both species but also distinct expression profiles of phospholipases A2 (PLA2) and some snake venom metalloproteinases (SVMPs) and snake venom serine proteinases (SVSPs) isoforms. We found evidence of modular expression regulation of several toxin isoforms implicated in venom divergence and observed correlated expression of several transcription factors. We did not find strong evidence for miRNAs shaping interspecific divergence of the venom phenotypes, but we identified a subset of toxin isoforms whose final expression may be fine-tuned by specific miRNAs. Sequence analysis on orthologous toxins showed a high rate of substitutions between PLA2s, which indicates that these toxins may be under strong positive selection or represent paralogous toxins in these species. Our results support other recent studies in suggesting that gene regulation is a principal mode of venom evolution across recent timescales, especially among species with conserved ecotypes.

Tiglon enables accurate transcriptome assembly via integrating mappings of different aligners

Full-length transcript reconstruction has a pivotal role in RNA-seq data analysis. In this research, we present a new genome-guided transcriptome assembly algorithm, namely Tiglon, which integrates multiple alignments of different mapping tools and builds the labeled splice graphs, followed by a label-based dynamic path-searching strategy to reconstruct the transcripts. We evaluate Tiglon on a simulated dataset and 12 real datasets under the Hisat2 and Star mappings. The results indicate that the integrating techniques of Tiglon exhibit great superiority over the state-of-the-art assemblers, including StringTie2 and Scallop, depending on Hisat2 alignments, Star alignments, or the merged alignments of both. Especially, Tiglon is significantly powerful in recovering lowly expressed transcripts.

•

Tiglon is designed for integrating multiple alignments to assemble transcripts

•

Integrating alignments of different aligners is helpful for transcriptome assembly

•

Tiglon proposes a new graph model called the labeled splice graph

•

Our experiments demonstrate that Tiglon outperforms the leading assemblers

A transcriptomic reevaluation of the accessory olfactory organ in Bichir (Polypterus senegalus)

Fish possess one olfactory organ called the olfactory epithelium (OE), by which various chemical substances are detected. On the other hand, tetrapods possess two independent olfactory organs called the main olfactory epithelium (MOE) and vomeronasal organ (VNO), each of which mainly detects general odorants and pheromones, respectively. Traditionally, the VNO, so-called concentrations of vomeronasal neurons, was believed to have originated in tetrapods. However, recent studies have identified a primordial VNO in lungfish, implying that the origin of the VNO was earlier than traditionally expected. In this study, we examined the presence/absence of the VNO in the olfactory organ of bichir (Polypterus senegalus), which is the most ancestral group of extant bony vertebrates. In particular, we conducted a transcriptomic evaluation of the accessory olfactory organ (AOO), which is anatomically separated from the main olfactory organ (MOO) in bichir. As a result, several landmark genes specific to the VNO and MOE in tetrapods were both expressed in the MOO and AOO, suggesting that these organs were not functionally distinct in terms of pheromone and odorant detection. Instead, differentially expressed gene (DEG) analysis showed that DEGs in AOO were enriched in genes for cilia movement, implying its additional and specific function in efficient water uptake into the nasal cavity other than chemosensing. This transcriptomic study provides novel insight into the long-standing question of AOO function in bichir and suggests that VNO originated in the lineage of lobe-finned fish during vertebrate evolution.

A simple guide to de novo transcriptome assembly and annotation

A transcriptome constructed from short-read RNA sequencing (RNA-seq) is an easily attainable proxy catalog of protein-coding genes when genome assembly is unnecessary, expensive or difficult. In the absence of a sequenced genome to guide the reconstruction process, the transcriptome must be assembled de novo using only the information available in the RNA-seq reads. Subsequently, the sequences must be annotated in order to identify sequence-intrinsic and evolutionary features in them (for example, protein-coding regions). Although straightforward at first glance, de novo transcriptome assembly and annotation can quickly prove to be challenging undertakings. In addition to familiarizing themselves with the conceptual and technical intricacies of the tasks at hand and the numerous pre- and post-processing steps involved, those interested must also grapple with an overwhelmingly large choice of tools. The lack of standardized workflows, fast pace of development of new tools and techniques and paucity of authoritative literature have served to exacerbate the difficulty of the task even further. Here, we present a comprehensive overview of de novo transcriptome assembly and annotation. We discuss the procedures involved, including pre- and post-processing steps, and present a compendium of corresponding tools.