Homologous Positions Research Articles

Functional similarities and differences among subunits of the nascent polypeptide-associated complex (NAC) of Saccharomyces cerevisiae.

Protein factors bind ribosomes near the tunnel exit, facilitating protein trafficking and folding. In eukaryotes, the heterodimeric nascent polypeptide-associated complex (NAC) is the most abundant - equimolar to ribosomes. Saccharomyces cerevisiae has a minor β-type subunit (Nacβ2) in addition to abundant Nacβ1, and therefore two NAC heterodimers, α/β1 and α/β12. The additional beta NAC gene arose at the time of the whole genome duplication that occurred in the S. cerevisiae lineage. Nacβ2 has been implicated in regulating the fate of mRNA encoding ribosomal protein Rpl4 during translation via its interaction with the Caf130 subunit of the regulatory CCR4-Not complex. We found that Nacβ2 residues just C-terminal to the globular domain are required for its interaction with Caf130 and its negative effect on growth of cells lacking Acl4, the specialized chaperone for Rpl4. Substitution of these Nacβ2 residues at homologous positions in Nacβ1 results in a chimeric protein that interacts with Caf130 and slows the growth of ∆acl4 cells lacking Nacβ2. Furthermore, alteration of residues in the N-terminus of Nacβ2 or chimeric Nacβ1 previously shown to affect ribosome binding overcomes the growth defect of ∆acl4. Our results are consistent with a model in which Nacβ2's ribosome association per se, or its precise positioning, is necessary for productive recruitment of CCR4-Not via its interaction with the Caf130 subunit to drive Rpl4 mRNA degradation.

Peptide Nucleic Acid-Mediated Regulation of CRISPR-Cas9 Specificity.

Although CRISPR-Cas9 gene therapies have proven to be a powerful tool across many applications, improvements are necessary to increase the specificity of this technology. Cas9 cutting in off-target sites remains an issue that limits CRISPR's application in human-based therapies. Treatment of autosomal dominant diseases also remains a challenge when mutant alleles differ from the wild-type sequence by only one base pair. Here, we utilize synthetic peptide nucleic acids (PNAs) that bind selected spacer sequences in the guide RNA (gRNA) to increase Cas9 specificity up to 10-fold. We interrogate variations in PNA length, binding position, and degree of homology with the gRNA. Our findings reveal that PNAs bound in the region distal to the protospacer adjacent motif (PAM) site effectively enhance specificity in both on-target/off-target and allele-specific scenarios. In addition, we demonstrate that introducing deliberate mismatches between PNAs bound in the PAM-proximal region of the gRNA can modulate Cas9 activity in an allele-specific manner. These advancements hold promise for addressing current limitations and expanding the therapeutic potential of CRISPR technology.

TRIM71 mutations cause a neurodevelopmental syndrome featuring ventriculomegaly and hydrocephalus.

Congenital hydrocephalus (CH), characterized by cerebral ventriculomegaly, is one of the most common reasons for pediatric brain surgery. Recent studies have implicated lin-41 (lineage variant 41)/TRIM71 (tripartite motif 71) as a candidate CH risk gene, however, TRIM71 variants have not been systematically examined in a large patient cohort or conclusively linked with an OMIM syndrome. Through cross-sectional analysis of the largest assembled cohort of patients with cerebral ventriculomegaly, including neurosurgically-treated CH (totaling 2,697 parent-proband trios and 8,091 total exomes), we identified 13 protein-altering de novo variants (DNVs) in TRIM71 in unrelated children exhibiting variable ventriculomegaly, CH, developmental delay, dysmorphic features, and other structural brain defects including corpus callosum dysgenesis and white matter hypoplasia. Eight unrelated patients were found to harbor arginine variants, including two recurrent missense DNVs, at homologous positions in RPXGV motifs of different NHL domains. Seven additional patients with rare, damaging, unphased or transmitted variants of uncertain significance were also identified. NHL-domain variants of TRIM71 exhibited impaired binding to the canonical TRIM71 target CDKN1A; other variants failed to direct the subcellular localization of TRIM71 to processing bodies. Single-cell transcriptomic analysis of human embryos revealed expression of TRIM71 in early first-trimester neural stem cells of the brain. These data show TRIM71 is essential for human brain morphogenesis and that TRIM71 mutations cause a novel neurodevelopmental syndrome featuring ventriculomegaly and CH.

