Reversion Mosaicism in Primary Immunodeficiency Diseases

Reversion mosaicism has been reported in an increasing number of genetic disorders including primary immunodeficiency diseases. Several mechanisms can mediate somatic reversion of inherited mutations. Back mutations restore wild-type sequences, whereas second-site mutations result in compensatory changes. In addition, intragenic recombination, chromosomal deletions, and copy-neutral loss of heterozygosity have been demonstrated in mosaic individuals. Revertant cells that have regained wild-type function may be associated with milder disease phenotypes in some immunodeficient patients with reversion mosaicism. Revertant cells can also be responsible for immune dysregulation. Studies identifying a large variety of genetic changes in the same individual further support a frequent occurrence of reversion mosaicism in primary immunodeficiency diseases. This phenomenon also provides unique opportunities to evaluate the biological effects of restored gene expression in different cell lineages. In this paper, we review the recent findings of reversion mosaicism in primary immunodeficiency diseases and discuss its clinical implications.

Multi-Omics Analyses Reveal the Regulatory Network and the Function of ZmUGTs in Maize Defense Response

Maize is one of the major crops in the world; however, diseases caused by various pathogens seriously affect its yield and quality. The maize Rp1-D21 mutant (mt) caused by the intragenic recombination between two nucleotide-binding, leucine-rich repeat (NLR) proteins, exhibits autoactive hypersensitive response (HR). In this study, we integrated transcriptomic and metabolomic analyses to identify differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) in Rp1-D21 mt compared to the wild type (WT). Genes involved in pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) and effector-triggered immunity (ETI) were enriched among the DEGs. The salicylic acid (SA) pathway and the phenylpropanoid biosynthesis pathway were induced at both the transcriptional and metabolic levels. The DAMs identified included lipids, flavones, and phenolic acids, including 2,5-DHBA O-hexoside, the production of which is catalyzed by uridinediphosphate (UDP)-dependent glycosyltransferase (UGT). Four maize UGTs (ZmUGTs) homologous genes were among the DEGs. Functional analysis by transient co-expression in Nicotiana benthamiana showed that ZmUGT9250 and ZmUGT5174, but not ZmUGT9256 and ZmUGT8707, partially suppressed the HR triggered by Rp1-D21 or its N-terminal coiled-coil signaling domain (CCD21). None of the four ZmUGTs interacted physically with CCD21 in yeast two-hybrid or co-immunoprecipitation assays. We discuss the possibility that ZmUGTs might be involved in defense response by regulating SA homeostasis.

Multiple Mechanisms Drive the Evolutionary Adaptation of Phytophthora infestans Effector Avr1 to Host Resistance

Effectors, a group of small proteins secreted by pathogens, play a central role in antagonistic interactions between plant hosts and pathogens. The evolution of effector genes threatens plant disease management and sustainable food production, but population genetic analyses to understand evolutionary mechanisms of effector genes are limited compared to molecular and functional studies. Here we investigated the evolution of the Avr1 effector gene from 111 Phytophthora infestans isolates collected from six areas covering three potato cropping regions in China using a population genetic approach. High genetic variation of the effector gene resulted from diverse mechanisms including base substitution, pre-termination, intragenic recombination and diversifying selection. Nearly 80% of the 111 sequences had a point mutation in the 512th nucleotide (T512G), which generated a pre-termination stop codon truncating 38 amino acids in the C-terminal, suggesting that the C-terminal may not be essential to ecological and biological functions of P. infestans. A significant correlation between the frequency of Avr1 sequences with the pre-termination and annual mean temperature in the collection sites suggests that thermal heterogeneity might be one of contributors to the diversifying selection, although biological and biochemical mechanisms of the likely thermal adaptation are not known currently. Our results highlight the risk of rapid adaptation of P. infestans and possibly other pathogens as well to host resistance, and the application of eco-evolutionary principles is necessary for sustainable disease management in agricultural ecosystems.

Virulence effector SidJ evolution in Legionella pneumophila is driven by positive selection and intragenic recombination

