Amino Acid Sequence Divergence Research Articles

Whole-genome sequences from wild-type and laboratory-evolved strains define the alleleome and establish its hallmarks

The genomic diversity across strains of a species forms the genetic basis for differences in their behavior. A large-scale assessment of sequence variation has been made possible by the growing availability of strain-specific whole-genome sequences (WGS) and with the advent of large-scale databases of laboratory-acquired mutations. We define the Escherichia coli "alleleome" through a genome-scale assessment of amino acid (AA) sequence diversity in open reading frames across 2,661 WGS from wild-type strains. We observe a highly conserved alleleome enriched in mutations unlikely to affect protein function. In contrast, 33,000 mutations acquired in laboratory evolution experiments result in more severe AA substitutions that are rarely achieved by natural selection. Large-scale assessment of the alleleome establishes a method for the quantification of bacterial allelic diversity, reveals opportunities for synthetic biology to explore novel sequence space, and offers insights into the constraints governing evolution.

Characterization of TRAF genes and their responses to Vibrio anguillarum challenge in Argopecten scallops

The tumor necrosis factor receptor-related factor (TRAF) family has been reported to be involved in many immune pathways, such as TNFR, TLR, NLR, and RLR in animals. However, little is known about the roles of TRAF genes in the innate immune of Argopecten scallops. In this study, we first identified five TRAF genes, including TRAF2, TRAF3, TRAF4, TRAF6 and TRAF7, but not TRAF1 and TRAF5, from both the bay scallop A. irradians (Air) and the Peruvian scallop A. purpuratus (Apu). The phylogenetic analysis showed that the TRAF genes in Argopecten scallops (AiTRAF) belong to the branch of molluscan TRAF family, which lacks TRAF1 and TRAF5. Since TRAF6 is a key bridge factor in the tumor necrosis factor superfamily and plays an important role in innate and adaptive immunity, we cloned the ORFs of the TRAF6 gene in both A. irradians and A. purpuratus, as well as in two reciprocal hybrids (Aip for the hybrid Air × Apu and Api for the hybrid Apu × Air). Differences in conformational and post-translational modification resulted from the variation in amino acid sequences may cause differences in activity among them. Analysis of conserved motifs and protein structural domains revealed that AiTRAF contains typical structural domains similar to those of other mollusks and has the same conserved motifs. Tissue expression of TRAF in Argopecten scallops challenged by Vibrio anguillarum was examined by qRT-PCR. The results showed that AiTRAF were higher in gill and hepatopancreas. When challenged by Vibrio anguillarum, the expression of AiTRAF was significantly increased compared with the control group, indicating that AiTRAF may play an important role in the immunity of scallops. In addition, the expression of TRAF was higher in Api and Aip than in Air when challenged by Vibrio anguillarum, suggesting that TRAF may have contributed to the high resistance of Api and Aip to Vibrio anguillarum. The results of this study may provide new insights into the evolution and function of TRAF genes in bivalves and ultimately benefit scallop breeding.

Proteolytic activity and heat resistance of the protease AprX from Pseudomonas in relation to genotypic characteristics

AprX is an alkaline metalloprotease produced by Pseudomonas spp. and encoded by its initial gene of the aprX-lipA operon. The intrinsic diversity among Pseudomonas spp. regarding their proteolytic activity is the main challenge for the development of accurate methods for spoilage prediction of ultra-high temperature (UHT) treated milk in the dairy industry. In the present study, 56 Pseudomonas strains were characterized by assessing their proteolytic activity in milk before and after lab-scale UHT treatment. From these, 24 strains were selected based on their proteolytic activity for whole genome sequencing (WGS) to identify common genotypic characteristics that correlated with the observed variations in proteolytic activity. Four groups (A1, A2, B and N) were determined based on operon aprX-lipA sequence similarities. These alignment groups were observed to significantly influence the proteolytic activity of the strains, with an average proteolytic activity of A1 > A2 > B > N. The lab-scale UHT treatment did not significantly influence their proteolytic activity, indicating a high thermal stability of proteases among strains. Amino acid sequence variation of biologically-relevant motifs in the AprX sequence, namely the Zn2+-binding motif at the catalytic domain and the C-terminal type I secretion signaling mechanism, were found to be highly conserved within alignment groups. These motifs could serve as future potential genetic biomarkers for determination of alignment groups and thereby strain spoilage potential.

