Donor DNA designing
This study discovered four target sites in the 2065 bp-long zebrafish
lepb gene sequence for the sin hyperactivated recombinase-based RGR platform. Two target sites in the direct strand of the
lepb gene were utilized to locate the target locus to be mutated at the 1738
th position, in exon 1. A guanine (G) to thymine (T) substitution was designed to create a nonsense mutation, enabling a premature upstream stop codon (GAG to TAG). The corresponding amino acid relation with the target locus revealed glutamate (E97) in the 168 aa-long zebrafish
lepb protein. The truncated
lepb protein was created by substituting glutamate (‘E97’) with ‘#’ (Stop codon), resulting in a shortened protein size from 168 aa to 96 aa.
The wild and truncated
lepb of
D. rerio were considered membrane proteins due to their transmembrane region (residues 2-24, KSSMIFCLLISSLVAVSISRPTA) of the protein (Fig 1a and 1b). The wild
lepb contains four cysteine residues (C8, C112, C148 and C157) (Fig 1a) and the C8 only remained in the truncated
lepb (Fig1b). An intramolecular disulfide bridge formed by two cysteine residues (especially from the CD loop and at the C-terminal end) is crucial for the 3D structure, stability and bioactivity of the
lep protein
(Bakshi et al., 2022). Like other vertebrates
(Cai et al., 2024), wild
lepb of zebrafish also has a disulfide bridge forming two cysteine residues, such as C112 in the CD loop and C157 at the C-terminal end, which are absent in the truncated
lepb protein, resulting in reduced stability, improper folding and impaired bioactivity of the
lepb protein.
The wild
lepb contains 52.98% alpha helix and 3.57% beta-turn, while the truncated
lepb has 51.04% alpha helix and 2.08% beta-turn. The alpha helix and beta sheets are crucial for the secondary and tertiary structures of the protein (
Rajeshwar and Abraham 2024). The truncation results in a reduction of amino acids in these critical regions, affecting the protein’s secondary and tertiary structures. Key features such as a conserved disulfide bridge and hydrophobic amino acids are crucial for the stable four-helix packing required for proper protein folding and function across vertebrates
(Londraville et al., 2017; Al Mughram et al., 2023). Zebrafish wild
lepb also consists of four anti-parallel alpha helices denoted as helix-A (residues 25-47, PEDRIRIIARTTISRIKKIKDEH), helix-B (residues 63-80, PIDGLSSVLSYLSYLQLR), helix-C (residues 86-109, AQHLQQVQIDLETLLRTLEELAVS) and helix-D (residues 133-153, NYLHLLELQRFLEKLCLNIDK), linked with two long AB (residues 48-62, FQMSPEIDFGPDIDN) and CD (residues 110-132, QGCPLPNPETPVHKEETAFPVTS) loops and one short BC (residues 81-85, LHVPP) loop (Fig 1a) like
lep of other vertebrates
(Bakshi et al., 2022; Casado et al., 2023).
Like other vertebrates’
lep (Greco et al., 2021), the zebrafish
lepb interacts with its receptor (
lepr) via three binding sites. Site-I forms the hexameric
lep-
lepr complex through residues of the C-terminus of helix-D and AB loop, which are essential for receptor activation
(Peelman et al., 2006). Site-II contains residues of helix-A and helix-C, facilitating the high-affinity interaction with the CRH2 domain of
lepr (Iserentant et al., 2005), while site-III connects the N-terminus of helix-D with the Ig-like domain of
lepr (Rock et al., 1996), enabling the conformational changes to activate the
lepr. However, the truncated
lepb fails to activate
lepr due to the lack of partial helix-C, the complete CD loop, helix-D and the C-terminal end (Fig 1b), which disrupts complex formation and receptor dimerization, thereby inhibiting signalling pathways (JAK2/STAT3, JAK2/STAT5, PI3K/IRS/AKT and SHP2/ERK pathways). Especially, the anorexigenic action of
lepb will be restricted, allowing fish to take a feed effectively and gain more body weight. Therefore, the fish growth rate will be increased by enhanced FCR and will quickly attain the marketable size compared to non-edited fish (
Del Vecchio et al., 2021).
Additionally, the 3D models of wild
lepb and truncated
lepb proteins of
D. rerio were constructed using the SWISS-MODEL server with the A0A671LBE5.1.A template further confirmed the structural modifications resulting from the mutation. The global model quality estimate (GMQE) scores were 0.80 and 0.82 for wild
lepb and truncated
lepb proteins, respectively, indicating satisfactory alignment quality, ranging from 0 to 1
(Waterhouse et al., 2024). High GMQE scores supported the model validation. The Qualitative Model Energy ANalysis (QMEAN) scores were -0.84 for wild
lepb and -1.17 for truncated
lepb, indicating a slight decline in model quality after truncation. However, scores were near zero, exhibiting strong alignment with experimental structures
(Benkert et al., 2011). In MolProbity analysis, the residues in the favoured region of wild
lepb and truncated
lepb proteins of
D. rerio were 90.30% and 94.68%, respectively. The 1.82% outliers were in wild
lepb, while truncated
lepb had none. Both models showed over 90% of residues in the favoured region and less than 2% in the outlier region, signifying reliable model quality
(Sobolev et al., 2020). The molprobity scores
(Williams et al., 2018) were 1.12 for wild
lepb and 1.34 for truncated
lepb with the latter suggesting decreased quality due to structural alterations. The Ramachandran plots ensured exceptional quality, with most residues in the favoured (green) and allowed (light green) regions. The ProSA-web Z-scores were -7.02 for wild
lepb and -3.53 for truncated
lepb, further confirming the proteins’ validity, as negative Z-scores
(Elengoe et al., 2014; Singh et al., 2018) indicate good quality. The various analyses revealed that the nonsense mutation at the 1738
th nucleotide base (G) (locus) in the zebrafish
lepb can effectively knock out its expression.
Construction of Donor DNA
To construct the donor DNA (pAS67), gDNA was isolated from zebrafish tissue and pectoral fin and the
lepb gene fragment (Fr-20) was amplified through PCR with a targeted amplicon length of 1808 bp. The Fr-20 contains the potential locus flanked by RGR target sites. After gel elution, the purified Fr-20 was cloned into the pJET1.2 blunt vector (plasmid pAS66) and transformed into the competent DH5α strain,
E. coli cells, where six well-separated colonies were screened using Colony PCR to confirm the existence of Fr-20. The recombinant plasmid DNAs (pAS66 A and B) were isolated from the two putative transformants and Sanger DNA sequencing verified the presence of insert-DNA (Fr-20) with a potential locus, guanine (G) (Fig 2a). Subsequently, Donor DNA (pAS67) containing mutagenic amplified product (Fr-20*) was constructed using the recombinant plasmid (pAS66) as a template according to the Q5 SDM kit protocol, resulting in a premature upstream stop codon (GAG to TAG) mutation by substituting the guanine (G) with thymine (T) at the 2016
th nucleotide position of recombinant plasmid pAS66. The mutagenic recombinant plasmids (pAS67 A and B) were isolated from the two putative mutated recombinants. The success of the mutagenesis was confirmed through sequencing, ensuring a guanine (G) substituted precisely by thymine (T) (Fig 2b) at the expected nucleotide position in the mutagenic donor DNA plasmids (pAS67 A and B). Successful donor DNA construction was achieved through site-directed mutagenesis, with sequencing verifying the intended mutation, thymine (T) substituting guanine (G) at the target site. The sequencing chromatograms showed distinct peaks confirming the presence of the mutation at the target site, with no additional alterations detected in the surrounding regions, confirming the specificity and accuracy of the mutagenesis process (Fig 2b). Additionally, no off-target mutations were observed, validating the precision of the experimental approach. It supports that the anorexigenic function of the zebrafish
lepb can effectively be knocked out by using the generated donor DNA in this study.
Generally, zebrafish
lep knocked out mutants exhibit increased feed intake and an obese phenotype
(Audira et al., 2018; Hu et al., 2022). In future
in vivo validation studies, the genotyping, multi-omics integration, body length-weight measurement, 3D locomotor activity, aggressiveness test, circadian rhythm test, passive avoidance test (Short-term memory test), colour preferences assay, predator avoidance test, biochemical assays (to determine the appetite-controlling hormone, oxidative, neurotransmitter, anti-oxidative capacity, lipid peroxidation, stress hormone, DNA damage and inflammation markers) and proximate analysis can be used further to examine the changes occurred in the
lepb knocked out mutants physiological, phenotypical and behavioural outcomes
(Audira et al., 2018; Del Vecchio et al., 2021).