Coronavirus M Protein Trafficking in Epithelial Cells Utilizes a Myosin Vb Splice Variant and Rab10.

The membrane (M) glycoprotein of coronaviruses (CoVs) serves as the nidus for virion assembly. Using a yeast two-hybrid screen, we identified the interaction of the cytosolic tail of Murine Hepatitis Virus (MHV-CoV) M protein with Myosin Vb (MYO5B), specifically with the alternative splice variant of cellular MYO5B including exon D (MYO5B+D), which mediates interaction with Rab10. When co-expressed in human lung epithelial A549 and canine kidney epithelial MDCK cells, MYO5B+D co-localized with the MHV-CoV M protein, as well as with the M proteins from Porcine Epidemic Diarrhea Virus (PEDV-CoV), Middle East Respiratory Syndrome (MERS-CoV) and Severe Acute Respiratory Syndrome 2 (SARS-CoV-2). Co-expressed M proteins and MYO5B+D co-localized with endogenous Rab10 and Rab11a. We identified point mutations in MHV-CoV M that blocked the interaction with MYO5B+D in yeast 2-hybrid assays. One of these point mutations (E121K) was previously shown to block MHV-CoV virion assembly and its interaction with MYO5B+D. The E to K mutation at homologous positions in PEDV-CoV, MERS-CoV and SARS-CoV-2 M proteins also blocked colocalization with MYO5B+D. The knockdown of Rab10 blocked the co-localization of M proteins with MYO5B+D and was rescued by re-expression of CFP-Rab10. Our results suggest that CoV M proteins traffic through Rab10-containing systems, in association with MYO5B+D.

Structural role of K224 in taniborbactam inhibition of NDM-1.

Taniborbactam (TAN; VNRX-5133) is a novel bicyclic boronic acid β-lactamase inhibitor (BLI) being developed in combination with cefepime (FEP). TAN inhibits both serine and some metallo-β-lactamases. Previously, the substitution R228L in VIM-24 was shown to increase activity against oxyimino-cephalosporins like FEP and ceftazidime (CAZ). We hypothesized that substitutions at K224, the homologous position in NDM-1, could impact FEP/TAN resistance. To evaluate this, a library of codon-optimized NDM K224X clones for minimum inhibitory concentration (MIC) measurements was constructed; steady-state kinetics and molecular docking simulations were next performed. Surprisingly, our investigation revealed that the addition of TAN restored FEP susceptibility only for NDM-1, as the MICs for the other 19 K224X variants remained comparable to those of FEP alone. Moreover, compared to NDM-1, all K224X variants displayed significantly lower MICs for imipenem, tebipenem, and cefiderocol (32-, 133-, and 33-fold lower, respectively). In contrast, susceptibility to CAZ was mostly unaffected. Kinetic assays with the K224I variant, the only variant with hydrolytic activity to FEP comparable to NDM-1, confirmed that the inhibitory capacity of TAN was modestly compromised (IC50 0.01 µM vs 0.14 µM for NDM-1). Lastly, structural modeling and docking simulations of TAN in NDM-1 and in the K224I variant revealed that the hydrogen bond between TAN's carboxylate with K224 is essential for the productive binding of TAN to the NDM-1 active site. In addition to the report of NDM-9 (E149K) as FEP/TAN resistant, this study demonstrates the fundamental role of single amino acid substitutions in the inhibition of NDM-1 by TAN.

Noncovalent binding of carbofuran to acetylcholinesterase from Homo sapiens, Danio rerio, Apis mellifera and Caenorhabditis elegans: Homology modelling, molecular docking and dynamics, and quantum biochemistry description

