A rare CTSC mutation in Papillon-Lefèvre Syndrome results in abolished serine protease activity and reduced NET formation but otherwise normal neutrophil function

Papillon-Lefèvre Syndrome (PLS) is an autosomal recessive monogenic disease caused by loss-of-function mutations in the CTSC gene, thus preventing the synthesis of the protease Cathepsin C (CTSC) in a proteolytically active form. CTSC is responsible for the activation of the pro-forms of the neutrophil serine proteases (NSPs; Elastase, Proteinase 3 and Cathepsin G), suggesting its involvement in a variety of neutrophil functions. In PLS neutrophils, the lack of CTSC protease activity leads to inactivity of the NSPs. Clinically, PLS is characterized by an early, typically pre-pubertal, onset of severe periodontal pathology and palmoplantar hyperkeratosis. However, PLS is not considered an immune deficiency as patients do not typically suffer from recurrent and severe (bacterial and fungal) infections. In this study we investigated an unusual CTSC mutation in two siblings with PLS, a 503A>G substitution in exon 4 of the CTSC gene, expected to result in an amino acid replacement from tyrosine to cysteine at position 168 of the CTSC protein. Both patients bearing this mutation presented with pronounced periodontal pathology. The characteristics and functions of neutrophils from patients homozygous for the 503A>G CTSC mutation were compared to another previously described PLS mutation (755A>T), and a small cohort of healthy volunteers. Neutrophil lysates from patients with the 503A>G substitution lacked CTSC protein and did not display any CTSC or NSP activity, yet neutrophil counts, morphology, priming, chemotaxis, radical production, and regulation of apoptosis were without any overt signs of alteration. However, NET formation upon PMA-stimulation was found to be severely depressed, but not abolished, in PLS neutrophils.

A Small Step, a Giant Leap: Somatic Hypermutation of a Single Amino Acid Leads to Anti-La Autoreactivity

The anti-La mab 312B, which was established by hybridoma technology from human-La transgenic mice after adoptive transfer of anti-human La T cells, immunoprecipitates both native eukaryotic human and murine La protein. Therefore, it represents a true anti-La autoantibody. During maturation, the anti-La mab 312B acquired somatic hypermutations (SHMs) which resulted in the replacement of four aa in the complementarity determining regions (CDR) and seven aa in the framework regions. The recombinant derivative of the anti-La mab 312B in which all the SHMs were corrected to the germline sequence failed to recognize the La antigen. We therefore wanted to learn which SHM(s) is (are) responsible for anti-La autoreactivity. Humanization of the 312B ab by grafting its CDR regions to a human Ig backbone confirms that the CDR sequences are mainly responsible for anti-La autoreactivity. Finally, we identified that a single amino acid replacement (D > Y) in the germline sequence of the CDR3 region of the heavy chain of the anti-La mab 312B is sufficient for anti-La autoreactivity.

In Silico and In Vivo Analysis of Amino Acid Substitutions That Cause Laminopathies

Mutations in the LMNA gene cause diseases called laminopathies. LMNA encodes lamins A and C, intermediate filaments with multiple roles at the nuclear envelope. LMNA mutations are frequently single base changes that cause diverse disease phenotypes affecting muscles, nerves, and fat. Disease-associated amino acid substitutions were mapped in silico onto three-dimensional structures of lamin A/C, revealing no apparent genotype–phenotype connections. In silico analyses revealed that seven of nine predicted partner protein binding pockets in the Ig-like fold domain correspond to sites of disease-associated amino acid substitutions. Different amino acid substitutions at the same position within lamin A/C cause distinct diseases, raising the question of whether the nature of the amino acid replacement or genetic background differences contribute to disease phenotypes. Substitutions at R249 in the rod domain cause muscular dystrophies with varying severity. To address this variability, we modeled R249Q and R249W in Drosophila Lamin C, an orthologue of LMNA. Larval body wall muscles expressing mutant Lamin C caused abnormal nuclear morphology and premature death. When expressed in indirect flight muscles, R249W caused a greater number of adults with wing posturing defects than R249Q, consistent with observations that R249W and R249Q cause distinct muscular dystrophies, with R249W more severe. In this case, the nature of the amino acid replacement appears to dictate muscle disease severity. Together, our findings illustrate the utility of Drosophila for predicting muscle disease severity and pathogenicity of variants of unknown significance.

