InterProScan: protein domains identifier

Genomics of Secondarily Temperate Adaptation in the Only Non-Antarctic Icefish

White-blooded Antarctic icefishes, a family within the adaptive radiation of Antarctic notothenioid fishes, are an example of extreme biological specialization to both the chronic cold of the Southern Ocean and life without hemoglobin. As a result, icefishes display derived physiology that limits them to the cold and highly oxygenated Antarctic waters. Against these constraints, remarkably one species, the pike icefish Champsocephalus esox, successfully colonized temperate South American waters. To study the genetic mechanisms underlying secondarily temperate adaptation in icefishes, we generated chromosome-level genome assemblies of both C. esox and its Antarctic sister species, Champsocephalus gunnari. The C. esox genome is similar in structure and organization to that of its Antarctic congener; however, we observe evidence of chromosomal rearrangements coinciding with regions of elevated genetic divergence in pike icefish populations. We also find several key biological pathways under selection, including genes related to mitochondria and vision, highlighting candidates behind temperate adaptation in C. esox. Substantial antifreeze glycoprotein (AFGP) pseudogenization has occurred in the pike icefish, likely due to relaxed selection following ancestral escape from Antarctica. The canonical AFGP locus organization is conserved in C. esox and C. gunnari, but both show a translocation of two AFGP copies to a separate locus, previously unobserved in cryonotothenioids. Altogether, the study of this secondarily temperate species provides an insight into the mechanisms underlying adaptation to ecologically disparate environments in this otherwise highly specialized group.

A chromosome-level genome assembly of the beet armyworm Spodoptera exigua

BACKGROUND

The beet armyworm Spodoptera exigua is a polyphagous caterpillar that causes serious damage to many species of crops and vegetables. To gain insight into how this polyphagous insect differs from less harmful oligophagous species, we generated a chromosome-level assembly and compared it to closely related species with the same or different feeding habits.

RESULTS

Based on Illumina and Pacific Biosciences data and Hi-C technology, 425.6 Mb of genome sequences were anchored and oriented into 31 linkage groups, with an N50 length of 14.8 Mb. A total of 24,649 gene models were predicted, of which 97.4% were identified in the genome assembly. Chemosensory genes are vital for locating food: of the four main families, odorant-binding proteins, chemosensory proteins and olfactory receptors showed little difference, whereas gustatory receptors are greatly expanded in S. exigua. Examination of other polyphagous insects confirmed this difference from oligophagous congeners and further identified the bitter receptor subfamily as being particularly affected.

CONCLUSION

Our high-quality genome sequence for beet armyworm identified a key expansion of the bitter gustatory receptor subfamily in this and other pests that differs crucially from more benign relatives and offers insight into the biology and possible future means of control for these economically important insects.

Comparative bioinformatics analysis and abiotic stress responses of expansin proteins in Cucurbitaceae members: watermelon and melon

Watermelon and melon are members of the Cucurbitaceae family including economically significant crops in the world. The expansin protein family, which is one of the members of the cell wall, breaks down the non-covalent bonds between cell wall polysaccharides, causing pressure-dependent cell expansion. Comparative bioinformatics and molecular characterization analysis of the expansin protein family were carried out in the watermelon (Citrullus lanatus) and melon (Cucumis melo) plants in the study. Gene expression levels of expansin family members were analyzed in leaf and root tissues of watermelon and melon under ABA, drought, heat, cold, and salt stress conditions by quantitative real-time PCR analysis. After comprehensive searches, 40 expansin proteins (22 ClaEXPA, 14 ClaEXPLA, and 4 ClaEXPB) in watermelon and 43 expansin proteins (19 CmEXPA, 15 CmEXPLA, 3 CmEXPB, and 6 CmEXPLB) in melon were identified. The greatest orthologous genes were identified with soybean expansin genes for watermelon and melon. However, the latest divergence time between orthologous genes was determined with poplar expansin genes for watermelon and melon expansin genes. ClaEXPA-04, ClaEXPA-09, ClaEXPB-01, ClaEXPB-03, and ClaEXPLA-13 genes in watermelon and CmEXPA-12, CmEXPA-10, and CmEXPLA-01 genes in melon can be involved in tissue development and abiotic stress response of the plant. The current study combining bioinformatics and experimental analysis can provide a detailed characterization of the expansin superfamily which has roles in growth and reaction to the stress of the plant. The study ensures detailed data for future studies examining gene functions including the roles in plant growth and stress conditions.

Multi-omics reveal differentiation and maintenance of dimorphic flowers in an alpine plant on the Qinghai-Tibet Plateau

Dimorphic flowers growing on a single individual plant play a critical role in extreme adaption and reproductive assurance in plants and have high ecological and evolutionary significance. However, the omics bases underlying such a differentiation and maintenance remain largely unknown. We aimed to investigate this through genomic, transcriptome and metabolomic analyses of dimorphic flowers in an alpine biennial, Sinoswertia tetraptera (Gentianaceae). A high-quality chromosome-level genome sequence (903 Mb) was first assembled for S. tetraptera with 31,359 protein-coding genes annotated. Two rounds of recent independent whole-genome duplication (WGD) were revealed. Numerous genes from the recent species-specific WGD were found to be differentially expressed in the two types of flowers, and this may have helped contribute to the origin of this innovative trait. The genes with contrasting expressions between flowers were related to biosynthesis of hormones, floral pigments (carotenoids and flavonoids) and iridoid compounds, which are involved in both flower development and colour. Metabolomic analyses similarly suggested differential concentrations of these chemicals in the two types of flowers. The expression interactions between multiple genes may together lead to contrasting morphology and chemical concentration and open versus closed pollination of the dimorphic flowers in this species for reproductive assurance.

Identification of potential proteins translated from circular RNA splice variants

Circular RNAs (circRNAs) are covalently closed RNA molecules generated from precursor RNAs by the head-to-tail backsplicing of exons. Hundreds of studies demonstrated that circRNAs are ubiquitously expressed and regulate cellular events by modulating microRNA (miRNA) and RNA-binding protein (RBP) activities. A few circRNAs are also known to translate into functional polypeptides regulating cellular physiology. All these functions primarily depend on the full-length sequence of the circRNAs. CircRNA backsplice junction sequence is the key to identifying circRNAs and their full-length mature sequence. However, some multi-exonic circRNAs exist in different isoforms sharing identical backsplice junction sequences and are termed circRNA splice variants. Here, we analyzed the previously published HeLa cell RNA-seq datasets to identify circRNA splice variants using the de novo module of the CIRCexplorer2 circRNA annotation pipeline. A subset of circRNAs with splice variants was validated by the circRNA-rolling circle amplification (circRNA-RCA) method. Interestingly, several validated circRNAs were predicted to translate into proteins by the riboCIRC database. Furthermore, polyribosome fractionation followed by quantitative PCR confirmed the association of a subset of circRNAs with polyribosome supporting their protein-coding potential. Finally, bioinformatics analysis of proteins derived from splice variants of circCORO1C and circASPH suggested altered protein sequences and structures that could affect their physiological functions. Together, our study identified novel circRNA splice variants and their potential translation into protein isoforms which may regulate various physiological processes.

cotH Genes Are Necessary for Normal Spore Formation and Virulence in Mucor lusitanicus

Mucormycosis is an invasive fungal infection caused by certain members of the fungal order of Mucorales. The species most frequently identified as the etiological agents of mucormycosis belong to the genera Rhizopus, Lichtheimia, and Mucor. The frequency of systemic mucormycosis has been increasing, mainly because of increasing numbers of susceptible patients. Furthermore, Mucorales display intrinsic resistance to the majority of routinely used antifungal agents (e.g., echinocandins and short-tailed azoles), which limits the number of possible therapeutic options. All the above-mentioned issues urge the improvement of molecular identification methods and the discovery of new antifungal targets and strategies. Spore coat proteins (CotH) constitute a kinase family present in many pathogenic bacteria and fungi and participate in the spore formation in these organisms. Moreover, some of them can act as virulence factors being receptors of the human GRP78 protein during Rhizopus delemar-induced mucormycosis. We identified 17 cotH-like genes in the Mucor lusitanicus genome database. Successful disruption of five cotH genes in Mucor was performed using the CRISPR-Cas9 system. The CotH3 and CotH4 proteins play a role in adaptation to different temperatures as well as in developing the cell wall structure. We also show CotH4 protein is involved in spore wall formation by affecting the total chitin content and, thus, the composition of the spore wall. The role of CotH3 and CotH4 proteins in virulence was confirmed in two invertebrate models and a diabetic ketoacidosis (DKA) mouse model. IMPORTANCE Current treatment options for mucormycosis are inadequate, resulting in high mortality rates, especially among immunosuppressed patients. The development of novel therapies for mucormycosis has been hampered by lack of understanding of the pathogenetic mechanisms. The importance of the cell surface CotH proteins in the pathogenesis of Rhizopus-mediated mucormycosis has been recently described. However, the contribution of this family of proteins to the virulence of other mucoralean fungi and their functionality in vital processes remain undefined. Through the use of the CRISPR-Case9 gene disruption system, we demonstrate the importance of several of the CotH proteins to the virulence of Mucor lusitanicus by using three infection models. We also report on the importance of one of these proteins, CotH4, to spore wall formation by affecting chitin content. Therefore, our studies extend the importance of CotH proteins to Mucor and identify the mechanism by which one of the CotH proteins contributes to the development of a normal fungal cell wall, thereby indicating that this family of proteins can be targeted for future development of novel therapeutic strategies of mucormycosis.