Although various regulatory agencies have banned or severely restricted the use of carbofuran (CAR), recent reports indicate the presence of CAR residues in both cultivated and wild areas. This pesticide is a potent inhibitor of acetylcholinesterase (AChE), which acts by preventing the hydrolysis of acetylcholine (ACh). Given the critical role of AChE::ACh in the proper functioning of the nervous system, we thought it appropriate to investigate the binding of CAR to AChEs from Homo sapiens, Danio rerio, Apis mellifera, and Caenorhabditis elegans using homology modelling, molecular docking, molecular dynamics, and quantum biochemistry. Molecular docking and dynamics results indicated peculiar structural behavior in each AChE::CAR system. Quantum biochemistry results showed similar affinities for all complexes, confirming the description of carbofuran as a broad-spectrum pesticide, and have a limited correlation with IC50 values. We found the following decreasing affinity order of AChE species: H. sapiens > A. mellifera > C. elegans > D. rerio. The computational results suggest that CAR occupies different pockets in the AChEs studied. In addition, our results showed that CAR binds to hsAChE and ceAChE in a very similar manner: it has high affinities for the same subsites in both species and forms hydrogen bonds with residues (hsTYR124 and ceTRP107) occupying homologous positions in the peripheral site. This suggests that this nematode is a potential model to evaluate the toxicity of carbamates, even though the sequence identity between them is only 41 %. Interestingly, we also observed that the catalytic histidines of drAChE and amAChE exhibited favorable contacts with carbofuran, suggesting that the non-covalent binding of carbofuran to these proteins may promote faster carbamylation rates than the binding modes to human and worm acetylcholinesterases. Our computational results provide a better understanding of the binding mechanisms in these complexes, as well as new insights into the mechanism of carbamylation.

Comprehensive Analysis of CircRNA Expression Profiles in Multiple Tissues of Pigs.

Circular RNAs (circRNAs) are a class of non-coding RNAs with diverse functions, and previous studies have reported that circRNAs are involved in the growth and development of pigs. However, studies about porcine circRNAs over the past few years have focused on a limited number of tissues. Based on 215 publicly available RNA sequencing (RNA-seq) samples, we conducted a comprehensive analysis of circRNAs in nine pig tissues, namely, the gallbladder, heart, liver, longissimus dorsi, lung, ovary, pituitary, skeletal muscle, and spleen. Here, we identified a total of 82,528 circRNAs and discovered 3818 novel circRNAs that were not reported in the CircAtlas database. Moreover, we obtained 492 housekeeping circRNAs and 3489 tissue-specific circRNAs. The housekeeping circRNAs were enriched in signaling pathways regulating basic biological tissue activities, such as chromatin remodeling, nuclear-transcribed mRNA catabolic process, and protein methylation. The tissue-specific circRNAs were enriched in signaling pathways related to tissue-specific functions, such as muscle system process in skeletal muscle, cilium organization in pituitary, and cortical cytoskeleton in ovary. Through weighted gene co-expression network analysis, we identified 14 modules comprising 1377 hub circRNAs. Additionally, we explored circRNA-miRNA-mRNA networks to elucidate the interaction relationships between tissue-specific circRNAs and tissue-specific genes. Furthermore, our conservation analysis revealed that 19.29% of circRNAs in pigs shared homologous positions with their counterparts in humans. In summary, this extensive profiling of housekeeping, tissue-specific, and co-expressed circRNAs provides valuable insights into understanding the molecular mechanisms of pig transcriptional expression, ultimately deepening our understanding of genetic and biological processes.

Informal arts engagement programs in disadvantaged schools: student aspirations and creative limits

ABSTRACT It is widely argued that the arts have a range of cultural, economic, and educational benefits. However, under state austerity arts curricula are devalued in favour of industry skills. To address this gap in arts education, a new type of student focussed informal arts engagement program has emerged. This article draws on a qualitative study of disadvantaged Australian secondary students’ experience of an arts engagement program and explores their experiences through Bourdieusian concepts. We observe how students’ homologous position allowed an immediate appreciation of the arts and note how their habitus frame the arts as ‘work’ and as a technical accomplishment. An illusio in the arts as career emerged from these understandings. We suggest informal arts programs act as a collective gift within a weak cycle of reciprocity, but without expanded in-school opportunities fall short in offering students new ways of understanding the place and value of the arts.

