Advancing human disease research with fish evolutionary mutant models.

Model organism research is essential to understand disease mechanisms. However, laboratory-induced genetic models can lack genetic variation and often fail to mimic the spectrum of disease severity. Evolutionary mutant models (EMMs) are species with evolved phenotypes that mimic human disease. EMMs complement traditional laboratory models by providing unique avenues to study gene-by-environment interactions, modular mutations in noncoding regions, and their evolved compensations. EMMs have improved our understanding of complex diseases, including cancer, diabetes, and aging, and illuminated mechanisms in many organs. Rapid advancements of sequencing and genome-editing technologies have catapulted the utility of EMMs, particularly in fish. Fish are the most diverse group of vertebrates, exhibiting a kaleidoscope of specialized phenotypes, many that would be pathogenic in humans but are adaptive in the species' specialized habitat. Importantly, evolved compensations can suggest avenues for novel disease therapies. This review summarizes current research using fish EMMs to advance our understanding of human disease.

Aldh2 is a lineage-specific metabolic gatekeeper in melanocyte stem cells.

Melanocyte stem cells (McSCs) in zebrafish serve as an on-demand source of melanocytes during growth and regeneration, but metabolic programs associated with their activation and regenerative processes are not well known. Here, using live imaging coupled with scRNA-sequencing, we discovered that, during regeneration, quiescent McSCs activate a dormant embryonic neural crest transcriptional program followed by an aldehyde dehydrogenase (Aldh) 2 metabolic switch to generate progeny. Unexpectedly, although ALDH2 is well known for its aldehyde-clearing mechanisms, we find that, in regenerating McSCs, Aldh2 activity is required to generate formate - the one-carbon (1C) building block for nucleotide biosynthesis - through formaldehyde metabolism. Consequently, we find that disrupting the 1C cycle with low doses of methotrexate causes melanocyte regeneration defects. In the absence of Aldh2, we find that purines are the metabolic end product sufficient for activated McSCs to generate progeny. Together, our work reveals McSCs undergo a two-step cell state transition during regeneration, and that the reaction products of Aldh2 enzymes have tissue-specific stem cell functions that meet metabolic demands in regeneration.

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis.

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when the effects of homozygosity for recessive mutations are studied in non-isogenic backgrounds, genes located proximal to the mutation on the same chromosome often appear over-represented among those genes identified as differentially expressed (DE). One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect functional responses to the mutation but, instead, result from an unequal distribution of expression quantitative trait loci (eQTLs) between sample groups of mutant or wild-type genotypes. This is problematic because eQTL expression differences are difficult to distinguish from genes that are DE due to functional responses to a mutation. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) are also observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between those sample groups subjected to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting-based approaches. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. Our method for DAR analysis requires only RNA-sequencing data, facilitating its application across new and existing datasets.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae.

BACKGROUND: The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. RESULTS: We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary. CONCLUSIONS: Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Cold-Driven Hemoglobin Evolution in Antarctic Notothenioid Fishes Prior to Hemoglobin Gene Loss in White-Blooded Icefishes.

Expression of multiple hemoglobin isoforms with differing physiochemical properties likely helps species adapt to different environmental and physiological conditions. Antarctic notothenioid fishes inhabit the icy Southern Ocean and display fewer hemoglobin isoforms, each with less affinity for oxygen than temperate relatives. Reduced hemoglobin multiplicity was proposed to result from relaxed selective pressure in the cold, thermally stable, and highly oxygenated Antarctic waters. These conditions also permitted the survival and diversification of white-blooded icefishes, the only vertebrates living without hemoglobin. To understand hemoglobin evolution during adaptation to freezing water, we analyzed hemoglobin genes from 36 notothenioid genome assemblies. Results showed that adaptation to frigid conditions shaped hemoglobin gene evolution by episodic diversifying selection concomitant with cold adaptation and by pervasive evolution in Antarctic notothenioids compared to temperate relatives, likely a continuing adaptation to Antarctic conditions. Analysis of hemoglobin gene expression in adult hematopoietic organs in various temperate and Antarctic species further revealed a switch in hemoglobin gene expression underlying hemoglobin multiplicity reduction in Antarctic fish, leading to a single hemoglobin isoform in adult plunderfishes and dragonfishes, the sister groups to icefishes. The predicted high hemoglobin multiplicity in Antarctic fish embryos based on transcriptomic data, however, raises questions about the molecular bases and physiological implications of diverse hemoglobin isoforms in embryos compared to adults. This analysis supports the hypothesis that the last common icefish ancestor was vulnerable to detrimental mutations affecting the single ancestral expressed alpha- and beta-globin gene pair, potentially predisposing their subsequent loss.

