Transcription and expression of Plasmodium falciparum histidine-rich proteins in different stages and strains : implications for rapid diagnostic tests

Background: Although rapid diagnostic tests (RDTs) for Plasmodium falciparum infection that target histidine rich protein 2 (PfHRP2) are generally sensitive, their performance has been reported to be variable. One possible explanation for variable test performance is differences in expression level of PfHRP in different parasite isolates. Methods: Total RNA and protein were extracted from synchronised cultures of 7 P. falciparum lines over 5 time points of the life cycle, and from synchronised ring stages of 10 falciparum lines. Using quantitative real-time polymerase chain reaction, Western blot analysis and ELISA we investigated variations in the transcription and protein levels of pfhrp2, pfhrp3 and PfHRP respectively in the different parasite lines, over the parasite intraerythrocytic life cycle. Results: Transcription of pfhrp2 and pfhrp3 in different parasite lines over the parasite life cycle was observed to vary relative to the control parasite K1. In some parasite lines very low transcription of these genes was observed. The peak transcription was observed in ring-stage parasites. Pfhrp2 transcription was observed to be consistently higher than pfhrp3 transcription within parasite lines. The intraerythrocytic lifecycle stage at which the peak level of protein was present varied across strains. Total protein levels were more constant relative to total mRNA transcription, however a maximum 24 fold difference in expression at ring-stage parasites relative to the K1 strain was observed. Conclusions: The levels of transcription of pfhrp2 and pfhrp3, and protein expression of PfHRP varied between different P. falciparum strains. This variation may impact on the detection sensitivity of PfHRP2-detecting RDTs.

Plasmodium falciparum parasites lacking histidine-rich protein 2 and 3 : a review and recommendations for accurate reporting

Malaria rapid diagnostic tests (RDTs) play a critical role in malaria case management, surveillance and case investigations. Test performance is largely determined by design and quality characteristics, such as detection sensitivity, specificity, and thermal stability. However, parasite characteristics such as variable or absent expression of antigens targeted by RDTs can also affect RDT performance. Plasmodium falciparum parasites lacking the PfHRP2 protein, the most common target antigen for detection of P. falciparum, have been reported in some regions. Therefore, accurately mapping the presence and prevalence of P. falciparum parasites lacking pfhrp2 would be an important step so that RDTs targeting alternative antigens, or microscopy, can be preferentially selected for use in such regions. Herein the available evidence and molecular basis for identifying malaria parasites lacking PfHRP2 is reviewed, and a set of recommended procedures to apply for future investigations for parasites lacking PfHRP2, is proposed.

Modelling the dynamics of Plasmodium falciparum histidine-rich protein 2 in human malaria to better understand malaria rapid diagnostic test performance

BACKGROUND: Effective diagnosis of malaria is a major component of case management. Rapid diagnostic tests (RDTs) based on Plasmodium falciparumhistidine-rich protein 2 (PfHRP2) are popular for diagnosis of this most virulent malaria infection. However, concerns have been raised about the longevity of the PfHRP2 antigenaemia following curative treatment in endemic regions. METHODS: A model of PfHRP2 production and decay was developed to mimic the kinetics of PfHRP2 antigenaemia during infections. Data from two human infection studies was used to fit the model, and to investigate PfHRP2 kinetics. Four malaria RDTs were assessed in the laboratory to determine the minimum detectable concentration of PfHRP2. RESULTS: Fitting of the PfHRP2 dynamics model indicated that in malaria naive hosts, P. falciparum parasites of the 3D7 strain produce 1.4 x 10(-)(1)(3) g of PfHRP2 per parasite per replication cycle. The four RDTs had minimum detection thresholds between 6.9 and 27.8 ng/mL. Combining these detection thresholds with the kinetics of PfHRP2, it is predicted that as few as 8 parasites/muL may be required to maintain a positive RDT in a chronic infection. CONCLUSIONS: The results of the model indicate that good quality PfHRP2-based RDTs should be able to detect parasites on the first day of symptoms, and that the persistence of the antigen will cause the tests to remain positive for at least seven days after treatment. The duration of a positive test result following curative treatment is dependent on the duration and density of parasitaemia prior to treatment and the presence and affinity of anti-PfHRP2 antibodies.