The structured coalescent in the context of gene copy number variation

The Structured Coalescent was introduced to describe the coalescent process in spatially subdivided populations with migration. Here, we re-interpret migration routes of individuals in the original model as “migration routes” of single genes in tandemly arranged gene arrays. A gene copy may change its position within the array via unequal recombination. Hence, in a coalescent framework, two copies sampled from two chromosomes may coalesce only if they are at exactly homologous positions. Otherwise, one or multiple recombination events have to occur before they can coalesce, thereby increasing mean coalescence time and expected genetic diversity among the copies in a gene array.We explicitly calculate the transition probabilities on these routes backward in time. We simulate the structured coalescent with migration and coalescence rates informed by the unequal recombination process of gene copies. With this novel interpretation of population structure models we determine coalescence times and expected genetic diversity in samples of orthologous and paralogous copies from a gene family. As a case study, we discuss the site frequency spectrum of a small gene family in the two scenarios of high and of no gene copy number variation among individuals. These examples underline the significance of our model, since standard test-statistics may lead to misinterpretations when analyzing sequence data of multi-copy genes due to their different expected genetic diversity.

Bourdieusian Boundary-Making, Social Networks, and Capital Conversion: Inequality among International Degree Holders in Hong Kong

Sociological research richly documents the many ways through which education becomes a form of convertible capital, but focuses less on the cultural schemas that graduates possess and use to respond to disruptions of capital conversion processes. Using the case of international degree holders in Hong Kong, this article draws upon Bourdieu’s theory of practice to interrogate the cultural schemas that valorize international degrees when their conversion pathways to economic capital are subjectively perceived to weaken. This article unearths the role of social networks in embedding cultural schemas and their effects on relations within the field: when faced with diminishing economic returns, international degree holders hold fast to their schemas vis-à-vis fellow international graduates and reconceptualize their degrees as symbolic capital to cope with the loss by enacting symbolic violence against domestic degree holders. Class boundaries are ultimately entrenched when international degree graduates valorize their cultural capital gains and legitimate their economic capital losses. Doing so compromises their class interests by forcing themselves into an interstitial position between different fields: though they occupy dominating homologous positions in the cultural field, they choose to overlook their dominated homologous positions in the economic field.

Modulation of plant acetyl-CoA synthetase activity by post-translational lysine acetylation.

Acetyl-CoA synthetase (ACS) is one of several enzymes that generate the key metabolic intermediate, acetyl-CoA. In microbes and mammals ACS activity is regulated by the post-translational acetylation of a key lysine residue. ACS in plant cells is part of a two-enzyme system that maintains acetate homeostasis, but its post-translational regulation is unknown. This study demonstrates that the plant ACS activity can be regulated by the acetylation of a specific lysine residue that is positioned in a homologous position as the microbial and mammalian ACS sequences that regulates ACS activity, occurring in the middle of a conserved motif, near the carboxyl-end of the protein. The inhibitory effect of the acetylation of residue Lys-622 of the Arabidopsis ACS was demonstrated by site-directed mutagenesis of this residue, including its genetic substitution with the non-canonical N-ε-acetyl-lysine residue. This latter modification lowered the catalytic efficiency of the enzyme by a factor of more than 500-fold. Michaelis-Menten kinetic analysis of the mutant enzyme indicates that this acetylation affects the first half-reaction of the ACS catalyzed reaction, namely, the formation of the acetyl adenylate enzyme intermediate. The post-translational acetylation of the plant ACS could affect acetate flux in the plastids and overall acetate homeostasis.

Ivacaftor-Mediated Potentiation of ABCB4 Missense Mutations Affecting Critical Motifs of the NBDs: Repositioning Perspectives for Hepatobiliary Diseases.