Variability of genes encoding nonstructural proteins of rotavirus А (Reoviridae: Rotavirus: Rotavirus A) genotype G9P[8] during the period of dominance in the territory of Nizhny Novgorod (central part of Russia) (2011-2020)

In Russia, rotavirus A is the main cause of severe viral gastroenteritis in young children. The molecular features that allow a rotavirus of a particular genotype to gain an evolutionary advantage remain unclear, therefore, the study of the genetic diversity of rotaviruses based on genes encoding nonstructural proteins (NSPs) responsible for the reproduction of the virus in the cell is an urgent task. To study the genetic diversity of rotaviruses of genotype G9P[8], which dominated Nizhny Novgorod in 20112020, based on genes encoding nonstructural proteins. Rotavirus-positive samples were subjected to PCR-genotyping and sequencing of NSP1 NSP5 genes. Phylogenetic analysis was carried out in the MEGA X program. In the period 20112020, G9P[8] rotaviruses with four variants of the NSP2 gene were co-circulating in Nizhny Novgorod. New alleles were noted in 2012 (N1-a-III), 2016 (N1-a-IV) and in 2019 (N1-a-II). The appearance of new variants of other genes occurred in 2014 (E1-3, NSP4), 2018 (T1-a3-III, NSP3) and in 2019 (A1-b-II, NSP1). NSP2 gene had the most variable amino acid sequence (16 substitutions), 2 to 7 substitutions were observed in NSP1, NSP3 and NSP4, NSP5 was conservative. The results obtained are consistent with the literature data and indicate the participation of NSP genes in maintaining the heterogeneity of the rotavirus population. Until 2018, the genetic diversity of rotaviruses in Nizhny Novgorod was determined by the circulation of strains carrying several alleles of the NSP2 gene and conservative genes NSP1, NSP3NSP5. By the end of the study period, new variants of the genotype G9P[8] were formed in the population, carrying previously unknown combinations of alleles of nonstructural genes.

Evidence of a Set of Core-Function Genes in 16 Bacillus Podoviral Genomes with Considerable Genomic Diversity.

Bacteriophage genomes represent an enormous level of genetic diversity and provide considerable potential to acquire new insights about viral genome evolution. In this study, the genome sequences of sixteen Bacillus-infecting bacteriophages were explored through comparative genomics approaches to reveal shared and unique characteristics. These bacteriophages are in the Salasmaviridae family with small (18,548-27,206 bp) double-stranded DNA genomes encoding 25-46 predicted open reading frames. We observe extensive nucleotide and amino acid sequence divergence among a set of core-function genes that present clear synteny. We identify two examples of sequence directed recombination within essential genes, as well as explore the expansion of gene content in these genomes through the introduction of novel open reading frames. Together, these findings highlight the complex evolutionary relationships of phage genomes that include old, common origins as well as new components introduced through mosaicism.

The Impact of Geographical Variation in Plasmodium knowlesi Apical Membrane Protein 1 (PkAMA-1) on Invasion Dynamics of P. knowlesi.

Plasmodium knowlesi has emerged as an important zoonotic parasite that causes persistent symptomatic malaria in humans. The signs and symptoms of malaria are attributed to the blood stages of the parasites, which start from the invasion of erythrocytes by the blood stage merozoites. The apical membrane protein 1 (AMA-1) plays an important role in the invasion. In this study, we constructed and expressed recombinant PkAMA-1 domain II (PkAMA-1-DII) representing the predominant haplotypes from Peninsular Malaysia and Malaysian Borneo and raised specific antibodies against the recombinant proteins in rabbits. Despite the minor amino acid sequence variation, antibodies raised against haplotypes from Peninsular Malaysia and Malaysian Borneo demonstrated different invasion inhibition (46.81% and 39.45%, respectively) to P. knowlesi A1-H.1, a reference strain derived from Peninsular Malaysia. Here, we demonstrated how a minor variation in a conserved parasite protein could cast a significant impact on parasite invasion biology, suggesting a complex host-switching of P. knowlesi from different locations. This may challenge the implementation of a standardized One Health approach against the transmission of knowlesi malaria.

Construction of Peptide Library in Mammalian Cells by dsDNA-Based Strategy.