Identification of optimal epitopes for Plasmodium falciparum rapid diagnostic tests that target histidine-rich proteins 2 and 3

Rapid diagnostic tests (RDTs) represent important tools to diagnose malaria infection. To improve understanding of the variable performance of RDTs that detect the major target in Plasmodium falciparum, namely, histidine-rich protein 2 (HRP2), and to inform the design of better tests, we undertook detailed mapping of the epitopes recognized by eight HRP-specific monoclonal antibodies (MAbs). To investigate the geographic skewing of this polymorphic protein, we analyzed the distribution of these epitopes in parasites from geographically diverse areas. To identify an ideal amino acid motif for a MAb to target in HRP2 and in the related protein HRP3, we used a purpose-designed script to perform bioinformatic analysis of 448 distinct gene sequences from pfhrp2 and from 99 sequences from the closely related gene pfhrp3. The frequency and distribution of these motifs were also compared to the MAb epitopes. Heat stability testing of MAbs immobilized on nitrocellulose membranes was also performed. Results of these experiments enabled the identification of MAbs with the most desirable characteristics for inclusion in RDTs, including copy number and coverage of target epitopes, geographic skewing, heat stability, and match with the most abundant amino acid motifs identified. This study therefore informs the selection of MAbs to include in malaria RDTs as well as in the generation of improved MAbs that should improve the performance of HRP-detecting malaria RDTs.

Global sequence variation in the histidine-rich proteins 2 and 3 of Plasmodium falciparum: implications for the performance of malaria rapid diagnostic tests

Background Accurate diagnosis is essential for prompt and appropriate treatment of malaria. While rapid diagnostic tests (RDTs) offer great potential to improve malaria diagnosis, the sensitivity of RDTs has been reported to be highly variable. One possible factor contributing to variable test performance is the diversity of parasite antigens. This is of particular concern for Plasmodium falciparum histidine-rich protein 2 (PfHRP2)-detecting RDTs since PfHRP2 has been reported to be highly variable in isolates of the Asia-Pacific region. Methods The pfhrp2 exon 2 fragment from 458 isolates of P. falciparum collected from 38 countries was amplified and sequenced. For a subset of 80 isolates, the exon 2 fragment of histidine-rich protein 3 (pfhrp3) was also amplified and sequenced. DNA sequence and statistical analysis of the variation observed in these genes was conducted. The potential impact of the pfhrp2 variation on RDT detection rates was examined by analysing the relationship between sequence characteristics of this gene and the results of the WHO product testing of malaria RDTs: Round 1 (2008), for 34 PfHRP2-detecting RDTs. Results Sequence analysis revealed extensive variations in the number and arrangement of various repeats encoded by the genes in parasite populations world-wide. However, no statistically robust correlation between gene structure and RDT detection rate for P. falciparum parasites at 200 parasites per microlitre was identified. Conclusions The results suggest that despite extreme sequence variation, diversity of PfHRP2 does not appear to be a major cause of RDT sensitivity variation.

Circulating antibodies against Plasmodium falciparum histidine-rich proteins 2 interfere with antigen detection by rapid diagnostic tests