Diversity of the Antimicrobial Peptide Genes in Collembola

Multidrug-resistant bacteria are a current health crisis threatening the world's population, and scientists are looking for new drugs to combat them. Antimicrobial peptides (AMPs), which are part of the organism's innate immune system, are a promising new drug class as they can disrupt bacterial cell membranes. This study explored antimicrobial peptide genes in collembola, a non-insect hexapod lineage that has survived in microbe-rich habitats for millions of years, and their antimicrobial peptides have not been thoroughly investigated. We used in silico analysis (homology-based gene identification, physicochemical and antimicrobial property prediction) to identify AMP genes from the genomes and transcriptomes of five collembola representing three main suborders: Entomobryomorpha (Orchesella cincta, Sinella curviseta), Poduromorpha (Holacanthella duospinosa, Anurida maritima), and Symphypleona (Sminthurus viridis). We identified 45 genes belonging to five AMP families, including (a) cysteine-rich peptides: diapausin, defensin, and Alo; (b) linear α-helical peptide without cysteine: cecropin; (c) glycine-rich peptide: diptericin. Frequent gene gains and losses were observed in their evolution. Based on the functions of their orthologs in insects, these AMPs potentially have broad activity against bacteria, fungi, and viruses. This study provides candidate collembolan AMPs for further functional analysis that could lead to medicinal use.

Citrus genomic resources unravel putative genetic determinants of Huanglongbing pathogenicity

Citrus HLB caused by Candidatus Liberibacter asiaticus is a pathogen-triggered immune disease. Here, we identified putative genetic determinants of HLB pathogenicity by integrating citrus genomic resources to characterize the pan-genome of accessions that differ in their response to HLB. Genome-wide association mapping and analysis of allele-specific expression between susceptible, tolerant, and resistant accessions further refined candidates underlying the response to HLB. We first developed a phased diploid assembly of Citrus sinensis 'Newhall' genome and produced resequencing data for 91 citrus accessions that differ in their response to HLB. These data were combined with previous resequencing data from 356 accessions for genome-wide association mapping of the HLB response. Genes determinants for HLB pathogenicity were associated with host immune response, ROS production, and antioxidants. Overall, this study has provided a significant resource of citrus genomic data and identified candidate genes to be further explored to understand the genetic determinants of HLB pathogenicity.

A novel binary pesticidal protein from Chryseobacterium arthrosphaerae controls western corn rootworm by a different mode of action to existing commercial pesticidal proteins

The western corn rootworm (WCR) Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae) remains one of the economically most important pests of maize (Zea mays) due to its adaptive capabilities to pest management options. This includes the ability to develop resistance to some of the commercial pesticidal proteins originating from different strains of Bacillus thuringiensis. Although urgently needed, the discovery of new, environmentally safe agents with new modes of action is a challenge. In this study we report the discovery of a new family of binary pesticidal proteins isolated from several Chryseobacterium species. These novel binary proteins, referred to as GDI0005A and GDI0006A, produced as recombinant proteins, prevent growth and increase mortality of WCR larvae, as does the bacteria. These effects were found both in susceptible and resistant WCR colonies to Cry3Bb1 and Cry34Ab1/Cry35Ab1 (reassigned Gpp34Ab1/Tpp35Ab1). This suggests GDI0005A and GDI0006A may not share the same binding sites as those commercially deployed proteins and thereby possess a new mode of action. This paves the way towards the development of novel biological or biotechnological management solutions urgently needed against rootworms.

Early antitermination in the atypical coliphage mEp021 mediated by the Gp17 protein

The coliphage mEp021 belongs to a phage group with a unique immunity repressor, and its life cycle requires the host factor Nus. mEp021 has been classified as non-lambdoid based on its specific characteristics. The mEp021 genome carries a gene encoding an N_λ-like antiterminator protein, termed Gp17, and three nut sites (nutL, nutR1, and nutR2). Analysis of plasmid constructs containing these nut sites, a transcription terminator, and a GFP reporter gene showed high levels of fluorescence when Gp17 was expressed, but not in its absence. Like lambdoid N proteins, Gp17 has an arginine-rich motif (ARM), and mutations in its arginine codons inhibit its function. In infection assays using the mutant phage mEp021ΔGp17::Kan (where gp17 has been deleted), gene transcripts located downstream of transcription terminators were obtained only when Gp17 was expressed. In contrast to phage lambda, mEp021 virus particle production was partially restored (>1/3 relative to wild type) when nus mutants (nusA1, nusB5, nusC60, and nusE71) were infected with mEp021 and Gp17 was overexpressed. Our results suggest that RNA polymerase reads through the third nut site (nutR2), which is more than 7.9 kbp downstream of nutR1.

Genome-Wide Identification of G Protein-Coupled Receptors in Ciliated Eukaryotes

G protein-coupled receptors (GPCRs) are the largest family of transmembrane receptors and play important roles in many physiological processes. As a representative group of protozoa, ciliates represent the highest stage of eukaryotic cell differentiation and evolution in terms of their reproductive mode, two-state karyotype, and extremely diverse cytogenesis patterns. GPCRs have been poorly reported in ciliates. In this study, we identified 492 GPCRs in 24 ciliates. Using the existing classification system for animals, GPCRs in ciliates can be assigned to four families, including families A, B, E, and F. Most (377 members) belong to family A. The number of GPCRs is extremely different in different ciliates; the Heterotrichea ciliates usually have more GPCRs than other ciliates. Parasitic or symbiotic ciliates usually have only a few GPCRs. Gene/genome duplication events seem to play important roles in the expansion of the GPCR superfamily in ciliates. GPCRs in ciliates displayed seven typical domain organizations. GPCRs in an ortholog group are common and conserved in all ciliates. The gene expression analysis of the members in this conserved ortholog group in the model ciliate, Tetrahymena thermophila, suggested that these GPCRs play important roles in the life cycle of ciliates. In summary, this study provides the first comprehensive genome-wide identification of GPCRs in ciliates, improving our understanding of the evolution and function of GPCR in ciliates.

GUN Mutants: New Weapons To Unravel Ascospore Germination Regulation in the Model Fungus Podospora anserina

In Podospora anserina as in many other Ascomycetes, ascospore germination is a regulated process that requires the breaking of dormancy. Despite its importance in survival and dispersal, ascospore germination in filamentous fungi has been poorly investigated, and little is known about its regulation and genetic control. We have designed a positive genetic screen that led to the isolation of mutants showing uncontrolled germination, the GUN (Germination UNcontrolled) mutants. Here, we report on the characterization of the gun1SG (Spontaneous Germination) mutant. We show that gun1SG is mutated in Pa_6_1340, the ortholog of Magnaporthe oryzae Pth2, which encodes a carnitine-acetyltransferase (CAT) involved in the shuttling of acetyl coenzyme A between peroxisomes and mitochondria and which is required for appressorium development. Bioinformatic analysis revealed that the mutated residue (I441) is highly conserved among Fungi and that the mutation has a deleterious impact on the protein function. We show that GUN1 is essential for ascospore germination and that the protein is localized both in mitochondria and in peroxisomes. Finally, epistasis studies allowed us to place GUN1 together with the PaMpk2 MAPK pathway upstream of the PaNox2/PaPls1 complex in the regulation of ascospore germination. In addition, we show that GUN1 plays a role in appressorium functioning. The pivotal role of GUN1, the ortholog of Pth2, in ascospore germination and in appressorium functioning reinforces the idea of a common genetic regulation governing both appressorium development and melanized ascospore germination. Furthermore, we characterize the second CAT encoded in P. anserina genome, Pa_3_7660/GUP1, and we show that the function of both CATs is conserved in P. anserina. IMPORTANCE The regulation of ascospore germination in filamentous fungi has been poorly investigated so far. To unravel new genes involved in this regulation pathway, we conducted a genetic screen in Podospora anserina, and we isolated 57 mutants affected in ascospore germination. Here, we describe the Germination UNcontrolled One (gun1SG) mutant, and we characterize the gene affected. GUN1 is a peroxisomal/mitochondrial carnitine-acetyltransferase required for acetyl coenzyme A shuttling between both organelles, and we show that GUN1 is a pleiotropic gene also involved in appressorium functioning similarly to its ortholog, the pathogenesis factor Pth2, in the plant pathogen Magnaporthe oryzae. Given the similarities in the regulation of appressorium development and ascospore germination, we speculate that discovering new genes controlling ascospore germination in P. anserina may lead to the discovery of new pathogenesis factors in pathogenic fungi. The characterization of GUN1, the ortholog of M. oryzae Pth2, represents a proof of concept.

