A honey bee symbiont buffers larvae against nutritional stress through lysine supplementation

Honey bees, the worlds most significant agricultural pollinator, have suffered dramatic losses in the last few decades. These losses are largely due to the synergistic effects of multiple stressors, the most pervasive of which is limited nutrition. The effects of poor nutrition are most damaging in the developing larvae of honey bees, who mature into workers unable to meet the needs of their colony. It is therefore essential that we better understand the nutritional landscape experienced by honey bee larvae. In this study, we characterize the metabolic capabilities of a honey bee larvae-associated bacterium, Bombella apis (formerly Parasaccharibacter apium), and its effects on the nutritional resilience of larvae. We found that B. apis is the only bacterium associated with larvae that can withstand the antimicrobial larval diet. Further, we found that B. apis can synthesize all essential amino acids and significantly alters the amino acid content of synthetic larval diet, largely by increasing the essential amino acid lysine. Analyses of gene gain/loss across the phylogeny suggest that two distinct cationic amino acid transporters were gained by B. apis ancestors, and the transporter LysE is conserved across all sequenced strains of B. apis. This result suggests that amino acid export is a key feature conserved within the Bombella clade. Finally, we tested the impact of B. apis on developing honey bee larvae subjected to nutritional stress and found that larvae supplemented with B. apis are bolstered against mass reduction despite limited nutrition. Together, these data suggest an important role of B. apis as a nutritional mutualist of honey bee larvae.

Comparative genome analysis of plant ascomycete fungal pathogens with different lifestyles reveals distinctive virulence strategies