Measuring Phylogenetic Information of Incomplete Sequence Data

Widely used approaches for extracting phylogenetic information from aligned sets of molecular sequences rely upon probabilistic models of nucleotide substitution or amino-acid replacement. The phylogenetic information that can be extracted depends on the number of columns in the sequence alignment and will be decreased when the alignment contains gaps due to insertion or deletion events. Motivated by the measurement of information loss, we suggest assessment of the Effective Sequence Length (ESL) of an aligned data set. The ESL can differ from the actual number of columns in a sequence alignment because of the presence of alignment gaps. Furthermore, the estimation of phylogenetic information is affected by model misspecification. Inevitably, the actual process of molecular evolution differs from the probabilistic models employed to describe this process. This disparity means the amount of phylogenetic information in an actual sequence alignment will differ from the amount in a simulated data set of equal size, which motivated us to develop a new test for model adequacy. Via theory and empirical data analysis, we show how to disentangle the effects of gaps and model misspecification. By comparing the Fisher information of actual and simulated sequences, we identify which alignment sites and tree branches are most affected by gaps and model misspecification.

Evolutionary Profile for (Host and Viral) MLKL Indicates Its Activities as a Battlefront for Extensive Counteradaptation

Pathogen infection triggers host innate defenses which may result in the activation of regulated cell death (RCD) pathways such as apoptosis. Given a vital role in immunity, apoptotic effectors are often counteracted by pathogen-encoded antagonists. Mounting evidence indicates that programmed necrosis, which is mediated by the RIPK3/MLKL axis and termed necroptosis, evolved as a countermeasure to pathogen-mediated inhibition of apoptosis. Yet, it is unclear whether components of this emerging RCD pathway display signatures associated with pathogen conflict that are rare in combination but common to key host defense factors, namely, rapid evolution, viral homolog (virolog), and cytokine induction. We leveraged evolutionary sequence analysis that examines rates of amino acid replacement, which revealed: 1) strong and recurrent signatures of positive selection for primate and bat RIPK3 and MLKL, and 2) elevated rates of amino acid substitution on multiple RIPK3/MLKL surfaces suggestive of past antagonism with multiple, distinct pathogen-encoded inhibitors. Furthermore, our phylogenomics analysis across poxvirus genomes illuminated volatile patterns of evolution for a recently described MLKL viral homolog. Specifically, poxviral MLKLs have undergone numerous gene replacements mediated by duplication and deletion events. In addition, MLKL protein expression is stimulated by interferons in human and mouse cells. Thus, MLKL displays all three hallmarks of pivotal immune factors of which only a handful of factors like OAS1 exhibit. These data support the hypothesis that over evolutionary time MLKL functions—which may include execution of necroptosis—have served as a major determinant of infection outcomes despite gene loss in some host genomes.

The Algerian Chapter of SARS-CoV-2 Pandemic: An Evolutionary, Genetic, and Epidemiological Prospect

To explore the SARS-CoV-2 pandemic in Algeria, a dataset comprising ninety-five genomes originating from SARS-CoV-2 sampled from Algeria and other countries worldwide, from 24 December 2019, through 4 March 2021, was thoroughly examined. While performing a multi-component analysis regarding the Algerian outbreak, the toolkit of phylogenetic, phylogeographic, haplotype, and genomic analysis were effectively implemented. We estimated the Time to the Most Recent Common Ancestor (TMRCA) in reference to the Algerian pandemic and highlighted the multiple introductions of the disease and the missing data depicted in the transmission loop. In addition, we emphasized the significant role played by local and international travels in disease dissemination. Most importantly, we unveiled mutational patterns, the effect of unique mutations on corresponding proteins, and the relatedness regarding the Algerian sequences to other sequences worldwide. Our results revealed individual amino-acid replacements such as the deleterious replacement A23T in the orf3a gene in Algeria_EPI_ISL_418241. Additionally, a connection between Algeria_EPI_ISL_420037 and sequences originating from the USA was observed through a USA characteristic amino-acid replacement T1004I in the nsp3 gene, found in the aforementioned Algerian sequence. Similarly, successful tracing could be established, such as Algeria/G37318-8849/2020|EPI_ISL_766863, which was imported from Saudi Arabia during the pilgrimage. Lastly, we assessed the Algerian mitigation measures regarding disease containment using statistical analyses.

