Incautious design of shRNAs for stable overexpression of miRNAs could result in generation of undesired isomiRs.

shRNA-mediated strategy of miRNA overexpression based on RNA Polymerase III (Pol III) expression cassettes is widely used for miRNA functional studies. For some miRNAs, e.g., encoded in the genome as a part of a polycistronic miRNA cluster, it is most likely the only way for their individual stable overexpression. Here we have revealed that expression of miRNAs longer than 19 nt (e.g. 23 nt in length hsa-miR-93-5p) using such approach could be accompanied by undesired predominant generation of 5' end miRNA isoforms (5'-isomiRs). Extra U residues (up to five) added by Pol III at the 3' end of the transcribed shRNA during transcription termination could cause a shift in the Dicer cleavage position of the shRNA. This results in the formation of 5'-isomiRs, which have a significantly altered seed region compared to the initially encoded canonical hsa-miR-93-5p. We demonstrated that the commonly used qPCR method is insensitive to the formation of 5'-isomiRs and cannot be used to confirm miRNA overexpression. However, the predominant expression of 5'-isomiRs without three or four first nucleotides instead of the canonical isoform could be disclosed based on miRNA-Seq analysis. Moreover, mRNA sequencing data showed that the 5'-isomiRs of hsa-miR-93-5p presumably regulate their own mRNA targets. Thus, omitting miRNA-Seq analysis may lead to erroneous conclusions regarding revealed mRNA targets and possible molecular mechanisms in which studied miRNA is involved. Overall, the presented results show that structures of shRNAs for stable overexpression of miRNAs requires careful design to avoid generation of undesired 5'-isomiRs.

IGFBP6 regulates extracellular vesicles formation via cholesterol abundance in MDA-MB-231 cells.

Breast cancer recurrence is associated with the growth of disseminated cancer cells that separate from the primary tumor before surgical treatment and hormonal therapy and form a metastatic niche in distant organs. We previously demonstrated that IGFBP6 expression is associated with the risk of early relapse of luminal breast cancer. Knockdown of IGFBP6 in MDA-MB-231 breast cancer cells increased their invasiveness, proliferation, and metastatic potential. In addition, the knockdown of IGFBP6 leads to impaired lipid metabolism. In this study, we demonstrated that the knockdown of the IGFBP6 gene, a highly selective inhibitor of IGF-II, led to a significant decline in the number of secreted extracellular vesicles (EVs) and altered cholesterol metabolism in MDA-MB-231 cells. Knockdown of IGFBP6 led to a decrease in the essential proteins responsible for the biogenesis of cholesterol LDLR and LSS, which reduced the amount by more than 13 times. In addition, the knockdown of IGFBP6 led to a possible change in the profile of adhesion molecules on the surface of EVs. The expression of L1CAM, IGSF3, EpCAM, CD24, and CD44 decreased, and the expression of EGFR increased. We can conclude that the negative prognostic value of low expression of this gene could be associated with increased activity of IGF2 in tumor-associated fibroblasts due to low secretion of IGFBP6 by tumor cells. In addition, changing the profile of adhesion molecules on the surface of tumor EVs may contribute to the more efficient formation of metastatic niches.

α-Latrotoxin Tetramers Spontaneously Form Two-Dimensional Crystals in Solution and Coordinated Multi-Pore Assemblies in Biological Membranes.

α-Latrotoxin (α-LTX) was found to form two-dimensional (2D) monolayer arrays in solution at relatively low concentrations (0.1 mg/mL), with the toxin tetramer constituting a unit cell. The crystals were imaged using cryogenic electron microscopy (cryoEM), and image analysis yielded a ~12 Å projection map. At this resolution, no major conformational changes between the crystalline and solution states of α-LTX tetramers were observed. Electrophysiological studies showed that, under the conditions of crystallization, α-LTX simultaneously formed multiple channels in biological membranes that displayed coordinated gating. Two types of channels with conductance levels of 120 and 208 pS were identified. Furthermore, we observed two distinct tetramer conformations of tetramers both when observed as monodisperse single particles and within the 2D crystals, with pore diameters of 11 and 13.5 Å, suggestive of a flickering pore in the middle of the tetramer, which may correspond to the two states of toxin channels with different conductance levels. We discuss the structural changes that occur in α-LTX tetramers in solution and propose a mechanism of α-LTX insertion into the membrane. The propensity of α-LTX tetramers to form 2D crystals may explain many features of α-LTX toxicology and suggest that other pore-forming toxins may also form arrays of channels to exert maximal toxic effect.

