The Pfam protein families database

Unraveling the Diverse Roles of Neglected Genes Containing Domains of Unknown Function (DUFs): Progress and Perspective

Domain of unknown function (DUF) is a general term for many uncharacterized domains with two distinct features: relatively conservative amino acid sequence and unknown function of the domain. In the Pfam 35.0 database, 4795 (24%) gene families belong to the DUF type, yet, their functions remain to be explored. This review summarizes the characteristics of the DUF protein families and their functions in regulating plant growth and development, generating responses to biotic and abiotic stress, and other regulatory roles in plant life. Though very limited information is available about these proteins yet, by taking advantage of emerging omics and bioinformatic tools, functional studies of DUF proteins could be utilized in future molecular studies.

Genome-Wide Analysis of SIMILAR TO RCD ONE (SRO) Family Revealed Their Roles in Abiotic Stress in Poplar

SIMILAR TO RCD ONE (SRO) gene family is a small plant-specific gene family responsible for growth, development, and stress responses. In particular, it plays a vital role in responding to abiotic stresses such as salt, drought, and heavy metals. Poplar SROs are rarely reported to date. In this study, a total of nine SRO genes were identified from Populus simonii × Populus nigra, which are more similar to dicotyledon SRO members. According to phylogenetic analysis, the nine PtSROs can be divided into two groups, and the members in the same cluster have a similar structure. There were some cis-regulatory elements related to abiotic stress response and hormone-induced factors identified in the promoter regions of PtSROs members. Subcellular localization and transcriptional activation activity of PtSRO members revealed a consistent expression profile of the genes with similar structural profiles. In addition, both RT-qPCR and RNA-Seq results indicated that PtSRO members responded to PEG-6000, NaCl, and ABA stress in the roots and leaves of Populus simonii × Populus nigra. The PtSRO genes displayed different expression patterns and peaked at different time points in the two tissues, which was more significant in the leaves. Among them, PtSRO1c and PtSRO2c were more prominent in response to abiotic stress. Furthermore, protein interaction prediction showed that the nine PtSROs might interact with a broad range of transcription factors (TFs) involved in stress responses. In conclusion, the study provides a solid basis for functional analysis of the SRO gene family in abiotic stress responses in poplar.

Identification and Comparative Expression Profiles of Candidate Olfactory Receptors in the Transcriptomes of the Important Egg Parasitoid Wasp Anastatus japonicus Ashmead (Hymenoptera: Eupelmidae)

Anastatus japonicus Ashmead is an egg parasitoid wasp important for the biological control of fruit crop pests. The olfaction of parasitoids is crucial to searching for host pests in fruit crops. In this study, we sequenced and analyzed the antennal and abdominal transcriptomes of A. japonicus to better understand the olfactory mechanisms in this species. A total of 201 putative olfactory receptor genes were identified, including 184 odorant receptors (ORs) and 17 ionotropic receptors (IRs). Then, we assayed the tissue-specific and sex-biased expression profiles of those genes based on the transcriptional levels. In total, 165 ORs and 15 IRs had upregulated expression in the antennae. The expression levels of 133 ORs, including odorant receptor co-receptor (AjapORco), and 10 IRs, including AjapIR8a, were significantly different between the female and male antennae. Our results provide valuable information for further studies on the molecular mechanisms of the olfactory system in A. japonicus.

Transcriptome sequencing of Cocos nucifera leaves in response to Rhynchophorus ferrugineus infestation

Red palm weevil (RPW, Rhynchophorus ferrugineus) is an invasive pest of palms. In China, coconut (Cocos nucifera) production is being significantly affected by the RPW attack. To develop a long-term RPW control strategy, host-plant resistance is the most sustainable option. In this regard, the availability of transcriptome sequencing data from RPW-infected coconut plants can be highly useful. Therefore, the present study assessed coconut leaf physiological responses and transcriptional changes after different days of RPW attack i.e., 5, 10, 15, 20, and 25 days after infestation (DAI). A comparison of physiological data indicated that populations with the higher number of RPW insects i.e., population C (15 males +21 females) and D (20 males +28 females) triggered higher antioxidant enzyme activities. We used this data to study the transcriptomic responses on 5 and 20 DAI. Of the 38,432 detected transcripts, 3,984, 1,981, 3,925, and 2,257 were differentially expressed in CK (control/no RPW)_vs._C (5 DAI), CK_vs._D (5 DAI), CK_vs._C (20 DAI), and CK_vs._D (20 DAI), respectively. These transcripts were enriched in plant-pathogen interaction, phenylpropanoid/flavonoid biosynthesis, amino sugar and nucleotide sugar metabolism, plant hormone signal transduction, mitogen-activated protein kinase, and reactive oxygen scavenging pathway. We discuss these results and present several candidate genes to be manipulated for developing a sustainable strategy to control RPW attack regarding host-plant resistance. Furthermore, these findings provide a basis for developing effective early and late RPW attack detection strategies.

Hybridization and introgression are prevalent in Southern European Erysimum (Brassicaceae) species

BACKGROUND AND AIMS

Hybridization is a common and important force in plant evolution. One of its outcomes is introgression - the transfer of small genomic regions from one taxon to another by hybridization and repeated backcrossing. This process is believed to be common in glacial refugia, where range expansions and contractions can lead to cycles of sympatry and isolation, creating conditions for extensive hybridization and introgression. Polyploidization is another genome-wide process with a major influence on plant evolution. Both hybridization and polyploidization can have complex effects on plant evolution. However, these effects are often difficult to understand in recently evolved species complexes.

METHODS

We combined flow cytometry, analyses of transcriptomic sequences and pollen tube growth assays to investigate the consequences of polyploidization, hybridization and introgression on the recent evolution of several Erysimum (Brassicaceae) species from the South of the Iberian Peninsula, a well-known glacial refugium. This species complex differentiated in the last 2 million years, and its evolution has been hypothesized to be determined mainly by polyploidization, interspecific hybridization and introgression.

KEY RESULTS

Our results support a scenario of widespread hybridization involving both extant and 'ghost' taxa. Several taxa studied here, most notably those with purple corollas, are polyploids, probably of allopolyploid origin. Moreover, hybridization in this group might be an ongoing phenomenon, as pre-zygotic barriers appeared weak in many cases.

CONCLUSIONS

The evolution of Erysimum spp. has been determined by hybridization to a large extent. Species with purple (polyploids) and yellow flowers (mostly diploid) exhibit a strong signature of introgression in their genomes, indicating that hybridization occurred regardless of colour and across ploidy levels. Although the adaptive value of such genomic exchanges remains unclear, our results demonstrate the significance of hybridization for plant diversification, which should be taken into account when studying plant evolution.

Novel insights on genes and pathways involved in Pinus elliottii response to resinosis

Pinus elliottii, an important coniferous timber species, has recently become one of the most popular sources of resin in China. Resinosis is a common disease that may negatively affect pine tree growth and production. In this study, we used single-molecule real-time sequencing and Illumina RNA sequencing to generate an accurate transcriptome for P. elliottii. The transcriptome included 90,026 transcripts, 5160 long non-coding RNAs and 7710 transcription factors. We then analyzed RNA-sequencing, small RNA-sequencing and degradome data to identify genes, miRNAs and key miRNA-target pairs involved in response to resinosis in P. elliottii. We identified 1305 genes and 1151 miRNAs exhibiting significant differential expression in response to resinosis. According to the degradome sequencing analysis, 318 differentially expressed transcripts were targets of 14 differentially expressed miRNAs. Our study has provided resources for further functional characterization of genes and miRNAs involved in resinosis in P. elliottii, which should aid the future disease-resistance breeding of this species.