A complex genetic architecture in zebrafish relatives Danio quagga and D. kyathit underlies development of stripes and spots.

Vertebrate pigmentation is a fundamentally important, multifaceted phenotype. Zebrafish, Danio rerio, has been a valuable model for understanding genetics and development of pigment pattern formation due to its genetic and experimental tractability, advantages that are shared across several Danio species having a striking array of pigment patterns. Here, we use the sister species D. quagga and D. kyathit, with stripes and spots, respectively, to understand how natural genetic variation impacts phenotypes at cellular and organismal levels. We first show that D. quagga and D. kyathit phenotypes resemble those of wild-type D. rerio and several single locus mutants of D. rerio, respectively, in a morphospace defined by pattern variation along dorsoventral and anteroposterior axes. We then identify differences in patterning at the cellular level between D. quagga and D. kyathit by repeated daily imaging during pattern development and quantitative comparisons of adult phenotypes, revealing that patterns are similar initially but diverge ontogenetically. To assess the genetic architecture of these differences, we employ reduced-representation sequencing of second-generation hybrids. Despite the similarity of D. quagga to D. rerio, and D. kyathit to some D. rerio mutants, our analyses reveal a complex genetic basis for differences between D. quagga and D. kyathit, with several quantitative trait loci contributing to variation in overall pattern and cellular phenotypes, epistatic interactions between loci, and abundant segregating variation within species. Our findings provide a window into the evolutionary genetics of pattern-forming mechanisms in Danio and highlight the complexity of differences that can arise even between sister species. Further studies of natural genetic diversity underlying pattern variation in D. quagga and D. kyathit should provide insights complementary to those from zebrafish mutant phenotypes and more distant species comparisons.

A Y-linked anti-Müllerian hormone type-II receptor is the sex-determining gene in ayu, Plecoglossus altivelis.

Whole-genome duplication and genome compaction are thought to have played important roles in teleost fish evolution. Ayu (or sweetfish), Plecoglossus altivelis, belongs to the superorder Stomiati, order Osmeriformes. Stomiati is phylogenetically classified as sister taxa of Neoteleostei. Thus, ayu holds an important position in the fish tree of life. Although ayu is economically important for the food industry and recreational fishing in Japan, few genomic resources are available for this species. To address this problem, we produced a draft genome sequence of ayu by whole-genome shotgun sequencing and constructed linkage maps using a genotyping-by-sequencing approach. Syntenic analyses of ayu and other teleost fish provided information about chromosomal rearrangements during the divergence of Stomiati, Protacanthopterygii and Neoteleostei. The size of the ayu genome indicates that genome compaction occurred after the divergence of the family Osmeridae. Ayu has an XX/XY sex-determination system for which we identified sex-associated loci by a genome-wide association study by genotyping-by-sequencing and whole-genome resequencing using wild populations. Genome-wide association mapping using wild ayu populations revealed three sex-linked scaffolds (total, 2.03 Mb). Comparison of whole-genome resequencing mapping coverage between males and females identified male-specific regions in sex-linked scaffolds. A duplicate copy of the anti-Müllerian hormone type-II receptor gene (amhr2bY) was found within these male-specific regions, distinct from the autosomal copy of amhr2. Expression of the Y-linked amhr2 gene was male-specific in sox9b-positive somatic cells surrounding germ cells in undifferentiated gonads, whereas autosomal amhr2 transcripts were detected in somatic cells in sexually undifferentiated gonads of both genetic males and females. Loss-of-function mutation for amhr2bY induced male to female sex reversal. Taken together with the known role of Amh and Amhr2 in sex differentiation, these results indicate that the paralog of amhr2 on the ayu Y chromosome determines genetic sex, and the male-specific amh-amhr2 pathway is critical for testicular differentiation in ayu.