Mathematical model explains differences in Omicron and Delta SARS-CoV-2 dynamics in Caco-2 and Calu-3 cells.

Within-host infection dynamics of Omicron dramatically differs from previous variants of SARS-CoV-2. However, little is still known about which parameters of virus-cell interplay contribute to the observed attenuated replication and pathogenicity of Omicron. Mathematical models, often expressed as systems of differential equations, are frequently employed to study the infection dynamics of various viruses. Adopting such models for results of in vitro experiments can be beneficial in a number of aspects, such as model simplification (e.g., the absence of adaptive immune response and innate immunity cells), better measurement accuracy, and the possibility to measure additional data types in comparison with in vivo case. In this study, we consider a refinement of our previously developed and validated model based on a system of integro-differential equations. We fit the model to the experimental data of Omicron and Delta infections in Caco-2 (human intestinal epithelium model) and Calu-3 (lung epithelium model) cell lines. The data include known information on initial conditions, infectious virus titers, and intracellular viral RNA measurements at several time points post-infection. The model accurately explains the experimental data for both variants in both cell lines using only three variant- and cell-line-specific parameters. Namely, the cell entry rate is significantly lower for Omicron, and Omicron triggers a stronger cytokine production rate (i.e., innate immune response) in infected cells, ultimately making uninfected cells resistant to the virus. Notably, differences in only a single parameter (e.g., cell entry rate) are insufficient to obtain a reliable model fit for the experimental data.

Viral and host small RNA transcriptome analysis of SARS-CoV-1 and SARS-CoV-2-infected human cells reveals novel viral short RNAs.

RNA viruses have been shown to express various short RNAs, some of which have regulatory roles during replication, transcription, and translation of viral genomes. However, short viral RNAs generated from SARS-CoV-1 and SARS-CoV-2 genomic RNAs remained largely unexplored, possibly due limitations of the widely used library preparation methods for small RNA deep sequencing and corresponding data processing. By analyzing publicly available small RNA sequencing datasets, we observed that human Calu-3 cells infected by SARS-CoV-1 or SARS-CoV-2 accumulate multiple previously unreported short viral RNAs. In addition, we verified the presence of the five most abundant SARS-CoV-2 short viral RNAs in SARS-CoV-2-infected human lung adenocarcinoma cells by quantitative PCR. Interestingly, the copy number of the observed SARS-CoV-2 short viral RNAs dramatically exceeded the expression of previously reported viral microRNAs in the same cells. We hypothesize that the reported SARS-CoV-2 short viral RNAs could serve as biomarkers for early infection stages due to their high abundance. Furthermore, unlike SARS-CoV-1, the SARS-CoV-2 infection induced significant (Benjamini-Hochberg-corrected p-value <0.05) deregulation of Y-RNA, transfer RNA, vault RNA, as well as more than 300 endogenous short RNAs that aligned predominantly to human protein-coding and long noncoding RNA transcripts. In particular, more than 20-fold upregulation of reads derived from Y-RNA (and several transfer RNAs) have been documented in RNA-seq datasets from SARS-CoV-2 infected cells. Finally, a significant proportion of short RNAs derived from full-length viral genomes also aligned to various human genome (hg38) sequences, suggesting opportunities to investigate regulatory roles of short viral RNAs during infection. Further characterization of the small RNA landscape of both viral and host genomes is clearly warranted to improve our understanding of molecular events related to infection and to design more efficient strategies for therapeutic interventions as well as early diagnosis.