Genome-wide identification and characterization of exapted transposable elements in the large genome of sunflower (Helianthus annuus L.)

Transposable elements (TEs) are an important source of genome variability, playing many roles in the evolution of eukaryotic species. Besides well-known phenomena, TEs may undergo the exaptation process and generate the so-called exapted transposable element genes (ETEs). Here we present a genome-wide survey of ETEs in the large genome of sunflower (Helianthus annuus L.), in which the massive amount of TEs, provides a significant source for exaptation. A library of sunflower TEs was used to build TE-specific Hidden Markov Model profiles, to search for all available sunflower gene products. In doing so, 20 016 putative ETEs were identified and further investigated for the characteristics that distinguish TEs from genes, leading to the validation of 3530 ETEs. The analysis of ETEs transcription patterns under different stress conditions showed a differential regulation triggered by treatments mimicking biotic and abiotic stress; furthermore, the distribution of functional domains of differentially regulated ETEs revealed a relevant presence of domains involved in many aspects of cellular functions. A comparative genomic investigation was performed including species representative of Asterids and appropriate outgroups: the bulk of ETEs that resulted were specific to the sunflower, while few ETEs presented orthologues in the genome of all analyzed species, making the hypothesis of a conserved function. This study highlights the crucial role played by exaptation, actively contributing to species evolution.

SNP-based high-density linkage map construction and QTL mapping of black spot disease resistance in Chinese sand pear

Black spot disease (PBS) caused by Alternaria alternata is an economic disease of pear (Pyrus pyrifolia Nakai). Developing cultivars with durable PBS resistance traits is an important research objective for improving pear germplasm. The Deshengxiang is a popular pear variety in China and resistant to PBS. This study aimed to detect quantitative trait loci (QTL) associated with PBS resistance trait in pear and determine closely linked molecular markers by specific locus amplified fragment sequencing (SLAF-seq). F1 population resulting from a cross between "Deshengxiang" (female) and "Guiguan," a susceptible (male) variety, was developed and evaluated in 2016 and 2017. SLAF technology was used to discover SNPs in the F1 individuals and subsequently a high-density genetic linkage map for PBS resistance was constructed which contained 17,604 SNP markers. Based on the linkage map, the markers were distributed into 17 linkage groups, spanning 1548.48 cM, with a mean marker distance of 0.09 cM, representing the densest genetic map of the genus Pyrus. QTL analysis of PBS resistance identified a locus strongly related to PBS resistance at 77.68 ~ 112.99 cM on linkage group 15, which was further narrowed down to 93.79 ~ 112.99 cM. Two markers, Marker94293 and Marker94206, located at 97.47 and 102.93 cM, were closely associated with PBS resistance, with a Δ (SNP index) value of 0.46. Co-localization of QTL interval, bioinformatics analysis, and functional annotation revealed PBS putative candidate genes. Overall, the high-density pear linkage map is a suitable reference for mapping PBS resistance trait, QTL, and genes identified in this study contribute information that could be useful for PBS improvement in pear.

Anaerobic methanotroph 'Candidatus Methanoperedens nitroreducens' has a pleomorphic life cycle

'Candidatus Methanoperedens' are anaerobic methanotrophic (ANME) archaea with global importance to methane cycling. Here meta-omics and fluorescence in situ hybridization (FISH) were applied to characterize a bioreactor dominated by 'Candidatus Methanoperedens nitroreducens' performing anaerobic methane oxidation coupled to nitrate reduction. Unexpectedly, FISH revealed the stable co-existence of two 'Ca. M. nitroreducens' morphotypes: the archetypal coccobacilli microcolonies and previously unreported planktonic rods. Metagenomic analysis showed that the 'Ca. M. nitroreducens' morphotypes were genomically identical but had distinct gene expression profiles for proteins associated with carbon metabolism, motility and cell division. In addition, a third distinct phenotype was observed, with some coccobacilli 'Ca. M. nitroreducens' storing carbon as polyhydroxyalkanoates. The phenotypic variation of 'Ca. M. nitroreducens' probably aids their survival and dispersal in the face of sub-optimal environmental conditions. These findings further demonstrate the remarkable ability of members of the 'Ca. Methanoperedens' to adapt to their environment.

Chromosome-level genome assembly reveals potential epigenetic mechanisms of the thermal tolerance in the oriental fruit fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis (Hendel), has very strong ecological adaptability and phenotypic plasticity. Here, the genome of B. dorsalis was assembled into 549.45 Mb sequences with a contig N50 length of 12.81 Mb. Among, 95.67 % assembled genome sequences were anchored on six chromosomes with an N50 length of 94.63 Mb. According to the basic characteristics of the sex chromosomes of Tephritidae, the X chromosome of B. dorsalis was identified. Significant gene expansions were detected in several important gene families related to adaptability. In particular, we annotated 50 histone modification enzymes (HMEs) in this genome. A comparative transcriptome analysis indicated that 12 HME genes were differentially expressed in two thermo-tolerant strains (heat and cold). Interestingly, four and seven of the 12 HME genes responded to heat shock or cold hardening, respectively. These evidences suggested that the histone modification as an epigenetic modification may be involved in the thermal tolerance of B. dorsalis, but with different regulation mechanisms in thermal acclimation and hardening. The high quality genome of B. dorsalis provides an invaluable resource for further functional genomic study. Moreover, comparative genomic analysis will shed insights on revealing the mechanisms of adaptive evolution in this fly.

RNA-Seq and genetic diversity analysis of faba bean (Vicia faba L.) varieties in China

Background

Faba bean (Vicia faba L) is one of the most important legumes in the world. However, there is relatively little genomic information available for this species owing to its large genome. The lack of data impedes the discovery of molecular markers and subsequent genetic research in faba bean. The objective of this study was to analyze the faba bean transcriptome, and to develop simple sequence repeat (SSR) markers to determine the genetic diversity of 226 faba bean varieties derived from different regions in China.

Methods

Faba bean varieties with different phenotype were used in transcriptome analysis. The functions of the unigenes were analyzed using various database. SSR markers were developed and the polymorphic markers were selected to conduct genetic diversity analysis.

Results

A total of 92.43 Gb of sequencing data was obtained in this study, and 133,487 unigene sequences with a total length of 178,152,541 bp were assembled. A total of 5,200 SSR markers were developed on the basis of RNA-Seq analysis. Then, 200 SSR markers were used to evaluate polymorphisms. In total, 103 (51.5%) SSR markers showed significant and repeatable bands between different faba bean varieties. Clustering analysis revealed that 226 faba bean materials were divided into five groups. Genetic diversity analysis revealed that the relationship between different faba beans in China was related, especially in the same region. These results provided a valuable data resource for annotating genes to different categories and developing SSR markers.

Multi-locus genome-wide association studies reveal genomic regions and putative candidate genes associated with leaf spot diseases in African groundnut (Arachis hypogaea L.) germplasm

Early leaf spot (ELS) and late leaf spot (LLS) diseases are the two most destructive groundnut diseases in Ghana resulting in ≤ 70% yield losses which is controlled largely by chemical method. To develop leaf spot resistant varieties, the present study was undertaken to identify single nucleotide polymorphism (SNP) markers and putative candidate genes underlying both ELS and LLS. In this study, six multi-locus models of genome-wide association study were conducted with the best linear unbiased predictor obtained from 294 African groundnut germplasm screened for ELS and LLS as well as image-based indices of leaf spot diseases severity in 2020 and 2021 and 8,772 high-quality SNPs from a 48 K SNP array Axiom platform. Ninety-seven SNPs associated with ELS, LLS and five image-based indices across the chromosomes in the 2 two sub-genomes. From these, twenty-nine unique SNPs were detected by at least two models for one or more traits across 16 chromosomes with explained phenotypic variation ranging from 0.01 - 62.76%, with exception of chromosome (Chr) 08 (Chr08), Chr10, Chr11, and Chr19. Seventeen potential candidate genes were predicted at ± 300 kbp of the stable/prominent SNP positions (12 and 5, down- and upstream, respectively). The results from this study provide a basis for understanding the genetic architecture of ELS and LLS diseases in African groundnut germplasm, and the associated SNPs and predicted candidate genes would be valuable for breeding leaf spot diseases resistant varieties upon further validation.

