To catch a hijacker: abundance, evolution and genetic diversity of P4-like bacteriophage satellites

Bacteriophages (phages) are bacterial parasites that can themselves be parasitized by phage satellites. The molecular mechanisms used by satellites to hijack phages are sometimes understood in great detail, but the origins, abundance, distribution and composition of these elements are poorly known. Here, we show that P4-like elements are present in more than 30% of the genomes of Enterobacterales, and in almost half of those of Escherichia coli , sometimes in multiple distinct copies. We identified over 1000 P4-like elements with very conserved genetic organization of the core genome and a few hotspots with highly variable genes. These elements are never found in plasmids and have very little homology to known phages, suggesting an independent evolutionary origin. Instead, they are scattered across chromosomes, possibly because their integrases are often exchanged with other elements. The rooted phylogenies of hijacking functions are correlated and suggest longstanding coevolution. They also reveal broad host ranges in P4-like elements, as almost identical elements can be found in distinct bacterial genera. Our results show that P4-like phage satellites constitute a very distinct, widespread and ancient family of mobile genetic elements. They pave the way for studying the molecular evolution of antagonistic interactions between phages and their satellites. This article is part of the theme issue ‘The secret lives of microbial mobile genetic elements’.

Computational Inference, Validation, and Analysis of 5’UTR-Leader Sequences of Alleles of Immunoglobulin Heavy Chain Variable Genes

Upstream and downstream sequences of immunoglobulin genes may affect the expression of such genes. However, these sequences are rarely studied or characterized in most studies of immunoglobulin repertoires. Inference from large, rearranged immunoglobulin transcriptome data sets offers an opportunity to define the upstream regions (5’-untranslated regions and leader sequences). We have now established a new data pre-processing procedure to eliminate artifacts caused by a 5’-RACE library generation process, reanalyzed a previously studied data set defining human immunoglobulin heavy chain genes, and identified novel upstream regions, as well as previously identified upstream regions that may have been identified in error. Upstream sequences were also identified for a set of previously uncharacterized germline gene alleles. Several novel upstream region variants were validated, for instance by their segregation to a single haplotype in heterozygotic subjects. SNPs representing several sequence variants were identified from population data. Finally, based on the outcomes of the analysis, we define a set of testable hypotheses with respect to the placement of particular alleles in complex IGHV locus haplotypes, and discuss the evolutionary relatedness of particular heavy chain variable genes based on sequences of their upstream regions.

Single-cell RNA sequencing analysis of mouse follicular somatic cells

Within the development of ovarian follicle, in addition to cell proliferation and differentiation, sophisticated cell–cell cross talks are established among follicular somatic cells such as granulosa cells (GCs) and theca cells. To systematically reveal the cell differentiation and signal transductions in follicular somatic cells, we collected the mouse follicular somatic cells from secondary to ovulatory stage, and analyzed the single cell transcriptomes. Having data filtered and screened, we found 6883 high variable genes in 4888 single cells. Then follicular somatic cells were clustered into 26 cell clusters, including 18 GC clusters, 4 theca endocrine cell (TEC) clusters, and 4 other somatic cell clusters, which include immune cells and Acta2 positive theca externa cells. From our data, we found there was metabolic reprogramming happened during GC differentiation. We also found both Cyp19a1 and Cyp11a1 could be expressed in TECs. We analyzed the expression patterns of genes associated with cell–cell interactions such as steroid hormone receptor genes, insulin signaling genes, and cytokine/transformation growth factor beta associated genes in all cell clusters. Lastly, we clustered the highly variable genes into 300 gene clusters, which could be used to search new genes involved in follicle development. These transcriptomes of follicular somatic cells provide us potential clues to reveal how mammals regulating follicle development and could help us find targets to improve oocyte quality for women with low fertility.

Stabilizing Genetically Unstable Simple Sequence Repeats in the Campylobacter jejuni Genome by Multiplex Genome Editing: a Reliable Approach for Delineating Multiple Phase-Variable Genes

Campylobacter jejuni is the leading bacterial cause of foodborne gastroenteritis in developed countries and occasionally progresses to the autoimmune disease Guillain-Barré syndrome. A relatively large number of hypermutable simple sequence repeat (SSR) tracts in the C. jejuni genome markedly decreases its phenotypic stability through reversible changes in the ON or OFF expression states of the genes in which they reside, a phenomenon called phase variation.

Phase Variation During Host Colonization and Invasion by Campylobacter jejuni and Other Campylobacter Species

Phase variation (PV) is a phenomenon common to a variety of bacterial species for niche adaption and survival in challenging environments. Among Campylobacter species, PV depends on the presence of intergenic and intragenic hypermutable G/C homopolymeric tracts. The presence of phase-variable genes is of especial interest for species that cause foodborne or zoonotic infections in humans. PV influences the formation and the structure of the lipooligosaccharide, flagella, and capsule in Campylobacter species. PV of components of these molecules is potentially important during invasion of host tissues, spread within hosts and transmission between hosts. Motility is a critical phenotype that is potentially modulated by PV. Variation in the status of the phase-variable genes has been observed to occur during colonization in chickens and mouse infection models. Interestingly, PV is also involved in bacterial survival of attack by bacteriophages even during chicken colonization. This review aims to explore and discuss observations of PV during model and natural infections by Campylobacter species and how PV may affect strategies for fighting infections by this foodborne pathogen.

Genome Scan for Variable Genes Involved in Environmental Adaptations of Nubian Ibex