While different display technologies, represented by phage display, have been widely used in drug discovery, they still can hardly achieve function-based peptide screening, which in most cases is performed in mammalian cells. And most attempts to screen functional peptides with mammalian platforms utilized plasmids to store coding information. Our previous work established double-stranded DNAs (dsDNAs) as innovative biological parts to implement AND-gate genetic circuits in mammalian cells. In the current study, we employ dsDNAs with terminal NNK degenerate codons to implement AND-gate genetic circuits and generate peptide libraries in mammalian cells. This dsDNA-based AND-gate (DBAG) peptide library construction strategy is easy to perform, requiring only PCR reaction and cell transfection. High-throughput sequencing (HTS) and single-cell sequencing results revealed both peptide length and amino acid sequence diversity of DBAG peptide libraries. Moreover, as a feasibility test of this strategy, we identified an MDM2-interacting peptide by applying the DBAG peptide library to a mammalian cell-based two-hybrid system. Our work establishes dsDNAs with terminal degenerate codons as biological parts to build peptide libraries in mammalian cells, which may have great application potential in the future.

Interaction of C-terminal p53 isoforms depends strongly upon DNA sequence and topology

The p53 protein is a key tumor suppressor and the most commonly mutated and down-regulated protein in human tumors. It functions mainly through interaction with DNA, and p53 acts as a transcription factor that recognizes the so-called p53 target sites on the promoters of various genes. P53 has been shown to exist as many isoforms, including three C-terminal isoforms that are produced by alternative splicing. Because the C-terminal domain is responsible for sequence-nonspecific binding and regulation of p53 binding, we have analyzed DNA recognition by these C-terminal isoforms. Using atomic force microscopy, we show for the first time that all C-terminal isoforms recognize superhelical DNA. It is particularly noteworthy that a sequence-specific p53 consensus binding site is bound by p53α and β isoforms with similar affinities, whilst p53α shows higher binding to a quadruplex sequence than both p53β and p53γ, and p53γ loses preferential binding to both the consensus binding sequence and the quadruplex-forming sequence. These results show the important role of the variable p53 C-terminal amino acid sequences for DNA recognition.

493. Defining background shared antibody sequences between unrelated healthy individuals (public clonotypes) to support future studies on specific infectious disease related conditions.

Abstract Background Public clonotypes, antibodies against specific antigens in unrelated individuals that have genetic similarities, have been shown in a variety of infections, including SARS-CoV-2 and HIV. Likely, there are shared antibody responses between individuals for many infections. To explore antibody responses that would coincide with specific infectious diseases that may set off chronic illnesses, such as Multiple Sclerosis or Alzheimer's disease, defining the background shared clonotypes is needed to differentiate disease from normal background public clonotype responses. Methods Heavy chain variable sequences were retrieved from public biorepositories (Bioproject PRJNA486667) composed of 43 healthy persons, and two groups of HIV infected persons; 114 with broadly neutralizing antibodies and 91 without broadly neutralizing antibodies. We utilized the Immcantation package of software run on our SUNY Buffalo computational cluster. After PRESTo annotation, duplicate sequences were collapsed and sequences of only single counts were removed. Clonal groups were determined using ChangeO requiring IGHV, IGHJ, and CDR3 amino acid sequence to be perfectly matched. Figures and statistics were generated with immcantation, excel, and graphpad prism 8. Results 244850 heavy chain sequences from 43 healthy controls were compared for exact matches to predicted germline variable segment and CDR3 amino acid sequence and identified 0.23% as public clonotypes. Comparison to 205 HIV + individuals (a total of 1.4 million comparative sequences) showed that 2.35% of heavy chain sequences were seen in more than one individual. Generally, public clonotypes had shorter CDR3s (peak of 9 amino acids). VH 3-9, 3-30 and 4-34 were the most commonly used variable segments in public clonotypes. Common exact match CDR3 sequences using a variety of variable sequences, including an 11 amino acid CDR3 sequence motif, were also discovered. Conclusion This early work has identified several public clonotypes that are shared among subjects who are HIV positive and otherwise healthy people. Defining the sequences commonly seen between individuals can assist in specifying antibody responses specific to disease states from larger sequence databases. Disclosures All Authors: No reported disclosures.