Computational Methods and Deep Learning for Elucidating Protein Interaction Networks

Protein interactions play a critical role in all biological processes, but experimental identification of protein interactions is a time- and resource-intensive process. The advances in next-generation sequencing and multi-omics technologies have greatly benefited large-scale predictions of protein interactions using machine learning methods. A wide range of tools have been developed to predict protein-protein, protein-nucleic acid, and protein-drug interactions. Here, we discuss the applications, methods, and challenges faced when employing the various prediction methods. We also briefly describe ways to overcome the challenges and prospective future developments in the field of protein interaction biology.

Toward generalizable prediction of antibody thermostability using machine learning on sequence and structure features

Over the last three decades, the appeal for monoclonal antibodies (mAbs) as therapeutics has been steadily increasing as evident with FDA's recent landmark approval of the 100th mAb. Unlike mAbs that bind to single targets, multispecific biologics (msAbs) have garnered particular interest owing to the advantage of engaging distinct targets. One important modular component of msAbs is the single-chain variable fragment (scFv). Despite the exquisite specificity and affinity of these scFv modules, their relatively poor thermostability often hampers their development as a potential therapeutic drug. In recent years, engineering antibody sequences to enhance their stability by mutations has gained considerable momentum. As experimental methods for antibody engineering are time-intensive, laborious and expensive, computational methods serve as a fast and inexpensive alternative to conventional routes. In this work, we show two machine learning approaches - one with pre-trained language models (PTLM) capturing functional effects of sequence variation, and second, a supervised convolutional neural network (CNN) trained with Rosetta energetic features - to better classify thermostable scFv variants from sequence. Both of these models are trained over temperature-specific data (TS50 measurements) derived from multiple libraries of scFv sequences. On out-of-distribution (refers to the fact that the out-of-distribution sequnes are blind to the algorithm) sequences, we show that a sufficiently simple CNN model performs better than general pre-trained language models trained on diverse protein sequences (average Spearman correlation coefficient, ρ , of 0.4 as opposed to 0.15). On the other hand, an antibody-specific language model performs comparatively better than the CNN model on the same task (ρ = 0.52). Further, we demonstrate that for an independent mAb with available thermal melting temperatures for 20 experimentally characterized thermostable mutations, these models trained on TS50 data could identify 18 residue positions and 5 identical amino-acid mutations showing remarkable generalizability. Our results suggest that such models can be broadly applicable for improving the biological characteristics of antibodies. Further, transferring such models for alternative physicochemical properties of scFvs can have potential applications in optimizing large-scale production and delivery of mAbs or bsAbs.

A high-quality genome compendium of the human gut microbiome of Inner Mongolians

Metagenome-based resources have revealed the diversity and function of the human gut microbiome, but further understanding is limited by insufficient genome quality and a lack of samples from typically understudied populations. Here we used hybrid long-read PromethION and short-read HiSeq sequencing to characterize the faecal microbiota of 60 Inner Mongolian individuals (n = 180 samples over three time points) who were part of a probiotic yogurt intervention trial. We present the Inner Mongolian Gut Genome catalogue, comprising 802 closed and 5,927 high-quality metagenome-assembled genomes. This approach achieved high genome continuity and substantially increased the resolution of genomic elements, including ribosomal RNA operons, metabolic gene clusters, prophages and insertion sequences. Particularly, we report the ribosomal RNA operon copy numbers for uncultured species, over 12,000 previously undescribed gut prophages and the distribution of insertion sequence elements across gut bacteria. Overall, these data provide a high-quality, large-scale resource for studying the human gut microbiota.

The inside scoop: Comparative genomics of two intranuclear bacteria, "Candidatus Berkiella cookevillensis" and "Candidatus Berkiella aquae"

"Candidatus Berkiella cookevillensis" (strain CC99) and "Candidatus Berkiella aquae" (strain HT99), belonging to the Coxiellaceae family, are gram-negative bacteria isolated from amoebae in biofilms present in human-constructed water systems. Both bacteria are obligately intracellular, requiring host cells for growth and replication. The intracellular bacteria-containing vacuoles of both bacteria closely associate with or enter the nuclei of their host cells. In this study, we analyzed the genome sequences of CC99 and HT99 to better understand their biology and intracellular lifestyles. The CC99 genome has a size of 2.9Mb (37.9% GC) and contains 2,651 protein-encoding genes (PEGs) while the HT99 genome has a size of 3.6Mb (39.4% GC) and contains 3,238 PEGs. Both bacteria encode high proportions of hypothetical proteins (CC99: 46.5%; HT99: 51.3%). The central metabolic pathways of both bacteria appear largely intact. Genes for enzymes involved in the glycolytic pathway, the non-oxidative branch of the phosphate pathway, the tricarboxylic acid pathway, and the respiratory chain were present. Both bacteria, however, are missing genes for the synthesis of several amino acids, suggesting reliance on their host for amino acids and intermediates. Genes for type I and type IV (dot/icm) secretion systems as well as type IV pili were identified in both bacteria. Moreover, both bacteria contain genes encoding large numbers of putative effector proteins, including several with eukaryotic-like domains such as, ankyrin repeats, tetratricopeptide repeats, and leucine-rich repeats, characteristic of other intracellular bacteria.

Temporal Analysis of Gene Expression and Isoform Switching in Brown Bears (Ursus arctos)

Hibernation in brown bears is an annual process involving multiple physiologically distinct seasons-hibernation, active, and hyperphagia. While recent studies have characterized broad patterns of differential gene regulation and isoform usage between hibernation and active seasons, patterns of gene and isoform expression during hyperphagia remain relatively poorly understood. The hyperphagia stage occurs between active and hibernation seasons and involves the accumulation of large fat reserves in preparation for hibernation. Here, we use time-series analyses of gene expression and isoform usage to interrogate transcriptomic regulation associated with all three seasons. We identify a large number of genes with significant differential isoform usage (DIU) across seasons and show that these patterns of isoform usage are largely tissue-specific. We also show that DIU and differential gene-level expression responses are generally non-overlapping, with only a small subset of multi-isoform genes showing evidence of both gene-level expression changes and changes in isoform usage across seasons. Additionally, we investigate nuanced regulation of candidate genes involved in the insulin signaling pathway and find evidence of hyperphagia-specific gene expression and isoform regulation that may enhance fat accumulation during hyperphagia. Our findings highlight the value of using temporal analyses of both gene- and isoform-level gene expression when interrogating complex physiological phenotypes and provide new insight into the mechanisms underlying seasonal changes in bear physiology.

Phylogenetic, Structural and Functional Evolution of the LHC Gene Family in Plant Species