Spontaneous metastasis xenograft models link CD44 isoform 4 to angiogenesis, hypoxia, EMT and mitochondria-related pathways in colorectal cancer.

Hematogenous metastasis limits the survival of colorectal cancer (CRC) patients. Here, we illuminated the roles of CD44 isoforms in this process. Isoforms 3 and 4 were predominantly expressed in CRC patients. CD44 isoform 4 indicated poor outcome and correlated with epithelial-mesenchymal transition (EMT) and decreased oxidative phosphorylation (OxPhos) in patients; opposite associations were found for isoform 3. Pan-CD44 knockdown (kd) independently impaired primary tumor formation and abrogated distant metastasis in CRC xenografts. The xenograft tumors mainly expressed the clinically relevant CD44 isoforms 3 and 4. Both isoforms were enhanced in the paranecrotic, hypoxic tumor regions but were generally absent in lung metastases. Upon CD44 kd, tumor angiogenesis was increased in the paranecrotic areas, accompanied by reduced hypoxia-inducible factor-1α and CEACAM5 but increased E-cadherin expression. Mitochondrial genes and proteins were induced upon pan-CD44 kd, as were OxPhos genes. Hypoxia increased VEGF release from tumor spheres, particularly upon CD44 kd. Genes affected upon CD44 kd in xenografts specifically overlapped concordantly with genes correlating with CD44 isoform 4 (but not isoform 3) in patients, validating the clinical relevance of the used model and highlighting the metastasis-promoting role of CD44 isoform 4.

Challenges in characterization of transcriptomes of extracellular vesicles and non-vesicular extracellular RNA carriers.

Since its original discovery over a decade ago, extracellular RNA (exRNA) has been found in all biological fluids. Furthermore, extracellular microRNA has been shown to be involved in communication between various cell types. Importantly, the exRNA is protected from RNases degradation by certain carriers including membrane vesicles and non-vesicular protein nanoparticles. Each type of carrier has its unique exRNA profile, which may vary depending on cell type and physiological conditions. To clarify putative mechanisms of intercellular communication mediated by exRNA, the RNA profile of each carrier has to be characterized. While current methods of biofluids fractionation are continuously improving, they fail to completely separate exRNA carriers. Likewise, most popular library preparation approaches for RNA sequencing do not allow obtaining exhaustive and unbiased data on exRNA transcriptome. In this mini review we discuss ongoing progress in the field of exRNA, with the focus on exRNA carriers, analyze the key methodological challenges and provide recommendations on how the latter could be overcome.

RNA-binding proteins regulating the CD44 alternative splicing.

Alternative splicing is often deregulated in cancer, and cancer-specific isoform switches are part of the oncogenic transformation of cells. Accumulating evidence indicates that isoforms of the multifunctional cell-surface glycoprotein CD44 play different roles in cancer cells as compared to normal cells. In particular, the shift of CD44 isoforms is required for epithelial to mesenchymal transition (EMT) and is crucial for the maintenance of pluripotency in normal human cells and the acquisition of cancer stem cells phenotype for malignant cells. The growing and seemingly promising use of splicing inhibitors for treating cancer and other pathologies gives hope for the prospect of using such an approach to regulate CD44 alternative splicing. This review integrates current knowledge about regulating CD44 alternative splicing by RNA-binding proteins.

Metabolic effects of vasopressin in pathophysiology of diabetic kidney disease.

The diabetic kidney disease (DKD) is the major cause of the chronic kidney disease (CKD). Enhanced plasma vasopressin (VP) levels have been associated with the pathophysiology of DKD and CKD. Stimulation of VP release in DKD is caused by glucose-dependent reset of the osmostat leading to secondary pathophysiologic effects mediated by distinct VP receptor types. VP is a stress hormone exhibiting the antidiuretic action in the kidney along with broad adaptive effects in other organs. Excessive activation of the vasopressin type 2 (V2) receptor in the kidney leads to glomerular hyperfiltration and nephron loss, whereas stimulation of vasopressin V1a or V1b receptors in the liver, pancreas, and adrenal glands promotes catabolic metabolism for energy mobilization, enhancing glucose production and aggravating DKD. Increasing availability of selective VP receptor antagonists opens new therapeutic windows separating the renal and extra-renal VP effects for the concrete applications. Improved understanding of these paradigms is mandatory for further drug design and translational implementation. The present concise review focuses on metabolic effects of VP affecting DKD pathophysiology.