Uncovering the involvement of DoDELLA1-interacting proteins in development by characterizing the DoDELLA gene family in Dendrobium officinale

BACKGROUND

Gibberellins (GAs) are widely involved in plant growth and development. DELLA proteins are key regulators of plant development and a negative regulatory factor of GA. Dendrobium officinale is a valuable traditional Chinese medicine, but little is known about D. officinale DELLA proteins. Assessing the function of D. officinale DELLA proteins would provide an understanding of their roles in this orchid's development.

RESULTS

In this study, the D. officinale DELLA gene family was identified. The function of DoDELLA1 was analyzed in detail. qRT-PCR analysis showed that the expression levels of all DoDELLA genes were significantly up-regulated in multiple shoots and GA3-treated leaves. DoDELLA1 and DoDELLA3 were significantly up-regulated in response to salt stress but were significantly down-regulated under drought stress. DoDELLA1 was localized in the nucleus. A strong interaction was observed between DoDELLA1 and DoMYB39 or DoMYB308, but a weak interaction with DoWAT1.

CONCLUSIONS

In D. officinale, a developmental regulatory network involves a close link between DELLA and other key proteins in this orchid's life cycle. DELLA plays a crucial role in D. officinale development.

Genome-wide identification and expression analysis of the HD2 protein family and its response to drought and salt stress in Gossypium species

Histone deacetylase 2 (HD2) proteins play an important role in the regulation of gene expression. This helps with the growth and development of plants and also plays a crucial role in responses to biotic and abiotic stress es. HD2s comprise a C₂H₂-type Zn²⁺ finger at their C-terminal and an HD2 label, deacetylation and phosphorylation sites, and NLS motifs at their N-terminal. In this study, a total of 27 HD2 members were identified, using Hidden Markov model profiles, in two diploid cotton genomes (Gossypium raimondii and Gossypium arboretum) and two tetraploid cotton genomes (Gossypium hirsutum and Gossypium barbadense). These cotton HD2 members were classified into 10 major phylogenetic groups (I-X), of which group III was found to be the largest with 13 cotton HD2 members. An evolutionary investigation showed that the expansion of HD2 members primarily occurred as a result of segmental duplication in paralogous gene pairs. Further qRT-PCR validation of nine putative genes using RNA-Seq data suggested that GhHDT3D.2 exhibits significantly higher levels of expression at 12h, 24h, 48h, and 72h of exposure to both drought and salt stress conditions compared to a control measure at 0h. Furthermore, gene ontology, pathways, and co-expression network study of GhHDT3D.2 gene affirmed their significance in drought and salt stress responses.

Comparative Modeling and Analysis of Extremophilic D-Ala-D-Ala Carboxypeptidases

Understanding the molecular adaptations of organisms to extreme environments requires a comparative analysis of protein structure, function, and dynamics across species found in different environmental conditions. Computational studies can be particularly useful in this pursuit, allowing exploratory studies of large numbers of proteins under different thermal and chemical conditions that would be infeasible to carry out experimentally. Here, we perform such a study of the MEROPS family S11, S12, and S13 proteases from psychophilic, mesophilic, and thermophilic bacteria. Using a combination of protein structure prediction, atomistic molecular dynamics, and trajectory analysis, we examine both conserved features and trends across thermal groups. Our findings suggest a number of hypotheses for experimental investigation.

Wound healing and regeneration in the reef building coral Acropora millepora

Branching scleractinian corals are niche-constructing organisms, providing continuously-growing, structural foundation for spectacularly biodiverse coral reef ecosystems. A large part of their success lies in the ability to quickly regenerate following mechanical damage. Even now, when the corals undergo great decline due to anthropogenic weather and storm extremes, it is surprising how little is known about molecular mechanisms governing regeneration in these iconic organisms. In this study, we used RNA-seq to identify genes involved in the regeneration of Acropora millepora, starting with the initial wound closure up to complete rebuilding of lost structures. Many of the differentially expressed genes we found in the wound healing steps are homologues of genes known to be involved in wound healing and regeneration of bilaterian and other cnidarian species, prominently including multiple components of FGF and Wnt signalling pathways. Comparison between genes involved in wound healing and continuous growth of the colony demonstrates both similarity and distinctiveness of the genetic programmes controlling these processes. A striking example is specific expression of c-Fos, a transcription factor with conserved role in early injury response, during the earliest stages of wound healing of A. millepora. By comparing results obtained in diverse experimental conditions including a closed-loop, recirculating aquarium and a flow-through system of marine station, we have demonstrated feasibility of using zooxanthellate scleractinian corals as experimental models in fundamental biology research, including studies of regeneration.

Identification of AHL Synthase in Desulfovibrio vulgaris Hildenborough Using an In-Silico Methodology

Sulfate-reducing bacteria (SRB) are anaerobic bacteria that form biofilm and induce corrosion on various material surfaces. The quorum sensing (QS) system that employs acyl homoserine lactone (AHL)-type QS molecules primarily govern biofilm formation. Studies on SRB have reported the presence of AHL, but no AHL synthase have been annotated in SRB so far. In this computational study, we used a combination of data mining, multiple sequence alignment (MSA), homology modeling and docking to decode a putative AHL synthase in the model SRB, Desulfovibrio vulgaris Hildenborough (DvH). Through data mining, we shortlisted 111 AHL synthase genes. Conserved domain analysis of 111 AHL synthase genes generated a consensus sequence. Subsequent MSA of the consensus sequence with DvH genome indicated that DVU_2486 (previously uncharacterized protein from acetyltransferase family) is the gene encoding for AHL synthase. Homology modeling revealed the existence of seven α-helices and six β sheets in the DvH AHL synthase. The amalgamated study of hydrophobicity, binding energy, and tunnels and cavities revealed that Leu99, Trp104, Arg139, Trp97, and Tyr36 are the crucial amino acids that govern the catalytic center of this putative synthase. Identifying AHL synthase in DvH would provide more comprehensive knowledge on QS mechanism and help design strategies to control biofilm formation.

Two haplotype-resolved, gap-free genome assemblies for Actinidia latifolia and Actinidia chinensis shed light on the regulatory mechanisms of vitamin C and sucrose metabolism in kiwifruit

Kiwifruit is a recently domesticated horticultural fruit crop with substantial economic and nutritional value, especially because of the high content of vitamin C in its fruit. In this study, we de novo assembled two telomere-to-telomere kiwifruit genomes from Actinidia chinensis var. 'Donghong' (DH) and Actinidia latifolia 'Kuoye' (KY), with total lengths of 608 327 852 and 640 561 626 bp for 29 chromosomes, respectively. With a burst of structural variants involving inversion, translocations, and duplications within 8.39 million years, the metabolite content of DH and KY exhibited differences in saccharides, lignans, and vitamins. A regulatory ERF098 transcription factor family has expanded in KY and Actinidia eriantha, both of which have ultra-high vitamin C content. With each assembly phased into two complete haplotypes, we identified allelic variations between two sets of haplotypes, leading to protein sequence variations in 26 494 and 27 773 gene loci and allele-specific expression of 4687 and 12 238 homozygous gene pairs. Synchronized metabolome and transcriptome changes during DH fruit development revealed the same dynamic patterns in expression levels and metabolite contents; free fatty acids and flavonols accumulated in the early stages, but sugar substances and amino acids accumulated in the late stages. The AcSWEET9b gene that exhibits allelic dominance was further identified to positively correlate with high sucrose content in fruit. Compared with wild varieties and other Actinidia species, AcSWEET9b promoters were selected in red-flesh kiwifruits that have increased fruit sucrose content, providing a possible explanation on why red-flesh kiwifruits are sweeter. Collectively, these two gap-free kiwifruit genomes provide a valuable genetic resource for investigating domestication mechanisms and genome-based breeding of kiwifruit.