Light-harvesting chlorophyll a/b-binding (LHC) superfamily proteins play a vital role in photosynthesis. Although the physiological and biochemical functions of LHC genes have been well-characterized, the structural evolution and functional differentiation of the products need to be further studied. In this paper, we report the genome-wide identification and phylogenetic analysis of LHC genes in photosynthetic organisms. A total of 1222 non-redundant members of the LHC family were identified from 42 species. According to the phylogenetic clustering of their homologues with Arabidopsis thaliana, they can be divided into four subfamilies. In the subsequent evolution of land plants, a whole-genome replication (WGD) event was the driving force for the evolution and expansion of the LHC superfamily, with its copy numbers rapidly increasing in angiosperms. The selection pressure of photosystem II sub-unit S (PsbS) and ferrochelatase (FCII) families were higher than other subfamilies. In addition, the transcriptional expression profiles of LHC gene family members in different tissues and their expression patterns under exogenous abiotic stress conditions significantly differed, and the LHC genes are highly expressed in mature leaves, which is consistent with the conclusion that LHC is mainly involved in the capture and transmission of light energy in photosynthesis. According to the expression pattern and copy number of LHC genes in land plants, we propose different evolutionary trajectories in this gene family. This study provides a basis for understanding the molecular evolutionary characteristics and evolution patterns of plant LHCs.

Identification of hub genes associated with follicle development in multiple births sheep by WGCNA

Sheep exhibit a distinct estrous cycle that includes four different phases: proestrus, estrus, late estrus, and luteal phase. As the estrous cycle repeats, follicular development regularly alternates. We thus investigated ovarian transcriptome data from each of the four phases using weighted gene co-expression network analysis (WGCNA) to identify modules, pathways, and genes essential to follicle growth and development. We clustered mRNA and long non-coding RNA (lncRNA) into different modules by WGCNA, and calculated correlation coefficients between genes and Stages of the estrous cycle. Co-expression of the black module (cor = 0.81, P<0.001) and the yellow module (cor = 0.61, P<0.04) was found to be critical for follicle growth and development. A total of 2066 genes comprising the black and yellow modules was used for functional enrichment. The results reveal that these genes are mainly enriched in Cell cycle, PI3K-Akt signaling pathway, Oocyte meiosis, Apoptosis, and other important signaling pathways. We also identified seven hub genes (BUB1B, MAD2L1, ASPM, HSD3B1, WDHD1, CENPA, and MXI1) that may play a role in follicle development. Our study may provide several important new markers allowing in depth exploration of the genetic basis for multiparous reproduction in sheep.

Metagenomics and metatranscriptomics reveal broadly distributed, active, novel methanotrophs in the Gulf of Mexico hypoxic zone and in the marine water column

The northern Gulf of Mexico (nGOM) hypoxic zone is a shallow water environment where methane, a potent greenhouse gas, fluxes from sediments to bottom water and remains trapped due to summertime stratification. When the water column is destratified, an active planktonic methanotrophic community could mitigate the efflux of methane, which accumulates to high concentrations, to the atmosphere. To investigate the possibility of such a biofilter in the nGOM hypoxic zone we performed metagenome assembly, and metagenomic and metatranscriptomic read mapping. Methane monooxygenase (pmoA) was an abundant transcript, yet few canonical methanotrophs have been reported in this environment, suggesting a role for non-canonical methanotrophs. To determine the identity of these methanotrophs, we reconstructed six novel metagenome-assembled genomes (MAGs) in the Planctomycetota, Verrucomicrobiota and one putative Latescibacterota, each with at least one pmoA gene copy. Based on ribosomal protein phylogeny, closely related microbes (mostly from Tara Oceans) and isolate genomes were selected and co-analyzed with the nGOM MAGs. Gene annotation and read mapping suggested that there is a large, diverse and unrecognized community of active aerobic methanotrophs in the nGOM hypoxic zone and in the global ocean that could mitigate methane flux to the atmosphere.

A shotgun approach to explore the bacterial diversity and a brief insight into the glycoside hydrolases of Samiti lake located in the Eastern Himalayas

BACKGROUND

The Himalayas have always been an enigma and, being biodiversity hotspots, are considered extremely important from an ecological point of view. Recent advances in studies regarding high-altitude lakes have garnered relevant importance as these habitats could harbor potential psychrophilic and psychrotrophic microbes with bio-prospective applications. Contemplating the above scenario, the present study has been undertaken to understand the diversity and the functional capacities of the microbes thriving in this lake.

RESULTS

In our present study on Samiti Lake, the abundance of Proteobacteria as the major phylum was seen in both the soil and water samples. Incase of the ABSLW (water) and ABS1 (soil) sample, 148,066 and 239,754 predicted genes, were taken for functional analysis. The KEGG analysis showed that ABSLW and ABS1 had 122,911 and 160,268, genes assigned to KO terms respectively. Whereas in case of COG functional analysis, 104,334 and 130,191 genes were assigned to different COG classes for ABSLW and ABS1 respectively. Further, on studying the glycoside hydrolases, an abundance of GH13, GH2, GH3, GH43, and GH23 in both the soil and water samples were seen.

CONCLUSION

Our study has provided a comprehensive report about the bacterial diversity and functional capacities of microbes thriving in Samiti Lake.  It has also thrown some light on the occurrence of glycoside hydrolases in this region, as they have numerous biotechnological applications in different sectors.

Toxic effects and transcriptional responses in zebrafish liver cells following perfluorooctanoic acid exposure

As a typical type of persistent organic pollutant, perfluorooctanoic acid (PFOA) is pervasive in the environment. Multiple studies have found that PFOA has hepatotoxicity, but the mechanism remains poorly understood. In this study, the toxic effects of different concentrations of PFOA on zebrafish liver cells were systematically assessed by recording cell survival, ultrastructural observations, and transcriptome analyses. The results showed that the inhibition of cell viability and the massive accumulation of autophagic vacuoles were observed at 400 µM PFOA, while transcriptomic changes occurred with treatments of 1 and 400 µM PFOA. The transcription levels of 1055 (977 up- and 78 down-regulated genes) and 520 (446 up- and 74 down-regulated genes) genes were significantly changed after treatment with 1 and 400 µM PFOA, respectively. Based on Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis, significant expression changes were observed in autophagy, tight junction, signal transduction, immune system, endocrine system, and metabolism-related pathways, indicating that such processes were greatly affected by PFOA exposure. The findings of this study will provide a scientific basis for the toxic effects and potential toxic mechanisms of PFOA on zebrafish, and provide information for ecological risk assessments.

A common transcriptional mechanism involving R-loop and RNA abasic site regulates an enhancer RNA of APOE

RNA is modified by hundreds of chemical reactions and folds into innumerable shapes. However, the regulatory role of RNA sequence and structure and how dysregulation leads to diseases remain largely unknown. Here, we uncovered a mechanism where RNA abasic sites in R-loops regulate transcription by pausing RNA polymerase II. We found an enhancer RNA, AANCR, that regulates the transcription and expression of apolipoprotein E (APOE). In some human cells such as fibroblasts, AANCR is folded into an R-loop and modified by N-glycosidic cleavage; in this form, AANCR is a partially transcribed nonfunctional enhancer and APOE is not expressed. In contrast, in other cell types including hepatocytes and under stress, AANCR does not form a stable R-loop as its sequence is not modified, so it is transcribed into a full-length enhancer that promotes APOE expression. DNA sequence variants in AANCR are associated significantly with APOE expression and Alzheimer's Disease, thus AANCR is a modifier of Alzheimer's Disease. Besides AANCR, thousands of noncoding RNAs are regulated by abasic sites in R-loops. Together our data reveal the essentiality of the folding and modification of RNA in cellular regulation and demonstrate that dysregulation underlies common complex diseases such as Alzheimer's disease.

A high-quality genome assembly and annotation of Quercus acutissima Carruth

INTRODUCTION

Quercus acutissima is an economic and ecological tree species often used for afforestation of arid and semi-arid lands and is considered as an excellent tree for soil and water conservation.