FishmiRNA: An Evolutionarily Supported MicroRNA Annotation and Expression Database for Ray-Finned Fishes.

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression involved in countless biological processes and are widely studied across metazoans. Although miRNA research continues to grow, the large community of fish miRNA researchers lacks exhaustive resources consistent among species. To fill this gap, we developed FishmiRNA, an evolutionarily supported miRNA annotation and expression database for ray-finned fishes: www.fishmirna.org. The self-explanatory database contains detailed, manually curated miRNA annotations with orthology relationships rigorously established by sequence similarity and conserved syntenies, and expression data provided for each detected mature miRNA. In just few clicks, users can download the annotation and expression database in several convenient formats either in its entirety or a subset. Simple filters and Blast search options also permit the simultaneous exploration and visual comparison of expression data for up to any ten mature miRNAs across species and organs. FishmiRNA was specifically designed for ease of use to reach a wide audience.

Coordinated patterning of zebrafish caudal fin symmetry by a central and two peripheral organizers.

BACKGROUND: Caudal fin symmetry characterizes teleosts and likely contributes to their evolutionary success. However, the coordinated development and patterning of skeletal elements establishing external symmetry remains incompletely understood. We explore the spatiotemporal emergence of caudal skeletal elements in zebrafish to consider evolutionary and developmental origins of caudal fin symmetry. RESULTS: Transgenic reporters and skeletal staining reveal that the hypural diastema-defining gap between hypurals 2 and 3 forms early and separates progenitors of two plates of connective tissue. Two sets of central principal rays (CPRs) synchronously, sequentially, and symmetrically emerge around the diastema. The two dorsal- and ventral-most rays (peripheral principal rays, PPRs) arise independently and earlier than adjacent CPRs. Muscle and tendon markers reveal that different muscles attach to CPR and PPR sets. CONCLUSIONS: We propose that caudal fin symmetry originates from a central organizer that establishes the hypural diastema and bidirectionally patterns surrounding tissue into two plates of connective tissue and two mirrored sets of CPRs. Further, two peripheral organizers unidirectionally specify PPRs, forming a symmetric "composite" fin derived from three fields. Distinct CPR and PPR ontogenies may represent developmental modules conferring ray identities, muscle connections, and biomechanical properties. Our model contextualizes mechanistic studies of teleost fin morphological variation.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs.

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3'-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.

Evolution of the nitric oxide synthase family in vertebrates and novel insights in gill development.

Nitric oxide (NO) is an ancestral key signalling molecule essential for life and has enormous versatility in biological systems, including cardiovascular homeostasis, neurotransmission and immunity. Although our knowledge of NO synthases (Nos), the enzymes that synthesize NO in vivo, is substantial, the origin of a large and diversified repertoire of nos gene orthologues in fishes with respect to tetrapods remains a puzzle. The recent identification of nos3 in the ray-finned fish spotted gar, which was considered lost in this lineage, changed this perspective. This finding prompted us to explore nos gene evolution, surveying vertebrate species representing key evolutionary nodes. This study provides noteworthy findings: first, nos2 experienced several lineage-specific gene duplications and losses. Second, nos3 was found to be lost independently in two different teleost lineages, Elopomorpha and Clupeocephala. Third, the expression of at least one nos paralogue in the gills of developing shark, bichir, sturgeon, and gar, but not in lamprey, suggests that nos expression in this organ may have arisen in the last common ancestor of gnathostomes. These results provide a framework for continuing research on nos genes' roles, highlighting subfunctionalization and reciprocal loss of function that occurred in different lineages during vertebrate genome duplications.