Evolution of the WRKY66 Gene Family and Its Mutations Generated by the CRISPR/Cas9 System Increase the Sensitivity to Salt Stress in Arabidopsis

Group Ⅲ WRKY transcription factors (TFs) play pivotal roles in responding to the diverse abiotic stress and secondary metabolism of plants. However, the evolution and function of WRKY66 remains unclear. Here, WRKY66 homologs were traced back to the origin of terrestrial plants and found to have been subjected to both motifs' gain and loss, and purifying selection. A phylogenetic analysis showed that 145 WRKY66 genes could be divided into three main clades (Clade A-C). The substitution rate tests indicated that the WRKY66 lineage was significantly different from others. A sequence analysis displayed that the WRKY66 homologs had conserved WRKY and C2HC motifs with higher proportions of crucial amino acid residues in the average abundance. The AtWRKY66 is a nuclear protein, salt- and ABA- inducible transcription activator. Simultaneously, under salt stress and ABA treatments, the superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activities, as well as the seed germination rates of Atwrky66-knockdown plants generated by the clustered, regularly interspaced, short palindromic repeats/CRISPR-associated 9 (CRISPR/Cas9) system, were all lower than those of wild type (WT) plants, but the relative electrolyte leakage (REL) was higher, indicating the increased sensitivities of the knockdown plants to the salt stress and ABA treatments. Moreover, RNA-seq and qRT-PCR analyses revealed that several regulatory genes in the ABA-mediated signaling pathway involved in stress response of the knockdown plants were significantly regulated, being evidenced by the more moderate expressions of the genes. Therefore, the AtWRKY66 likely acts as a positive regulator in the salt stress response, which may be involved in an ABA-mediated signaling pathway.

Pangenomic analysis of Wolbachia provides insight into the evolution of host adaptation and cytoplasmic incompatibility factor genes

Introduction

The genus Wolbachia provides a typical example of intracellular bacteria that infect the germline of arthropods and filarial nematodes worldwide. Their importance as biological regulators of invertebrates, so it is particularly important to study the evolution, divergence and host adaptation of these bacteria at the genome-wide level.

Methods

Here, we used publicly available Wolbachia genomes to reconstruct their evolutionary history and explore their adaptation under host selection.

Results

Our findings indicate that segmental and single-gene duplications, such as DNA methylase, bZIP transcription factor, heat shock protein 90, in single monophyletic Wolbachia lineages (including supergroups A and B) may be responsible for improving the ability to adapt to a broad host range in arthropod-infecting strains. In contrast to A strains, high genetic diversity and rapidly evolving gene families occur in B strains, which may promote the ability of supergroup B strains to adapt to new hosts and their large-scale spreading. In addition, we hypothesize that there might have been two independent horizontal transfer events of cif genes in two sublineages of supergroup A strains. Interestingly, during the independent evolution of supergroup A and B strains, the rapid evolution of cif genes in supergroup B strains resulted in the loss of their functional domain, reflected in a possible decrease in the proportion of induced cytoplasmic incompatibility (CI) strains.

Discussion

This present study highlights for reconstructing of evolutionary history, addressing host adaptation-related evolution and exploring the origin and divergence of CI genes in each Wolbachia supergroup. Our results thus not only provide a basis for further exploring the evolutionary history of Wolbachia adaptation under host selection but also reveal a new research direction for studying the molecular regulation of Wolbachia- induced cytoplasmic incompatibility.

Computational models for prediction of protein-protein interaction in rice and Magnaporthe grisea

INTRODUCTION

Plant-microbe interactions play a vital role in the development of strategies to manage pathogen-induced destructive diseases that cause enormous crop losses every year. Rice blast is one of the severe diseases to rice Oryza sativa (O. sativa) due to Magnaporthe grisea (M. grisea) fungus. Protein-protein interaction (PPI) between rice and fungus plays a key role in causing rice blast disease.

METHODS

In this paper, four genomic information-based models such as (i) the interolog, (ii) the domain, (iii) the gene ontology, and (iv) the phylogenetic-based model are developed for predicting the interaction between O. sativa and M. grisea in a whole-genome scale.

RESULTS AND DISCUSSION

A total of 59,430 interacting pairs between 1,801 rice proteins and 135 blast fungus proteins are obtained from the four models. Furthermore, a machine learning model is developed to assess the predicted interactions. Using composition-based amino acid composition (AAC) and conjoint triad (CT) features, an accuracy of 88% and 89% is achieved, respectively. When tested on the experimental dataset, the CT feature provides the highest accuracy of 95%. Furthermore, the specificity of the model is verified with other pathogen-host datasets where less accuracy is obtained, which confirmed that the model is specific to O. sativa and M. grisea. Understanding the molecular processes behind rice resistance to blast fungus begins with the identification of PPIs, and these predicted PPIs will be useful for drug design in the plant science community.

Integrated Transcriptome and Proteome Analysis Reveals that the Antimicrobial Griseofulvin Targets Didymella segeticola Beta-Tubulin to Control Tea Leaf Spot

Because effective control measures are lacking, tea leaf spot caused by Didymella segeticola results in huge tea (Camellia sinensis) production losses on tea plantations in Guizhou Province, southwestern China. Screening for natural antimicrobial agents with higher control effects against this pathogen and studying their modes of action may contribute to disease management. Here, Penicillium griseofulvum-derived antimicrobial griseofulvin (GSF) can inhibit the hyphal growth of D. segeticola strain GZSQ-4, with a half-maximal effective concentration of 0.37 μg/ml in vitro and a higher curative efficacy at a lower dose of 25 μg/ml for detached tea twigs. GSF induces deformed and slightly curly hyphae with enlarged ends, with protoplasts agglutinated in the hyphae, and higher numbers of hyphal protuberances. GSF alters hyphal morphology and the subcellular structure's order. The integrated transcriptome and proteome data revealed that the transport of materials in cells, cellular movement, and mitosis were modulated by GSF. Molecular docking indicated that beta-tubulin was the most potent target of GSF, with a binding free energy of -13.59 kcal/mol, and microscale thermophoresis indicated that the dissociation constant (Kd) value of GSF binding to beta-tubulin 1, compared with beta-tubulin 2, was significantly lower. Thus, GSF potentially targets beta-tubulin 1 to disturb the chromosomal separation and fungal mitosis, thereby inhibiting hyphal growth.

Functional Characterization of Potato UBC13-UEV1s Genes Required for Ubiquitin Lys63 Chain to Polyubiquitination

Ubiquitin-conjugating enzymes (E2s/UBC) are components of the ubiquitin proteasome system (UPS), and the ubiquitin-conjugating enzyme variant (UEV) is one of E2s (ubiquitin-conjugating enzymes, UBC) subfamily. The UEVs and UBC13 play an auxiliary role in mediating Lys63-linked polyUb chain assembly, which is correlated with target protein non-proteolytic functions, such as DNA repair or response to stress. However, the collaborative mechanism of StUBC13 (homologue of AtUBC13) and StUEVs (the UEVs in potato) involved in potato are not fully understood understood. Here, we identified two StUBC13 and seven StUEVs from potato genome. We analyzed protein motif and conserved domain, gene structure, phylogenetic features, cis-acting elements of StUBC13 and StUEVs. Subsequently, we screened StUBC13 partners protein and verified interaction between StUBC13 and StUEVs using yeast two-hybrid, split luciferase complementation (SLC) and bimolecular fluorescence complementation (BiFC) approach. The expression profile and qRT-PCR analysis suggested that StUBC13 and StUEVs gene exhibited a tissue-specific expression and were induced by different stress. Overall, this investigative study provides a comprehensive reference and view for further functional research on StUBC13 and StUEV1s in potato.

Phylogenomic analysis of Wolbachia genomes from the Darwin Tree of Life biodiversity genomics project