METHODS

Here, we combined PacBio long reads, Hi-C, and Illumina short reads to assemble Q. acutissima genome.

RESULTS

We generated a 957.1 Mb genome with a contig N50 of 1.2 Mb and scaffold N50 of 77.0 Mb. The repetitive sequences constituted 55.63% of the genome, among which long terminal repeats were the majority and accounted for 23.07% of the genome. Ab initio, homology-based and RNA sequence-based gene prediction identified 29,889 protein-coding genes, of which 82.6% could be functionally annotated. Phylogenetic analysis showed that Q. acutissima and Q. variabilis were differentiated around 3.6 million years ago, and showed no evidence of species-specific whole genome duplication.

CONCLUSION

The assembled and annotated high-quality Q. acutissima genome not only promises to accelerate the species molecular biology studies and breeding, but also promotes genome level evolutionary studies.

Combined Transcriptome and Metabolome Profiling Provide Insights into Cold Responses in Rapeseed (Brassica napus L.) Genotypes with Contrasting Cold-Stress Sensitivity

Low temperature is a major environmental factor, which limits rapeseed (Brassica napus L.) growth, development, and productivity. So far, the physiological and molecular mechanisms of rapeseed responses to cold stress are not fully understood. Here, we explored the transcriptome and metabolome profiles of two rapeseed genotypes with contrasting cold responses, i.e., XY15 (cold-sensitive) and GX74 (cold-tolerant). The global metabolome profiling detected 545 metabolites in siliques of both genotypes before (CK) and after cold-stress treatment (LW). The contents of several sugar metabolites were affected by cold stress with the most accumulated saccharides being 3-dehydro-L-threonic acid, D-xylonic acid, inositol, D-mannose, D-fructose, D-glucose, and L-glucose. A total of 1943 and 5239 differentially expressed genes were identified from the transcriptome sequencing in XY15CK_vs_XY15LW and GX74CK_vs_GX74LW, respectively. We observed that genes enriched in sugar metabolism and biosynthesis-related pathways, photosynthesis, reactive oxygen species scavenging, phytohormone, and MAPK signaling were highly expressed in GX74LW. In addition, several genes associated with cold-tolerance-related pathways, e.g., the CBF-COR pathway and MAPK signaling, were specifically expressed in GX74LW. Contrarily, genes in the above-mentioned pathways were mostly downregulated in XY15LW. Thus, our results indicate the involvement of these pathways in the differential cold-stress responses in XY15 and GX74.

Shared Patterns of Gene Expression and Protein Evolution Associated with Adaptation to Desert Environments in Rodents

Desert specialization has arisen multiple times across rodents and is often associated with a suite of convergent phenotypes, including modification of the kidneys to mitigate water loss. However, the extent to which phenotypic convergence in desert rodents is mirrored at the molecular level is unknown. Here, we sequenced kidney mRNA and assembled transcriptomes for three pairs of rodent species to search for shared differences in gene expression and amino acid sequence associated with adaptation to deserts. We conducted phylogenetically independent comparisons between a desert specialist and a non-desert relative in three families representing ∼70 million years of evolution. Overall, patterns of gene expression faithfully recapitulated the phylogeny of these six taxa providing a strong evolutionary signal in levels of mRNA abundance. We also found that 8.6% of all genes showed shared patterns of expression divergence between desert and non-desert taxa, much of which likely reflects convergent evolution, and representing more than expected by chance under a model of independent gene evolution. In addition to these shared changes, we observed many species-pair-specific changes in gene expression indicating that instances of adaptation to deserts include a combination of unique and shared changes. Patterns of protein evolution revealed a small number of genes showing evidence of positive selection, the majority of which did not show shared changes in gene expression. Overall, our results suggest that convergent changes in gene regulation play an important role in the complex trait of desert adaptation in rodents.

A comprehensive Bioconductor ecosystem for the design of CRISPR guide RNAs across nucleases and technologies

