Protein Sequence Alignment Analysis Research Articles

PKHD1L1, a gene involved in the stereocilia coat, causes autosomal recessive nonsyndromic hearing loss

Identification of genes associated with nonsyndromic hearing loss is a crucial endeavor given the substantial number of individuals who remain without a diagnosis after even the most advanced genetic testing. PKHD1L1 was established as necessary for the formation of the cochlear hair-cell stereociliary coat and causes hearing loss in mice and zebrafish when mutated. We sought to determine if biallelic variants in PKHD1L1 also cause hearing loss in humans. Exome sequencing was performed on DNA of four families segregating autosomal recessive nonsyndromic sensorineural hearing loss. Compound heterozygous p.[(Gly129Ser)];p.[(Gly1314Val)] and p.[(Gly605Arg)];p[(Leu2818TyrfsTer5)], homozygous missense p.(His2479Gln) and nonsense p.(Arg3381Ter) variants were identified in PKHD1L1 that were predicted to be damaging using in silico pathogenicity prediction methods. In vitro functional analysis of two missense variants was performed using purified recombinant PKHD1L1 protein fragments. We then evaluated protein thermodynamic stability with and without the missense variants found in one of the families and performed a minigene splicing assay for another variant. In silico molecular modeling using AlphaFold2 and protein sequence alignment analysis were carried out to further explore potential variant effects on structure. In vitro functional assessment indicated that both engineered PKHD1L1 p.(Gly129Ser) and p.(Gly1314Val) mutant constructs significantly reduced the folding and structural stabilities of the expressed protein fragments, providing further evidence to support pathogenicity of these variants. Minigene assay of the c.1813G>A p.(Gly605Arg) variant, located at the boundary of exon 17, revealed exon skipping leading to an in-frame deletion of 48 amino acids. In silico molecular modeling exposed key structural features that might suggest PKHD1L1 protein destabilization. Multiple lines of evidence collectively associate PKHD1L1 with nonsyndromic mild–moderate to severe sensorineural hearing loss. PKHD1L1 testing in individuals with mild–moderate hearing loss may identify further affected families.

Alignment of human aquaporin 4 and ß-amyloid proteins may indicate involvement of ß-amyloid in brain water homeostasis and prevention of brain edema.

Alignment of human aquaporin 4 and ß-amyloid proteins may indicate involvement of ß-amyloid in brain water homeostasis and prevention of brain edema.

Low grain weight, a new allele of BRITTLE CULM12, affects grain size through regulating GW7 expression in rice.

Grain weight is a major determinant in rice yield, which is tightly associated with grain size. However, the underlying molecular mechanisms that control this trait remain unclear. Here, we report a rice (Oryza sativa) mutant, low grain weight (lgw), which shows that reduced grain length is caused by decreased cell elongation and proliferation. Map-based cloning revealed that all mutant phenotypes resulted from a nine-base pair (bp) deletion in LGW, which encodes the kinesin-like protein BRITTLE CULM12 (BC12). Protein sequence alignment analysis revealed that the mutation site was located at the nuclear localization signal (NLS) of LGW/BC12, resulting in the lgw protein not being located in the nucleus. LGW is preferentially expressed in both culms and roots, as well as in the early developing panicles. Overexpression of LGW increased the grain length, indicating that LGW is a positive regulator for regulating grain length. In addition, LGW/BC12 is directly bound to the promoter of GW7 and activates its expression. Elevating the GW7 expression levels in lgw plants rescued the small grain size phenotype. We conclude that LGW regulates grain development by directly binding to the GW7 promoter and activating its expression. Our findings revealed that LGW plays an important role in regulating grain size, and manipulation of this gene provides a new strategy for regulating grain weight in rice.

Expression Pattern and Functional Characterization of PISTILLATA Ortholog Associated With the Formation of Petaloid Sepals in Double-Flower Eriobotrya japonica (Rosaceae).

Double-flower Eriobotrya japonica, of which one phenotype is homeotic transformation of sepals into petals, is a new germplasm for revealing the molecular mechanisms underlying the floral organ transformation. Herein, we analyzed the sequence, expression pattern and functional characterization of EjPI, which encoded a B-class floral homeotic protein referred to as PISTILLATA ortholog, from genetically cognate single-flower and double-flower E. japonica. Phylogenetic analysis suggested that the EjPI gene was assigned to the rosids PI/GLO lineage. Analysis of protein sequence alignments showed that EjPI has typical domains of M, I, K, and C, and includes a distinctive PI motif at the C-terminal region. Compared with asterids PI/GLO lineage, the K1 and K3 subdomains of EjPI both contain a single amino acid difference. Subcellular localization of EjPI was determined to be in the nucleus. Expression pattern analysis revealed that EjPI expressed not only in petals, filament, and anther in single-flower E. japonica, but also in petaloid sepals in double-flower E. japonica. Meanwhile, there were high correlation between EjPI transcript level and petaloid area within a sepal. Furthermore, 35S::EjPI transgenic wild-type Arabidopsis caused the homeotic transformation of the first whorl sepals into petaloid sepals. Ectopic expression of EjPI in transgenic pi-1 mutant Arabidopsis rescued normal petals and stamens. These results suggest expression pattern of EjPI is associated with the formation of petaloid sepal. Our study provides the potential application of EjPI for biotechnical engineering to create petaloid sepals or regulate floral organ identity in angiosperms.