Background Rapid diagnostic tests (RDTs) for detection of Plasmodium falciparum infection that target P. falciparum histidine-rich protein 2 (PfHRP2), a protein that circulates in the blood of patients infected with this species of malaria, are widely used to guide case management. Understanding determinants of PfHRP2 availability in circulation is therefore essential to understanding the performance of PfHRP2-detecting RDTs. Methods The possibility that pre-formed host anti-PfHRP2 antibodies may block target antigen detection, thereby causing false negative test results was investigated in this study. Results Anti-PfHRP2 antibodies were detected in 19/75 (25%) of plasma samples collected from patients with acute malaria from Cambodia, Nigeria and the Philippines, as well as in 3/28 (10.7%) asymptomatic Solomon Islands residents. Pre-incubation of plasma samples from subjects with high-titre anti-PfHRP2 antibodies with soluble PfHRP2 blocked the detection of the target antigen on two of the three brands of RDTs tested, leading to false negative results. Pre-incubation of the plasma with intact parasitized erythrocytes resulted in a reduction of band intensity at the highest parasite density, and a reduction of lower detection threshold by ten-fold on all three brands of RDTs tested. Conclusions These observations indicate possible reduced sensitivity for diagnosis of P. falciparum malaria using PfHRP2-detecting RDTs among people with high levels of specific antibodies and low density infection, as well as possible interference with tests configured to detect soluble PfHRP2 in saliva or urine samples. Further investigations are required to assess the impact of pre-formed anti-PfHRP2 antibodies on RDT performance in different transmission settings.

Biochemical characterization of Aprataxin, the protein deficient in Ataxia with Oculomotor Apraxia type 1