The success of CRISPR-mediated gene perturbation studies is highly dependent on the quality of gRNAs, and several tools have been developed to enable optimal gRNA design. However, these tools are not all adaptable to the latest CRISPR modalities or nucleases, nor do they offer comprehensive annotation methods for advanced CRISPR applications. Here, we present a new ecosystem of R packages, called crisprVerse, that enables efficient gRNA design and annotation for a multitude of CRISPR technologies. This includes CRISPR knockout (CRISPRko), CRISPR activation (CRISPRa), CRISPR interference (CRISPRi), CRISPR base editing (CRISPRbe) and CRISPR knockdown (CRISPRkd). The core package, crisprDesign, offers a user-friendly and unified interface to add off-target annotations, rich gene and SNP annotations, and on- and off-target activity scores. These functionalities are enabled for any RNA- or DNA-targeting nucleases, including Cas9, Cas12, and Cas13. The crisprVerse ecosystem is open-source and deployed through the Bioconductor project ( https://github.com/crisprVerse ).

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Identification of Essential Proteins Based on Local Random Walk and Adaptive Multi-View Multi-Label Learning

Accumulating evidences have indicated that essential proteins play vital roles in human physiological process. In recent years, although researches on prediction of essential proteins have been developing rapidly, there are as well various limitations such as unsatisfactory data suitability, low accuracy of predictive results and so on. In this manuscript, a novel method called RWAMVL was proposed to predict essential proteins based on the Random Walk and the Adaptive Multi-View multi-label Learning. In RWAMVL, considering that the inherent noise is ubiquitous in existing datasets of known protein-protein interactions (PPIs), a variety of different features including biological features of proteins and topological features of PPI networks were obtained by adopting adaptive multi-view multi-label learning first. And then, an improved random walk method was designed to detect essential proteins based on these different features. Finally, in order to verify the predictive performance of RWAMVL, intensive experiments were done to compare it with multiple state-of-the-art predictive methods under different expeditionary frameworks. And as a result, RWAMVL was proven that it can achieve better prediction accuracy than all those competitive methods, which demonstrated as well that RWAMVL may be a potential tool for prediction of key proteins in the future.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Comparative genomics and pangenomics of vancomycin-resistant and susceptible Enterococcus faecium from Irish hospitals

Introduction. Enterococcus faecium has emerged as an important nosocomial pathogen, which is increasingly difficult to treat due to the genetic acquisition of vancomycin resistance. Ireland has a recalcitrant vancomycin-resistant bloodstream infection rate compared to other developed countries.

Hypothesis/Gap statement. Vancomycin resistance rates persist amongst E. faecium isolates from Irish hospitals. The evolutionary genomics governing these trends have not been fully elucidated.

Methodology. A set of 28 vancomycin-resistant isolates was sequenced to construct a dataset alongside 61 other publicly available Irish genomes. This dataset was extensively analysed using in silico methodologies (comparative genomics, pangenomics, phylogenetics, genotypics and comparative functional analyses) to uncover distinct evolutionary, coevolutionary and clinically relevant population trends.

Results. These results suggest that a stable (in terms of genome size, GC% and number of genes), yet genetically diverse population (in terms of gene content) of E. faecium persists in Ireland with acquired resistance arising via plasmid acquisition (vanA) or, to a lesser extent, chromosomal recombination (vanB). Population analysis revealed five clusters with one cluster partitioned into four clades which transcend isolation dates. Pangenomic and recombination analyses revealed an open (whole genome and chromosomal specific) pangenome illustrating a rampant evolutionary pattern. Comparative resistomics and virulomics uncovered distinct chromosomal and mobilomal propensity for multidrug resistance, widespread chromosomal point-mutation-mediated resistance and chromosomally harboured arsenals of virulence factors. Interestingly, a potential difference in biofilm formation strategies was highlighted by coevolutionary analysis, suggesting differential biofilm genotypes between vanA and vanB isolates.

Conclusions. These results highlight the evolutionary history of Irish E. faecium isolates and may provide insight into underlying infection dynamics in a clinical setting. Due to the apparent ease of vancomycin resistance acquisition over time, susceptible E. faecium should be concurrently reduced in Irish hospitals to mitigate potential resistant infections.

The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB

BACKGROUND

In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens. Pyrenophora tritici-repentis (Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements.

RESULTS

A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr's adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr's effectors: ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb 'Starship' transposon (dubbed 'Horizon') with a clearly defined target site and target site duplications. 'Horizon' was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer of ToxA) was nested within this newly identified Starship. Additionally, ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative 'Starship' (dubbed 'Icarus') in a ToxB-producing isolate. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and 'Icarus' were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical inheritance. Finally, the genome architecture of Ptr was defined as 'one-compartment' based on calculated gene distances and evolutionary rates.

CONCLUSIONS

These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr.

Genome sequencing and comparative analysis of Ficus benghalensis and Ficus religiosa species reveal evolutionary mechanisms of longevity

Ficus benghalensis and Ficus religiosa are large woody trees well known for their long lifespan, ecological and traditional significance, and medicinal properties. To understand the genomic and evolutionary aspects of these characteristics, the whole genomes of these Ficus species were sequenced using 10x Genomics linked reads and Oxford Nanopore long reads. The draft genomes of F. benghalensis and F. religiosa comprised of 392.89 Mbp and 332.97 Mbp, respectively. We established the genome-wide phylogenetic positions of the two Ficus species with respect to 50 other Angiosperm species. Comparative evolutionary analyses with other phylogenetically closer Eudicot species revealed adaptive evolution in genes involved in key cellular mechanisms associated with prolonged survival including phytohormones signaling, senescence, disease resistance, and abiotic stress tolerance, which provide genomic insights into the mechanisms conferring longevity and suggest that longevity is a multifaceted phenomenon. This study also provides clues on the existence of CAM pathway in these Ficus species.

Multiomics Profiling of Toxins in the Venom of the Amazonian Spider Acanthoscurria juruenicola

Acanthoscurria juruenicola is an Amazonian spider described for the first time almost a century ago. However, little is known about their venom composition. Here, we present a multiomics characterization of A. juruenicola venom by a combination of transcriptomics, proteomics, and peptidomics approaches. Transcriptomics of female venom glands resulted in 93,979 unique assembled mRNA transcript encoding proteins. A total of 92 proteins were identified in the venom by mass spectrometry, including 14 mature cysteine-rich peptides (CRPs). Quantitative analysis showed that CRPs, cysteine-rich secretory proteins, metalloproteases, carbonic anhydrases, and hyaluronidase comprise >90% of the venom proteome. Relative quantification of venom toxins was performed by DIA and DDA, revealing converging profiles of female and male specimens by both methods. Biochemical assays confirmed the presence of active hyaluronidases, phospholipases, and proteases in the venom. Moreover, the venom promoted in vivo paralytic activities in crickets, consistent with the high concentration of CRPs. Overall, we report a comprehensive analysis of the arsenal of toxins of A. juruenicola and highlight their potential biotechnological and pharmacological applications. Mass spectrometry data were deposited to the ProteomeXchange Consortium via the PRIDE repository with the dataset identifier PXD013149 and via the MassIVE repository with the dataset identifier MSV000087777.

mTOR Contributes to the Proteome Diversity through Transcriptome-Wide Alternative Splicing

The mammalian target of rapamycin (mTOR) pathway is crucial in energy metabolism and cell proliferation. Previously, we reported transcriptome-wide 3′-untranslated region (UTR) shortening by alternative polyadenylation upon mTOR activation and its impact on the proteome. Here, we further interrogated the mTOR-activated transcriptome and found that hyperactivation of mTOR promotes transcriptome-wide exon skipping/exclusion, producing short isoform transcripts from genes. This widespread exon skipping confers multifarious regulations in the mTOR-controlled functional proteomics: AS in coding regions widely affects the protein length and functional domains. They also alter the half-life of proteins and affect the regulatory post-translational modifications. Among the RNA processing factors differentially regulated by mTOR signaling, we found that SRSF3 mechanistically facilitates exon skipping in the mTOR-activated transcriptome. This study reveals a role of mTOR in AS regulation and demonstrates that widespread AS is a multifaceted modulator of the mTOR-regulated functional proteome.

Ribovirus classification by a polymerase barcode sequence

RNA viruses encoding a polymerase gene (riboviruses) dominate the known eukaryotic virome. High-throughput sequencing is revealing a wealth of new riboviruses known only from sequence, precluding classification by traditional taxonomic methods. Sequence classification is often based on polymerase sequences, but standardised methods to support this approach are currently lacking. To address this need, we describe the polymerase palmprint, a segment of the palm sub-domain robustly delineated by well-conserved catalytic motifs. We present an algorithm, Palmscan, which identifies palmprints in nucleotide and amino acid sequences; PALMdb, a collection of palmprints derived from public sequence databases; and palmID, a public website implementing palmprint identification, search, and annotation. Together, these methods demonstrate a proof-of-concept workflow for high-throughput characterisation of RNA viruses, paving the path for the continued rapid growth in RNA virus discovery anticipated in the coming decade.

Antibacterial T6SS effectors with a VRR-Nuc domain are structure-specific nucleases

The type VI secretion system (T6SS) secretes antibacterial effectors into target competitors. Salmonella spp. encode five phylogenetically distinct T6SSs. Here, we characterize the function of the SPI-22 T6SS of Salmonella bongori showing that it has antibacterial activity and identify a group of antibacterial T6SS effectors (TseV1-4) containing an N-terminal PAAR-like domain and a C-terminal VRR-Nuc domain encoded next to cognate immunity proteins with a DUF3396 domain (TsiV1-4). TseV2 and TseV3 are toxic when expressed in Escherichia coli and bacterial competition assays confirm that TseV2 and TseV3 are secreted by the SPI-22 T6SS. Phylogenetic analysis reveals that TseV1-4 are evolutionarily related to enzymes involved in DNA repair. TseV3 recognizes specific DNA structures and preferentially cleave splayed arms, generating DNA double-strand breaks and inducing the SOS response in target cells. The crystal structure of the TseV3:TsiV3 complex reveals that the immunity protein likely blocks the effector interaction with the DNA substrate. These results expand our knowledge on the function of Salmonella pathogenicity islands, the evolution of toxins used in biological conflicts, and the endogenous mechanisms regulating the activity of these toxins.

Different skeletal protein toolkits achieve similar structure and performance in the tropical coral Stylophora pistillata and the temperate Oculina patagonica

Stony corals (order: Scleractinia) differ in growth form and structure. While stony corals have gained the ability to form their aragonite skeleton once in their evolution, the suite of proteins involved in skeletogenesis is different for different coral species. This led to the conclusion that the organic portion of their skeleton can undergo rapid evolutionary changes by independently evolving new biomineralization-related proteins. Here, we used liquid chromatography-tandem mass spectrometry to sequence skeletogenic proteins extracted from the encrusting temperate coral Oculina patagonica. We compare it to the previously published skeletal proteome of the branching subtropical corals Stylophora pistillata as both are regarded as highly resilient to environmental changes. We further characterized the skeletal organic matrix (OM) composition of both taxa and tested their effects on the mineral formation using a series of overgrowth experiments on calcite seeds. We found that each species utilizes a different set of proteins containing different amino acid compositions and achieve a different morphology modification capacity on calcite overgrowth. Our results further support the hypothesis that the different coral taxa utilize a species-specific protein set comprised of independent gene co-option to construct their own unique organic matrix framework. While the protein set differs between species, the specific predicted roles of the whole set appear to underline similar functional roles. They include assisting in forming the extracellular matrix, nucleation of the mineral and cell signaling. Nevertheless, the different composition might be the reason for the varying organization of the mineral growth in the presence of a particular skeletal OM, ultimately forming their distinct morphologies.

Genome-Wide Investigation and Expression Analysis of the Nitraria sibirica Pall. CIPK Gene Family

The calcineurin B-like-interacting protein kinase (CIPK) protein family plays a key role in the plant calcium ion-mediated signal transduction pathway, which regulates a plant's response to abiotic stress. Nitraria sibirica pall. (N. sibirica) is a halophyte with a strong tolerance for high salt environments, yet how it is able to deal with salt stress on a molecular level is still unknown. Due to their function as described in other plant species, CIPK genes are prime candidates for a role in salt stress signaling in N. sibirica. In this study, we identified and analyzed the phylogenetic makeup and gene expression of the N. sibirica CIPK gene family. A total of 14 CIPKs were identified from the N. sibirica genome and were clustered into seven groups based on their phylogeny. The promoters of NsCIPK genes contained multiple elements involved in hormonal and stress response. Synteny analysis identified a total of three pairs of synteny relationships between NsCIPK genes. Each gene showed its own specific expression pattern across different tissues, with the overall expression of CIPK6 being the lowest, and that of CIPK20 being the highest. Almost all CIPK genes tended to respond to salt, drought, and cold stress, but with different sensitivity levels. In this study, we have provided a general description of the NsCIPK gene family and its expression, which will be of great significance for further understanding of the NsCIPK gene family function.

Disease resistance in coral is mediated by distinct adaptive and plastic gene expression profiles

Infectious diseases are an increasing threat to coral reefs, resulting in altered community structure and hindering the functional contributions of disease-susceptible species. We exposed seven reef-building coral species from the Caribbean to white plague disease and determined processes involved in (i) lesion progression, (ii) within-species gene expression plasticity, and (iii) expression-level adaptation among species that lead to differences in disease risk. Gene expression networks enriched in immune genes and cytoskeletal arrangement processes were correlated to lesion progression rates. Whether or not a coral developed a lesion was mediated by plasticity in genes involved in extracellular matrix maintenance, autophagy, and apoptosis, while resistant coral species had constitutively higher expression of intracellular protein trafficking. This study offers insight into the process involved in lesion progression and within- and between-species dynamics that lead to differences in disease risk that is evident on current Caribbean reefs.

Identification of candidate sex-specific genomic regions in male and female Asian arowana genomes

BACKGROUND

Asian arowana, Scleropages formosus, is one of the most expensive aquarium fish species worldwide. Its sex, however, cannot be distinguished clearly at any development stage, which impedes captive breeding and species protection for this endangered aquarium fish.

RESULTS

To discover molecular clues to the sex of Asian arowana, we sequenced 26.5 Gb of PacBio HiFi reads and 179.2 Gb of Hi-C reads for 1 male fish and also sequenced 106.5 Gb of Illumina reads, 36.0 Gb of PacBio Sequel reads, and 80.7 Gb of Hi-C reads for 1 female individual. The final male and female genome assemblies were approximately 756.8 Mb and 781.5 Mb in length and contained 25,262 and 25,328 protein-coding genes, respectively. We also resequenced the genomes of 15 male and 15 female individuals with approximately 722.1 Gb of Illumina reads. A genome-wide association study identified several potentially divergent regions between male and female individuals. In these regions, cd48 and cfap52 could be candidate genes for sex determination of Asian arowana. We also found some structural variations in few chromosomes between male and female individuals.

CONCLUSION

We provided an improved reference genome assembly of female arowana and generated the first sequenced genome of 1 male individual. These valuable genetic resources and resequencing data may improve global aquarium fish research.

Evolutionary relationship analysis of STARD gene family and VvSTARD5 improves tolerance of salt stress in transgenic tomatoes

The steroidogenic acute regulatory protein-related lipid transfer domain (STARD) forms a protein that can bind membrane-derived phospholipid second messengers and plasma membranes. Although it has been reported in many plants, the evolutionary relationship of the STARD gene family has not been systematically analyzed, and functions of the HD-START and HD-START-MEKHLA domain subgroup genes under hormone and abiotic stress are also unclear in grapes. This study identified and analyzed 23 VvSTARD genes, which were distinctly divided into five subgroups according to five conserved domain types. The analyses of codon preference, selective pressure, and synteny relationship revealed that grape had higher homology with Arabidopsis compared with rice. Interestingly, the expression levels of VvSTARD genes in subgroups 1, 2, and 3 exhibited significant upregulation under NaCl treatment at 24 h, but VvSTARD genes in subgroups 4 and 5 were upregulated under methyl jasmonate (MeJA) treatment at 24 h. The subcellular localization showed that VvSTARD5 was localized in the nucleus. Additionally, under NaCl treatment at 24 h, there were an obvious decrease in the relative electrical leakages and the content of malondialdehyde (MDA), while the relative expression level of VvSTARD5 and content of proline were obviously enhanced in three transgenic lines. Therefore, the overexpression of VvSTARD5 greatly increased the salt tolerance of transgenic tomatoes. Collectively, this study preliminarily explores the comprehensive function of the STARD gene family in grapes and verifies the function of VvSTARD5 in response to salt.

Complete genome sequence and phylogenetic analysis of medicinal plant Abrus cantoniensis for evolutionary research and germplasm utilization

Abrus cantoniensis Hance, a native medicinal plant in southern China, is officially recorded in the Chinese Pharmacopoeia. Here, we presented the first high-quality genome in Abrus genus, A. cantoniensis genome, as well as the detailed genomic information. The assembled genome size was 381.27 Mb with a scaffold N50 of 18.95 Mb, and 98.97% of the assembled sequences were anchored on 11 pseudochromosomes. The A. cantoniensis genome comprised 25,058 protein-coding genes and 45.12% of the assemblies were repetitive sequences. Comparative genome analysis suggested that chromosome translocation and inversion played an important role in the differentiation of Abrus. In addition, 24 toxin-related genes were identified, which formed two tandem gene clusters on chromosomes 2 and 3. The chromosome-level genome of A. cantoniensis obtained in this work provides a valuable resource for understanding the evolution, active ingredient biosynthesis, and genetic improvement for A. cantoniensis and Abrus species.

Vitamin B12-auxotrophy in dinoflagellates caused by incomplete or absent cobalamin-independent methionine synthase genes (metE)

Dinoflagellates are responsible for most marine harmful algal blooms (HABs) and play vital roles in many ocean processes. More than 90% of dinoflagellates are vitamin B12 auxotrophs and that B12 availability can control dinoflagellate HABs, yet the genetic basis of B12 auxotrophy in dinoflagellates in the framework of the ecology of dinoflagellates and particularly HABs, which was the objective of this work. Here, we investigated the presence, phylogeny, and transcription of two methionine synthase genes (B12-dependent metH and B12-independent metE) via searching and assembling transcripts and genes from transcriptomic and genomic databases, cloning 38 cDNA isoforms of the two genes from 14 strains of dinoflagellates, measuring the expression at different scenarios of B12, and comprehensive phylogenetic analyses of more than 100 organisms. We found that 1) metH was present in all 58 dinoflagellates accessible and metE was present in 40 of 58 species, 2) all metE genes lacked N-terminal domains, 3) metE of dinoflagellates were phylogenetically distinct from other known metE genes, and 4) expression of metH in dinoflagellates was responsive to exogenous B12 levels while expression of metE was not responding as that of genuine metE genes. We conclude that most, hypothetically all, dinoflagellates have either non-functional metE genes lacking N-terminal domain for most species, or do not possess metE for other species, which provides the genetic basis for the widespread nature of B12 auxotrophy in dinoflagellates. The work elucidated a fundamental aspect of the nutritional ecology of dinoflagellates.