A benchmark study of sequence alignment methods for protein clustering

BackgroundProtein sequence alignment analyses have become a crucial step for many bioinformatics studies during the past decades. Multiple sequence alignment (MSA) and pair-wise sequence alignment (PSA) are two major approaches in sequence alignment. Former benchmark studies revealed drawbacks of MSA methods on nucleotide sequence alignments. To test whether similar drawbacks also influence protein sequence alignment analyses, we propose a new benchmark framework for protein clustering based on cluster validity. This new framework directly reflects the biological ground truth of the application scenarios that adopt sequence alignments, and evaluates the alignment quality according to the achievement of the biological goal, rather than the comparison on sequence level only, which averts the biases introduced by alignment scores or manual alignment templates. Compared with former studies, we calculate the cluster validity score based on sequence distances instead of clustering results. This strategy could avoid the influence brought by different clustering methods thus make results more dependable.ResultsResults showed that PSA methods performed better than MSA methods on most of the BAliBASE benchmark datasets. Analyses on the 80 re-sampled benchmark datasets constructed by randomly choosing 90% of each dataset 10 times showed similar results.ConclusionsThese results validated that the drawbacks of MSA methods revealed in nucleotide level also existed in protein sequence alignment analyses and affect the accuracy of results.

Power law tails in phylogenetic systems

Covariance analysis of protein sequence alignments uses coevolving pairs of sequence positions to predict features of protein structure and function. However, current methods ignore the phylogenetic relationships between sequences, potentially corrupting the identification of covarying positions. Here, we use random matrix theory to demonstrate the existence of a power law tail that distinguishes the spectrum of covariance caused by phylogeny from that caused by structural interactions. The power law is essentially independent of the phylogenetic tree topology, depending on just two parameters-the sequence length and the average branch length. We demonstrate that these power law tails are ubiquitous in the large protein sequence alignments used to predict contacts in 3D structure, as predicted by our theory. This suggests that to decouple phylogenetic effects from the interactions between sequence distal sites that control biological function, it is necessary to remove or down-weight the eigenvectors of the covariance matrix with largest eigenvalues. We confirm that truncating these eigenvectors improves contact prediction.

Sequence alignment analysis of proteins involved in platelet-endothelial cell interaction identifies molecular incompatibilities between Homo sapiens and Sus scrofa

Background: Platelets play a vital role in acute humoral xenograft rejection (AHXR) in pig-to-primate xenotransplantation, presenting as microvascular thrombosis in the graft and/or consumptive coagulopathy in the recipient. Adhesion and aggregation of primate platelets to the activated vascular endothelial cells through sequential binding of ligands on endothelial cells and subendothelial matrix ultimately trigger a complex biological process of prothrombotic signaling cascades. Increasing evidence suggests that the molecular incompatibilities in effector molecules across species may partially contribute to dysregulated microvascular thrombosis in xenografts. Method: We selected candidate genes from the NCBI database with keywords: platelet-endothelium interaction，platelet adhension, platelet aggregation, and subendothelial matrix ligands. Pair-wise amino acid alignments were made using the Emboss Needle method. Emboss needle created optimal global alignment of the amino acid sequences of human genes and pig genes using ClustalW2.Results: Most of the proteins involved in platelet-EC interaction in Homo sapiens share high sequence similarity with their homologues in Sus scrofa. Cytokines that potentially induce endothelial damage (such as CD40L, TNF-α) were highly conserved between Homo sapiens and Sus scrofa. Some endothelium-derived cytokines (such as IL-8, CCL2, CCL5) that can induce platelet activation or enhance aggregation share high sequence similarity between Homo sapiens and sus scrofa. Some regulators that potentially transduce inhibitory signaling to control platelet activation or complement activation have relatively poor sequence identity between Homo sapiens and Sus scrofa, and some even lack their homologues in Sus scrofa.Conclusion: The characteristic of sequence similarity of proteins involved in platelet-EC interaction indicated the molecular incompatibilities between humans (Homo sapiens) and pigs (Sus scrofa). This study provides a clue for explanation of excessive platelet activity in pig-to-primate xenotransplantation.

The Phylogenetic Analysis and Identification of a Novel Remorin Member from Rice(<i>Oryza sativa</i> L) by Proteomics under Salt Stress