The signature of SARS-CoV-2 evolution reflects selective pressures within human guts.

In somatic cells, microRNAs (miRNAs) bind to the genomes of RNA viruses and influence their translation and replication. In London and Berlin samples represented in GISAID database, we traced severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages and divided these sequenced in two groups, "Ancestral variants" and "Omicrons," and analyzed them through the prism of the tissue-specific binding between host miRNAs and viral messenger RNAs. We demonstrate a significant number of miRNA-binding sites in the NSP4 region of the SARS-CoV-2 genome, with evidence of evolutionary pressure within this region exerted by human intestinal miRNAs. Notably, in infected cells, NSP4 promotes the formation of double-membrane vesicles, which serve as the scaffolds for replication-transcriptional complexes and protect viral RNA from intracellular destruction. In 3 years of selection, the loss of many miRNA-binding sites in general and those within the NSP4 in particular has shaped the SARS-CoV-2 genomes. With that, the descendants of the BA.2 variants were promoted as dominant strains, which define current momentum of the pandemics.

Immunomodulatory Effect and Biological Significance of ß-Glucans.

ß-glucan, one of the homopolysaccharides composed of D-glucose, exists widely in cereals and microorganisms and possesses various biological activities, including anti-inflammatory, antioxidant, and anti-tumor properties. More recently, there has been mounting proof that ß-glucan functions as a physiologically active "biological response modulator (BRM)", promoting dendritic cell maturation, cytokine secretion, and regulating adaptive immune responses-all of which are directly connected with ß-glucan-regulated glucan receptors. This review focuses on the sources, structures, immune regulation, and receptor recognition mechanisms of ß-glucan.

Development of a novel mathematical model that explains SARS-CoV-2 infection dynamics in Caco-2 cells.

Mathematical modeling is widely used to study within-host viral dynamics. However, to the best of our knowledge, for the case of SARS-CoV-2 such analyses were mainly conducted with the use of viral load data and for the wild type (WT) variant of the virus. In addition, only few studies analyzed models for in vitro data, which are less noisy and more reproducible. In this work we collected multiple data types for SARS-CoV-2-infected Caco-2 cell lines, including infectious virus titers, measurements of intracellular viral RNA, cell viability data and percentage of infected cells for the WT and Delta variants. We showed that standard models cannot explain some key observations given the absence of cytopathic effect in human cell lines. We propose a novel mathematical model for in vitro SARS-CoV-2 dynamics, which included explicit modeling of intracellular events such as exhaustion of cellular resources required for virus production. The model also explicitly considers innate immune response. The proposed model accurately explained experimental data. Attenuated replication of the Delta variant in Caco-2 cells could be explained by our model on the basis of just two parameters: decreased cell entry rate and increased cytokine production rate.

ELOVL5 and IGFBP6 genes modulate sensitivity of breast cancer cells to ferroptosis.