Background Pathogens have evolved diverse lifestyles and adopted pivotal new roles in both natural ecosystems and human environments. However, the molecular mechanisms underlying their adaptation to new lifestyles are obscure. Comparative genomics was adopted to determine distinct strategies of plant ascomycete fungal pathogens with different lifestyles and to elucidate their distinctive virulence strategies. Results We found that plant ascomycete biotrophs exhibited lower gene gain and loss events and loss of CAZyme-encoding genes involved in plant cell wall degradation and biosynthesis gene clusters for the production of secondary metabolites in the genome. Comparison with the candidate effectome detected distinctive variations between plant biotrophic pathogens and other groups (including human, necrotrophic and hemibiotrophic pathogens). The results revealed the biotroph-specific and lifestyle-conserved candidate effector families. These data have been configured in web-based genome browser applications for public display (http://lab.malab.cn/soft/PFPG). This resource allows researchers to profile the genome, proteome, secretome and effectome of plant fungal pathogens. Conclusions Our findings demonstrated different genome evolution strategies of plant fungal pathogens with different lifestyles and explored their lifestyle-conserved and specific candidate effectors. It will provide a new basis for discovering the novel effectors and their pathogenic mechanisms.

Comparative Genomics and Gene Pool Analysis Reveal the Decrease of Genome Diversity and Gene Number in Rice Blast Fungi by Stable Adaption with Rice

Magnaporthe oryzae caused huge losses in rice and wheat production worldwide. Comparing to long-term co-evolution history with rice, wheat-infecting isolates were new-emerging. To reveal the genetic differences between rice and wheat blast on global genomic scale, 109 whole-genome sequences of M. oryzae from rice, wheat, and other hosts were reanalyzed in this study. We found that the rice lineage had gone through stronger selective sweep and fewer conserved genes than those of Triticum and Lolium lineages, which indicated that rice blast fungi adapted to rice by gene loss and rapid evolution of specific loci. Furthermore, 228 genes associated with host adaptation of M. oryzae were found by presence/absence variation (PAV) analyses. The functional annotation of these genes found that the fine turning of genes gain/loss involved with transport and transcription factor, thiol metabolism, and nucleotide metabolism respectively are major mechanisms for rice adaption. This result implies that genetic base of specific host plant may lead to gene gain/loss variation of pathogens, so as to enhance their adaptability to host. Further characterization of these specific loci and their roles in adaption and evaluation of the fungi may eventually lead to understanding of interaction mechanism and develop new strategies of the disease management.

Rapid methicillin resistance diversification in Staphylococcus epidermidis colonizing human neonates

Early in life, infants are colonized with multiple bacterial strains whose differences in gene content can have important health consequences. Metagenomics-based approaches have revealed gene content differences between different strains co-colonizing newborns, but less is known about the rate, mechanism, and phenotypic consequences of gene content diversification within strains. Here, focusing on Staphylococcus epidermidis, we whole-genome sequence and phenotype more than 600 isolates from newborns. Within days of birth, infants are co-colonized with a highly personalized repertoire of S. epidermidis strains, which are spread across the newborn body. Comparing the genomes of multiple isolates of each strain, we find very little evidence of adaptive evolution via single-nucleotide polymorphisms. By contrast, we observe gene content differences even between otherwise genetically identical cells, including variation of the clinically important methicillin resistance gene, mecA, suggesting rapid gene gain and loss events at rates higher than point mutations. Mapping the genomic architecture of structural variants by long-read Nanopore sequencing, we find that deleted regions were always flanked by direct repeats, consistent with site-specific recombination. However, we find that even within a single genetic background, recombination occurs at multiple, often non-canonical repeats, leading to the rapid evolution of patient-specific diverse structural variants in the SCCmec island and to differences in antibiotic resistance.

Polysaccharide-Bacteria Interactions From the Lens of Evolutionary Ecology

Microbes have the unique ability to break down the complex polysaccharides that make up the bulk of organic matter, initiating a cascade of events that leads to their recycling. Traditionally, the rate of organic matter degradation is perceived to be limited by the chemical and physical structure of polymers. Recent advances in microbial ecology, however, suggest that polysaccharide persistence can result from non-linear growth dynamics created by the coexistence of alternate degradation strategies, metabolic roles as well as by ecological interactions between microbes. This complex “landscape” of degradation strategies and interspecific interactions present in natural microbial communities appears to be far from evolutionarily stable, as frequent gene gain and loss reshape enzymatic repertoires and metabolic roles. In this perspective, we discuss six challenges at the heart of this problem, ranging from the evolution of genetic repertoires, phenotypic heterogeneity in clonal populations, the development of a trait-based ecology, and the impact of metabolic interactions and microbial cooperation on degradation rates. We aim to reframe some of the key questions in the study of polysaccharide-bacteria interactions in the context of eco-evolutionary dynamics, highlighting possible research directions that, if pursued, would advance our understanding of polysaccharide degraders at the interface between biochemistry, ecology and evolution.

Large-scale gene gains and losses molded the NLR defense arsenal during the Cucurbita evolution

Main conclusion Genome-wide annotation reveals that the gene birth–death process of the Cucurbita R family is associated with a species-specific diversification of TNL and CNL protein classes. Abstract The Cucurbitaceae family includes nearly 1000 plant species known universally as cucurbits. Cucurbita genus includes many economically important worldwide crops vulnerable to more than 200 pathogens. Therefore, the identification of pathogen-recognition genes is of utmost importance for this genus. The major class of plant-resistance (R) genes encodes nucleotide-binding site and leucine-rich repeat (NLR) proteins, and is divided into three sub-classes namely, TIR-NB-LRR (TNL), CC-NB-LRR (CNL) and RPW8-NB-LRR (RNL). Although the characterization of the NLR gene family has been carried out in important Cucurbita species, this information is still linked to the availability of sequenced genomes. In this study, we analyzed 40 de novo transcriptomes and 5 genome assemblies, which were explored to investigate the Cucurbita expressed-NLR (eNLR) and NLR repertoires using an ad hoc gene annotation approach. Over 1850 NLR-encoding genes were identified, finely characterized and compared to 96 well-characterized plant R-genes. The maximum likelihood analyses revealed an unusual diversification of CNL/TNL genes and a strong RNL conservation. Indeed, several gene gain and loss events have shaped the Cucurbita NLR family. Finally, to provide a first validation step Cucurbita, eNLRs were explored by real-time PCR analysis. The NLR repertories of the 12 Cucurbita species presented in this paper will be useful to discover novel R-genes.

Zanthoxylum-specific whole genome duplication and recent activity of transposable elements in the highly repetitive paleotetraploid Z. bungeanum genome

Zanthoxylum bungeanum is an important spice and medicinal plant that is unique for its accumulation of abundant secondary metabolites, which create a characteristic aroma and tingling sensation in the mouth. Owing to the high proportion of repetitive sequences, high heterozygosity, and increased chromosome number of Z. bungeanum, the assembly of its chromosomal pseudomolecules is extremely challenging. Here, we present a genome sequence for Z. bungeanum, with a dramatically expanded size of 4.23 Gb, assembled into 68 chromosomes. This genome is approximately tenfold larger than that of its close relative Citrus sinensis. After the divergence of Zanthoxylum and Citrus, the lineage-specific whole-genome duplication event η-WGD approximately 26.8 million years ago (MYA) and the recent transposable element (TE) burst ~6.41 MYA account for the substantial genome expansion in Z. bungeanum. The independent Zanthoxylum-specific WGD event was followed by numerous fusion/fission events that shaped the genomic architecture. Integrative genomic and transcriptomic analyses suggested that prominent species-specific gene family expansions and changes in gene expression have shaped the biosynthesis of sanshools, terpenoids, and anthocyanins, which contribute to the special flavor and appearance of Z. bungeanum. In summary, the reference genome provides a valuable model for studying the impact of WGDs with recent TE activity on gene gain and loss and genome reconstruction and provides resources to accelerate Zanthoxylum improvement.

Genome-wide analysis of BpDof genes and the tolerance to drought stress in birch (Betula platyphylla)

Background DNA binding with one finger (Dof) proteins are plant-specific transcription factors playing vital roles in developmental processes and stress responses in plants. Nevertheless, the characterizations, expression patterns, and functions of the Dof family under drought stress (a key determinant of plant physiology and metabolic homeostasis) in woody plants remain unclear. Methods The birch (Betula platyphylla var. mandshuric) genome and plant TFDB database were used to identify Dof gene family members in birch plants. ClustalW2 of BioEdit v7.2.1, MEGA v7.0, ExPASy ProtParam tool, Subloc, TMHMM v2.0, GSDS v2.0, MEME, TBtools, KaKs Calculator v2.0, and PlantCARE were respectively used to align the BpDof sequences, build a phylogenetic tree, identify the physicochemical properties, analyze the chromosomal distribution and synteny, and identify the cis-elements in the promoter regions of the 26 BpDof genes. Additionally, the birch seedlings were exposed to PEG6000-simulated drought stress, and the expression patterns of the BpDof genes in different tissues were analyzed by qRT-PCR. The histochemical staining and the evaluation of physiological indexes were performed to assess the plant tolerance to drought with transient overexpression of BpDof4, BpDof11, and BpDof17 genes. SPSS software and ANOVA were used to conduct all statistical analyses and determine statistically significant differences between results. Results A total of 26 BpDof genes were identified in birch via whole-genome analysis. The conserved Dof domain with a C(x)2C(x)21C(x)2C zinc finger motif was present in all BpDof proteins. These birch BpDofs were classified into four groups (A to D) according to the phylogenetic analysis of Arabidopsis thaliana Dof genes. BpDof proteins within the same group mostly possessed similar motifs, as detected by conserved motif analysis. The exon–intron analysis revealed that the structures of BpDof genes differed, indicating probable gene gain and lose during the BpDof evolution. The chromosomal distribution and synteny analysis showed that the 26 BpDofs were unevenly distributed on 14 chromosomes, and seven duplication events among six chromosomes were found. Cis-acting elements were abundant in the promoter regions of the 26 BpDof genes. qRT-PCR revealed that the expression of the 26 BpDof genes was differentially regulated by drought stress among roots, stems, and leaves. Most BpDof genes responded to drought stress, and BpDof4, BpDof11, and BpDof17­ were significantly up-regulated. Therefore, plants overexpressing these three genes were generated to investigate drought stress tolerance. The BpDof4-, BpDof11-, and BpDof17­-overexpressing plants showed promoted reactive oxygen species (ROS) scavenging capabilities and less severe cell damage, suggesting that they conferred enhanced drought tolerance in birch. This study provided an in-depth insight into the structure, evolution, expression, and function of the Dof gene family in plants.

Feedback regulation of Notch signaling and myogenesis connected by MyoD–Dll1 axis

Muscle precursor cells known as myoblasts are essential for muscle development and regeneration. Notch signaling is an ancient intercellular communication mechanism that plays prominent roles in controlling the myogenic program of myoblasts. Currently whether and how the myogenic cues feedback to refine Notch activities in these cells are largely unknown. Here, by mouse and human gene gain/loss-of-function studies, we report that MyoD directly turns on the expression of Notch-ligand gene Dll1 which activates Notch pathway to prevent precautious differentiation in neighboring myoblasts, while autonomously inhibits Notch to facilitate a myogenic program in Dll1 expressing cells. Mechanistically, we studied cis-regulatory DNA motifs underlying the MyoD–Dll1–Notch axis in vivo by characterizing myogenesis of a novel E-box deficient mouse model, as well as in human cells through CRISPR-mediated interference. These results uncovered the crucial transcriptional mechanism that mediates the reciprocal controls of Notch and myogenesis.