Effector proteins translocated by the Dot/Icm type IV secretion system determine the virulence of Legionella pneumophila (L. pneumophila). Among these effectors, members of the SidE family (SidEs) regulate several cellular processes through a unique phosphoribosyl ubiquitination mechanism mediated by another effector, SidJ. Host-cell calmodulin (CaM) activates SidJ to glutamylate the SidEs of ubiquitin (Ub) ligases and to make a balanced Ub ligase activity. Given the central role of SidJ in this regulatory process, studying the nature of evolution of sidJ is important to understand the virulence of L. pneumophila and the interaction between the bacteria and its hosts. By studying sidJ from a large number of L. pneumophila strains and using various molecular evolution algorithms, we demonstrated that intragenic recombination drove the evolution of sidJ and contributed to sidJ diversification. Additionally, we showed that four codons of sidJ which are located in the N-terminal (NTD) (codons 58 and 200) and C-terminal (CTD) (codons 868 and 869) domains, but not in the kinase domain (KD) had been subjected to strong positive selection pressure, and variable mutation profiles of these codons were identified. Protein structural modeling of SidJ provided possible explanations for these mutations. Codons 868 and 869 mutations might engage in regulating the interactions of SidJ with CaM through hydrogen bonds and affect the CaM docking to SidJ. Mutation in codon 58 of SidJ might affect the distribution of main-chain atoms that are associated with the interaction with CaM. In contrast, mutations in codon 200 might influence the α-helix stability in the NTD. These mutations might be important to balance Ub ligase activity for different L. pneumophila hosts. This study first reported that intragenic recombination and positive Darwinian selection both shaped the genetic plasticity of sidJ, contributing to a deeper understanding of the adaptive mechanisms of this intracellular bacterium to different hosts.

A Case of Intragenic Recombination Dramatically Impacting the Phage WO Genetic Diversity in Gall Wasps

The phage WO was characterized in Wolbachia, a strictly intracellular bacterium causing several reproductive alterations in its arthropod hosts. This study aimed to screen the presence of Wolbachia and phage WO in 15 gall wasp species from six provinces of southern China to investigate their diversity and prevalence patterns. A high incidence of Wolbachia infection was determined in the gall wasp species, with an infection rate of 86.7% (13/15). Moreover, seven species had double or multiple infections. All Wolbachia-infected gall wasp species were found to harbor phage WO. The gall wasp species infected with a single Wolbachia strain were found to harbor a single phage WO type. On the contrary, almost all species with double or multiple Wolbachia infections harbored a high level of phage WO diversity (ranging from three to 27 types). Six horizontal transfer events of phage WO in Wolbachia were found to be associated with gall wasps, which shared identical orf7 sequences among their respective accomplices. The transfer potentially took place through gall inducers and associated inquilines infected with or without Wolbachia. Furthermore, 10 putative recombination events were identified from Andricus hakonensis and Andricus sp2, which harbored multiple phage WO types, suggesting that intragenic recombination was the important evolutionary force, which effectively promoted the high level of phage WO diversity associated with gall wasps.

Horizontal gene transfer and recombination analysis of SARS-CoV-2 genes helps discover its close relatives and shed light on its origin

Background The SARS-CoV-2 pandemic is one of the greatest global medical and social challenges that have emerged in recent history. Human coronavirus strains discovered during previous SARS outbreaks have been hypothesized to pass from bats to humans using intermediate hosts, e.g. civets for SARS-CoV and camels for MERS-CoV. The discovery of an intermediate host of SARS-CoV-2 and the identification of specific mechanism of its emergence in humans are topics of primary evolutionary importance. In this study we investigate the evolutionary patterns of 11 main genes of SARS-CoV-2. Previous studies suggested that the genome of SARS-CoV-2 is highly similar to the horseshoe bat coronavirus RaTG13 for most of the genes and to some Malayan pangolin coronavirus (CoV) strains for the receptor binding (RB) domain of the spike protein. Results We provide a detailed list of statistically significant horizontal gene transfer and recombination events (both intergenic and intragenic) inferred for each of 11 main genes of the SARS-CoV-2 genome. Our analysis reveals that two continuous regions of genes S and N of SARS-CoV-2 may result from intragenic recombination between RaTG13 and Guangdong (GD) Pangolin CoVs. Statistically significant gene transfer-recombination events between RaTG13 and GD Pangolin CoV have been identified in region [1215–1425] of gene S and region [534–727] of gene N. Moreover, some statistically significant recombination events between the ancestors of SARS-CoV-2, RaTG13, GD Pangolin CoV and bat CoV ZC45-ZXC21 coronaviruses have been identified in genes ORF1ab, S, ORF3a, ORF7a, ORF8 and N. Furthermore, topology-based clustering of gene trees inferred for 25 CoV organisms revealed a three-way evolution of coronavirus genes, with gene phylogenies of ORF1ab, S and N forming the first cluster, gene phylogenies of ORF3a, E, M, ORF6, ORF7a, ORF7b and ORF8 forming the second cluster, and phylogeny of gene ORF10 forming the third cluster. Conclusions The results of our horizontal gene transfer and recombination analysis suggest that SARS-CoV-2 could not only be a chimera virus resulting from recombination of the bat RaTG13 and Guangdong pangolin coronaviruses but also a close relative of the bat CoV ZC45 and ZXC21 strains. They also indicate that a GD pangolin may be an intermediate host of this dangerous virus.