ABCB4 (ATP-binding cassette subfamily B member 4) is a hepatocanalicular floppase involved in biliary phosphatidylcholine (PC) secretion. Variations in the ABCB4 gene give rise to several biliary diseases, including progressive familial intrahepatic cholestasis type 3 (PFIC3), an autosomal recessive disease that can be lethal in the absence of liver transplantation. In this study, we investigated the effect and potential rescue of ten ABCB4 missense variations in NBD1:NBD2 homologous positions (Y403H/Y1043H, K435M/K1075M, E558K/E1200A, D564G/D1206G and H589Y/H1231Y) all localized at the conserved and functionally critical motifs of ABC transporters, six of which are mutated in patients. By combining structure analysis and in vitro studies, we found that all ten mutants were normally processed and localized at the canalicular membrane of HepG2 cells, but showed dramatically impaired PC transport activity that was significantly rescued by treatment with the clinically approved CFTR potentiator ivacaftor. Our results provide evidence that functional ABCB4 mutations are rescued by ivacaftor, paving the way for the repositioning of this potentiator for the treatment of selected patients with PFIC3 caused by mutations in the ATP-binding sites of ABCB4.

Multiplexed LNP-mRNA vaccination against pathogenic coronavirus species.

Although COVID-19 vaccines have been developed, multiple pathogenic coronavirus species exist, urging on development of multispecies coronavirus vaccines. Here we develop prototype lipid nanoparticle (LNP)-mRNA vaccine candidates against SARS-CoV-2 Delta, SARS-CoV, and MERS-CoV, and we test how multiplexing LNP-mRNAs can induce effective immune responses in animal models. Triplex and duplex LNP-mRNA vaccinations induce antigen-specific antibody responses against SARS-CoV-2, SARS-CoV, and MERS-CoV. Single-cell RNA sequencing profiles the global systemic immune repertoires and respective transcriptome signatures of vaccinated animals, revealing a systemic increase in activated B cells and differential gene expression across major adaptive immune cells. Sequential vaccination shows potent antibody responses against all three species, significantly stronger than simultaneous vaccination in mixture. These data demonstrate the feasibility, antibody responses, and single-cell immune profiles of multispecies coronavirus vaccination. The direct comparison between simultaneous and sequential vaccination offers insights into optimization of vaccination schedules to provide broad and potent antibody immunity against three major pathogenic coronavirus species.

Shared Transcriptional Machinery at Homologous Alleles Leads to Reduced Transcription in Early Drosophila Embryos.

The mechanism by which transcriptional machinery is recruited to enhancers and promoters to regulate gene expression is one of the most challenging and extensively studied questions in modern biology. We explored the possibility that interallelic interactions between two homologous alleles might affect gene regulation. Using an MS2- and PP7-based, allele-specific live imaging assay, we visualized de novo transcripts of a reporter gene in hemizygous and homozygous Drosophila embryos. Surprisingly, each homozygous allele produced fewer RNAs than the corresponding hemizygous allele, suggesting the possibility of allelic competition in homozygotes. However, the competition was not observed when the enhancer-promoter interaction was weakened by placing the reporter construct in a different chromosome location or by moving the enhancer further away from the promoter. Moreover, the reporter gene showed reduced transcriptional activity when a partial transcription unit (either an enhancer or reporter gene only) was in the homologous position. We propose that the transcriptional machinery that binds both the enhancer and promoter regions, such as RNA Pol II or preinitiation complexes, may be responsible for the allelic competition. We showed that the degree of allelic interference increased over developmental time as more Pol II was needed to activate zygotic genes. Such allelic competition was observed for an endogenous gene as well. Our study provides new insights into the role of 3D interallelic interactions in gene regulation.

Global, highly specific and fast filtering of alignment seeds