Neurodegenerative disorders are heterogenous in nature and include a range of ataxias with oculomotor apraxia, which are characterised by a wide variety of neurological and ophthalmological features. This family includes recessive and dominant disorders. A subfamily of autosomal recessive cerebellar ataxias are characterised by defects in the cellular response to DNA damage. These include the well characterised disorders Ataxia-Telangiectasia (A-T) and Ataxia-Telangiectasia Like Disorder (A-TLD) as well as the recently identified diseases Spinocerebellar ataxia with axonal neuropathy Type 1 (SCAN1), Ataxia with Oculomotor Apraxia Type 2 (AOA2), as well as the subject of this thesis, Ataxia with Oculomotor Apraxia Type 1 (AOA1). AOA1 is caused by mutations in the APTX gene, which is located at chromosomal locus 9p13. This gene codes for the 342 amino acid protein Aprataxin. Mutations in APTX cause destabilization of Aprataxin, thus AOA1 is a result of Aprataxin deficiency. Aprataxin has three functional domains, an N-terminal Forkhead Associated (FHA) phosphoprotein interaction domain, a central Histidine Triad (HIT) nucleotide hydrolase domain and a C-terminal C2H2 zinc finger. Aprataxins FHA domain has homology to FHA domain of the DNA repair protein 5’ polynucleotide kinase 3’ phosphatase (PNKP). PNKP interacts with a range of DNA repair proteins via its FHA domain and plays a critical role in processing damaged DNA termini. The presence of this domain with a nucleotide hydrolase domain and a DNA binding motif implicated that Aprataxin may be involved in DNA repair and that AOA1 may be caused by a DNA repair deficit. This was substantiated by the interaction of Aprataxin with proteins involved in the repair of both single and double strand DNA breaks (XRay Cross-Complementing 1, XRCC4 and Poly-ADP Ribose Polymerase-1) and the hypersensitivity of AOA1 patient cell lines to single and double strand break inducing agents. At the commencement of this study little was known about the in vitro and in vivo properties of Aprataxin. Initially this study focused on generation of recombinant Aprataxin proteins to facilitate examination of the in vitro properties of Aprataxin. Using recombinant Aprataxin proteins I found that Aprataxin binds to double stranded DNA. Consistent with a role for Aprataxin as a DNA repair enzyme, this binding is not sequence specific. I also report that the HIT domain of Aprataxin hydrolyses adenosine derivatives and interestingly found that this activity is competitively inhibited by DNA. This provided initial evidence that DNA binds to the HIT domain of Aprataxin. The interaction of DNA with the nucleotide hydrolase domain of Aprataxin provided initial evidence that Aprataxin may be a DNA-processing factor. Following these studies, Aprataxin was found to hydrolyse 5’adenylated DNA, which can be generated by unscheduled ligation at DNA breaks with non-standard termini. I found that cell extracts from AOA1 patients do not have DNA-adenylate hydrolase activity indicating that Aprataxin is the only DNA-adenylate hydrolase in mammalian cells. I further characterised this activity by examining the contribution of the zinc finger and FHA domains to DNA-adenylate hydrolysis by the HIT domain. I found that deletion of the zinc finger ablated the activity of the HIT domain against adenylated DNA, indicating that the zinc finger may be required for the formation of a stable enzyme-substrate complex. Deletion of the FHA domain stimulated DNA-adenylate hydrolysis, which indicated that the activity of the HIT domain may be regulated by the FHA domain. Given that the FHA domain is involved in protein-protein interactions I propose that the activity of Aprataxins HIT domain may be regulated by proteins which interact with its FHA domain. We examined this possibility by measuring the DNA-adenylate hydrolase activity of extracts from cells deficient for the Aprataxin-interacting DNA repair proteins XRCC1 and PARP-1. XRCC1 deficiency did not affect Aprataxin activity but I found that Aprataxin is destabilized in the absence of PARP-1, resulting in a deficiency of DNA-adenylate hydrolase activity in PARP-1 knockout cells. This implies a critical role for PARP-1 in the stabilization of Aprataxin. Conversely I found that PARP-1 is destabilized in the absence of Aprataxin. PARP-1 is a central player in a number of DNA repair mechanisms and this implies that not only do AOA1 cells lack Aprataxin, they may also have defects in PARP-1 dependant cellular functions. Based on this I identified a defect in a PARP-1 dependant DNA repair mechanism in AOA1 cells. Additionally, I identified elevated levels of oxidized DNA in AOA1 cells, which is indicative of a defect in Base Excision Repair (BER). I attribute this to the reduced level of the BER protein Apurinic Endonuclease 1 (APE1) I identified in Aprataxin deficient cells. This study has identified and characterised multiple DNA repair defects in AOA1 cells, indicating that Aprataxin deficiency has far-reaching cellular consequences. Consistent with the literature, I show that Aprataxin is a nuclear protein with nucleoplasmic and nucleolar distribution. Previous studies have shown that Aprataxin interacts with the nucleolar rRNA processing factor nucleolin and that AOA1 cells appear to have a mild defect in rRNA synthesis. Given the nucleolar localization of Aprataxin I examined the protein-protein interactions of Aprataxin and found that Aprataxin interacts with a number of rRNA transcription and processing factors. Based on this and the nucleolar localization of Aprataxin I proposed that Aprataxin may have an alternative role in the nucleolus. I therefore examined the transcriptional activity of Aprataxin deficient cells using nucleotide analogue incorporation. I found that AOA1 cells do not display a defect in basal levels of RNA synthesis, however they display defective transcriptional responses to DNA damage. In summary, this thesis demonstrates that Aprataxin is a DNA repair enzyme responsible for the repair of adenylated DNA termini and that it is required for stabilization of at least two other DNA repair proteins. Thus not only do AOA1 cells have no Aprataxin protein or activity, they have additional deficiencies in PolyADP Ribose Polymerase-1 and Apurinic Endonuclease 1 dependant DNA repair mechanisms. I additionally demonstrate DNA-damage inducible transcriptional defects in AOA1 cells, indicating that Aprataxin deficiency confers a broad range of cellular defects and highlighting the complexity of the cellular response to DNA damage and the multiple defects which result from Aprataxin deficiency. My detailed characterization of the cellular consequences of Aprataxin deficiency provides an important contribution to our understanding of interlinking DNA repair processes.

Kallikrein-related peptidase 4 activation of protease-activated receptor family members and association with prostate cancer

