Functional and molecular features of the calmodulin-interacting protein IQCG required for haematopoiesis in zebrafish

Genome-wide association study reveals the underlying regulatory mechanisms of red blood traits in Anadara granosa

BACKGROUND

Anadara granosa, commonly known as the blood clam, exhibits the unusual characteristic of having red blood among invertebrates. There is significant individual variation in blood color intensity among blood clams; individuals with vibrant red blood are deemed healthier and exhibit stronger stress resistance. However, the molecular basis underlying these red blood traits (RBTs) remains poorly understood.

RESULTS

In this study, we performed genome-wide association studies (GWAS) in a population of 300 A. granosa individuals, focusing on RBTs as measured by hemoglobin concentration (HC), total hemocyte count (THC), and heme concentration (HEME). Our analysis identified 18 single nucleotide polymorphisms (SNPs) correlated with RBTs, subsequently selected 117 candidate genes within a 100 kb flanking region of these SNPs, potentially involved in the RBTs of A. granosa. Moreover, we discovered two haplotype blocks specifically associated with THC and HEME. Further analysis revealed eight genes (Septin7, Hox5, Cbfa2t3, Avpr1b, Hhex, Eif2ak3, Glrk, and Rpl35a) that significantly influence RBTs. Notably, a heterozygous A/T mutation in the 3'UTR of Cbfa2t3 was found to promote blood cell proliferation. These genes suggest that the hematopoietic function plays a significant role in the variability of RBTs in A. granosa.

CONCLUSIONS

Our findings reveal a conservation of the regulatory mechanisms of RBTs between blood clams and vertebrates. The results not only provide a scientific basis for selective breeding in blood clams, but also offer deeper insights into the evolutionary mechanisms of RBTs in invertebrates.

Deficiency of IQCH causes male infertility in humans and mice

IQ motif-containing proteins can be recognized by calmodulin (CaM) and are essential for many biological processes. However, the role of IQ motif-containing proteins in spermatogenesis is largely unknown. In this study, we identified a loss-of-function mutation in the novel gene IQ motif-containing H (IQCH) in a Chinese family with male infertility characterized by a cracked flagellar axoneme and abnormal mitochondrial structure. To verify the function of IQCH, Iqch knockout (KO) mice were generated via CRISPR-Cas9 technology. As expected, the Iqch KO male mice exhibited impaired fertility, which was related to deficient acrosome activity and abnormal structures of the axoneme and mitochondria, mirroring the patient phenotypes. Mechanistically, IQCH can bind to CaM and subsequently regulate the expression of RNA-binding proteins (especially HNRPAB), which are indispensable for spermatogenesis. Overall, this study revealed the function of IQCH, expanded the role of IQ motif-containing proteins in reproductive processes, and provided important guidance for genetic counseling and genetic diagnosis of male infertility.

Identification of IQCH as a calmodulin-associated protein required for sperm motility in humans

Sperm fertilization ability mainly relies on proper sperm progression through the female genital tract and capacitation, which involves phosphorylation signaling pathways triggered by calcium and bicarbonate. We performed exome sequencing of an infertile asthenozoospermic patient and identified truncating variants in MAP7D3, encoding a microtubule-associated protein, and IQCH, encoding a protein of unknown function with enzymatic and signaling features. We demonstrate the deleterious impact of both variants on sperm transcripts and proteins from the patient. We show that, in vitro, patient spermatozoa could not induce the phosphorylation cascades associated with capacitation. We also provide evidence for IQCH association with calmodulin, a well-established calcium-binding protein that regulates the calmodulin kinase. Notably, we describe IQCH spatial distribution around the sperm axoneme, supporting its function within flagella. Overall, our work highlights the cumulative pathological impact of gene mutations and identifies IQCH as a key protein required for sperm motility and capacitation.

IQUB deficiency causes male infertility by affecting the activity of p-ERK1/2/RSPH3

STUDY QUESTION

Can new genetic factors responsible for male infertility be identified, especially for those characterized by asthenospermia despite normal sperm morphology?

SUMMARY ANSWER

We identified the novel pathogenetic gene IQ motif and ubiquitin-like domain-containing (IQUB) as responsible for male infertility characterized by asthenospermia, involving sperm radial spoke defects.

WHAT IS KNOWN ALREADY

