The fibroblast growth factor receptors (FGFRs) share a canonical structure comprising an extracellular region with three immunoglobulin-like (Ig-like) domains, a stretch of acidic residues between the first and second Ig-like domains termed the acidic box, a transmembrane helix, and an intracellular domain with tyrosine kinase activity (
Mohammadi et al. 2005). Binding of FGFs and heparan sulfate to the FGFRs induces dimerization (
Mohammadi et al. 2005). An active complex of two FGFs, two FGFRs and one or two heparan sulfate chains results in trans-phosphorylation of specific tyrosine kinase residues in the intra cellular domain of the FGFR, recruiting adapter proteins and enabling downstream signaling by diverse pathways including the PLCγ pathway, P-I-3 Kinase-AKT/PKB pathway and the Ras-MAP kinase pathway (
Dailey et al. 2005).
Fgfrl1 is the most diverged of the Fgfrs (
Wiedemann and Trueb 2000;
Kim et al. 2001). Although it is structurally related to the other Fgfrs, with the NH
2-terminal of the protein having three Ig-like domains, six cysteines, an acidic box and one transmembrane domain, Fgfrl1 lacks both the histidine–alanine–valine motif that is critical for binding other FGFR molecules and the tyrosine kinase domain that is required for intracellular signaling (
Wiedemann and Trueb 2001;
Kim et al. 2001;
Sleeman et al. 2001). Fgfrl1 is increased at sites of cell–cell contact and the protein is hypothesized to promote cell adhesion (
Rieckmann et al. 2008). Fgfrl1 can therefore be considered similar to a nectin, a class of Ig-like, transmembrane cell adhesion molecule required for contact inhibition of cell movement and proliferation (
Takai et al. 2008). Finally, there is some evidence that Fgfrl1 may act as a decoy receptor to coordinate the distribution of free FGF ligand and inhibit FGF signaling, as Fgfrl1 can inhibit the FGF-stimulated growth of MG63 cells (
Trueb et al. 2003) and inhibit the FGF-stimulated activity of an FGF-inducible responsive element in a luciferase experiment (
Rieckmann et al. 2009).
Our cDNA arrays examined comparative gene expression from the diaphragms of sex-matched,
Fgfrl1 homozygous null mice and wildtype littermates at E18.5. The murine diaphragm forms between E11.5 and E13.5, (
Leinwand et al. 2002), but E18.5 was chosen because of the ease of diaphragmatic dissection at this later stage of gestation. We also believe that an examination of this time period is still useful, as at least some genes that have been implicated in diaphragm development (for example,
Slit3 and
Fog2), are expressed at this later stage of development in addition to earlier time periods. The expression of
Fgfrl1 was reduced in the
Fgfrl1 homozygous null mice compared to wildtype littermates (;
p = 0.004), validating the array methodology. Only seven genes besides
Fgfrl1 showed significant decreases in expression level in the mutant diaphragms when compared with wildtype littermates:
Tpm3, Myl2, Lrtm1, Myh4, Myl3, Myh7 and
Hephl1. Tropomyosin-3 (
Tpm3; OMIM 191030) showed the greatest reduction in expression in the
Fgfrl1 mutant mice (). This gene encodes a component of the thin filaments of the sarcomere, and loss of function mutations in
TPM3 have been reported in patients with nemaline myopathy and congenital fiber type disproportion (
Laing et al. 1995;
Clarke et al. 2008;
Lehtokari et al. 2008). Interestingly, diaphragm hypoplasia is a recognized complication of severe nemaline myopathy (
Wada et al. 1996;
Danhaive et al. 2007). However, in patients with
TPM3 loss of function mutations, the diameter of the type I slow fibers is reduced, as
TPM3 encodes the slow isoform of skeletal muscle α-tropomyosin (
Lehtokari et al. 2008). As the
Fgfrl1 null mice had normal muscle fiber histology and the ratio of fiber subtypes was unaltered when compared with normal littermates (
Baertschi et al. 2007), this suggests that the reduction in
Tpm3 expression in the
Fgfrl1 null mice was a secondary phenomenon, attributable either to the deficiency of
Fgfrl1 and/or to the muscular hypoplasia of the diaphragm.
The expression arrays also demonstrated reduced expression for four other genes that encode sarcomeric proteins: Myosin, light chain-2 (
Myl2; OMIM 160781), Myosin, heavy chain-4 (
Myh4; OMIM 160742), Myosin, light chain-3 (
Myl3; OMIM 160790), and Myosin, heavy chain-7 (
Myh7, OMIM 160760).
Myl3 and
Myh7 both encode slow-twitch, type I muscle isoforms (
Jandreski et al. 1987), whereas
Myl2 regulates myosin ATPase activity in smooth muscle. Heterozygous, missense
Myl2 mutations have been described in hypertrophic cardiomyopathy (
Poetter et al. 1996;
Flavigny et al. 1998;
Richard et al. 2003). Mutations in
Myh7 have also been described in hypertrophic and dilated cardiomyopathies (
García-Castro et al. 2003;
Richard et al. 2003;
Hershberger et al. 2008;
Møller et al. (2009) and have been estimated to account for up to 40–50% of hypertrophic cardiomyopathies.
Myh7 mutations have also been reported in myosin storage myopathy and Laing myopathy (
Meredith et al. 2004), but diaphragm weakness is not a common component of these myopathies.
Myl3 mutations are rare in dilated and hypertrophic cardiomyopathies (
Richard et al. 2003;
Fokstuen et al. 2008;
Møller et al. (2009), whereas mutations in
Myh4 have not been reported in cardiac disease. We were not able to find any link between mutations in these genes and diaphragmatic defects.
Increased expression of
Tpm3 and
Myh7, both components of the slow-twitch fibers that predominate in the diaphragm, has previously been noted in wildtype diaphragms compared to wildtype hindlimbs in mice (
Porter et al. 2004). As these genes are expressed in the muscular component of the developing diaphragm, it is plausible that the differences in expression seen for the sarcomere genes could be related to the muscular hypoplasia of the diaphragm per se, independent from
Fgfrl1 expression, or the differences in expression could be related the muscular hypoplasia of the diaphragm in addition to the absence of
Fgfrl1 expression. For example, the reduced expression of
Tpm3 and the other sarcomere genes in the
Fgfrl1 null diaphragms may reflect reduced myoblast adhesion and secondary diaphragm thinning due to the postulated role of
Fgfrl1 in cell adhesion. However, at a minimum, our results demonstrate that the expression of these genes is important for the late stages of murine diaphragm development in utero.
Little is known regarding the functions of
Lrtm1 and
Hephl1, the remaining two differentially expressed genes from our array studies. Leucine-rich transmembrane protein 1 (
Lrtm1) has three coding exons that produce a 356 amino acid membrane protein of unknown function (
Kim et al. 2007). Transcripts have been isolated from many cDNAs from the pineal gland, lung and eye (Unigene;
http://www.ncbi.nlm.nih.gov/unigene), and there does not seem to be an obvious relationship with muscle hypoplasia. A BlastP analysis (
http://blast.ncbi.nlm.nih.gov/Blast.cgi) revealed that Lrtm1 shows strong homology to Slit3 (1e-19) and Slit2 (4e-16), suggesting that Lrtm1 and Slit3 could share a similar function. This is interesting, as
Slit3 homozygous null mice have diaphragmatic hernias resulting from a failure of the septum transversum, or central tendon, of the diaphragm, to separate from the liver (
Yuan et al. 2003). Slit3 was localized to the plasma membrane and mitochondria of a macrophage cell line, RAW264.7, and inhibition of
Slit3 reduced cell motility and the activity of two small GTPases known to be involved in actin cytoskeleton organization,
Cdc42 and
Rac (
Begum et al. 2004;
Tanno et al. 2007). Rac1 and Cdc42 have been shown to be essential for myoblast fusion in the diaphragm in conditional mutant mice (
Vasyutina et al. 2009). Further studies are needed to investigate whether
Lrtm1 is also involved in cell motility or myoblast fusion in the diaphragm and to define its relationship with
Fgfrl1. The Slit proteins have also been shown be important in axon guidance and cell migration (
Ma and Tessier-Lavigne 2007). The protein similarity between Lrtm1 and Slit3 is intriguing, as it suggests that there may be more than one gene with a role in cell migration that is important in diaphragm development, perhaps acting by providing guidance cues for muscle cell or muscle nerve cell migration, if Lrtm1 does indeed share a similar function to Slit3.
The last gene,
Hephl1, has greatest similarity to
Hephaestin (
Heph; OMIM 300167), a gene encoding a transmembrane protein with ferroxidase activity that can function as an iron transporter in intestinal cells (
Vulpe et al. 1999;
Griffiths et al. 2005).
Heph is mutated in sex-linked anemia (sla) mice that have a microcytic, hypochromic anemia with deficient export of iron from enterocytes (
Petrak and Vyoral 2005), but Heph shows no known link with human disease and we were unable to find any connection between Heph and diaphragm development.
We chose to re-sequence
FGFRL1 in patients with CDH because of reports of two independent mouse models of loss of
Fgfrl1 gene function that had muscular hypoplasia of distinct, distal regions of the diaphragm (
Baertschi et al. 2007;
Catela et al. 2009) and the finding of
FGFRL1 as a candidate gene for CDH associated with 4p deletions. Our re-sequencing studies did not identify any new sequence alterations in
FGFRL1 in CDH patients, nor did they demonstrate a significant difference in the frequency of SNPs between control and patient population groups (). Despite the relatively small number of genotyped patients, our interpretation of the results is that the CDH patients studied were representative of this population and that further sequencing would be unlikely to uncover any significant differences from normal controls in terms of the allele frequencies for these SNPs in
FGFRL1. In addition, several of the SNPs have relatively high heterozygosity scores (for example, see p.Pro20Pro and p.Pro362Gln; ), thus making biological significance possibly less likely due to the relatively common frequency of carriers. However, one CDH patient who had a 4p16.3 deletion that included
FGFRL1 was hemizygous for the minor allele for three of the six SNPs, raising the possibility that the CDH in this patient may have been related to reduced
FGFRL1 expression caused by his deletion and these SNPs, or by the deletion of
FGFRL1 and at least one other gene at 4p16.3.
Sequencing of
FGFRL1 has previously been performed in 55 patients with congenital skeletal malformations comprising 22 patients with non-syndromic craniosynostosis and 19 with short stature (
Rieckmann et al. 2008). In a female with a clinical diagnosis of Antley-Bixler syndrome (craniosynostosis, radio-ulnar synostosis and genital anomalies), a frameshift mutation in exon 6 of
FGFRL1 caused by an insertion of four nucleotides was identified (
Rieckmann et al. 2008). The mutant and wildtype protein both localized to the plasma membrane, but the mutant protein lacked a sorting signal and was not removed from the plasma membrane as quickly, enabling interactions with FGF ligands of a longer duration and thus suggesting a gain of function, in contrast to Fgfrl1 wildtype protein, which was removed from the plasma membrane more quickly (
Rieckmann et al. 2008). This patient was also re-sequenced for P450 oxidoreductase mutations, a gene known to cause Antley–Bixler syndrome, and a missense mutation, p. Ala287Pro, and a novel mutation at the 3′ splice junction of exon 8 changing the donor sequence AG
gtacca into AG
atacca (
Rieckmann et al. 2008) were identified. There are no other published re-sequencing studies of the
FGFRL1 gene. Our failure to detect mutations or significant polymorphisms in
FGFRL1 in CDH patients could be construed as further evidence that diaphragm hypoplasia can be a pathologically and genetically distinct process from CDH.