BackgroundAn important initial phase of arguably most homology search and alignment methods such as required for genome alignments is seed finding. The seed finding step is crucial to curb the runtime as potential alignments are restricted to and anchored at the sequence position pairs that constitute the seed. To identify seeds, it is good practice to use sets of spaced seed patterns, a method that locally compares two sequences and requires exact matches at certain positions only.ResultsWe introduce a new method for filtering alignment seeds that we call geometric hashing. Geometric hashing achieves a high specificity by combining non-local information from different seeds using a simple hash function that only requires a constant and small amount of additional time per spaced seed. Geometric hashing was tested on the task of finding homologous positions in the coding regions of human and mouse genome sequences. Thereby, the number of false positives was decreased about million-fold over sets of spaced seeds while maintaining a very high sensitivity.ConclusionsAn additional geometric hashing filtering phase could improve the run-time, accuracy or both of programs for various homology-search-and-align tasks.

Distinct functional roles for the M4 α-helix from each homologous subunit in the heteropentameric ligand-gated ion channel nAChR

The outermost lipid-exposed α-helix (M4) in each of the homologous α, β, δ, and γ/ε subunits of the muscle nicotinic acetylcholine receptor (nAChR) has previously been proposed to act as a lipid sensor. However, the mechanism by which this sensor would function is not clear. To explore how the M4 α-helix from each subunit in human adult muscle nAChR influences function, and thus explore its putative role in lipid sensing, we functionally characterized alanine mutations at every residue in αM4, βM4, δM4, and εM4, along with both alanine and deletion mutations in the post-M4 region of each subunit. Although no critical interactions involving residues on M4 or in post-M4 were identified, we found that numerous mutations at the M4–M1/M3 interface altered the agonist-induced response. In addition, homologous mutations in M4 in different subunits were found to have different effects on channel function. The functional effects of multiple mutations either along M4 in one subunit or at homologous positions of M4 in different subunits were also found to be additive. Finally, when characterized in both Xenopus oocytes and human embryonic kidney 293T cells, select αM4 mutations displayed cell-specific phenotypes, possibly because of the different membrane lipid environments. Collectively, our data suggest different functional roles for the M4 α-helix in each heteromeric nAChR subunit and predict that lipid sensing involving M4 occurs primarily through the cumulative interactions at the M4–M1/M3 interface, as opposed to the alteration of specific interactions that are critical to channel function.

Bi-alignments with affine gaps costs

BackgroundCommonly, sequence and structure elements are assumed to evolve congruently, such that homologous sequence positions correspond to homologous structural features. Assuming congruent evolution, alignments based on sequence and structure similarity can therefore optimize both similarities at the same time in a single alignment. To model incongruent evolution, where sequence and structural features diverge positionally, we recently introduced bi-alignments. This generalization of sequence and structure-based alignments is best understood as alignments of two distinct pairwise alignments of the same entities: one modeling sequence similarity, the other structural similarity.ResultsOptimal bi-alignments with affine gap costs (or affine shift cost) for two constituent alignments can be computed exactly in quartic space and time. Even bi-alignments with affine shift and gap cost, as well as bi-alignment with sub-additive gap cost are optimized efficiently. Affine gap-cost bi-alignment of large proteins (sim 930 aa) can be computed.ConclusionAffine cost bi-alignments are of practical interest to study shifts of protein sequences and protein structures relative to each other.AvailabilityThe affine cost bi-alignment algorithm has been implemented in Python 3 and Cython. It is available as free software from https://github.com/s-will/BiAlign/releases/tag/v0.3 and as bioconda package bialign.

A Myosin Vb Splice Variant Regulates Coronavirus M Protein Trafficking in Polarized Epithelial Cells