To date, only a few genes have been found to be responsible for asthenospermia with normal sperm morphology. Iqub, encoding the IQUB protein, is highly and specifically expressed in murine testes and interacts with the proteins radial spoke head 3 (RSPH3), CEP295 N-terminal like (CEP295NL or DDC8), glutathione S-transferase mu 1 (GSTM1) and outer dense fiber of sperm tails 1 (ODF1) in the yeast two-hybrid system.

STUDY DESIGN, SIZE, DURATION

The IQUB variant was identified by whole-exome sequencing in a cohort of 126 male infertility patients with typical asthenospermia recruited between 2015 and 2020. Knockout (KO) and knockin (KI) mouse models, scanning and transmission electron microscopy (TEM), and other functional assays were performed, between 2019 and 2021.

PARTICIPANTS/MATERIALS, SETTING, METHODS

The IQUB variant was identified by whole-exome sequencing and confirmed by Sanger sequencing. Iqub KO and KI mice were constructed to mimic the phenotype of the affected individual. After recapitulating the phenotype of human male infertility, scanning and TEM were performed to check the ultrastructure of the sperm. Western blot and co-immunoprecipitation were performed to clarify the pathological mechanism of the IQUB variant.

MAIN RESULTS AND THE ROLE OF CHANCE

We identified a homozygous nonsense IQUB variant (NM_001282855.2:c.942T> G(p.Tyr314*)) from an infertile male. Iqub KO and KI mice mimicked the infertility phenotype and confirmed IQUB to be the pathogenetic gene. Scanning and TEM showed that sperm of both the mouse models and the affected individual had radial spoke defects. The functional assay suggested that IQUB may recruit calmodulin in lower Ca2+ environments to facilitate the normal assembly of radial spokes by inhibiting the activity of RSPH3/p-ERK1/2 (a nontypical AKAP (A-Kinase Anchoring Protein) forming by RSPH3 and phosphorylation of extracellular signal-regulated kinase 1 and 2 (p-ERK1/2)).

LIMITATIONS, REASONS FOR CAUTION

Additional cases are needed to confirm the genetic contribution of IQUB variants to male infertility. In addition, because no IQUB antibody is available for immunofluorescence and the polyclonal antibody we generated was only effective in western blotting, immunostaining for IQUB was not performed in this study. Therefore, this study lacks direct in vivo proof to confirm the effect of the variant on IQUB protein level.

WIDER IMPLICATIONS OF THE FINDINGS

Our results suggest a causal relation between IQUB variants and male infertility owing to asthenospermia, and partly clarify the pathological mechanism of IQUB variants. This expands our knowledge of the genes involved in human sperm asthenospermia and potentially provides a new genetic marker for male infertility.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the National Key Research and Development Program of China (2021YFC2700100), the National Natural Science Foundation of China (32130029, 82171643, 81971450, 82001538, and 81971382) and the Guangdong Science and Technology Department Guangdong-Hong Kong-Macao Joint Innovation Project (2020A0505140003). There are no competing interests to declare.

TRIAL REGISTRATION NUMBER

N/A.

Further Insights on RNA Expression and Sperm Motility

Asthenozoospermia is one of the main causes of male infertility and it is characterized by reduced sperm motility. Several mutations in genes that code for structural or functional constituents of the sperm have already been identified as known causes of asthenozoospermia. In contrast, the role of sperm RNA in regulating sperm motility is still not fully understood. Consequently, here we aim to contribute to the knowledge regarding the expression of sperm RNA, and ultimately, to provide further insights into its relationship with sperm motility. We investigated the expression of a group of mRNAs by using real-time PCR (CATSPER3, CFAP44, CRHR1, HIP1, IQCG KRT34, LRRC6, QRICH2, RSPH6A, SPATA33 and TEKT2) and the highest score corresponding to the target miRNA for each mRNA in asthenozoospermic and normozoospermic individuals. We observed a reduced expression of all mRNAs and miRNAs in asthenozoospermic patients compared to controls, with a more accentuated reduction in patients with progressive sperm motility lower than 15%. Our work provides further insights regarding the role of RNA in regulating sperm motility. Further studies are required to determine how these genes and their corresponding miRNA act regarding sperm motility, particularly KRT34 and CRHR1, which have not previously been seen to play a significant role in regulating sperm motility.