Genome characterization of Hirame novirhabdovirus (HIRRV) isolate CNPo2015 and transcriptome analysis of Hirame natural embryo (HINAE) cells infected with CNPo2015

Hirame novirhabdovirus (HIRRV) is a fish rhabdovirus belonging to family Rhabdoviridae, genus Novirhabdovirus, which is highly contagious and virulent, and causes hemorrhagic disease in many fish species. In the present work, the whole genome sequence of HIRRV strain CNPo2015 that previously isolated from cultured flounders was obtained using high-throughput sequencing. It consists of 10,998 nucleotides and encodes six viral proteins arranged in order of 3′-N-P-M-G-NV-L-5′. Among Novirhabdovirus, L protein of CNPo2015 possessed the lowest amino acid sequence divergence with HIRRV isolate CA 9703 and HIRRV 080113, and the highest with Snakehead rhabdovirus. Furthermore, the immune response of Hirame natural embryo (HINAE) cell line to HIRRV infection was characterized by RNA-seq, and the results showed that 1976 differentially expressed genes (DEGs) including 1219 up-regulated and 727 down-regulated genes were identified in the HINAE cells infected with HIRRV at 48 h post infection (hpi). Several KEGG pathways were significantly enriched in the viral infected cells, such as cytokine-cytokine receptor interaction, JAK-STAT signaling pathway, cell cycle, apoptosis, RIG-I-like receptors signaling pathway and P13K-AKT signaling pathway. Post viral infection, the flow cytometric Annexin V/PI assay found that apoptotic rate of HINAE cells showed a slight increase within 3 days and then the early and late apoptotic rate were significantly increased to 41 ± 2.65% and 12.37 ± 2.61% at day 4, respectively. Meanwhile, qRT-PCR results also showed that six apoptosis-related genes (BCL2L1, CASPASE 3, CASPASE 10, FAS, AKT and CDK1) were significantly upregulated. This investigation has not only enriched our knowledge of sequence difference characteristics between CNPo2015 and other Novirhabdoviruses, but also provided a data basis for deeper understanding of immune responses in flounder cells post viral infection.

Diversity and molecular epidemiology of Ostreid herpesvirus 1 in farmed Crassostrea gigas in Australia: Geographic clusters and implications for “microvariants” in global mortality events

Since 2010, mass mortality events known as Pacific oyster mortality syndrome (POMS) have occurred in Crassostrea gigas in Australia associated with Ostreid herpesvirus 1. The virus was thought to be an OsHV-1 µVar or “microvariant”, i.e. one of the dominant variants associated with POMS in Europe, but there are few data to characterize the genotype in Australia. Consequently, the genetic identity and diversity of the virus was determined to understand the epidemiology of the disease in Australia. Samples were analysed from diseased C. gigas over five summer seasons between 2011 and 2016 in POMS-affected estuaries: Georges River in New South Wales (NSW), Hawkesbury River (NSW) and Pitt Water in Tasmania. Sequencing was attempted for six genomic regions. Numerous variants were identified among these regions (n = 100 isolates) while twelve variants were identified from concatenated nucleotide sequences (n = 61 isolates). Nucleotide diversity of the seven genotypes of C region among Australian isolates (Pi 0.99 × 10−3) was the lowest globally. All Australian isolates grouped in a cluster distinct from other OsHV-1 isolates worldwide. This is the first report that Australian outbreaks of POMS were associated with OsHV-1 distinct from OsHV-1 reference genotype, µVar and other microvariants from other countries. The findings illustrate that microvariants are not the only variants of OsHV-1 associated with mass mortality events in C. gigas. In addition, there was mutually exclusive spatial clustering of viral genomic and amino acid sequence variants between estuaries, and a possible association between genotype/amino acid sequence and the prevalence and severity of POMS, as this differed between these estuaries. The sequencing findings supported prior epidemiological evidence for environmental reservoirs of OsHV-1 for POMS outbreaks in Australia.

Complete genome sequence of a novel polerovirus infecting chickpea (Cicer arietinum L.).