The membrane (M) glycoprotein of coronaviruses (CoVs) serves as the nidus for virion assembly. Using a yeast two‐hybrid screen, we identified the interaction of the cytosolic tail of Murine Hepatitis Virus (MHV‐CoV) M protein with Myosin Vb (MYO5B) (Figure 1A). MHV‐CoV M interacts specifically with the alternative splice variant of cellular MYO5B including Exon D (MYO5B+D), that also mediates interaction with cellular Rab10. When co‐expressed in human lung epithelial A549 and canine kidney epithelial MDCK cells, MYO5B+D co‐localized with MHV‐CoV M protein, as well as with M proteins from porcine epidemic diarrhea Virus (PEDV‐CoV), Middle East Respiratory Syndrome (MERS‐CoV) and Severe Acute Respiratory Syndrome 2 (SARS‐CoV‐2) (Figure 2). M‐GFP chimeric proteins co‐expressed with mCherry‐MYO5B+D also co‐localized with endogenous Rab10 and Rab11a (Figure 1B,C). We identified point mutations in MHV‐CoV M that blocked the interaction with MYO5B+D in yeast 2‐hybrid assays (Figure 1B). One of these point mutations (E121K) was previously shown to block MHV‐CoV virion assembly, and it blocked interaction with MYO5B+D. The E to K mutation at homologous positions in PEDV‐CoV, MERS‐CoV and SARS‐CoV‐2 M proteins also blocked colocalization with MYO5B+D. Knockdown of Rab10 blocked the co‐localization of the M proteins with MYO5B+D (Figure 2). Re‐expression of Cerulein‐Rab10 in Rab10 KD cells re‐established colocalization between M proteins and MYO5B+D (Figure 2). Our results suggest that the interaction of CoV M proteins with MYO5B+D may play a role in regulating their trafficking through Rab10‐containing membrane systems in epithelial cells.

A Collaborative Classroom Investigation of the Evolution of SABATH Methyltransferase Substrate Preference Shifts over 120 My of Flowering Plant History.

Next-generation sequencing has resulted in an explosion of available data, much of which remains unstudied in terms of biochemical function; yet, experimental characterization of these sequences has the potential to provide unprecedented insight into the evolution of enzyme activity. One way to make inroads into the experimental study of the voluminous data available is to engage students by integrating teaching and research in a college classroom such that eventually hundreds or thousands of enzymes may be characterized. In this study, we capitalize on this potential to focus on SABATH methyltransferase enzymes that have been shown to methylate the important plant hormone, salicylic acid (SA), to form methyl salicylate. We analyze data from 76 enzymes of flowering plant species in 23 orders and 41 families to investigate how widely conserved substrate preference is for SA methyltransferase orthologs. We find a high degree of conservation of substrate preference for SA over the structurally similar metabolite, benzoic acid, with recent switches that appear to be associated with gene duplication and at least three cases of functional compensation by paralogous enzymes. The presence of Met in active site position 150 is a useful predictor of SA methylation preference in SABATH methyltransferases but enzymes with other residues in the homologous position show the same substrate preference. Although our dense and systematic sampling of SABATH enzymes across angiosperms has revealed novel insights, this is merely the “tip of the iceberg” since thousands of sequences remain uncharacterized in this enzyme family alone.

Introgression of clubroot resistant gene into Brassica oleracea L. from Brassica rapa based on homoeologous exchange.

Clubroot is a soil-borne disease in cabbage (Brassica oleracea L. var. capitata L.) caused by Plasmodiophora brassicae, which poses a great threat to cabbage production. However, clubroot resistance (CR) genes in Brassica rapa could be introduced into the cabbage via breeding to make it clubroot resistant. In this study, CR genes from B. rapa were introduced into the cabbage genome and the mechanism of gene introgression was explored. Two methods were used to create CR materials: (i) The fertility of CR Ogura CMS cabbage germplasms containing CRa was restored by using an Ogura CMS restorer. After cytoplasmic replacement and microspore culture, CRa-positive microspore individuals were obtained. (ii) Distant hybridization was performed between cabbage and B. rapa, which contained three CR genes (CRa, CRb, and Pb8.1). Finally, BC2 individuals containing all three CR genes were obtained. Inoculation results showed that both CRa-positive microspore individuals and BC2 individuals containing three CR genes were resistant to race 4 of P. brassicae. Sequencing results from CRa-positive microspore individuals with specific molecular markers and genome-wide association study (GWAS) showed penetration at the homologous position of the cabbage genome by a 3.42Mb CRa containing a fragment from B. rapa; indicating homoeologous exchange (HE) as the theoretical basis for the introgression of CR resistance. The successful introduction of CR into the cabbage genome in the present study can provide useful clues for creating introgression lines within other species of interest.