The Darwin Tree of Life (DToL) project aims to sequence all described terrestrial and aquatic eukaryotic species found in Britain and Ireland. Reference genome sequences are generated from single individuals for each target species. In addition to the target genome, sequenced samples often contain genetic material from microbiomes, endosymbionts, parasites, and other cobionts. Wolbachia endosymbiotic bacteria are found in a diversity of terrestrial arthropods and nematodes, with supergroups A and B the most common in insects. We identified and assembled 110 complete Wolbachia genomes from 93 host species spanning 92 families by filtering data from 368 insect species generated by the DToL project. From 15 infected species, we assembled more than one Wolbachia genome, including cases where individuals carried simultaneous supergroup A and B infections. Different insect orders had distinct patterns of infection, with Lepidopteran hosts mostly infected with supergroup B, while infections in Diptera and Hymenoptera were dominated by A-type Wolbachia. Other than these large-scale order-level associations, host and Wolbachia phylogenies revealed no (or very limited) cophylogeny. This points to the occurrence of frequent host switching events, including between insect orders, in the evolutionary history of the Wolbachia pandemic. While supergroup A and B genomes had distinct GC% and GC skew, and B genomes had a larger core gene set and tended to be longer, it was the abundance of copies of bacteriophage WO who was a strong determinant of Wolbachia genome size. Mining raw genome data generated for reference genome assemblies is a robust way of identifying and analysing cobiont genomes and giving greater ecological context for their hosts.

Olfactory and gustatory receptor genes in fig wasps: Evolutionary insights from comparative studies

The mechanisms of chemoreception in fig wasps (Hymenoptera, Agaonidae) are of primary importance in their co-evolutionary relationship with the fig trees they pollinate. We used transcriptome sequences of 25 fig wasps in six genera that allowed a comparative approach to the evolution of key molecular components of fig wasp chemoreception: their odorant (OR) and gustatory (GR) receptor genes. In total, we identified 311 ORs and 47 GRs, with each species recording from 5 to 30 OR genes and 1-4 GR genes. 304 OR genes clustered into 18 orthologous groups known to be sensitive to cuticular hydrocarbons (CHC), pheromones, acids, alcohols and a variety of floral scents such as cineole, Linalool, and Heptanone. 45 GR genes clustered into 4 orthologous groups that contain sweet, bitter, CO₂ and undocumented receptors. Gene sequences in most orthologous groups varied greatly among species, except for ORco (60.0% conserved) and sweet receptors (30.7% conserved). Strong purifying selection of both odorant and gustatory genes was detected, as shown by low ω values. Signatures of positive selection were detected in loci from both OR and GR orthologous groups. Fig wasps have relatively few olfactory and especially gustatory receptors, reflecting the natural history of the system. Amino acid sequences nonetheless vary significantly between species and are consistent with the phylogenetic relationships among fig wasps. The differences in ORs within some orthologous groups from the same species, but different hosts and from closely related species from one host can reach as low as 49.3% and 9.8% respectively, implying the ORs of fig wasps can evolve rapidly to novel ecological environments. Our results provide a starting point for understanding the molecular basis of the chemosensory systems of fig wasps.

Metalloproteome of human-infective RNA viruses: a study towards understanding the role of metal ions in virology

Metalloproteins and metal-based inhibitors have been shown to effectively combat infectious diseases, particularly those caused by RNA viruses. In this study, a diverse set of bioinformatics methods was employed to identify metal-binding proteins of human RNA viruses. Seventy-three viral proteins with a high probability of being metal-binding proteins were identified. These proteins included 40 zinc-, 47 magnesium- and 14 manganese-binding proteins belonging to 29 viral species and eight significant viral families, including Coronaviridae, Flaviviridae and Retroviridae. Further functional characterization has revealed that these proteins play a critical role in several viral processes, including viral replication, fusion and host viral entry. They fall under the essential categories of viral proteins, including polymerase and protease enzymes. Magnesium ion is abundantly predicted to interact with these viral enzymes, followed by zinc. In addition, this study also examined the evolutionary aspects of predicted viral metalloproteins, offering essential insights into the metal utilization patterns among different viral species. The analysis indicates that the metal utilization patterns are conserved within the functional classes of the proteins. In conclusion, the findings of this study provide significant knowledge on viral metalloproteins that can serve as a valuable foundation for future research in this area.

The Clostridium-infecting filamentous phage CAK1 genome analysis allows to define a new potential clade of Tubulavirales

What we know about Tubulavirales, i.e. filamentous phages, essentially comes from Gram-negative-infecting Inoviridae. However, metagenomics recently suggests filamentous phages are much more widespread and diverse. Here, we report the complete sequence and functional annotation of CAK1, a 6.6 kb filamentous phage that was shown to chronically infect Clostridium beijerinckii 30 years ago and only represents the second filamentous phage cultivated on a Gram-positive bacterium. CAK1 has a typical filamentous phage modular genome with no homologs in databases and we were interested to compare it with a pig gut filamentous phage metagenomics dataset that we previously assembled and for which many filamentous phages were predicted to infect Clostridium species by bioinformatics means. CAK1 is distantly related to nine of these sequences, two of which have been predicted as Clostridium-associated. In itself, this small cluster of CAK1-connected sequences sheds light on the diversity of filamentous phages that putatively infect Clostridium species, and probably many other Gram-positive genera.

Does integument arise de novo or from pre-existing structures? ── Insights from the key regulatory genes controlling integument development

The origin of seeds is one of the key innovations in land plant evolution. Ovules are the developmental precursors of seeds. The integument is the envelope structure surrounding the nucellus within the ovule and developing into the seed coat when ovules mature upon fertilization. The question of whether the integument arise de novo or evolve from elaboration of pre-existing structures has caused much debate. By exploring the origin and evolution of the key regulatory genes controlling integument development and their functions during both individual and historical developmental processes, we showed the widespread presence of the homologs of ANT, CUC, BEL1, SPL, C3HDZ, INO, ATS, and ETT in seedless plant genomes. All of these genes have undergone duplication-divergence events in their evolutionary history, with most of the descendant paralogous suffering motif gain and/or loss in the coding regions. Expression and functional characterization have shown that these genes are key components of the genetic program that patterns leaf-like lateral organs. Serial homology can thus be postulated between integuments and other lateral organs in terms of the shared master regulatory genes. Given that the genetic program patterning leaf-like lateral organs formed in seedless plants, and was reused during seed origin, the integument is unlikely to arise de novo but evolved from the stem segment-specific modification of pre-existing serially homologous structures. The master 'switches' trigging the modification to specify the integument identity remain unclear. We propose a successive transformation model of integument origin.

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect

BACKGROUND

Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown.

RESULTS

A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts.

CONCLUSIONS

Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

Genome Characterization of Bacteriophage KPP-1, a Novel Member in the Subfamily Vequintavirinae, and Use of Its Endolysin for the Lysis of Multidrug-Resistant Klebsiella variicola In Vitro

Multidrug-resistant members of the Klebsiella pneumoniae complex have become a threat to human lives and animals, including aquatic animals, owing to the limited choice of antimicrobial treatments. Bacteriophages are effective natural tools available to fight against multidrug-resistant bacteria. The bacteriophage KPP-1 was found to be strictly lytic against K. variicola, a multidrug-resistant isolate, producing clear plaques. The genome sequence analysis of KPP-1 revealed that it comprised 143,369 base pairs with 47% overall GC content. A total of 272 genes (forward 161, complementary 111) encode for 17 tRNAs and 255 open reading frames (ORFs). Among them, 32 ORFs could be functionally annotated using the National Center for Biotechnology Information (NCBI) Protein Basic Local Alignment Search Tool (BLASTp) algorithm while 223 were found to code for hypothetical proteins. Comparative genomic analysis revealed that the closest neighbor of KPP-1 can be found in the genus Mydovirus of the subfamily Vequintavirinae. KPP-1 not only markedly suppressed the growth of the host but also worked synergistically with ampicillin. Useful genes for pathogen control such as endolysin (locus tag: KPP_11591) were found to have activity against multidrug-resistant isolate of K. variicola. Further studies are necessary to develop a strategy to control the emerging pathogen K. variicola using bacteriophages such as KPP-1.

Transcriptional and Physiological Analysis Reveal New Insights into the Regulation of Fertilization (N, P, K) on the Growth and Synthesis of Medicinal Components of Dendrobium denneanum