The complete genome sequence of a novel polerovirus identified in chickpea (C. arietinum L.) is presented. Its sequence was assembled using small RNA sequencing and assembly (sRSA) and confirmed by RT-PCR, 5' and 3' RACE, and Sanger sequencing. According to the current ICTV sequence demarcation criterion of greater than 10% amino acid (aa) sequence divergence in all gene products when compared to other poleroviruses, the newly identified polerovirus should be classified as a member of a new species, and we propose the name "chickpea leafroll virus" (CpLRV) for this virus. The genome of CpLRV is 5,770 nucleotides (nt) long and is organized into seven open reading frames (ORFs), designated as ORF0, ORF1, ORF2, ORF3a, ORF3, ORF4, and ORF5, which code for putative P0, P1, P1-P2, P3a, P3, P4, and P3-P5 proteins, respectively. The 5' untranslated region (UTR) consists of 27 nt, starting with the conserved sequence 5'-ACAAAA-3', which is typical of poleroviruses, while the 3' UTR consists of 229 nt. Phylogenetic analysis based on the aa sequences of P0, P1, P1-P2, P3, P4, and P3-P5 showed that CpLRV clustered with members of the genus Polerovirus and is closely related to chickpea chlorotic stunt virus (CpCSV) and faba bean polerovirus 1 (FBPV1). Recombination analysis suggested that CpLRV is a recombinant of two unknown viruses that share the highest nucleotide sequence similarity with FBPV1 (76.9% identity) and suakwa aphid-borne yellows virus (SAbYV) (64.8% identity). The putative recombination event was identified in the 5' region of the CpLRV genome, the region that encodes proteins P0, P1, and P1-P2. This is the first report of a polerovirus infecting chickpea in Mexico.

Effective Amino Acid Sequencing of Intact Filgrastim by Multimodal Mass Spectrometry with Topdownr.

Therapeutic proteins, known as biologicals, are an important and growing class of drugs for treatment of a series of human ailments. Amino acid sequence variants of therapeutic proteins can affect their safety and efficacy. Top-down mass spectrometry is well suited for the sequence analysis of intact therapeutic proteins. Fine-tuning of tandem mass spectrometry (MS/MS) fragmentation conditions is essential for maximizing the amino acid sequence coverage but is often time-consuming. We used topdownr, an automated and integrated multimodal approach to systematically assess high mass accuracy MS/MS fragmentation parameters to characterize filgrastim, a 19 kDa recombinant human granulocyte colony-stimulating factor used in treating neutropenia. A total of 276 different MS/MS conditions were systematically tested, including the following parameters: protein charge state, HCD and CID collision energy, ETD reaction time, ETD supplemental activation, and UVPD activation time. Stringent and accurate evaluation and annotation of the MS/MS data was achieved by requiring a fragment ion mass error of 5 ppm, considering reproducible N- and C-terminal fragment ions only, and excluding internal fragment ion assignments. We report the first EThcD and UVPD MS/MS analysis of intact filgrastim, and these two techniques combined resulted in 98% amino acid sequence coverage. By combining all tested fragmentation modes, we obtained near-complete amino acid sequence coverage (99.4%) of intact filgrastim.

Golden 2-like transcription factor contributes to the major QTL against rice black-streaked dwarf virus disease.

A major resistance QTL was identified on chromosome 6 in rice variety Wuke; both overexpression and knockdown experiments confirmed that OsGLK1 is the candidate gene for association with Rice black-streaked dwarf virus disease. Rice black-streaked dwarf virus disease is one of the most destructive rice viral diseases in China and East Asia. Progress has been limited in RBSDVD resistance breeding due to inadequate knowledge on the underlying functional genes. In this study, a major QTL for RBSDV (rice black-streaked dwarf virus) independent of SBPH (small brown planthopper) resistance was mapped in a 1.8Mb interval on chromosome 6 by using an F2:3 population originated from resistant rice variety Wuke. Representative transcripts within this region were analysed and three genes showing amino acid sequence variation in functional domains were selected for transformation. Overexpression experiments showed that one gene exhibited significant enhanced resistance compared to control lines, encoding protein involving Myb domain and probable transcription factor Golden 2-like1 (GLK1). Furthermore, OsGLK1 knockdown rice lines were investigated and the resistance ability was significantly declined without this gene compared to the wild type. Taken together, both overexpression and knockdown experiments strongly suggested that OsGLK1 plays an important role for RBSDV resistance and contributes to the major QTL. The study paves the way for elucidating the molecular mechanism underlying RBSDVD resistance and the molecular markers associated with OsGLK1 may be used for marker-assisted selection.