Two areas of particular importance in prostate cancer progression are primary tumour development and metastasis. These processes involve a number of physiological events, the mediators of which are still being discovered and characterised. Serine proteases have been shown to play a major role in cancer invasion and metastasis. The recently discovered phenomenon of their activation of a receptor family known as the protease activated receptors (PARs) has extended their physiological role to that of signaling molecule. Several serine proteases are expressed by malignant prostate cancer cells, including members of the kallikreinrelated peptidase (KLK) serine protease family, and increasingly these are being shown to be associated with prostate cancer progression. KLK4 is highly expressed in the prostate and expression levels increase during prostate cancer progression. Critically, recent studies have implicated KLK4 in processes associated with cancer. For example, the ectopic over-expression of KLK4 in prostate cancer cell lines results in an increased ability of these cells to form colonies, proliferate and migrate. In addition, it has been demonstrated that KLK4 is a potential mediator of cellular interactions between prostate cancer cells and osteoblasts (bone forming cells). The ability of KLK4 to influence cellular behaviour is believed to be through the selective cleavage of specific substrates. Identification of relevant in vivo substrates of KLK4 is critical to understanding the pathophysiological roles of this enzyme. Significantly, recent reports have demonstrated that several members of the KLK family are able to activate PARs. The PARs are relatively new members of the seven transmembrane domain containing G protein coupled receptor (GPCR) family. PARs are activated through proteolytic cleavage of their N-terminus by serine proteases, the resulting nascent N-terminal binds intramolecularly to initiate receptor activation. PARs are involved in a number of patho-physiological processes, including vascular repair and inflammation, and a growing body of evidence suggests roles in cancer. While expression of PAR family members has been documented in several types of cancers, including prostate, the role of these GPCRs in prostate cancer development and progression is yet to be examined. Interestingly, several studies have suggested potential roles in cellular invasion through the induction of cytoskeletal reorganisation and expression of basement membrane-degrading enzymes. Accordingly, this program of research focussed on the activation of the PARs by the prostate cancer associated enzyme KLK4, cellular processing of activated PARs and the expression pattern of receptor and agonist in prostate cancer. For these studies KLK4 was purified from the conditioned media of stably transfected Sf9 insect cells expressing a construct containing the complete human KLK4 coding sequence in frame with a V5 epitope and poly-histidine encoding sequences. The first aspect of this study was the further characterisation of this recombinant zymogen form of KLK4. The recombinant KLK4 zymogen was demonstrated to be activatable by the metalloendopeptidase thermolysin and amino terminal sequencing indicated that thermolysin activated KLK4 had the predicted N-terminus of mature active KLK4 (31IINED). Critically, removal of the pro-region successfully generated a catalytically active enzyme, with comparable activity to a previously published recombinant KLK4 produced from S2 insect cells. The second aspect of this study was the activation of the PARs by KLK4 and the initiation of signal transduction. This study demonstrated that KLK4 can activate PAR-1 and PAR-2 to mobilise intracellular Ca2+, but failed to activate PAR-4. Further, KLK4 activated PAR-1 and PAR-2 over distinct concentration ranges, with KLK4 activation and mobilisation of Ca2+ demonstrating higher efficacy through PAR-2. Thus, the remainder of this study focussed on PAR-2. KLK4 was demonstrated to directly cleave a synthetic peptide that mimicked the PAR-2 Nterminal activation sequence. Further, KLK4 mediated Ca2+ mobilisation through PAR-2 was accompanied by the initiation of the extra-cellular regulated kinase (ERK) cascade. The specificity of intracellular signaling mediated through PAR-2 by KLK4 activation was demonstrated by siRNA mediated protein depletion, with a reduction in PAR-2 protein levels correlating to a reduction in KLK4 mediated Ca2+mobilisation and ERK phosphorylation. The third aspect of this study examined cellular processing of KLK4 activated PAR- 2 in a prostate cancer cell line. PAR-2 was demonstrated to be expressed by five prostate derived cell lines including the prostate cancer cell line PC-3. It was also demonstrated by flow cytometry and confocal microscopy analyses that activation of PC-3 cell surface PAR-2 by KLK4 leads to internalisation of this receptor in a time dependent manner. Critically, in vivo relevance of the interaction between KLK4 and PAR-2 was established by the observation of the co-expression of receptor and agonist in primary prostate cancer and prostate cancer bone lesion samples by immunohistochemical analysis. Based on the results of this study a number of exciting future studies have been proposed, including, delineating differences in KLK4 cellular signaling via PAR-1 and PAR-2 and the role of PAR-1 and PAR-2 activation by KLK4 in prostate cancer cells and bone cells in prostate cancer progression.