Genotype-phenotype association and variant characterization in Diamond-Blackfan anemia caused by pathogenic variants in RPL35A

Diamond Blackfan anemia (DBA) is predominantly an autosomal dominant inherited red cell aplasia primarily caused by pathogenic germline variants in ribosomal protein genes. DBA due to pathogenic RPL35A variants has been associated with large 3q29 deletions and phenotypes not common in DBA. We conducted a multi-institutional genotypephenotype study of 45 patients with DBA associated with pathogenic RPL35A germline variants and curated the variant data on 21 additional cases from the literature. Genotype-phenotype analyses were conducted comparing patients with large deletions versus all other pathogenic variants in RPL35A. Twenty-two of the 45 cases had large deletions in RPL35A. After adjusting for multiple tests, a statistically significant association was observed between patients with a large deletion and steroid-resistant anemia, neutropenia, craniofacial abnormalities, chronic gastrointestinal problems, and intellectual disabilities (P<0.01) compared with all other pathogenic variants. Non-large deletion pathogenic variants were spread across RPL35Awith no apparent hot spot and 56% of the individual family variants were observed more than once. In this, the largest known study of DBA patients with pathogenic RPL35A variants, we determined that patients with large deletions have a more severe phenotype that is clinically different from those with non-large deletion variants. Genes of interest also deleted in the 3q29 region that could be associated with some of these phenotypic features include LMLN and IQCG. Management of DBA due to large RPL35A deletions may be challenging due to complex problems and require comprehensive assessments by multiple specialists including immunological, gastrointestinal, and developmental evaluations to provide optimal multidisciplinary care.

Structural Aspects and Prediction of Calmodulin-Binding Proteins

Calmodulin (CaM) is an important intracellular protein that binds Ca²⁺ and functions as a critical second messenger involved in numerous biological activities through extensive interactions with proteins and peptides. CaM’s ability to adapt to binding targets with different structures is related to the flexible central helix separating the N- and C-terminal lobes, which allows for conformational changes between extended and collapsed forms of the protein. CaM-binding targets are most often identified using prediction algorithms that utilize sequence and structural data to predict regions of peptides and proteins that can interact with CaM. In this review, we provide an overview of different CaM-binding proteins, the motifs through which they interact with CaM, and shared properties that make them good binding partners for CaM. Additionally, we discuss the historical and current methods for predicting CaM binding, and the similarities and differences between these methods and their relative success at prediction. As new CaM-binding proteins are identified and classified, we will gain a broader understanding of the biological processes regulated through changes in Ca²⁺ concentration through interactions with CaM.

Structural Diversity in Calmodulin - Peptide Interactions

Calmodulin (CaM) is a highly conserved eukaryotic Ca2+ sensor protein that is able to bind a large variety of target sequences without a defined consensus sequence. The recognition of this diverse target set allows CaM to take part in the regulation of several vital cell functions. To fully understand the structural basis of the regulation functions of CaM, the investigation of complexes of CaM and its targets is essential. In this minireview we give an outline of the different types of CaM - peptide complexes with 3D structure determined, also providing an overview of recently determined structures. We discuss factors defining the orientations of peptides within the complexes, as well as roles of anchoring residues. The emphasis is on complexes where multiple binding modes were found.

IQ motif containing D (IQCD), a new acrosomal protein involved in the acrosome reaction and fertilisation

The acrosome is single, large, dense-core secretory granule overlying the nucleus of most mammalian spermatozoa. Its exocytosis, the acrosome reaction, is a crucial event during fertilisation. In this study we identified a new acrosome-associated gene, namely IQ motif containing D (IQCD), expressed nearly in multiple tissues with highest expression levels in the testis. In mouse testis, Iqcd transcript accumulated from Postnatal Day (PND) 1 to adulthood. However, expression of IQCD protein at the testicular development stage started primarily from PND 18 and increased in an age-dependent manner until plateauing in adulthood. IQCD was primarily accumulated in the acrosome area of round and elongating spermatids within seminiferous tubules of the testes during the late stage of spermiogenesis; this immunolocalisation pattern is similar in mice and humans. IQCD levels in spermatozoa were significantly lower in IVF patients with total fertilisation failure or a low fertilisation rate than in healthy men. Anti-IQCD antibody significantly inhibited the acrosome reaction and slightly reduced protein tyrosine phosphorylation levels in human spermatozoa, but specifically blocked murine IVF. IQCD interacted with mammalian homolog of C. elegans uncoordinated gene 13 (Munc13) in spermatozoa and may participate in acrosome exocytosis. In conclusion, this study identified a new acrosomal protein, namely IQCD, which is involved in fertilisation and the acrosome reaction.