Bone microstructure and bone mineral density are not systemically different in Antarctic icefishes and related Antarctic notothenioids.

Ancestors of the Antarctic icefishes (family Channichthyidae) were benthic and had no swim bladder, making it energetically expensive to rise from the ocean floor. To exploit the water column, benthopelagic icefishes were hypothesized to have evolved a skeleton with "reduced bone," which gross anatomical data supported. Here, we tested the hypothesis that changes to icefish bones also occurred below the level of gross anatomy. Histology and micro-CT imaging of representative craniofacial bones (i.e., ceratohyal, frontal, dentary, and articular) of extant Antarctic fish species specifically evaluated two features that might cause the appearance of "reduced bone": bone microstructure (e.g., bone volume fraction and structure linear density) and bone mineral density (BMD, or mass of mineral per volume of bone). Measures of bone microstructure were not consistently different in bones from the icefishes Chaenocephalus aceratus and Champsocephalus gunnari, compared to the related benthic notothenioids Notothenia coriiceps and Gobionotothen gibberifrons. Some quantitative measures, such as bone volume fraction and structure linear density, were significantly increased in some icefish bones compared to homologous bones of non-icefish. However, such differences were rare, and no microstructural measures were consistently different in icefishes across all bones and species analyzed. Furthermore, BMD was similar among homologous bones of icefish and non-icefish Antarctic notothenioids. In summary, "reduced bone" in icefishes was not due to systemic changes in bone microstructure or BMD, raising the prospect that "reduced bone" in icefish occurs only at the gross anatomic level (i.e., smaller or fewer bones). Given that icefishes exhibit delayed skeletal development compared to non-icefish Antarctic fishes, combining these phenotypic data with genomic data might clarify genetic changes driving skeletal heterochrony.

Environmentally-induced sex reversal in fish with chromosomal vs. polygenic sex determination.

Sex ratio depends on sex determination mechanisms and is a key demographic parameter determining population viability and resilience to natural and anthropogenic stressors. There is increasing evidence that the environment can alter sex ratio even in genetically sex-determined species (GSD), as elevated temperature can cause female-to-male sex reversal (neomales). Alarmingly, neomales are being discovered in natural populations of several fish, amphibian and reptile species worldwide. Understanding the basis of neomale development is important for conservation biology. Among GSD species, it is unknown whether those with chromosomal sex determination (CSD), the most common system, will better resist the influence of high temperature than those with polygenic sex determination (PSD). Here, we compared the effects of elevated temperature in two wild zebrafish strains, Nadia (NA) and Ekkwill (EKW), which have CSD with a ZZ/ZW system, against the AB laboratory strain, which has PSD. First, we uncovered novel sex genotypes and the results showed that, at control temperature, the masculinization rate roughly doubled with the addition of each Z chromosome, while some ZW and WW fish of the wild strains became neomales. Surprisingly, we found that at elevated temperatures WW fish were just as likely as ZW fish to become neomales and that all strains were equally susceptible to masculinization. These results demonstrate that the Z chromosome is not essential for male development and that the dose of W buffers masculinization at the control temperature but not at elevated temperature. Furthermore, at the elevated temperature the testes of neomales, but not of normal males, contained more spermatozoa than at the control temperature. Our results show in an unprecedented way that, in a global warming scenario, CSD species may not necessarily be better protected against the masculinizing effect of elevated temperature than PSD species, and reveal genotype-by-temperature interactions in male sex determination and spermatogenesis.

Circulating miRNA repertoire as a biomarker of metabolic and reproductive states in rainbow trout.

