FORGE Canada Consortium: Outcomes of a 2-Year National Rare-Disease Gene-Discovery Project

Inherited monogenic disease has an enormous impact on the well-being of children and their families. Over half of the children living with one of these conditions are without a molecular diagnosis because of the rarity of the disease, the marked clinical heterogeneity, and the reality that there are thousands of rare diseases for which causative mutations have yet to be identified. It is in this context that in 2010 a Canadian consortium was formed to rapidly identify mutations causing a wide spectrum of pediatric-onset rare diseases by using whole-exome sequencing. The FORGE (Finding of Rare Disease Genes) Canada Consortium brought together clinicians and scientists from 21 genetics centers and three science and technology innovation centers from across Canada. From nation-wide requests for proposals, 264 disorders were selected for study from the 371 submitted; disease-causing variants (including in 67 genes not previously associated with human disease; 41 of these have been genetically or functionally validated, and 26 are currently under study) were identified for 146 disorders over a 2-year period. Here, we present our experience with four strategies employed for gene discovery and discuss FORGE’s impact in a number of realms, from clinical diagnostics to the broadening of the phenotypic spectrum of many diseases to the biological insight gained into both disease states and normal human development. Lastly, on the basis of this experience, we discuss the way forward for rare-disease genetic discovery both in Canada and internationally.     

Susceptibility of apple clearwing moth larvae,Synanthedon myopaeformis (Lepidoptera: Sesiidae) to Beauveria bassiana and Metarhizium brunneum

Apple clearwing moth larvae, Synanthedon myopaeformis (Lepidoptera: Sesiidae) were found to be susceptible to infection by two entomopathogenic fungi: an indigenous fungus isolated from S. myopaeformis cadavers and identified as Metarhizium brunneum (Petch); and Beauveria bassiana isolate GHA. In laboratory bioassays, larvae exhibited dose related mortality after exposure to both the M. brunneum and Beauveria bassiana with 7 day LC₅₀'s of 2.9×10⁵ and 3.4×10⁵ spores/mL, respectively. Larval mortalities caused by the two isolates at 1×10⁶ spores/mL were not significantly different and 73% of the M. brunneum-treated, and 76% of the B. bassiana-treated larvae were dead 7 days post treatment, with LT₅₀'s of 5.5 and 5.1 days, respectively.     

Mutation in the Gene Encoding Ubiquitin Ligase LRSAM1 in Patients with Charcot-Marie-Tooth Disease

Charcot-Marie-Tooth disease (CMT) represents a family of related sensorimotor neuropathies. We studied a large family from a rural eastern Canadian community, with multiple individuals suffering from a condition clinically most similar to autosomal recessive axonal CMT, or AR-CMT2. Homozygosity mapping with high-density SNP genotyping of six affected individuals from the family excluded 23 known genes for various subtypes of CMT and instead identified a single homozygous region on chromosome 9, at 122,423,730–129,841,977 Mbp, shared identical by state in all six affected individuals. A homozygous pathogenic variant was identified in the gene encoding leucine rich repeat and sterile alpha motif 1 (LRSAM1) by direct DNA sequencing of genes within the region in affected DNA samples. The single nucleotide change mutates an intronic consensus acceptor splicing site from AG to AA. Direct analysis of RNA from patient blood demonstrated aberrant splicing of the affected exon, causing an obligatory frameshift and premature truncation of the protein. Western blotting of immortalized cells from a homozygous patient showed complete absence of detectable protein, consistent with the splice site defect. LRSAM1 plays a role in membrane vesicle fusion during viral maturation and for proper adhesion of neuronal cells in culture. Other ubiquitin ligases play documented roles in neurodegenerative diseases. LRSAM1 is a strong candidate for the causal gene for the genetic disorder in our kindred.     

Colletotrichum gloeosporioides s.l. associated with Theobroma cacao and other plants in Panama: multilocus phylogenies distinguish host-associated pathogens from asymptomatic endophytes

Colletotrichum interacts with numerous plant species overtly as symptomatic pathogens and cryptically as asymptomatic endophytes. It is not known whether these contrasting ecological modes are optional strategies expressed by individual Colletotrichum species or whether a species' ecology is explicitly pathogenic or endophytic. We explored this question by inferring relationships among 77 C. gloeosporioides s.l. strains isolated from asymptomatic leaves and from anthracnose lesions on leaves and fruits of Theobroma cacao (cacao) and other plants from Panamá. ITS and 5'-tef1 were used to assess diversity and to delineate operational taxonomic units for multilocus phylogenetic analysis. The ITS and 5'-tef1 screens concordantly resolved four strongly supported lineages, clades A-D: Clade A includes the ex type of C. gloeosporioides, clade B includes the ex type ITS sequence of C. boninense, and clades C and D are unidentified. The ITS yielded limited resolution and support within all clades, in particular the C. gloeosporioides clade (A), the focal lineage dealt with in this study. In contrast the 5'-tef1 screen differentiated nine distinctive haplotype subgroups within the C. gloeosporioides clade that were concordant with phylogenetic terminals resolved in a five-locus nuclear phylogeny. Among these were two phylogenetic species associated with symptomatic infections specific to either cacao or mango and five phylogenetic species isolated principally as asymptomatic infections from cacao and other plant hosts. We formally describe two new species, C. tropicale and C. ignotum, that are frequent asymptomatic associates of cacao and other Neotropical plant species, and epitypify C. theobromicola, which is associated with foliar and fruit anthracnose lesions of cacao. Asymptomatic Colletotrichum strains isolated from cacao plants grown in China included six distinct C. gloeosporioides clade taxa, only one of which is known to occur in the Neotropics.     

Biphasic perimicrovillar membrane production following feeding by previously starved Dysdercus peruvianus (Hemiptera: Pyrrhocoridae)