Introduction: Relapse of breast cancer is one of the key obstacles to successful treatment. Previously we have shown that low expression of ELOVL5 and IGFBP6 genes in breast cancer tissue corresponded to poor prognosis. ELOVL5 participates directly in the elongation of polyunsaturated fatty acids (PUFAs) that are considered to play an important role in cancer cell metabolism. Thus, in this work we studied the changes in lipid metabolism in breast cancer cells with reduced expression of either ELOVL5 or IGFBP6 gene. Methods: MDA-MB-231 cells with a stable knockdown of either ELOVL5 or IGFBP6 gene were used in this study. Transcriptomic and proteomic analysis as well as RT-PCR were utilized to assess gene expression. Content of individual fatty acids in the cells was measured with HPLC-MS. HPLC was used for analysis of the kinetics of PUFAs uptake. Cell viability was measured with MTS assay. Flow cytometry was used to measure activation of apoptosis. Fluorescent microscopy was utilized to assess accumulation of ROS and formation of lipid droplets. Glutathione peroxidase activity was measured with a colorimetric assay. Results: We found that the knockdown of IGFBP6 gene led to significant changes in the profile of fatty acids in the cells and in the expression of many genes associated with lipid metabolism. As some PUFAs are known to inhibit proliferation and cause death of cancer cells, we also tested the response of the cells to single PUFAs and to combinations of docosahexaenoic acid (DHA, a n-3 PUFA) with standard chemotherapeutic drugs. Our data suggest that external PUFAs cause cell death by activation of ferroptosis, an iron-dependent mechanism of cell death with excessive lipid peroxidation. Moreover, both knockdowns increased cells' sensitivity to ferroptosis, probably due to a significant decrease in the activity of the antioxidant enzyme GPX4. Addition of DHA to commonly used chemotherapeutic drugs enhanced their effect significantly, especially for the cells with low expression of IGFBP6 gene. Discussion: The results of this study suggest that addition of PUFAs to the treatment regimen for the patients with low expression of IGFBP6 and ELOVL5 genes can be potentially beneficial and is worth testing in a clinically relevant setting.

HLA-A*01:01 allele diminishing in COVID-19 patients population associated with non-structural epitope abundance in CD8+ T-cell repertoire.

In mid-2021, the SARS-CoV-2 Delta variant caused the third wave of the COVID-19 pandemic in several countries worldwide. The pivotal studies were aimed at studying changes in the efficiency of neutralizing antibodies to the spike protein. However, much less attention was paid to the T-cell response and the presentation of virus peptides by MHC-I molecules. In this study, we compared the features of the HLA-I genotype in symptomatic patients with COVID-19 in the first and third waves of the pandemic. As a result, we could identify the diminishing of carriers of the HLA-A*01:01 allele in the third wave and demonstrate the unique properties of this allele. Thus, HLA-A*01:01-binding immunoprevalent epitopes are mostly derived from ORF1ab. A set of epitopes from ORF1ab was tested, and their high immunogenicity was confirmed. Moreover, analysis of the results of single-cell phenotyping of T-cells in recovered patients showed that the predominant phenotype in HLA-A*01:01 carriers is central memory T-cells. The predominance of T-lymphocytes of this phenotype may contribute to forming long-term T-cell immunity in carriers of this allele. Our results can be the basis for highly effective vaccines based on ORF1ab peptides.

Differential co-expression network analysis with DCoNA reveals isomiR targeting aberrations in prostate cancer.

MOTIVATION: One of the standard methods of high-throughput RNA sequencing analysis is differential expression. However, it does not detect changes in molecular regulation. In contrast to the standard differential expression analysis, differential co-expression one aims to detect pairs or clusters whose mutual expression changes between two conditions. RESULTS: We developed Differential Co-expression Network Analysis (DCoNA)-an open-source statistical tool that allows one to identify pair interactions, which correlation significantly changes between two conditions. Comparing DCoNA with the state-of-the-art analog, we showed that DCoNA is a faster, more accurate and less memory-consuming tool. We applied DCoNA to prostate mRNA/miRNA-seq data collected from The Cancer Genome Atlas (TCGA) and compared predicted regulatory interactions of miRNA isoforms (isomiRs) and their target mRNAs between normal and cancer samples. As a result, almost all highly expressed isomiRs lost negative correlation with their targets in prostate cancer samples compared to ones without the pathology. One exception to this trend was the canonical isomiR of hsa-miR-93-5p acquiring cancer-specific targets. Further analysis showed that cancer aggressiveness simultaneously increased with the expression level of this isomiR in both TCGA primary tumor samples and 153 blood plasma samples of P. Hertsen Moscow Oncology Research Institute patients' cohort analyzed by miRNA microarrays. AVAILABILITY AND IMPLEMENTATION: Source code and documentation of DCoNA are available at https://github.com/zhiyanov/DCoNA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A novel approach for a joint analysis of isomiR and mRNA expression data reveals features of isomiR targeting in breast cancer.