The Nubian ibex (Capra nubiana) is a wild goat species that inhabits the Sahara and Arabian deserts and is adapted to extreme ambient temperatures, intense solar radiation, and scarcity of food and water resources. To investigate desert adaptation, we explored the possible role of copy number variations (CNVs) in the evolution of Capra species with a specific focus on the environment of Capra nubiana. CNVs are structural genomic variations that have been implicated in phenotypic differences between species and could play a role in species adaptation. CNVs were inferred from Capra nubiana sequence data relative to the domestic goat reference genome using read-depth approach. We identified 191 CNVs overlapping with protein-coding genes mainly involved in biological processes such as innate immune response, xenobiotic metabolisms, and energy metabolisms. We found copy number variable genes involved in defense response to viral infections (Cluster of Differentiation 48, UL16 binding protein 3, Natural Killer Group 2D ligand 1-like, and Interferon-induced transmembrane protein 3), possibly suggesting their roles in Nubian ibex adaptations to viral infections. Additionally, we found copy number variable xenobiotic metabolism genes (carboxylesterase 1, Cytochrome P450 2D6, Glutathione S-transferase Mu 4, and UDP Glucuronosyltransferase-2B7), which are probably an adaptation of Nubian ibex to desert diets that are rich in plant secondary metabolites. Collectively, this study's results advance our understanding of CNVs and their possible roles in the adaptation of Nubian ibex to its environment. The copy number variable genes identified in Nubian ibex could be considered as subjects for further functional characterizations.

Computational inference, validation, and analysis of 5'UTR-leader sequences of alleles of immunoglobulin heavy chain variable genes

Upstream and downstream sequences of immunoglobulin genes may affect the expression of such genes. However, these sequences are rarely studied or characterized in most studies of immunoglobulin repertoires. Inference from large, rearranged immunoglobulin transcriptome data sets offers an opportunity to define the upstream regions (5'-untranslated regions and leader sequences). We have now established a new data pre-processing procedure to eliminate artifacts caused by a 5'-RACE library generation process, reanalyzed a previously studied data set defining human immunoglobulin heavy chain genes, and identified novel upstream regions, as well as previously identified upstream regions that may have been identified in error. Upstream sequences were also identified for a set of previously uncharacterized germline gene alleles. Several novel upstream region variants were validated, for instance by their segregation to a single haplotype in heterozygotic subjects. SNPs representing several sequence variants were identified from population data. Finally, based on the outcomes of the analysis, we define a set of testable hypotheses with respect to the placement of particular alleles in complex IGHV locus haplotypes, and discuss the evolutionary relatedness of particular heavy chain variable genes based on sequences of their upstream regions.

Sorghum Pan-Genome Explores the Functional Utility for Genomic-Assisted Breeding to Accelerate the Genetic Gain

Sorghum (Sorghum bicolor L.) is a staple food crops in the arid and rainfed production ecologies. Sorghum plays a critical role in resilient farming and is projected as a smart crop to overcome the food and nutritional insecurity in the developing world. The development and characterisation of the sorghum pan-genome will provide insight into genome diversity and functionality, supporting sorghum improvement. We built a sorghum pan-genome using reference genomes as well as 354 genetically diverse sorghum accessions belonging to different races. We explored the structural and functional characteristics of the pan-genome and explain its utility in supporting genetic gain. The newly-developed pan-genome has a total of 35,719 genes, a core genome of 16,821 genes and an average of 32,795 genes in each cultivar. The variable genes are enriched with environment responsive genes and classify the sorghum accessions according to their race. We show that 53% of genes display presence-absence variation, and some of these variable genes are predicted to be functionally associated with drought adaptation traits. Using more than two million SNPs from the pan-genome, association analysis identified 398 SNPs significantly associated with important agronomic traits, of which, 92 were in genes. Drought gene expression analysis identified 1,788 genes that are functionally linked to different conditions, of which 79 were absent from the reference genome assembly. This study provides comprehensive genomic diversity resources in sorghum which can be used in genome assisted crop improvement.

Stabilizing genetically unstable simple sequence repeats in the Campylobacter jejuni genome by multiplex genome editing: a reliable approach for delineating multiple phase-variable genes

Hypermutability of simple sequence repeats (SSR) through DNA slippage is a major mechanism of phase variation in Campylobacter jejuni . The presence of multiple SSR-mediated phase-variable genes encoding enzymes that modify surface structures, including capsular polysaccharide (CPS) and lipooligosaccharide (LOS), generates high levels of structural variants within bacterial populations, thereby promoting adaptation to selective pressures in host environments. Therefore, the phenotypic diversity generated by phase variation can limit the reproducibility of results with C. jejuni ; therefore, researchers need to genetically control the mutability of multiple SSRs. Here, we show that natural “cotransformation” is an effective method for C. jejuni genome editing. Cotransformation is a trait of naturally competent bacteria that causes uptake and integration of multiple different DNA fragments, which has been recently adapted to multiplex genome editing by natural transformation (MuGENT), a method for introducing multiple scarless mutations into the genomes of these bacteria. We found that the cotransformation frequencies of antibiotic resistance gene-marked DNA fragments and unmarked DNA fragments reached ~40% in C. jejuni . To examine the feasibility of MuGENT in C. jejuni , we “locked” either different polyG SSR tracts in strain NCTC11168 (which are located in the biosynthetic CPS and LOS gene clusters) into either the ON or OFF configurations by interrupting the continuous runs of G residues without changing the encoded amino acids. This approach, termed “MuGENT-SSR,” enabled the generation of all eight edits within 2 weeks and the identification of a phase-locked strain with a highly stable type of Penner serotyping, a CPS-based serotyping scheme. Furthermore, extensive genome editing of this strain by MuGENT-SSR identified a phase-variable gene that determines the Penner serotype of NCTC11168. Thus, MuGENT-SSR provides a platform for genetic and phenotypic engineering of genetically unstable C. jejuni , making it a reliable approach for elucidating the mechanisms underlying phase-variable expression of specific phenotypes.