BACKGROUND: Circulating miRNAs (c-miRNAs) are found in most, if not all, biological fluids and are becoming well-established non-invasive biomarkers of many human pathologies. However, their features in non-pathological contexts and whether their expression profiles reflect normal life history events have received little attention, especially in non-mammalian species. The aim of the present study was to investigate the potential of c-miRNAs to serve as biomarkers of reproductive and metabolic states in fish. RESULTS: The blood plasma was sampled throughout the reproductive cycle of female rainbow trout subjected to two different feeding regimes that triggered contrasting metabolic states. In addition, ovarian fluid was sampled at ovulation, and all samples were subjected to small RNA-seq analysis, leading to the establishment of a comprehensive miRNA repertoire (i.e., miRNAome) and enabling subsequent comparative analyses to a panel of RNA-seq libraries from a wide variety of tissues and organs. We showed that biological fluid miRNAomes are complex and encompass a high proportion of the overall rainbow trout miRNAome. While sharing a high proportion of common miRNAs, the blood plasma and ovarian fluid miRNAomes exhibited strong fluid-specific signatures. We further revealed that the blood plasma miRNAome significantly changed depending on metabolic and reproductive states. We subsequently identified three evolutionarily conserved muscle-specific miRNAs or myomiRs (miR-1-1/2-3p, miR-133a-1/2-3p, and miR-206-3p) that accumulated in the blood plasma in response to high feeding rates, making these myomiRs strong candidate biomarkers of active myogenesis. We also identified miR-202-5p as a candidate biomarker for reproductive success that could be used to predict ovulation and/or egg quality. CONCLUSIONS: Together, these promising results reveal the high potential of c-miRNAs, including evolutionarily conserved myomiRs, as physiologically relevant biomarker candidates and pave the way for the use of c-miRNAs for non-invasive phenotyping in various fish species.

Effect of Genetic Diagnosis on Patients with Previously Undiagnosed Disease.

BACKGROUND: Many patients remain without a diagnosis despite extensive medical evaluation. The Undiagnosed Diseases Network (UDN) was established to apply a multidisciplinary model in the evaluation of the most challenging cases and to identify the biologic characteristics of newly discovered diseases. The UDN, which is funded by the National Institutes of Health, was formed in 2014 as a network of seven clinical sites, two sequencing cores, and a coordinating center. Later, a central biorepository, a metabolomics core, and a model organisms screening center were added. METHODS: We evaluated patients who were referred to the UDN over a period of 20 months. The patients were required to have an undiagnosed condition despite thorough evaluation by a health care provider. We determined the rate of diagnosis among patients who subsequently had a complete evaluation, and we observed the effect of diagnosis on medical care. RESULTS: A total of 1519 patients (53% female) were referred to the UDN, of whom 601 (40%) were accepted for evaluation. Of the accepted patients, 192 (32%) had previously undergone exome sequencing. Symptoms were neurologic in 40% of the applicants, musculoskeletal in 10%, immunologic in 7%, gastrointestinal in 7%, and rheumatologic in 6%. Of the 382 patients who had a complete evaluation, 132 received a diagnosis, yielding a rate of diagnosis of 35%. A total of 15 diagnoses (11%) were made by clinical review alone, and 98 (74%) were made by exome or genome sequencing. Of the diagnoses, 21% led to recommendations regarding changes in therapy, 37% led to changes in diagnostic testing, and 36% led to variant-specific genetic counseling. We defined 31 new syndromes. CONCLUSIONS: The UDN established a diagnosis in 132 of the 382 patients who had a complete evaluation, yielding a rate of diagnosis of 35%. (Funded by the National Institutes of Health Common Fund.).

Unification of miRNA and isomiR research: the mirGFF3 format and the mirtop API.