Effect of Ca2+ on the promiscuous target-protein binding of calmodulin

Calmodulin (CaM) is a calcium sensing protein that regulates the function of a large number of proteins, thus playing a crucial part in many cell signaling pathways. CaM has the ability to bind more than 300 different target peptides in a Ca2+-dependent manner, mainly through the exposure of hydrophobic residues. How CaM can bind a large number of targets while retaining some selectivity is a fascinating open question. Here, we explore the mechanism of CaM selective promiscuity for selected target proteins. Analyzing enhanced sampling molecular dynamics simulations of Ca2+-bound and Ca2+-free CaM via spectral clustering has allowed us to identify distinct conformational states, characterized by interhelical angles, secondary structure determinants and the solvent exposure of specific residues. We searched for indicators of conformational selection by mapping solvent exposure of residues in these conformational states to contacts in structures of CaM/target peptide complexes. We thereby identified CaM states involved in various binding classes arranged along a depth binding gradient. Binding Ca2+ modifies the accessible hydrophobic surface of the two lobes and allows for deeper binding. Apo CaM indeed shows shallow binding involving predominantly polar and charged residues. Furthermore, binding to the C-terminal lobe of CaM appears selective and involves specific conformational states that can facilitate deep binding to target proteins, while binding to the N-terminal lobe appears to happen through a more flexible mechanism. Thus the long-ranged electrostatic interactions of the charged residues of the N-terminal lobe of CaM may initiate binding, while the short-ranged interactions of hydrophobic residues in the C-terminal lobe of CaM may account for selectivity. This work furthers our understanding of the mechanism of CaM binding and selectivity to different target proteins and paves the way towards a comprehensive model of CaM selectivity.

Widespread alternative exon usage in clinically distinct subtypes of Invasive Ductal Carcinoma

Cancer cells can have different patterns of exon usage of individual genes when compared to normal tissue, suggesting that alternative splicing may play a role in shaping the tumor phenotype. The discovery and identification of gene variants has increased dramatically with the introduction of RNA-sequencing technology, which enables whole transcriptome analysis of known, as well as novel isoforms. Here we report alternative splicing and transcriptional events among subtypes of invasive ductal carcinoma in The Cancer Genome Atlas (TCGA) Breast Invasive Carcinoma (BRCA) cohort. Alternative exon usage was widespread, and although common events were shared among three subtypes, ER+ HER2−, ER− HER2−, and HER2+, many events on the exon level were subtype specific. Additional RNA-seq analysis was carried out in an independent cohort of 43 ER+ HER2− and ER− HER2− primary breast tumors, confirming many of the exon events identified in the TCGA cohort. Alternative splicing and transcriptional events detected in five genes, MYO6, EPB41L1, TPD52, IQCG, and ACOX2 were validated by qRT-PCR in a third cohort of 40 ER+ HER2− and ER− HER2− patients, showing that these events were truly subtype specific.

Characterization of a Sea Urchin IQ Motif Containing Protein D as a Coactivator of Nuclear Receptors

Nuclear receptor (NR) interacting proteins, such as coactivators and corepressors, play a crucial role in specifying the transcriptional activity of the receptor. However, little is known about the functional features of the NR coregulators in marine invertebrates. Using the yeast two-hybrid screening method, a sea urchin oocyte cDNA library was screened for proteins that interact with the ligand-binding domain of human RXRα (hRXRα) as the bait protein in the presence of 9-cis retinoic acid. Here, we describe IQ motif containing protein D (IQCD) as an RXR-interacting coactivator. The open reading frame of Strongylocentrotus nudus IQCD (SnIQCD) cDNA contains 1464 bp encoding a protein of 487 amino acids. SnIQCD and the vertebrate IQCDs contain well-conserved C-terminal IQ motifs and coiled-coil domains. The interactions between RXRα and IQCD were confirmed by an immunoprecipitation assay and a mammal two-hybrid assay. RXRα preferentially interacted with the C-terminal half including IQ motif than the N-terminal half of SnIQCD. The coactivator interacting LXXLL motif in SnIQCD is not directly involved in the interaction with RXRα. SnIQCD overexpression increased the basal RXR transactivation of a RXR-responsive reporter gene. Furthermore, SnIQCD enhanced the transcriptional activity of RXR heterodimeric partners such as RAR, PPAR, and the steroid hormone receptor family members from mammals, teleost fish, and sea urchin. Taken together, we suggest that IQCD orthologs are able to function as transcriptional coactivators cooperating with NRs.