Dendrobium denneanum is an important medicinal and ornamental plant. Its ornamental and medicinal values are affected by its vegetative growth conditions and chemical composition accumulation. This study adopted an orthogonal experimental design to treat D. denneanum with nine different levels of nitrogen (N), potassium (K), and phosphorus (P). The morphological indicators of the plant were positively correlated with the nitrogen concentration. The polysaccharide content was the highest at 1500 mg·L^-1 nitrogen and 3000 mg·L^-1 phosphorous and was 26.84% greater than the control. The flavonoid content increased by 36.2% at 500 mg·L^-1 nitrogen, 2000 mg·L^-1 phosphorous, and 300 mg·L^-1 potassium. Principal component score analysis showed that nitrogen had the most significant impact on the various indicators of D. denneanum, followed by phosphorus and potassium. The comprehensive score showed that the T9 treatment (N: 1500 mg·L^-1, P: 3000 mg·L^-1, K: 500 mg·L^-1) had the strongest effect on D. denneanum. Transcriptional analysis showed that compared with the control, the T9 treatment led to 2277 differentially expressed genes (1230 upregulated and 1047 downregulated). This includes fifteen genes enriched in the MAPK signaling pathway, five genes in phenylpropanoid biosynthesis, and two genes in flavonoid biosynthesis. These genes may be involved in regulating plant growth and the biosynthesis of polysaccharides and flavonoids. This study provides guidance for the optimal use of N, P, and K in the cultivation of D. denneanum.

Label-free quantitative proteomics of arbuscular mycorrhizal Elaeagnus angustifolia seedlings provides insights into salt-stress tolerance mechanisms

INTRODUCTION

Soil salinization has become one of the most serious environmental issues globally. Excessive accumulation of soluble salts will adversely affect the survival, growth, and reproduction of plants. Elaeagnus angustifolia L., commonly known as oleaster or Russian olive, has the characteristics of tolerance to drought and salt. Arbuscular mycorrhizal (AM) fungi are considered to be bio-ameliorator of saline soils that can enhance the salt tolerance of the host plants. However, there is little information on the root proteomics of AM plants under salt stress.

METHODS

In this study, a label-free quantitative proteomics method was employed to identify the differentially abundant proteins in AM E. angustifolia seedlings under salt stress.

RESULTS

The results showed that a total of 170 proteins were significantly differentially regulated in E.angustifolia seedlings after AMF inoculation under salt stress. Mycorrhizal symbiosis helps the host plant E. angustifolia to respond positively to salt stress and enhances its salt tolerance by regulating the activities of some key proteins related to amino acid metabolism, lipid metabolism, and glutathione metabolism in root tissues.

CONCLUSION

Aspartate aminotransferase, dehydratase-enolase-phosphatase 1 (DEP1), phospholipases D, diacylglycerol kinase, glycerol-3-phosphate O-acyltransferases, and gamma-glutamyl transpeptidases may play important roles in mitigating the detrimental effect of salt stress on mycorrhizal E. angustifolia . In conclusion, these findings provide new insights into the salt-stress tolerance mechanisms of AM E. angustifolia seedlings and also clarify the role of AM fungi in the molecular regulation network of E. angustifolia under salt stress.

The putative metal-binding proteome of the Coronaviridae family

Metalloproteins are well-known for playing various physicochemical processes in all life forms, including viruses. Some life-threatening viruses (such as some members of the Coronaviridae family of viruses) are emerged and remerged frequently and are rapidly transmitted throughout the globe. This study aims to identify and characterize the metal-binding proteins (MBPs) of the Coronaviridae family of viruses and further provides insight into the MBP's role in sustaining and propagating viruses inside a host cell and in the outer environment. In this study, the available proteome of the Coronaviridae family was exploited. Identified potential MBPs were analyzed for their functional domains, structural aspects, and subcellular localization. We also demonstrate phylogenetic aspects of all predicted MBPs among other Coronaviridae family members to understand the evolutionary trend among their respective hosts. A total of 256 proteins from 51 different species of coronaviruses are predicted as MBPs. These MBPs perform various key roles in the replication and survival of viruses within the host cell. Cysteine, aspartic acid, threonine, and glutamine are key amino acid residues interacting with respective metal ions. Our observations also indicate that the metalloproteins of this family of viruses circulated and evolved in different hosts, which supports the zoonotic nature of coronaviruses. The comprehensive information on MBPs of the Coronaviridae family may be further helpful in designing novel therapeutic metalloprotein targets. Moreover, the study of viral MBPs can also help to understand the roles of MBPs in virus pathogenesis and virus-host interactions.

Annotation of the Turnera subulata (Passifloraceae) Draft Genome Reveals the S-Locus Evolved after the Divergence of Turneroideae from Passifloroideae in a Stepwise Manner

A majority of Turnera species (Passifloraceae) exhibit distyly, a reproductive system involving both self-incompatibility and reciprocal herkogamy. This system differs from self-incompatibility in Passiflora species. The genetic basis of distyly in Turnera is a supergene, restricted to the S-morph, and containing three S-genes. How supergenes and distyly evolved in Turnera, and the other Angiosperm families exhibiting distyly remain largely unknown. Unraveling the evolutionary origins in Turnera requires the generation of genomic resources and extensive phylogenetic analyses. Here, we present the annotated draft genome of the S-morph of distylous Turnera subulata. Our annotation allowed for phylogenetic analyses of the three S-genes' families across 56 plant species ranging from non-seed plants to eudicots. In addition to the phylogenetic analysis, we identified the three S-genes' closest paralogs in two species of Passiflora. Our analyses suggest that the S-locus evolved after the divergence of Passiflora and Turnera. Finally, to provide insights into the neofunctionalization of the S-genes, we compared expression patterns of the S-genes with close paralogs in Arabidopsis and Populus trichocarpa. The annotation of the T. subulata genome will provide a useful resource for future comparative work. Additionally, this work has provided insights into the convergent nature of distyly and the origin of supergenes.

Chromosome-level genome of black cutworm provides novel insights into polyphagy and seasonal migration in insects

BACKGROUND

The black cutworm, Agrotis ipsilon, is a serious global underground pest. Its distinct phenotypic traits, especially its polyphagy and ability to migrate long distances, contribute to its widening distribution and increasing difficulty of control. However, knowledge about these traits is still limited.

RESULTS

We generated a high-quality chromosome-level assembly of A. ipsilon using PacBio and Hi-C technology with a contig N50 length of ~ 6.7 Mb. Comparative genomic and transcriptomic analyses showed that detoxification-associated gene families were highly expanded and induced after insects fed on specific host plants. Knockout of genes that encoded two induced ABC transporters using CRISPR/Cas9 significantly reduced larval growth rate, consistent with their contribution to host adaptation. A comparative transcriptomic analysis between tethered-flight moths and migrating moths showed expression changes in the circadian rhythm gene AiCry2 involved in sensing photoperiod variations and may receipt magnetic fields accompanied by MagR and in genes that regulate the juvenile hormone pathway and energy metabolism, all involved in migration processes.

CONCLUSIONS

This study provides valuable genomic resources for elucidating the mechanisms involved in moth migration and developing innovative control strategies.

Marmesin and Marmelosin Interact with the Heparan Sulfatase-2 Active Site: Potential Mechanism for Phytochemicals from Bael Fruit Extract as Antitumor Therapeutics

Human heparan sulfatase-2 (HSULF-2) is an oncoprotein overexpressed in the surface of all types of tumor cells and its activity plays a critical role in cancer survival and progression. Our previous studies have shown that bael fruit extract, containing marmesin and marmelosin, inhibits the HSULF-2 activity and kills breast tumor cells, but the mechanism of these processes remains fairly known mainly because the HSULF-2's 3D structure is partially known. Herein, we aimed at providing an in silico molecular mechanism of the inhibition of human HSULF-2 by phytochemicals from bael fruit extract. Pharmacokinetic parameters of the main phytochemicals contained in the bael fruit extract, sequence-based 3D structure of human HSULF-2, and the interaction of bael fruit's phytochemicals with the enzyme active site was modeled, evaluated, and verified. Docking studies revealed marmesin and marmelosin as potential inhibitors with binding score -8.5 and -7.7 Kcal/mol; these results were validated using molecular dynamics simulations, which exhibited higher stability of the protein-ligand complexes. Taking together, with our earlier in vitro data, our computational analyses suggest that marmesin and marmelosin interact at the active site of HSULF-2 providing a potential mechanism for its inhibition and consequent antitumor activity by phytochemicals contained in the bael fruit extract.

Castro LF, Froufe E. The Crown Pearl V2: an improved genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758)