Horizontal gene transfer and recombination analysis of SARS-CoV-2 genes helps discover its close relatives and shed light on its origin

BackgroundThe SARS-CoV-2 pandemic is among the most dangerous infectious diseases that have emerged in recent history. Human CoV strains discovered during previous SARS outbreaks have been hypothesized to pass from bats to humans using intermediate hosts, e.g. civets for SARS-CoV and camels for MERS-CoV. The discovery of an intermediate host of SARS-CoV-2 and the identification of specific mechanism of its emergence in humans are topics of primary evolutionary importance. In this study we investigate the evolutionary patterns of 11 main genes of SARS-CoV-2. Previous studies suggested that the genome of SARS-CoV-2 is highly similar to the horseshoe bat coronavirus RaTG13 for most of the genes and to some Malayan pangolin coronavirus (CoV) strains for the receptor binding (RB) domain of the spike protein.ResultsWe provide a detailed list of statistically significant horizontal gene transfer and recombination events (both intergenic and intragenic) inferred for each of 11 main genes of the SARS-Cov-2 genome. Our analysis reveals that two continuous regions of genes S and N of SARS-CoV-2 may result from intragenic recombination between RaTG13 and Guangdong (GD) Pangolin CoVs. Statistically significant gene transfer-recombination events between RaTG13 and GD Pangolin CoV have been identified in region [1215-1425] of gene S and region [534-727] of gene N. Moreover, some significant recombination events between the ancestors of SARS-CoV-2, RaTG13, GD Pangolin CoV and bat CoV ZC45-ZXC21 coronaviruses have been identified in genes ORF1ab, S, ORF3a, ORF7a, ORF8 and N. Furthermore, topology-based clustering of gene trees inferred for 25 CoV organisms revealed a three-way evolution of coronavirus genes, with gene phylogenies of ORF1ab, S and N forming the first cluster, gene phylogenies of ORF3a, E, M, ORF6, ORF7a, ORF7b and ORF8 forming the second cluster, and phylogeny of gene ORF10 forming the third cluster.ConclusionsThe results of our horizontal gene transfer and recombination analysis suggest that SARS-Cov-2 could not only be a chimera resulting from recombination of the bat RaTG13 and Guangdong pangolin coronaviruses but also a close relative of the bat CoV ZC45 and ZXC21 strains. They also indicate that a GD pangolin may be an intermediate host of SARS-CoV-2.

Molecular Evolution of SARS-CoV-2 Structural Genes: Evidence of Positive Selection in the Spike Glycoprotein

Background: SARS-CoV-2 has caused a global pandemic since early 2020 and remains a serious public health issue worldwide. Four structural genes, envelope (E), membrane (M), nucleocapsid (N) and spike (S), play a key role in controlling entry into human cells and virion assembly of SARS-CoV-2. The evolution of these genes may determine infectivity of SARS-CoV-2, but thus far, little is known about them. Methods: We analyzed 3090 SARS-CoV-2 isolates from the GenBank database to determine the evolutionary patterns of the four structural genes by employing various molecular evolution algorithms. Results: Phylogenetic analyses showed that global SARS-CoV-2 isolates can be clustered into three to four major clades based upon protein sequence. Although intragenic recombination was not detected among different alleles, purifying selection has affected the evolution of these genes. By analyzing full genomic sequences of these alleles, our result revealed that codon 614 of the S glycoprotein has been subjected to a strong positive selection pressure, and a consistent D614G mutation was identified. Additionally, another potentially positive selection site at codon 5 in the signal sequence of the S protein was also identified with a consistent L5F mutation. The allele containing the D614G mutation has undergone significant expansion during SARS-CoV-2 transmission, implying a better adaptability of isolates with the mutation. Nevertheless, L5F allele expansion was found to be relatively restricted. The D614G mutation is located at subdomain 2 (SD2) of the C-terminal portion (CTP) of the S1 subunit. Protein structural modeling showed that the D614G mutation may cause the disruption of a salt bridge between S protein monomers and increase their flexibility, consequently promoting receptor binding domain (RBD) opening, virus attachment, and ultimately entry into host cells. Located at the signal sequence of S protein, the L5F mutation may facilitate protein folding, assembly, and secretion of the virus. Conclusions: This is the first reported evidence of positive Darwinian selection in the spike gene of SARS-CoV-2. This finding contributes to a broader understanding of the adaptive mechanisms of this virus, and provide insight for the development of novel therapeutic approaches, as well as the creation of effective vaccines, through targeting mutation sites.