A widely used procedure for selecting significant miRNA-mRNA or isomiR-mRNA pairs out of predicted interactions involves calculating the correlation between expression levels of miRNAs/isomiRs and mRNAs in a series of samples. In this manuscript, we aimed to assess the validity of this procedure by comparing isomiR-mRNA correlation profiles in sets of sequence-based predicted target mRNAs and non-target mRNAs (negative controls). Target prediction was carried out using RNA22 and TargetScan algorithms. Spearman's correlation analysis was conducted using miRNA and mRNA sequencing data of The Cancer Genome Atlas Breast Invasive Carcinoma (TCGA-BRCA) project. Luminal A, luminal B, basal-like breast cancer subtypes, and adjacent normal tissue samples were analyzed separately. Using the sets of putative targets and non-targets, we introduced adjusted isomiR targeting activity (ITA)-the number of negatively correlated potential isomiR targets adjusted by the background (estimated using non-target mRNAs). We found that for most isomiRs a significant negative correlation between isomiR-mRNA expression levels appeared more often in a set of predicted targets compared to the non-targets. This trend was detected for both classical seed region binding types (8mer, 7mer-m8, 7mer-A1, 6mer) predicted by TargetScan and the non-classical ones (G:U wobbles and up to one mismatch or unpaired nucleotide within seed sequence) predicted by RNA22. Adjusted ITA distributions were similar for target sites located in 3'-UTRs and coding mRNA sequences, while 5'-UTRs had much lower scores. Finally, we observed strong cancer subtype-specific patterns of isomiR activity, highlighting the differences between breast cancer molecular subtypes and normal tissues. Surprisingly, our target prediction- and correlation-based estimates of isomiR activities were practically non-correlated with the average isomiR expression levels neither in cancerous nor in normal samples.

On the peptide binding affinity changes in population-specific HLA repertoires to the SARS-CoV-2 variants Delta and Omicron.

OBJECTIVE: To investigate the changes of Spike protein-HLA binding affinity profiles between the Wuhan strain and two dominant variants, the Delta and the Omicron strains, among the Taiwanese, the British and the Russian populations. METHODS: The HLA frequencies and the HLA-peptide binding affinity profiles in the T-CoV database were combined to conduct the study. We focused on the public alleles in the three populations (HLA-A, HLA-B, HLA-C, HLA-DRB1, and/or HLA-DPA1/DPB1 alleles) and the altered peptides of the spike protein (compared to the Wuhan strain) in the Delta G/478K·V1 (B.1.617.2 + AY.1 + AY.2) and the Omicron (BA.1) strains. RESULTS: For the Delta strain, tight bindings of the altered peptides to the HLA alleles decrease in all three populations and almost vanish in the Taiwanese population. For the Omicron strain, tight bindings are mostly preserved for both HLA classes and in the Taiwanese and the British populations, with a slight reduction in HLA class II in the Taiwanese (1.4%), while the Russian population preserves a relatively high fraction of tight bindings for both HLA classes. CONCLUSION: We comprehensively reported the changes in the HLA-associated SARS-CoV-2 Spike protein peptide binding profiles among the Taiwanese, the British, and the Russian populations. Further studies are needed to understand the immunological mechanisms and the clinical value of our findings.

isomiRTar: a comprehensive portal of pan-cancer 5'-isomiR targeting.