Analysis of the Arabidopsis coilin mutant reveals a positive role of AtCOILIN in plant immunity

Biogenesis of ribonucleoproteins occurs in dynamic subnuclear compartments called Cajal bodies (CBs). COILIN is a critical scaffolding component essential for CB formation, composition, and activity. We recently showed that Arabidopsis (Arabidopsis thaliana) AtCOILIN is phosphorylated in response to bacterial elicitor treatment. Here, we further investigated the role of AtCOILIN in plant innate immunity. Atcoilin mutants are compromised in defense responses to bacterial pathogens. Besides confirming a role of AtCOILIN in alternative splicing (AS), Atcoilin showed differential expression of genes that are distinct from those of AS, including factors involved in RNA biogenesis, metabolism, plant immunity, and phytohormones. Atcoilin mutant plants have reduced levels of defense phytohormones. As expected, the mutant plants were more sensitive to the necrotrophic fungal pathogen Botrytis cinerea. Our findings reveal an important role for AtCOILIN in innate plant immunity.

Expression plasticity regulates intraspecific variation in the acclimatization potential of a reef-building coral

Phenotypic plasticity is an important ecological and evolutionary response for organisms experiencing environmental change, but the ubiquity of this capacity within coral species and across symbiont communities is unknown. We exposed ten genotypes of the reef-building coral Montipora capitata with divergent symbiont communities to four thermal pre-exposure profiles and quantified gene expression before stress testing 4 months later. Here we show two pre-exposure profiles significantly enhance thermal tolerance despite broadly different expression patterns and substantial variation in acclimatization potential based on coral genotype. There was no relationship between a genotype's basal thermal sensitivity and ability to acquire heat tolerance, including in corals harboring naturally tolerant symbionts, which illustrates the potential for additive improvements in coral response to climate change. These results represent durable improvements from short-term stress hardening of reef-building corals and substantial cryptic complexity in the capacity for plasticity.