An extension of Wang protein design model using Blosum62 substitution matrix

Humans life depends on the functionality of molecules in the body. One of these essential molecules is the protein that plays a vital role in our life, such that its malfunction can cause severe damages. Such roles make protein structure and its functionality necessary to understand. One of the problems that help us understand the relation between protein structure is the well-known protein design problem which attempts to find an amino acid sequence that can fold into a desired tertiary structure. However, despite having an acceptable accuracy in protein design, this accuracy is an identical percentage of amino acid retrieving. At the same time, it is well-known that amino acids can replace each other in evolution while the function and structure of protein stay the same. Thus the designed sequence does not have the opportunity to be close to the target in the evolutionary aspect. This paper presents an extension to Wang's deep learning model, which uses evolutionary information in the Blosum62 substitution matrix to take amino acid replacement probability into account while designing a sequence.

Identification and trajectory of growth of concerning SARS-CoV-2 variants in Houston, Texas, January through April 2021 based on 11,568 genome sequences

Genetic variants of the SARS-CoV-2 virus are of substantial concern because they can detrimentally alter the trajectory of the ongoing pandemic, and disease course in individual patients. Here we report genome sequences from 11,568 COVID-19 patients in the Houston Methodist healthcare system dispersed throughout the metroplex that were diagnosed from January 1, 2021 through April 30, 2021. This sample represents 94% of Houston Methodist cases and 4.6% of all reported cases in the metropolitan area during this period. The SARS-CoV-2 variant designated UK B.1.1.7 increased very rapidly, and now causes 75%-90% of all new cases in the Houston area. Five of the 2,543 B.1.1.7 genomes had an E484K change in spike protein. Compared with non-B.1.1.7 patients, individuals infected with B.1.1.7 had a significantly lower cycle threshold value (considered to be a proxy for higher virus load) and higher rate of hospitalization. Other variants (e.g., B.1.429, B.1.427, P.1, P.2, and R.1) also increased rapidly in frequency, although the magnitude was less than for B.1.1.7. We also identified 42 patients with a recently described R.1 variant that has an E484K amino acid replacement, and seven patients with the B.1.617 "India" variants. In the aggregate, our study shows the occurrence of a diverse array of concerning SARS-CoV-2 variants circulating in a major metropolitan area, documents B.1.1.7 as the major cause of new cases in Houston and heralds the arrival and spread of B.1.617 variants in the metroplex.

Modification of the transglucosylation properties of α-glucosidases from Aspergillus oryzae and Aspergillus sojae via a single critical amino acid replacement

We constructed enzyme variants of the α-glucosidases from Aspergillus oryzae (AoryAgdS) and Aspergillus sojae (AsojAgdL) by mutating the amino acid residue at position 450. AoryAgdS_H450R acquired the ability to produce considerable amounts of α-1,6-transglucosylation products, whereas AsojAgdL_R450H changed to produce more α-1,3- and α-1,4-transglucosylation products than α-1,6-products. The 450th amino acid residue is critical for the transglucosylation of these α-glucosidases.

Genomic signatures of a major adaptive event in the pathogenic fungus Melampsora larici-populina

Background: The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in non-model organisms. Host-pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Results: Using an array of genome scan methods, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype-phenotype statistical association and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although the selective sweep occurred only four years earlier, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus. Conclusions: Our results suggest that M. larici-populina underwent a soft selective sweep and possibly a prominent effect of outbreeding and recombination, which we speculate have increased the efficiency of selection.