Calcium triggers reversal of calmodulin on nested anti-parallel sites in the IQ motif of the neuronal voltage-dependent sodium channel NaV1.2

Several members of the voltage-gated sodium channel family are regulated by calmodulin (CaM) and ionic calcium. The neuronal voltage-gated sodium channel NaV1.2 contains binding sites for both apo (calcium-depleted) and calcium-saturated CaM. We have determined equilibrium dissociation constants for rat NaV1.2 IQ motif [IQRAYRRYLLK] binding to apo CaM (~3nM) and (Ca2+)4-CaM (~85nM), showing that apo CaM binding is favored by 30-fold. For both apo and (Ca2+)4-CaM, NMR demonstrated that NaV1.2 IQ motif peptide (NaV1.2IQp) exclusively made contacts with C-domain residues of CaM (CaMC). To understand how calcium triggers conformational change at the CaM-IQ interface, we determined a solution structure (2M5E.pdb) of (Ca2+)2-CaMC bound to NaV1.2IQp. The polarity of (Ca2+)2-CaMC relative to the IQ motif was opposite to that seen in apo CaMC-Nav1.2IQp (2KXW), revealing that CaMC recognizes nested, anti-parallel sites in Nav1.2IQp. Reversal of CaM may require transient release from the IQ motif during calcium binding, and facilitate a re-orientation of CaMN allowing interactions with non-IQ NaV1.2 residues or auxiliary regulatory proteins interacting in the vicinity of the IQ motif.

Ca2+-Induced Rigidity Change of the Myosin VIIa IQ Motif-Single α Helix Lever Arm Extension

Several unconventional myosins contain a highly charged single α helix (SAH) immediately following the calmodulin (CaM) binding IQ motifs, functioning to extend lever arms of these myosins. How such SAH is connected to the IQ motifs and whether the conformation of the IQ motifs-SAH segments are regulated by Ca²⁺ fluctuations are not known. Here, we demonstrate by solving its crystal structure that the predicted SAH of myosin VIIa (Myo7a) forms a stable SAH. The structure of Myo7a IQ5-SAH segment in complex with apo-CaM reveals that the SAH sequence can extend the length of the Myo7a lever arm. Although Ca²⁺-CaM remains bound to IQ5-SAH, the Ca²⁺-induced CaM binding mode change softens the conformation of the IQ5-SAH junction, revealing a Ca²⁺-induced lever arm flexibility change for Myo7a. We further demonstrate that the last IQ motif of several other myosins also binds to both apo- and Ca²⁺-CaM, suggesting a common Ca²⁺-induced conformational regulation mechanism.

Inhibition of the nuclear export of p65 and IQCG in leukemogenesis by NUP98-IQCG

NUP98 fuses with approximately 34 different partner genes via translocation in hematological malignancies. Transgenic or retrovirus-mediated bone marrow transplanted mouse models reveal the leukemogenesis of some NUP98-related fusion genes. We previously reported the fusion protein NUP98-IQ motif containing G (IQCG) in a myeloid/T lymphoid bi-phenoleukemia patient with t(3;11) and confirmed its leukemogenic ability. Herein, we demonstrated the association of NUP98-IQCG with CRM1, and found that NUP98-IQCG expression inhibits the CRM1-mediated nuclear export of p65 and enhances the transcriptional activity of nuclear factor-κB. Moreover, IQCG could be entrapped in the nucleus by NUP98-IQCG, and the fusion protein interacts with calmodulin via the IQ motif in a calcium-independent manner. Therefore, the inhibition of nuclear exports of p65 and IQCG might contribute to the leukemogenesis of NUP98-IQCG.