An androgen receptor mutation in the MDA-MB-453 cell line model of molecular apocrine breast cancer compromises receptor activity

Recent evidence indicates that the estrogen receptor-a-negative, androgen receptor (AR)- positive molecular apocrine subtype of breast cancer is driven by AR signaling. The MDA-MB-453 cell line is the prototypical model of this breast cancer subtype; its proliferation is stimulated by androgens such as 5a-dihydrotestosterone (DHT) but inhibited by the progestin medroxyprogesterone acetate (MPA) via AR-mediated mechanisms. We report here that the AR gene in MDAMB- 453 cells contains a G-T transversion in exon 7, resulting in a receptor variant with a glutamine to histidine substitution at amino acid 865 (Q865H) in the ligand binding domain. Compared with wild-type AR, the Q865H variant exhibited reduced sensitivity to DHT and MPA in transactivation assays in MDA-MB-453 and PC-3 cells but did not respond to non-androgenic ligands or receptor antagonists. Ligand binding, molecular modeling, mammalian two-hybrid and immunoblot assays revealed effects of the Q865H mutation on ligand dissociation, AR intramolecular interactions, and receptor stability. Microarray expression profiling demonstrated that DHT and MPA regulate distinct transcriptional programs in MDA-MB-453 cells. Gene Set Enrichment Analysis revealed that DHT- but not MPA-regulated genes were associated with estrogen-responsive transcriptomes from MCF-7 cells and the Wnt signaling pathway. These findings suggest that the divergent proliferative responses of MDA-MB-453 cells to DHT and MPA result from the different genetic programs elicited by these two ligands through the AR-Q865H variant. This work highlights the necessity to characterize additional models of molecular apocrine breast cancer to determine the precise role of AR signaling in this breast cancer subtype. Endocrine-Related Cancer (2012) 19 599–613

EUELC project : a multi-centre, multipurpose study to investigate early stage NSCLC, and to establish a biobank for ongoing collaboration

The European Early Lung Cancer (EUELC) project aims to determine if specific genetic alterations occurring in lung carcinogenesis are detectable in the respiratory epithelium. In order to pursue this objective, nonsmall cell lung cancer (NSCLC) patients with a very high risk of developing progressive lung cancer were recruited from 12 centres in eight European countries: France, Germany, southern Ireland, Italy, the Netherlands, Poland, Spain and the UK. In addition, NSCLC patients were followed up every 6 months for 36 months. A European Bronchial Tissue Bank was set up at the University of Liverpool (Liverpool, UK) to optimise the use of biological specimens. The molecular - pathological investigations were subdivided into specific work packages that were delivered by EUELC Partners. The work packages encompassed mutational analysis, genetic instability, methylation profiling, expression profiling utilising immunohistochemistry and chip-based technologies, as well as in-depth analysis of FHIT and RARβ genes, the telomerase catalytic subunit hTERT and genotyping of susceptibility genes in specific pathways. The EUELC project engendered a tremendous collaborative effort, and it enabled the EUELC Partners to establish protocols for assessing molecular biomarkers in early lung cancer with the view to using such biomarkers for early diagnosis and as intermediate end-points in future chemopreventive programmes. Copyright©ERS Journals Ltd 2009.

An improved method for undertaking limiting dilution assays for in vitro cloning of Plasmodium falciparum parasites