The development of perimicrovillar membranes (PMM) from midgut cells of starved and fed Dysdercus peruvianus was studied by using scanning electron microscopy (SEM), transmission electron microscopy (TEM) and assays for specific enzymatic markers of the perimicrovillar membranes (alpha-glucosidase), perimicrovillar space (aminopeptidase) and microvillar membranes (beta-glucosidase). High activities of these enzymes were observed 6h post-feeding and significant production of membranes was observed at 30 h post-feeding. In the gut cells of starved insects, the rough endoplasmic reticulum was organized in concentric bundles, with a greater number of mitochondria in the cellular apex. The presence of electron dense double-membrane vesicles and the production of PMM were not observed in this condition. Thirty hours post-feeding, a disorganization of the rough endoplasmic reticulum was observed, and it was possible to see double-membrane vesicles close to the cell apex. The membrane system formation was evident with a significant development of PMM in the midgut lumen. The luminal surface of the midgut during starvation and up to 48 h post-feeding was monitored using SEM. It was demonstrated that in the starved condition, the PMM was virtually absent from gut cells, except at the base of the microvilli. At 6h post-feeding, the microvilli were already completely covered with PMM, but with a maximum of PMM formation seen at 30 h post-feeding. Signals of PMM degradation were observed 48 h after pulse feeding.     

Normal levels of wild-type mitochondrial DNA maintain cytochrome c oxidase activity for two pathogenic mitochondrial DNA mutations but not for m.3243A-->G

Mitochondrial DNA (mtDNA) mutations are a common cause of human disease and accumulate as part of normal ageing and in common neurodegenerative disorders. Cells express a biochemical defect only when the proportion of mutated mtDNA exceeds a critical threshold, but it is not clear whether the actual cause of this defect is a loss of wild-type mtDNA, an excess of mutated mtDNA, or a combination of the two. Here, we show that segments of human skeletal muscle fibers harboring two pathogenic mtDNA mutations retain normal cytochrome c oxidase (COX) activity by maintaining a minimum amount of wild-type mtDNA. For these mutations, direct measurements of mutated and wild-type mtDNA molecules within the same skeletal muscle fiber are consistent with the "maintenance of wild type" hypothesis, which predicts that there is nonselective proliferation of mutated and wild-type mtDNA in response to the molecular defect. However, for the m.3243A-->G mutation, a superabundance of wild-type mtDNA was found in many muscle-fiber sections with negligible COX activity, indicating that the pathogenic mechanism for this particular mutation involves interference with the function of the wild-type mtDNA or wild-type gene products.     

Shifting distributions and speciation: species divergence during rapid climate change

Questions about how shifting distributions contribute to species diversification remain virtually without answer, even though rapid climate change during the Pleistocene clearly impacted genetic variation within many species. One factor that has prevented this question from being adequately addressed is the lack of precision associated with estimates of species divergence made from a single genetic locus and without incorporating processes that are biologically important as populations diverge. Analysis of DNA sequences from multiple variable loci in a coalescent framework that (i) corrects for gene divergence pre-dating speciation, and (ii) derives divergence-time estimates without making a priori assumptions about the processes underlying patterns of incomplete lineage sorting between species (i.e. allows for the possibility of gene flow during speciation), is critical to overcoming the inherent logistical and analytical difficulties of inferring the timing and mode of speciation during the dynamic Pleistocene. Estimates of species divergence that ignore these processes, use single locus data, or do both can dramatically overestimate species divergence. For example, using a coalescent approach with data from six loci, the divergence between two species of montane Melanoplus grasshoppers is estimated at between 200,000 and 300,000 years before present, far more recently than divergence estimates made using single-locus data or without the incorporation of population-level processes. Melanoplus grasshoppers radiated in the sky islands of the Rocky Mountains, and the analysis of divergence between these species suggests that the isolation of populations in multiple glacial refugia was an important factor in promoting speciation. Furthermore, the low estimates of gene flow between the species indicate that reproductive isolation must have evolved rapidly for the incipient species boundaries to be maintained through the subsequent glacial periods and shifts in species distributions.     

Tests of phenotypic and genetic concordance and their application to the conservation of Panamanian golden frogs (Anura, Bufonidae)

Evolutionarily significant units (ESUs) differ in the extent to which they capture, or even consider, adaptive variation, and most such designations are based solely on neutral genetic differences that may not capture variation relevant to species' adaptabilities to changing environmental conditions. While concordant patterns of divergence among data sets (i.e. neutral and potentially non-neutral characters) can strengthen ESU designations, determining whether such criteria are met for highly variable taxa is especially challenging. This study tests whether previously defined ESUs for endangered Panamanian golden frogs (Atelopus varius and Atelopus zeteki) exhibit concordant variation among multiple phenotypic traits and mitochondrial DNA sequences, and the extent to which such divergence corresponds to environmental differences. Multivariate analyses identify phenotypic and genetic differentiation consistent with proposed ESUs and support the status of A. varius and A. zeteki as separate species. Moreover, the significant association detected between ESU co-membership and genetic similarity, which remained strong after removing the effect of geographic distance, also indicates that genetic differences are not simply due to isolation by distance. Two phenotypic characters (body size and the extent of dorsal black patterning) that differ among ESUs also co-vary with environmental differences, suggesting that to the extent that these phenotypic differences are heritable, variation may be associated with adaptive divergence. Lastly, discriminant function analyses show that the frogs can be correctly assigned to ESUs based on simultaneous analysis of multiple characters. The study confirms the merit of conserving the previously proposed golden frog ESUs as well as demonstrates the utility and feasibility of combined analyses of ecological, morphological and genetic variation in evaluating ESUs, especially for highly variable taxa.