Contiguous assemblies are fundamental to deciphering the composition of extant genomes. In molluscs, this is considerably challenging owing to the large size of their genomes, heterozygosity, and widespread repetitive content. Consequently, long-read sequencing technologies are fundamental for high contiguity and quality. The first genome assembly of Margaritifera margaritifera (Linnaeus, 1758) (Mollusca: Bivalvia: Unionida), a culturally relevant, widespread, and highly threatened species of freshwater mussels, was recently generated. However, the resulting genome is highly fragmented since the assembly relied on short-read approaches. Here, an improved reference genome assembly was generated using a combination of PacBio CLR long reads and Illumina paired-end short reads. This genome assembly is 2.4 Gb long, organized into 1,700 scaffolds with a contig N50 length of 3.4 Mbp. The ab initio gene prediction resulted in 48,314 protein-coding genes. Our new assembly is a substantial improvement and an essential resource for studying this species' unique biological and evolutionary features, helping promote its conservation.

A Polyketide Synthetase Gene Cluster Is Responsible for Antibacterial Activity of Burkholderia contaminans MS14

Burkholderia contaminans MS14, isolated from a soil sample in Mississippi, is known for producing the novel antifungal compound occidiofungin. In addition, MS14 exhibits a broad range of antibacterial activities against common plant pathogens. Random mutagenesis and gene complementation indicate that four genes are required for antibacterial activity of strain MS14 against the fire blight pathogen Erwinia amylovora. With the aim of finding the biosynthetic gene cluster for the unknown antibacterial compound, we used RNA-seq to analyze the transcriptome of MS14 wild type and mutants lacking antibacterial activity. The twofold lower expressed genes in all mutants were studied, and a polyketide synthase (PKS) gene cluster was predicted to be directly involved in MS14 antibacterial activities. The nptII-resistance cassette and CRISPR-Cas9 systems were used to mutate the PKS gene cluster. Plate bioassays showed that either insertion or frame-shifting one of the PKS genes resulted in a loss of antibacterial activity. Considering that the antibacterial-defective mutants maintain the same antifungal activities as the wild-type strain, the results suggest that this PKS gene cluster is highly likely to be involved in or directly responsible for the production of MS14 antibacterial activity. Purification efforts revealed that the antibacterial activity of the compound synthesized by the gene cluster is sensitive to UV radiation. Nevertheless, these findings have provided more insights to understand the antibacterial activity of strain MS14.

Genomic profiling of climate adaptation in Aedes aegypti along an altitudinal gradient in Nepal indicates nongradual expansion of the disease vector

Driven by globalization, urbanization and climate change, the distribution range of invasive vector species has expanded to previously colder ecoregions. To reduce health-threatening impacts on humans, insect vectors are extensively studied. Population genomics can reveal the genomic basis of adaptation and help to identify emerging trends of vector expansion. By applying whole genome analyses and genotype-environment associations to populations of the main dengue vector Aedes aegypti, sampled along an altitudinal gradient in Nepal (200-1300 m), we identify putatively adaptive traits and describe the species' genomic footprint of climate adaptation to colder ecoregions. We found two differentiated clusters with significantly different allele frequencies in genes associated to climate adaptation between the highland population (1300 m) and all other lowland populations (≤800 m). We revealed nonsynonymous mutations in 13 of the candidate genes associated to either altitude, precipitation or cold tolerance and identified an isolation-by-environment differentiation pattern. Other than the expected gradual differentiation along the altitudinal gradient, our results reveal a distinct genomic differentiation of the highland population. Local high-altitude adaptation could be one explanation of the population's phenotypic cold tolerance. Carrying alleles relevant for survival under colder climate increases the likelihood of this highland population to a worldwide expansion into other colder ecoregions.

Comparison of Plasma Exosome Proteomes Between Obese and Non-Obese Patients with Type 2 Diabetes Mellitus

PURPOSE

Obesity is considered a promoter of type 2 diabetes mellitus (T2DM). However, the underlying mechanism remains unclear. This study aimed to identify plasma exosome differentially expressed proteins (DEPs) that are potentially involved in the development of obesity-related T2DM.

METHODS

Exosomes were isolated from the plasma of obese and non-obese T2DM patients (n = 10 for each group). A label-free quantitative mass spectrometry analysis was applied to identify plasma exosome DEPs in obese patients compared with non-obese patients, followed by bioinformatics analysis including GO annotation, KEGG analysis, subcellular localization prediction, transcription factor analysis, and protein-protein interaction (PPI) prediction.

RESULTS

We identified 2 significantly upregulated proteins (C9 and PON1) and 5 significantly downregulated proteins (HPX, A1BG, CFHR1, ANG, and CALM) in obese patients compared with those in non-obese patients. KEGG analysis demonstrated that the insulin signaling pathway was one of the pathways that significantly correlated with the DEPs. The DEPs were primarily localized in the extracellular space (5 out of 7). HMG-box and NF-Y beta might regulate the transcription of the DEPs. C9, PON1, HPX, and CFHR1 were present in the PPI network.

CONCLUSION

The plasma exosome DEPs are potentially responsible for the development of obesity-related T2DM possibly through the insulin signaling pathway and the interaction with other proteins. Our study may guide future research direction toward the pathogenesis of obesity-related T2DM.

Functional characterization of MANNOSE-BINDING LECTIN 1, a G-type lectin gene family member, in response to fungal pathogens of strawberry

The mannose-binding lectin gene MANNOSE-BINDING LECTIN 1 (MBL1) is a member of the G-type lectin family and is involved in defense in strawberry (Fragaria × ananassa). Genome-wide identification of the G-type lectin family was carried out in woodland strawberry, F. vesca, and 133 G-lectin genes were found. Their expression profiles were retrieved from available databases and indicated that many are actively expressed during plant development or interaction with pathogens. We selected MBL1 for further investigation and generated stable transgenic FaMBL1-overexpressing plants of F. ×ananassa to examine the role of this gene in defense. Plants were selected and evaluated for their contents of disease-related phytohormones and their reaction to biotic stresses, and this revealed that jasmonic acid decreased in the overexpressing lines compared with the wild-type (WT). Petioles of the overexpressing lines inoculated with Colletotrichum fioriniae had lower disease incidence than the WT, and leaves of these lines challenged by Botrytis cinerea showed significantly smaller lesion diameters than the WT and higher expression of CLASS II CHITINASE 2-1. Our results indicate that FaMBL1 plays important roles in strawberry response to fungal diseases caused by C. fioriniae and B. cinerea.

MeGATAs, functional generalists in interactions between cassava growth and development, and abiotic stresses

The proteins with DNA-binding preference to the consensus DNA sequence (A/T) GATA (A/G) belong to a GATA transcription factor family, with a wide array of biological processes in plants. Cassava (Manihot esculenta) is an important food crop with high production of starch in storage roots. Little was however known about cassava GATA domain-containing genes (MeGATAs). Thirty-six MeGATAs, MeGATA1 to MeGATA36, were found in this study. Some MeGATAs showed a collinear relationship with orthologous genes of Arabidopsis, poplar and potato, rice, maize and sorghum. Eight MeGATA-encoded proteins (MeGATAs) analysed were all localized in the nucleus. Some MeGATAs had potentials of binding ligands and/or enzyme activity. One pair of tandem-duplicated MeGATA17-MeGATA18 and 30 pairs of whole genome-duplicated MeGATAs were found. Fourteen MeGATAs showed low or no expression in the tissues. Nine analysed MeGATAs showed expression responses to abiotic stresses and exogenous phytohormones. Three groups of MeGATA protein interactions were found. Fifty-three miRNAs which can target 18 MeGATAs were identified. Eight MeGATAs were found to target other 292 cassava genes, which were directed to radial pattern formation and phyllome development by gene ontology enrichment, and autophagy by Kyoto Encyclopaedia of Genes and Genomes enrichment. These data suggest that MeGATAs are functional generalists in interactions between cassava growth and development, abiotic stresses and starch metabolism.

Gene abundance linked to climate zone: Parallel evolution of gene content along elevation gradients in lichenized fungi

Introduction

Intraspecific genomic variability affects a species' adaptive potential toward climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many of these species have remarkable environmental tolerances, and often form populations across different climate zones. Here, we combine comparative and population genomics to assess the presence and absence of genes in high and low elevation genomes of two lichenized fungi of the genus Umbilicaria.

Methods

The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked, which genes increase or decrease in abundance within populations along elevation.

Results

In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinase superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents and found that the number of this protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations.

Discussion

We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.

Identification and analysis of the secretome of plant pathogenic fungi reveals lifestyle adaptation