MOTIVATION: MicroRNAs (miRNAs) are small RNA molecules (∼22 nucleotide long) involved in post-transcriptional gene regulation. Advances in high-throughput sequencing technologies led to the discovery of isomiRs, which are miRNA sequence variants. While many miRNA-seq analysis tools exist, the diversity of output formats hinders accurate comparisons between tools and precludes data sharing and the development of common downstream analysis methods. RESULTS: To overcome this situation, we present here a community-based project, miRNA Transcriptomic Open Project (miRTOP) working towards the optimization of miRNA analyses. The aim of miRTOP is to promote the development of downstream isomiR analysis tools that are compatible with existing detection and quantification tools. Based on the existing GFF3 format, we first created a new standard format, mirGFF3, for the output of miRNA/isomiR detection and quantification results from small RNA-seq data. Additionally, we developed a command line Python tool, mirtop, to create and manage the mirGFF3 format. Currently, mirtop can convert into mirGFF3 the outputs of commonly used pipelines, such as seqbuster, isomiR-SEA, sRNAbench, Prost! as well as BAM files. Some tools have also incorporated the mirGFF3 format directly into their code, such as, miRge2.0, IsoMIRmap and OptimiR. Its open architecture enables any tool or pipeline to output or convert results into mirGFF3. Collectively, this isomiR categorization system, along with the accompanying mirGFF3 and mirtop API, provide a comprehensive solution for the standardization of miRNA and isomiR annotation, enabling data sharing, reporting, comparative analyses and benchmarking, while promoting the development of common miRNA methods focusing on downstream steps of miRNA detection, annotation and quantification. AVAILABILITY AND IMPLEMENTATION: https://github.com/miRTop/mirGFF3/ and https://github.com/miRTop/mirtop. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Heterozygous loss-of-function variants significantly expand the phenotypes associated with loss of GDF11.

PURPOSE: Growth differentiation factor 11 (GDF11) is a key signaling protein required for proper development of many organ systems. Only one prior study has associated an inherited GDF11 variant with a dominant human disease in a family with variable craniofacial and vertebral abnormalities. Here, we expand the phenotypic spectrum associated with GDF11 variants and document the nature of the variants. METHODS: We present a cohort of six probands with de novo and inherited nonsense/frameshift (4/6 patients) and missense (2/6) variants in GDF11. We generated gdf11 mutant zebrafish to model loss of gdf11 phenotypes and used an overexpression screen in Drosophila to test variant functionality. RESULTS: Patients with variants in GDF11 presented with craniofacial (5/6), vertebral (5/6), neurological (6/6), visual (4/6), cardiac (3/6), auditory (3/6), and connective tissue abnormalities (3/6). gdf11 mutant zebrafish show craniofacial abnormalities and body segmentation defects that match some patient phenotypes. Expression of the patients' variants in the fly showed that one nonsense variant in GDF11 is a severe loss-of-function (LOF) allele whereas the missense variants in our cohort are partial LOF variants. CONCLUSION: GDF11 is needed for human development, particularly neuronal development, and LOF GDF11 alleles can affect the development of numerous organs and tissues.

Transgene-mediated skeletal phenotypic variation in zebrafish.

When considering relationships between genotype and phenotype we frequently ignore the fact that the genome of a typical animal, notably including that of a fish and a human, harbours a huge amount of foreign DNA. Such DNA, in the form of transposable elements, can affect genome function in a major way, and transgene biology needs to be included in our understanding of the genome. Here we examine an unexpected phenotypic effect of the chromosomally integrated transgene fli1a-F-hsp70l:Gal4VP16 that serves as a model for transgene function generally. We examine larval fras1 mutant zebrafish (Danio rerio). Gal4VP16 is a potent transcriptional activator that is already well known for toxicity and mediating unusual transcriptional effects. In the presence of the transgene, phenotypes in the neural crest-derived craniofacial skeleton, notably fusions and shape changes associated with loss of function fras1 mutations, are made more severe, as we quantify by scoring phenotypic penetrance, the fraction of mutants expressing the trait. A very interesting feature is that the enhancements are highly specific for fras1 mutant phenotypes, occurring in the apparent absence of more widespread changes. Except for the features due to the fras1 mutation, the transgene-bearing larvae appear generally healthy and to be developing normally. The transgene behaves as a genetic partial dominant: a single copy is sufficient for the enhancements, yet, for some traits, two copies may exert a stronger effect. We made new strains bearing independent insertions of the fli1a-F-hsp70l:Gal4VP16 transgene in new locations in the genome, and observed increased severities of the same phenotypes as observed for the original insertion. This finding suggests that sequences within the transgene, for example Gal4VP16, are responsible for the enhancements, rather than the effect on neighbouring host sequences (such as an insertional mutation). The specificity and biological action underlying the traits are subjects of considerable interest for further investigation, as we discuss. Our findings show that work with transgenes needs to be undertaken with caution and attention to detail.