Genomic analysis reveals the role of integrative and conjugative elements in plant pathogenic bacteria

Background

ICEs are mobile genetic elements found integrated into bacterial chromosomes that can excise and be transferred to a new cell. They play an important role in horizontal gene transmission and carry accessory genes that may provide interesting phenotypes for the bacteria. Here, we seek to research the presence and the role of ICEs in 300 genomes of phytopathogenic bacteria with the greatest scientific and economic impact.

Results

Seventy-eight ICEs (45 distinct elements) were identified and characterized in chromosomes of Agrobacterium tumefaciens, Dickeya dadantii, and D. solani, Pectobacterium carotovorum and P. atrosepticum, Pseudomonas syringae, Ralstonia solanacearum Species Complex, and Xanthomonas campestris. Intriguingly, the co-occurrence of four ICEs was observed in some P. syringae strains. Moreover, we identified 31 novel elements, carrying 396 accessory genes with potential influence on virulence and fitness, such as genes coding for functions related to T3SS, cell wall degradation and resistance to heavy metals. We also present the analysis of previously reported data on the expression of cargo genes related to the virulence of P. atrosepticum ICEs, which evidences the role of these genes in the infection process of tobacco plants.

Conclusions

Altogether, this paper has highlighted the potential of ICEs to affect the pathogenicity and lifestyle of these phytopathogens and direct the spread of significant putative virulence genes in phytopathogenic bacteria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00275-1.

Proteomic and Transcriptomic Analyses in the Slipper Snail Crepidula fornicata Uncover Shell Matrix Genes Expressed During Adult and Larval Biomineralization

The gastropod shell is a composite composed of minerals and shell matrix proteins (SMPs). SMPs have been identified by proteomics in many molluscs, but few have been studied in detail. Open questions include (1) what gene regulatory networks regulate SMP expression, (2) what roles individual SMPs play in biomineralization, and (3) how the complement of SMPs changes over development. These questions are best addressed in a species in which gene perturbation studies are available; one such species is the slipper snail, Crepidula fornicata. Here, SEM and pXRD analysis demonstrated that the adult shell of C. fornicata exhibits crossed lamellar microstructure and is composed of aragonite. Using high-throughput proteomics we identified 185 SMPs occluded within the adult shell. Over half of the proteins in the shell proteome have known biomineralization domains, while at least 10% have no homologs in public databases. Differential gene expression analysis identified 20 SMP genes that are up-regulated in the shell-producing mantle tissue. Over half of these 20 SMPs are expressed during development with two, CfSMP1 and CfSMP2, expressed exclusively in the shell gland. Together, the description of the shell microstructure and a list of SMPs now sets the stage for studying the consequences of SMP gene knockdowns in molluscs.

A large-scale analysis of codon usage bias in 4868 bacterial genomes shows association of codon adaptation index with GC content, protein functional domains and bacterial phenotypes

Multiple synonymous codons code for the same amino acid, resulting in the degeneracy of the genetic code and in the preferred used of some codons called codon bias usage (CBU). We performed a large-scale analysis of codon usage bias analysing the distribution of the codon adaptation index (CAI) and the codon relative adaptiveness index (RA) in 4868 bacterial genomes. We found that CAI values differ significantly between protein functional domains and part of the protein outside domains and show how CAI, GC content and preferred usage of polymerase III alpha subunits are related. Additionally, we give evidence of the association between CAI and bacterial phenotypes.

The Genome of the Mimosoid Legume Prosopis cineraria, a Desert Tree

The mimosoid legumes are a clade of ~40 genera in the Caesalpinioideae subfamily of the Fabaceae that grow in tropical and subtropical regions. Unlike the better studied Papilionoideae, there are few genomic resources within this legume group. The tree Prosopis cineraria is native to the Near East and Indian subcontinent, where it thrives in very hot desert environments. To develop a tool to better understand desert plant adaptation mechanisms, we sequenced the P. cineraria genome to near-chromosomal assembly, with a total sequence length of ~691 Mb. We predicted 77,579 gene models (76,554 CDS, 361 rRNAs and 664 tRNAs) from the assembled genome, among them 55,325 (~72%) protein-coding genes that were functionally annotated. This genome was found to consist of over 58% repeat sequences, primarily long terminal repeats (LTR-)-retrotransposons. We find an expansion of terpenoid metabolism genes in P. cineraria and its relative Prosopis alba, but not in other legumes. We also observed an amplification of NBS-LRR disease-resistance genes correlated with LTR-associated retrotransposition, and identified 410 retrogenes with an active burst of chimeric retrogene creation that approximately occurred at the same time of divergence of P. cineraria from a common lineage with P. alba~23 Mya. These retrogenes include many biotic defense responses and abiotic stress stimulus responses, as well as the early Nodulin 93 gene. Nodulin 93 gene amplification is consistent with an adaptive response of the species to the low nitrogen in arid desert soil. Consistent with these results, our differentially expressed genes show a tissue specific expression of isoprenoid pathways in shoots, but not in roots, as well as important genes involved in abiotic salt stress in both tissues. Overall, the genome sequence of P. cineraria enriches our understanding of the genomic mechanisms of its disease resistance and abiotic stress tolerance. Thus, it is a very important step in crop and legume improvement.