Background Obtaining single parasite clones is required for many techniques in malaria research. Cloning by limiting dilution using microscopy-based assessment for parasite growth is an arduous and labor-intensive process. An alternative method for the detection of parasite growth in limiting dilution assays is using a commercial ELISA histidine-rich protein II (HRP2) detection kit. Methods Detection of parasite growth was undertaken using HRP2 ELISA and compared to thick film microscopy. An HRP2 protein standard was used to determine the detection threshold of the HRP2 ELISA assay, and a HRP2 release model was used to extrapolate the amount of parasite growth required for a positive result. Results The HRP2 ELISA was more sensitive than microscopy for detecting parasite growth. The minimum level of HRP2 protein detection of the ELISA was 0.11ng/ml. Modeling of HRP2 release determined that 2,116 parasites are required to complete a full erythrocytic cycle to produce sufficient HRP2 to be detected by the ELISA. Under standard culture conditions this number of parasites is likely to be reached between 8 to 14 days of culture. Conclusions This method provides an accurate and simple way for the detection of parasite growth in limiting dilution assays, reducing time and resources required in traditional methods. Furthermore the method uses spent culture media instead of the parasite-infected red blood cells, enabling culture to continue.

Affinity purification of plasmid DNA directly from crude bacterial cell lysates

We have shown previously that a sequence-specific DNA-binding protein based on the Lac repressor protein can isolate pre-purified DNA efficiently from simple buffer solution but our attempts to purify plasmids directly from crude starting materials were disappointing with unpractically low DNA yields. We have optimized tbe procedure and present a simple affinity methodology whereby plasmid DNA is purified directly by mixing two crude cell lysates, one cell lysate containing the plasmid and the other the protein affinity ligand, without the need for treatment by RNaseA. After IMAC chromatography, high purity supercoiled DNA is recovered in good yields of 100-150 μg plasmid per 200 mL shake flask culture. Moreover, the resulting DNA is free from linear or open-circular plasmid DNA, genomic DNA, RNA, and protein, to the limits of our detection. Furthermore, we show that lyophilized affinity ligand can be stored at room temperature and re-hydrated for use when required.

Comparative analysis and distribution of Omega-3 lcPUFA Biosynthesis genes in marine molluscs

Recent research has identified marine molluscs as an excellent source of omega-3 long-chain polyunsaturated fatty acids (lcPUFAs), based on their potential for endogenous synthesis of lcPUFAs. In this study we generated a representative list of fatty acyl desaturase (Fad) and elongation of very long-chain fatty acid (Elovl) genes from major orders of Phylum Mollusca, through the interrogation of transcriptome and genome sequences, and various publicly available databases. We have identified novel and uncharacterised Fad and Elovl sequences in the following species: Anadara trapezia, Nerita albicilla, Nerita melanotragus, Crassostrea gigas, Lottia gigantea, Aplysia californica, Loligo pealeii and Chlamys farreri. Based on alignments of translated protein sequences of Fad and Elovl genes, the haeme binding motif and histidine boxes of Fad proteins, and the histidine box and seventeen important amino acids in Elovl proteins, were highly conserved. Phylogenetic analysis of aligned reference sequences was used to reconstruct the evolutionary relationships for Fad and Elovl genes separately. Multiple, well resolved clades for both the Fad and Elovl sequences were observed, suggesting that repeated rounds of gene duplication best explain the distribution of Fad and Elovl proteins across the major orders of molluscs. For Elovl sequences, one clade contained the functionally characterised Elovl5 proteins, while another clade contained proteins hypothesised to have Elovl4 function. Additional well resolved clades consisted only of uncharacterised Elovl sequences. One clade from the Fad phylogeny contained only uncharacterised proteins, while the other clade contained functionally characterised delta-5 desaturase proteins. The discovery of an uncharacterised Fad clade is particularly interesting as these divergent proteins may have novel functions. Overall, this paper presents a number of novel Fad and Elovl genes suggesting that many mollusc groups possess most of the required enzymes for the synthesis of lcPUFAs.