Epigenetic factors Dnmt1 and Uhrf1 coordinate intestinal development.

Intestinal tract development is a coordinated process involving signaling among the progenitors and developing cells from all three germ layers. Development of endoderm-derived intestinal epithelium has been shown to depend on epigenetic modifications, but whether that is also the case for intestinal tract cell types from other germ layers remains unclear. We found that functional loss of a DNA methylation machinery component, ubiquitin-like protein containing PHD and RING finger domains 1 (uhrf1), leads to reduced numbers of ectoderm-derived enteric neurons and severe disruption of mesoderm-derived intestinal smooth muscle. Genetic chimeras revealed that Uhrf1 functions both cell-autonomously in enteric neuron precursors and cell-non-autonomously in surrounding intestinal cells, consistent with what is known about signaling interactions between these cell types that promote one another's development. Uhrf1 recruits the DNA methyltransferase Dnmt1 to unmethylated DNA during replication. Dnmt1 is also expressed in enteric neurons and smooth muscle progenitors. dnmt1 mutants have fewer enteric neurons and disrupted intestinal smooth muscle compared to wildtypes. Because dnmt1;uhrf1 double mutants have a similar phenotype to dnmt1 and uhrf1 single mutants, Dnmt1 and Uhrf1 must function together during enteric neuron and intestinal muscle development. This work shows that genes controlling epigenetic modifications are important to coordinate intestinal tract development, provides the first demonstration that these genes influence development of the ENS, and advances uhrf1 and dnmt1 as potential new Hirschsprung disease candidates.

Multiple independent chromosomal fusions accompanied the radiation of the Antarctic teleost genus Trematomus (Notothenioidei:Nototheniidae).

BACKGROUND: Chromosomal rearrangements are thought to be an important driving force underlying lineage diversification, but their link to speciation continues to be debated. Antarctic teleost fish of the family Nototheniidae (Notothenioidei) diversified in a changing environmental context, which led to ecological, morphological, and genetic differentiation among populations. In addition, extensive chromosomal repatterning accompanied species divergence in several clades. The most striking karyotypic changes involved the recent species radiation (about 10 My) of the genus Trematomus, with chromosomal pair numbers ranging between 29 and 12. These dramatic reductions in chromosome number resulted mostly from large-scale chromosome fusions. Multiple centric and/or tandem fusions have been hypothesized in at least seven of the twelve recognized Trematomus species. To reconstruct their evolutionary history, we employed comparative cytogenomics (BAC-FISH and chromosome painting) to reveal patterns of interspecific chromosomal orthologies across several notothenioid clades. RESULTS: We defined orthologous chromosomal segments of reference, termed Structural Units (SUs). SUs were identified in a total of 18 notothenioid species. We demonstrated for the first time that SUs were strongly conserved across every specimen examined, with chromosomal syntenies highlighting a paucity of intrachromosomal macro-rearrangements. Multiple independent fusions of these SUs were inferred in the Trematomus species, in contrast to the shared SU fusions in species of the sister lineage Notothenia. CONCLUSIONS: The SU segments were defined units of chromosomal rearrangement in the entire family Nototheiidae, which diverged from the other notothenioid families 20 My ago. Some of the identified chromosomal syntenies within the SUs were even conserved in their closest relatives, the family Eleginopsidae. Comparing the timing of acquisition of the fusions in the closely related genera Notothenia and Trematomus of the nototheniid species family, we conclude that they exhibit distinct chromosomal evolutionary histories, which may be relevant to different speciation scenarios.