Cytochrome P450 oxidoreductase variant A503V contributes to the increased CYP3A5 activity with tacrolimus in vitro

ABSTRACT Background Tacrolimus is a immunosuppressant drug in transplantation with a narrow therapeutic range and high pharmacokinetic variability. Clinical studies have revealed that POR*28 contributes enhanced tacrolimus clearance in CYP3A5 expressers. However, it remains unknown that how exactly the POR polymorphism could influence the metabolism of tacrolimus via CYP3A5 in vitro. Research design & methods A503V is an amino acid sequence variant encoded by POR*28. Wild-type (WT) and A503V POR, with WT CYP3A5 were expressed in recombinant HepG2 cells and reconstituted proteins. Michaelis constant (K m) and maximum velocity (V max) of CYP3A5 with tacrolimus as substrates were determined, and catalytic efficiency is expressed as V max/K m. Results Both WT and A503V POR down-regulated the CYP3A5 mRNA expression, and WT POR rather than A503V down-regulated the protein expression of CYP3A5 in recombinant HepG2 cells. Compared with WT POR, A503V increased metabolism of tacrolimus by CYP3A5 in both cellular and protein levels. Conclusion A503V can affect CYP3A5-catalyzed tacrolimus metabolism in vitro, which suggests that A503V has the potential to serve as a biomarker for tacrolimus treatment in transplantation recipients.

Molecular characterization and phylogenetic relationships among Rhynchophorus sp. haplotypes in Makkah Al-Mukarramah Region-KSA

The study aims at detecting and characterizing haplotypes of red palm weevil (RPW) Rhynchophorus sp. in the Western region of Saudi Arabia based on the cytochrome oxidase subunit I (COI) gene sequence. The results indicated the occurrence of 17 nucleotide substitutions, of which three were nonsynonymous (NS). These three NS substitutions resulted in the variation in amino acid sequence in three positions, out of 133. These amino acids are isoleucine/valine, glycine/serine, and arginine/histidine. Based on the chemical properties of the cytochrome C oxidase (COX) enzyme, it is likely that the change at the first position has no effect, while changes at the other two positions can affect the chemical properties of the enzyme. At the three-dimensional (3D) level, the first two positions exist at the border or inside loop 3–4 of the enzyme, while the third position exists inside loop 4–5. These two loops influence the folding pattern of the enzyme, thus, likely affecting the function of the enzyme. However, it is unlikely that variations in the three positions will affect the binding ability of the heme group, which promotes the action of the COX enzyme in the electron transport chain. Variations in chemical properties and 3D structure of COX enzyme might be an evolutionary process (positive selection) that promotes in-time and in-site adaptation to the insect. In conclusion, this study can be helpful in pest management programs and in tracing RPW geographic spread and migration in Saudi Arabia.

Genomic analysis of arginine vasopressin gene in riverine buffalo reveals its potential association with silent estrus behavior.

Poor estrus expression behavior causes suboptimal reproductive efficiency through poor conception rate. Various signaling pathways are involved in estrus expression but arginine vasopressin (AVP) gene with oxytocin predominantly regulates estrus behavior. This study aimed to perform genomic characterization and evolutionary dynamics of AVP gene through association testing of the novel polymorphic loci and comparative genomic analysis to explore the potential effect of AVP gene on estrus behavior of Nili-Ravi buffaloes. 198 Nili-Ravi buffaloes were screened for the quest of novel polymorphism in the AVP gene. In exon-1, five polymorphic sites were detected including deletion of two (c.47delA and c.57delA) nucleotides that caused drastic variation in subsequent amino acid sequence due to frame shift including functional short peptide of nine residues. The 3-D structure revealed a loss of transmembrane loop between 16 and 31 residues in Nili-Ravi buffalo AVP protein sequence, suggesting that missing loop apparently reduced the gene functionality in Nili-Ravi buffalo by inhibiting cellular reactions and muting the animal estrus cyclicity. Three polymorphisms detected in AVP gene were significantly associated with silent estrus (P <0.05). The comparative genomic analysis revealed that AVP gene is present on chromosome 14 having one conserved motif (Neurohypophysial) in buffalo. This study suggested the potential use of polymorphic sites as promising genetic markers for selection of buffaloes with pronounced estrus expression.