Rice root proteomics research identified to a new remorin1 protein,named OsREM1,which significantly increased (approximately 1.7-fold) compared with control under salt stress by Two-dimensional electrophoresis. It is found that the protein molecular weight is 29 kD, isoelectric point is 4.54. The bioinformatics analysis showed that: (1) The protein contains 195 amino acids, whose theoretical molecular weight is 21.3 kD, whose theoretical isoelectric point is 5.06. There are some differences with the experimental results, which may be the result of post-translational modification, such as glycosylation, phosphorylation; (2) The different crop remorin protein sequence alignment analysis found that the same rate of the amino acids is more than 50%, indicated that remorins are highly conserved protein; (3) Phylogenetic analysis showed OsREM1 had a low sequence identity with potao and tomato REM1, which is the agreement with traditional classification.

Extracellular domains play different roles in gap junction formation and docking compatibility

GJ (gap junction) channels mediate direct intercellular communication and play an important role in many physiological processes. Six connexins oligomerize to form a hemichannel and two hemichannels dock together end-to-end to form a GJ channel. Connexin extracellular domains (E1 and E2) have been shown to be important for the docking, but the molecular mechanisms behind the docking and formation of GJ channels are not clear. Recent developments in atomic GJ structure and functional studies on a series of connexin mutants revealed that E1 and E2 are likely to play different roles in the docking. Non-covalent interactions at the docking interface, including hydrogen bonds, are predicted to form between interdocked extracellular domains. Protein sequence alignment analysis on the docking compatible/incompatible connexins indicate that the E1 domain is important for the formation of the GJ channel and the E2 domain is important in the docking compatibility in heterotypic channels. Interestingly, the hydrogen-bond forming or equivalent residues in both E1 and E2 domains are mutational hot spots for connexin-linked human diseases. Understanding the molecular mechanisms of GJ docking can assist us to develop novel strategies in rescuing the disease-linked connexin mutants.

Protein Sequence Alignment Analysis by Local Covariation: Coevolution Statistics Detect Benchmark Alignment Errors

The use of sequence alignments to understand protein families is ubiquitous in molecular biology. High quality alignments are difficult to build and protein alignment remains one of the largest open problems in computational biology. Misalignments can lead to inferential errors about protein structure, folding, function, phylogeny, and residue importance. Identifying alignment errors is difficult because alignments are built and validated on the same primary criteria: sequence conservation. Local covariation identifies systematic misalignments and is independent of conservation. We demonstrate an alignment curation tool, LoCo, that integrates local covariation scores with the Jalview alignment editor. Using LoCo, we illustrate how local covariation is capable of identifying alignment errors due to the reduction of positional independence in the region of misalignment. We highlight three alignments from the benchmark database, BAliBASE 3, that contain regions of high local covariation, and investigate the causes to illustrate these types of scenarios. Two alignments contain sequential and structural shifts that cause elevated local covariation. Realignment of these misaligned segments reduces local covariation; these alternative alignments are supported with structural evidence. We also show that local covariation identifies active site residues in a validated alignment of paralogous structures. Loco is available at https://sourceforge.net/projects/locoprotein/files/

Aspergillus nidulans ArfB Plays a Role in Endocytosis and Polarized Growth

Filamentous fungi undergo polarized growth throughout most of their life cycles. The Spitzenkörper is an apical organelle composed primarily of vesicles that is unique to filamentous fungi and is likely to act as a vesicle supply center for tip growth. Vesicle assembly and trafficking are therefore important for hyphal growth. ADP ribosylation factors (Arfs), a group of small GTPase proteins, play an important role in nucleating vesicle assembly. Little is known about the role of Arfs in filamentous hyphal growth. We found that Aspergillus nidulans is predicted to encode six Arf family proteins. Analysis of protein sequence alignments suggests that A. nidulans ArfB shares similarity with ARF6 of Homo sapiens and Arf3p of Saccharomyces cerevisiae. An arfB null allele (arfB disrupted by a transposon [arfB::Tn]) was characterized by extended isotropic growth of germinating conidia followed by cell lysis or multiple, random germ tube emergence, consistent with a failure to establish polarity. The mutant germ tubes and hyphae that do form initially meander abnormally off of the axis of polarity and frequently exhibit dichotomous branching at cell apices, consistent with a defect in polarity maintenance. FM4-64 staining of the arfB::Tn strain revealed that another phenotypic characteristic seen for arfB::Tn is a reduction and delay in endocytosis. ArfB is myristoylated at its N terminus. Green fluorescent protein-tagged ArfB (ArfB::GFP) localizes to the plasma membrane and endomembranes and mutation (ArfB(G2A)::GFP) of the N-terminal myristoylation motif disperses the protein to the cytoplasm rather than to the membranes. These results demonstrate that ArfB functions in endocytosis to play important roles in polarity establishment during isotropic growth and polarity maintenance during hyphal extension.