Inaccurate cleavage of pri- and pre-miRNA hairpins by Drosha and Dicer results in the generation of miRNA isoforms known as isomiRs. isomiRs with 5'-end variations (5'-isomiRs) create a new dimension in miRNA research since they have different seed regions and distinct targetomes. We developed isomiRTar (https://isomirtar.hse.ru)-a comprehensive portal that allows one to analyze expression profiles and targeting activity of 5'-isomiRs in cancer. Using the Cancer Genome Atlas sequencing data, we compiled the list of 1022 5'-isomiRs expressed in 9282 tumor samples across 31 cancer types. Sequences of these isomiRs were used to predict target genes with miRDB and TargetScan. The putative interactions were then subjected to the co-expression analysis in each cancer type to identify isomiR-target pairs supported by significant negative correlations. Downstream analysis of the data deposited in isomiRTar revealed both cancer-specific and cancer-conserved 5'-isomiR expression landscapes. Pairs of isomiRs differing in one nucleotide shift from 5'-end had poorly overlapping targetomes with the median Jaccard index of 0.06. The analysis of colorectal cancer 5'-isomiR-mediated regulatory networks revealed promising candidate tumor suppressor isomiRs: hsa-miR-203a-3p-+1, hsa-miR-192-5p-+1 and hsa-miR-148a-3p-0. In summary, we believe that isomiRTar will help researchers find novel mechanisms of isomiR-mediated gene silencing in different types of cancer.

Specificity of viscumin revised. As probed with a printed glycan array.

Viscumin, a lectin used in anti-cancer therapy, was originally considered as ßGal recognizing protein; later, an ability to bind 6'-sialyl N-acetyllactosamine (6'SLN) terminated gangliosides was found. Here we probed viscumin with a printed glycan array (PGA) containing a large number of mammalian sulfated glycans, and found a strong binding to glycans with 6-O-SuGal moiety as lactose, N-acetyllactosamine (LN), di-N-acetyllactosamine (LacdiNAc), and even 6-O-SuGalNAcα (but not SiaTn). Also, the ability to bind some of αGal terminated glycans, including Galα1-3Galß1-4GlcNAc, was observed. Unexpectedly, only weak interaction was detected with parent neutral ß-galactosides including LN-LN-LN and branched (LN)2LN oligolactosamines; in the light of these data, one should not confidently classify viscumin as a ß-galactoside-binding lectin. Carrying out PGA in the presence of neutral or sulfated/sialylated glycan, together with sequential elution from lactose-sepharose and consideration of the protein structure, lead to the conclusion that two glycan-binding sites of viscumin have different specificities, one of which prefers charged sulfated and sialylated moieties.

ADGRL1 haploinsufficiency causes a variable spectrum of neurodevelopmental disorders in humans and alters synaptic activity and behavior in a mouse model.

ADGRL1 (latrophilin 1), a well-characterized adhesion G protein-coupled receptor, has been implicated in synaptic development, maturation, and activity. However, the role of ADGRL1 in human disease has been elusive. Here, we describe ten individuals with variable neurodevelopmental features including developmental delay, intellectual disability, attention deficit hyperactivity and autism spectrum disorders, and epilepsy, all heterozygous for variants in ADGRL1. In vitro, human ADGRL1 variants expressed in neuroblastoma cells showed faulty ligand-induced regulation of intracellular Ca2+ influx, consistent with haploinsufficiency. In vivo, Adgrl1 was knocked out in mice and studied on two genetic backgrounds. On a non-permissive background, mice carrying a heterozygous Adgrl1 null allele exhibited neurological and developmental abnormalities, while homozygous mice were non-viable. On a permissive background, knockout animals were also born at sub-Mendelian ratios, but many Adgrl1 null mice survived gestation and reached adulthood. Adgrl1-/- mice demonstrated stereotypic behaviors, sexual dysfunction, bimodal extremes of locomotion, augmented startle reflex, and attenuated pre-pulse inhibition, which responded to risperidone. Ex vivo synaptic preparations displayed increased spontaneous exocytosis of dopamine, acetylcholine, and glutamate, but Adgrl1-/- neurons formed synapses in vitro poorly. Overall, our findings demonstrate that ADGRL1 haploinsufficiency leads to consistent developmental, neurological, and behavioral abnormalities in mice and humans.