The secretory proteome plays an important role in the pathogenesis of phytopathogenic fungi. However, the relationship between the large-scale secretome of phytopathogenic fungi and their lifestyle is not fully understood. In the present study, the secretomes of 150 plant pathogenic fungi were predicted and the characteristics associated with different lifestyles were investigated. In total, 94,974 secreted proteins (SPs) were predicted from these fungi. The number of the SPs ranged from 64 to 1,662. Among these fungi, hemibiotrophic fungi had the highest number (average of 970) and proportion (7.1%) of SPs. Functional annotation showed that hemibiotrophic and necrotroph fungi, differ from biotrophic and symbiotic fungi, contained much more carbohydrate enzymes, especially polysaccharide lyases and carbohydrate esterases. Furthermore, the core and lifestyle-specific SPs orthogroups were identified. Twenty-seven core orthogroups contained 16% of the total SPs and their motif function annotation was represented by serine carboxypeptidase, carboxylesterase and asparaginase. In contrast, 97 lifestyle-specific orthogroups contained only 1% of the total SPs, with diverse functions such as PAN_AP in hemibiotroph-specific and flavin monooxygenases in necrotroph-specific. Moreover, obligate biotrophic fungi had the largest number of effectors (average of 150), followed by hemibiotrophic fungi (average of 120). Among these effectors, 4,155 had known functional annotation and pectin lyase had the highest proportion in the functionally annotated effectors. In addition, 32 sets of RNA-Seq data on pathogen-host interactions were collected and the expression levels of SPs were higher than that of non-SPs, and the expression level of effector genes was higher in biotrophic and hemibiotrophic fungi than in necrotrophic fungi, while secretase genes were highly expressed in necrotrophic fungi. Finally, the secretory activity of five predicted SPs from Setosphearia turcica was experimentally verified. In conclusion, our results provide a foundation for the study of pathogen-host interaction and help us to understand the fungal lifestyle adaptation.

Assembly and analytical validation of a metagenomic reference catalog of human gut microbiota based on co-barcoding sequencing

Human gut microbiota is associated with human health and disease, and is known to have the second-largest genome in the human body. The microbiota genome is important for their functions and metabolites; however, accurate genomic access to the microbiota of the human gut is hindered due to the difficulty of cultivating and the shortcomings of sequencing technology. Therefore, we applied the stLFR library construction method to assemble the microbiota genomes and demonstrated that assembly property outperformed standard metagenome sequencing. Using the assembled genomes as references, SNP, INDEL, and HGT gene analyses were performed. The results demonstrated significant differences in the number of SNPs and INDELs among different individuals. The individual displayed a unique species variation spectrum, and the similarity of strains within individuals decreased over time. In addition, the coverage depth analysis of the stLFR method shows that a sequencing depth of 60X is sufficient for SNP calling. HGT analysis revealed that the genes involved in replication, recombination and repair, mobilome prophages, and transposons were the most transferred genes among different bacterial species in individuals. A preliminary framework for human gut microbiome studies was established using the stLFR library construction method.

Identification and functional annotation of hypothetical proteins of uropathogenic Escherichia coli strain CFT073 towards designing antimicrobial drug targets

Urinary tract infections are a serious health concern worldwide, especially in developing countries. Escherichia coli strain CFT073 is a highly virulent pathogenic bacterial strain. CFT073 proteome contains 4897 proteins, out of which 992 have been classified as hypothetical proteins. Identification and characterization of hypothetical proteins can aid in the selection of targets for drug design. In this study, we studied the hypothetical proteins from the UPEC strain CFT073 using various computational tools. By NCBI-CDD, 376 protein sequences showed conserved domains. Based on the functional motifs in their primary sequences, we classified these 376 hypothetical proteins into 7 functional categories. Further KEGG database was used to find the roles of these hypothetical proteins in several pathways. Protein interaction network analysis of hypothetical proteins identified 53 proteins as highly interacting metabolic proteins. Virulence factor analysis of the proteins identified 8 proteins as virulent. We conducted a non-homology search for the identified proteins of UPEC in the available human proteome. We observed that 35 proteins are non-homologous to humans and hence could be selected for drug designing targets. Qualitative characterization of the selected 35 non-homologous hypothetical proteins including essentiality analysis and evaluation of druggability by similarity search against drug bank database was performed. Out of these 35 proteins, three-dimensional structures of six proteins (NP_752562.1, NP_756345.1, NP_754893.1, NP_756600.2, NP_755264.1 and NP_752994.1) could be successfully modelled. These new annotations can help to better understand disease mechanisms at the molecular level, as well as provide new targets for drug development against the UPEC strain CFT073.Communicated by Ramaswamy H. Sarma.

Ontology-driven analysis of marine metagenomics: what more can we learn from our data?

BACKGROUND

The proliferation of metagenomic sequencing technologies has enabled novel insights into the functional genomic potentials and taxonomic structure of microbial communities. However, cyberinfrastructure efforts to manage and enable the reproducible analysis of sequence data have not kept pace. Thus, there is increasing recognition of the need to make metagenomic data discoverable within machine-searchable frameworks compliant with the FAIR (Findability, Accessibility, Interoperability, and Reusability) principles for data stewardship. Although a variety of metagenomic web services exist, none currently leverage the hierarchically structured terminology encoded within common life science ontologies to programmatically discover data.

RESULTS

Here, we integrate large-scale marine metagenomic datasets with community-driven life science ontologies into a novel FAIR web service. This approach enables the retrieval of data discovered by intersecting the knowledge represented within ontologies against the functional genomic potential and taxonomic structure computed from marine sequencing data. Our findings highlight various microbial functional and taxonomic patterns relevant to the ecology of prokaryotes in various aquatic environments.

CONCLUSIONS

In this work, we present and evaluate a novel Semantic Web architecture that can be used to ask novel biological questions of existing marine metagenomic datasets. Finally, the FAIR ontology searchable data products provided by our API can be leveraged by future research efforts.

Chromosome-level genome and the identification of sex chromosomes in Uloborus diversus

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

Genome assembly of an Australian native grass species reveals a recent whole-genome duplication and biased gene retention of genes involved in stress response

BACKGROUND

The adaptive significance of polyploidy has been extensively debated, and chromosome-level genome assemblies of polyploids can provide insight into this. The Australian grass Bothriochloa decipiens belongs to the BCD clade, a group with a complex history of hybridization and polyploid. This is the first genome assembly and annotation of a species that belongs to this fascinating yet complex group.

FINDINGS

Using Illumina short reads, 10X Genomics linked reads, and Hi-C sequencing data, we assembled a highly contiguous genome of B. decipiens, with a total length of 1,218.22 Mb and scaffold N50 of 42.637 Mb. Comparative analysis revealed that the species experienced a relatively recent whole-genome duplication. We clustered the 20 major scaffolds, representing the 20 chromosomes, into the 2 subgenomes of the parental species using unique repeat signatures. We found evidence of biased fractionation and differences in the activity of transposable elements between the subgenomes prior to hybridization. Duplicates were enriched for genes involved in transcription and response to external stimuli, supporting a biased retention of duplicated genes following whole-genome duplication.

CONCLUSIONS

Our results support the hypotheses of a biased retention of duplicated genes following polyploidy and point to differences in repeat activity associated with subgenome dominance. B. decipiens is a widespread species with the ability to establish across many soil types, making it a prime candidate for climate change- resilient ecological restoration of Australian grasslands. This reference genome is a valuable resource for future population genomic research on Australian grasses.

Molecular mechanisms underlying hematophagia revealed by comparative analyses of leech genomes

BACKGROUND

Leeches have been used in traditional Chinese medicine since prehistoric times to treat a spectrum of ailments, but very little is known about their physiological, genetic, and evolutionary characteristics.

FINDINGS

We sequenced and assembled chromosome-level genomes of 3 leech species (bloodsucking Hirudo nipponia and Hirudinaria manillensis and nonbloodsucking Whitmania pigra). The dynamic population histories and genome-wide expression patterns of the 2 bloodsucking leech species were found to be similar. A combined analysis of the genomic and transcriptional data revealed that the bloodsucking leeches have a presumably enhanced auditory sense for prey location in relatively deep fresh water. The copy number of genes related to anticoagulation, analgesia, and anti-inflammation increased in the bloodsucking leeches, and their gene expressions responded dynamically to the bloodsucking process. Furthermore, the expanded FBN1 gene family may help in rapid body swelling of leeches after bloodsucking, and the expanded GLB3 gene family may be associated with long-term storage of prey blood in a leech's body.

CONCLUSIONS

The high-quality reference genomes and comprehensive datasets obtained in this study may facilitate innovations in the artificial culture and strain optimization of leeches.