Comparative analysis of human and bovine thyroglobulin structures

In biology, evolutionary conserved protein sequences show homologous physiological phenotypes in their structures and functions. If the protein has a vital function, its sequence is usually conserved across the species. However, in highly conserved protein there still remains small differences across the species. Upon protein–protein interaction (PPI), it is observed that the conserved proteins can have different binding partners that are considered to be caused by the small sequence variations in a specific domain. Thyroglobulin (TG) is the most commonly found protein in the thyroid gland of vertebrates and serves as the precursor of the thyroid hormones, tetraiodothyronine and triiodothyronine that are critical for growth, development and metabolism in vertebrates. In this study, we comparatively analyzed the sequences and structures of the highly conserved regions of TG from two different species in relation to their PPIs. In order to do so, we employed SIM for sequence alignment, STRING for PPI analysis and cryo-electron microscopy for 3D structural analysis. Our Cryo-EM model for TG of Bos taurus determined at 7.1 Å resolution fitted well with the previously published Cryo-EM model for Homo sapiens TG. By demonstrating overall structural homology between TGs from different species, we address that local amino acid sequence variation is sufficient to alter PPIs specific for the organism. We predict that our result will contribute to a deeper understanding in the evolutionary pattern applicable to many other proteins.

Generalizable strategy to analyze domains in the context of parent protein architecture: A CheW case study.

Domains are the three-dimensional building blocks of proteins. An individual domain can occur in a variety of domain architectures that perform unique functions and are subject to different evolutionary selective pressures. We describe an approach to evaluate the variability in amino acid sequences of a single domain across architectural contexts. The ability to distinguish different evolutionary outcomes of one protein domain can help determine whether existing knowledge about a specific domain will apply to an uncharacterized protein, lead to insights and hypotheses about function, and guide experimental priorities. We developed and tested our approach on CheW-like domains (PF01584), which mediate protein/protein interactions and are difficult to compare experimentally. CheW-like domains occur in CheW scaffolding proteins, CheA kinases, and CheV proteins that regulate bacterial chemotaxis. We analyzed 16 domain architectures that included 94% of all CheW-like domains found in nature. We identified six Classes of CheW-like domains with presumed functional differences. CheV and most CheW proteins contained Class 1 domains, whereas some CheW proteins contained Class 6 (~20%) or Class 2 (~1%) domains instead. Most CheA proteins contained Class 3 domains. CheA proteins with multiple Hpt domains contained Class 4 domains. CheA proteins with two CheW-like domains contained one Class 3 and one Class 5. We also created SimpLogo, an innovative method for visualizing amino acid composition across large sets of multiple sequence alignments of arbitrary length. SimpLogo offers substantial advantages over standard sequence logos for comparison and analysis of related protein sequences. The R package for SimpLogo is freely available.

Genome-Wide Association Studies Reveal Candidate Genes Related to Stem Diameter in Cucumber (Cucumis sativus L.).

The stem diameter, an important agronomic trait, affects cucumber growth and yield. However, no genes responsible for cucumber stem diameter have been identified yet. In this study, the stem diameter of 88 cucumber core germplasms were measured in spring 2020, autumn 2020 and autumn 2021, and a genome-wide association study (GWAS) was carried out based on the gene sequence and stem diameter of core germplasms. A total of eight loci (gSD1.1, gSD2.1, gSD3.1, gSD3.2, gSD4.1, gSD5.1, gSD5.2, and gSD6.1) significantly associated with cucumber stem diameter were detected. Of these, five loci (gSD1.1, gSD2.1, gSD3.1, gSD5.2, and gSD6.1) were repeatedly detected in two or more seasons and were considered as robust and reliable loci. Based on the linkage disequilibrium sequences of the associated SNP loci, 37 genes were selected. By further investigating the five loci via analyzing Arabidopsis homologous genes and gene haplotypes, five genes (CsaV3_1G028310, CsaV3_2G006960, CsaV3_3G009560, CsaV3_5G031320, and CsaV3_6G031260) showed variations in amino acid sequence between thick stem lines and thin stem lines. Expression pattern analyses of these genes also showed a significant difference between thick stem and thin stem lines. This study laid the foundation for gene cloning and molecular mechanism study of cucumber stem development.