Chromosomal mapping and candidate gene discovery of chicken developmental mutants and genome-wide variation analysis of MHC congenics

Orofacial Cleft and Mandibular Prognathism-Human Genetics and Animal Models

Many complex molecular interactions are involved in the process of craniofacial development. Consequently, the network is sensitive to genetic mutations that may result in congenital malformations of varying severity. The most common birth anomalies within the head and neck are orofacial clefts (OFCs) and prognathism. Orofacial clefts are disorders with a range of phenotypes such as the cleft of the lip with or without cleft palate and isolated form of cleft palate with unilateral and bilateral variations. They may occur as an isolated abnormality (nonsyndromic-NSCLP) or coexist with syndromic disorders. Another cause of malformations, prognathism or skeletal class III malocclusion, is characterized by the disproportionate overgrowth of the mandible with or without the hypoplasia of maxilla. Both syndromes may be caused by the presence of environmental factors, but the majority of them are hereditary. Several mutations are linked to those phenotypes. In this review, we summarize the current knowledge regarding the genetics of those phenotypes and describe genotype-phenotype correlations. We then present the animal models used to study these defects.

A Premature Stop Codon in RAF1 Is the Priority Candidate Causative Mutation of the Inherited Chicken Wingless-2 Developmental Syndrome

The chicken wingless-2 (wg-2) mutation is inherited in an autosomal recessive fashion, and the resulting phenotype in mutant (wg-2/wg-2) individuals is a developmental syndrome characterized by absent wings, truncated legs, craniofacial as well as skin and feather defects, and kidney malformations. Mapping and genotyping established that the mutation resides within 227 kilobases (kb) of chromosome 12 in a wg-2 congenic inbred line. A capture array was designed to target and sequence the candidate region along with flanking DNA in 24 birds from the line. Many point mutations and insertions or deletions were identified, and analysis of the linked variants indicated a point mutation predicted to cause a premature stop codon in the RAF1 gene. Expression studies were conducted inclusive of all genes in the candidate region. Interestingly, RAF1 transcription was elevated, yet the protein was absent in the mutants relative to normal individuals. RAF1 encodes a protein integral to the Ras/Raf/MAPK signaling pathway controlling cellular proliferation, and notably, human RASopathies are developmental syndromes caused by germline mutations in genes of this pathway. Our work indicates RAF1 as the priority candidate causative gene for wg-2 and provides a new animal model to study an important signaling pathway implicated in limb development, as well as RASopathies.

Two Proximally Close Priority Candidate Genes for diplopodia-1, an Autosomal Inherited Craniofacial-Limb Syndrome in the Chicken: MRE11 and GPR83

Next-generation sequencing (NGS) and expression technologies were utilized to investigate the genes and sequence elements in a 586 kb region of chicken chromosome 1 associated with the autosomal recessive diplopodia-1 (dp-1) mutation. This mutation shows a syndromic phenotype similar to known human developmental abnormalities (e.g., cleft palate, polydactyly, omphalocele [exposed viscera]). Toward our goal to ascertain the variant responsible, the entire 586 kb region was sequenced following utilization of a specifically designed capture array and to confirm/validate fine-mapping results. Bioinformatic analyses identified a total of 6142 sequence variants, which included SNPs, indels, and gaps. Of these, 778 SNPs, 146 micro-indels, and 581 gaps were unique to the UCD-Dp-1.003 inbred congenic line; those found within exons and splice sites were studied for contribution to the mutant phenotype. Upon further validation with additional mutant samples, a smaller subset (of variants [51]) remains linked to the mutation. Additionally, utilization of specific samples in the NGS technology was advantageous in that fine-mapping methodologies eliminated an additional 326 kb of sequence information on chromosome 1. Predicted and confirmed protein-coding genes within the smaller 260 kb region were assessed for their developmental expression patterns over several stages of early embryogenesis in regions/tissues of interest (e.g., digits, craniofacial region). Based on these results and known function in other vertebrates, 2 genes within 5 kb of each other, MRE11 and GPR83, are proposed as high-priority candidates for the dp-1 mutation.

Narrowing the wingless-2 mutation to a 227 kb candidate region on chicken chromosome 12

Wingless-2 (wg-2) is an autosomal recessive mutation in chicken that results in an embryonic lethal condition. Affected individuals exhibit a multisystem syndrome characterized by absent wings, truncated legs, and craniofacial, kidney, and feather malformations. Previously, work focused on phenotype description, establishing the autosomal recessive pattern of Mendelian inheritance and placing the mutation on an inbred genetic background to create the congenic line UCD Wingless-2.331. The research described in this paper employed the complementary tools of breeding, genetics, and genomics to map the chromosomal location of the mutation and successively narrow the size of the region for analysis of the causative element. Specifically, the wg-2 mutation was initially mapped to a 7 Mb region of chromosome 12 using an Illumina 3 K SNP array. Subsequent SNP genotyping and exon sequencing combined with analysis from improved genome assemblies narrowed the region of interest to a maximum size of 227 kb. Within this region, 3 validated and 3 predicted candidate genes are found, and these are described. The wg-2 mutation is a valuable resource to contribute to an improved understanding of the developmental pathways involved in chicken and avian limb development as well as serving as a model for human development, as the resulting syndrome shares features with human congenital disorders.

Utilizing the chicken as an animal model for human craniofacial ciliopathies

The chicken has been a particularly useful model for the study of craniofacial development and disease for over a century due to their relatively large size, accessibility, and amenability for classical bead implantation and transplant experiments. Several naturally occurring mutant lines with craniofacial anomalies also exist and have been heavily utilized by developmental biologist for several decades. Two of the most well known lines, talpid(2) (ta(2)) and talpid(3) (ta(3)), represent the first spontaneous mutants to have the causative genes identified. Despite having distinct genetic causes, both mutants have recently been identified as ciliopathic. Excitingly, both of these mutants have been classified as models for human craniofacial ciliopathies: Oral-facial-digital syndrome (ta(2)) and Joubert syndrome (ta(3)). Herein, we review and compare these two models of craniofacial disease and highlight what they have revealed about the molecular and cellular etiology of ciliopathies. Furthermore, we outline how applying classical avian experiments and new technological advances (transgenics and genome editing) with naturally occurring avian mutants can add a tremendous amount to what we currently know about craniofacial ciliopathies.

Sp6 and Sp8 transcription factors control AER formation and dorsal-ventral patterning in limb development

The formation and maintenance of the apical ectodermal ridge (AER) is critical for the outgrowth and patterning of the vertebrate limb. The induction of the AER is a complex process that relies on integrated interactions among the Fgf, Wnt, and Bmp signaling pathways that operate within the ectoderm and between the ectoderm and the mesoderm of the early limb bud. The transcription factors Sp6 and Sp8 are expressed in the limb ectoderm and AER during limb development. Sp6 mutant mice display a mild syndactyly phenotype while Sp8 mutants exhibit severe limb truncations. Both mutants show defects in AER maturation and in dorsal-ventral patterning. To gain further insights into the role Sp6 and Sp8 play in limb development, we have produced mice lacking both Sp6 and Sp8 activity in the limb ectoderm. Remarkably, the elimination or significant reduction in Sp6;Sp8 gene dosage leads to tetra-amelia; initial budding occurs, but neither Fgf8 nor En1 are activated. Mutants bearing a single functional allele of Sp8 (Sp6^−/−;Sp8^+/−) exhibit a split-hand/foot malformation phenotype with double dorsal digit tips probably due to an irregular and immature AER that is not maintained in the center of the bud and on the abnormal expansion of Wnt7a expression to the ventral ectoderm. Our data are compatible with Sp6 and Sp8 working together and in a dose-dependent manner as indispensable mediators of Wnt/βcatenin and Bmp signaling in the limb ectoderm. We suggest that the function of these factors links proximal-distal and dorsal-ventral patterning.

In this report we examined the functional roles of Sp6 and Sp8 during limb development using compound loss-of-function mutants. Sp6 and Sp8, two members of the Sp gene family, are expressed in the limb bud ectoderm and function downstream of WNT/βcatenin signaling for Fgf8 induction. The analysis of the allelic series shows that the progressive reduction in the dose of Sp6 and Sp8 gene products leads to predictable morphology, from syndactyly, to split hand/foot malformation, oligodactyly, truncation and finally amelia, indicating that these two factors act in a complementary manner. The molecular characterization of the mutant limbs reveal that Sp6/Sp8 are required in a dose-dependent manner for Fgf8 and En1 induction, thereby placing them as an important link between the induction of the AER and the establishment of dorsal-ventral patterning during limb development.

The cellular and molecular etiology of the craniofacial defects in the avian ciliopathic mutant talpid2

talpid(2) is an avian autosomal recessive mutant with a myriad of congenital malformations, including polydactyly and facial clefting. Although phenotypically similar to talpid(3), talpid(2) has a distinct facial phenotype and an unknown cellular, molecular and genetic basis. We set out to determine the etiology of the craniofacial phenotype of this mutant. We confirmed that primary cilia were disrupted in talpid(2) mutants. Molecularly, we found disruptions in Hedgehog signaling. Post-translational processing of GLI2 and GLI3 was aberrant in the developing facial prominences. Although both GLI2 and GLI3 processing were disrupted in talpid(2) mutants, only GLI3 activator levels were significantly altered in the nucleus. Through additional fine mapping and whole-genome sequencing, we determined that the talpid(2) phenotype was linked to a 1.4 Mb region on GGA1q that contained the gene encoding the ciliary protein C2CD3. We cloned the avian ortholog of C2CD3 and found its expression was ubiquitous, but most robust in the developing limbs and facial prominences. Furthermore, we found that C2CD3 is localized proximal to the ciliary axoneme and is important for docking the mother centriole to the ciliary vesicle and cell membrane. Finally, we identified a 19 bp deletion in talpid(2) C2CD3 that produces a premature stop codon, and thus a truncated protein, as the likely causal allele for the phenotype. Together, these data provide insight into the cellular, molecular and genetic etiology of the talpid(2) phenotype. Our data suggest that, although the talpid(2) and talpid(3) mutations affect a common ciliogenesis pathway, they are caused by mutations in different ciliary proteins that result in differences in craniofacial phenotype.

Characterization of insertional variation of porcine endogenous retroviruses in six different pig breeds

Pigs may need to be exploited as xenotransplantation donors due to the shortage of human organs, tissues and cells. Porcine endogenous retroviruses (PERVs) are a significant obstacle to xenotransplantation because they can infect human cells in vitro and have the potential for transmission of unexpected pathogens to humans. In this research, 101 pigs, including four commercial breeds (23 Berkshire, 13 Duroc, 22 Landrace and 14 Yorkshire pigs), one native breed (19 Korean native pigs) and one miniature breed (10 NIH miniature pigs) were used to investigate insertional variations for 11 PERV loci (three PERV-A, six PERV-B and two PERV-C). Over 60% of the pigs harbored one PERV-A (907F8) integration and five PERV-B (B3-3G, B3-7G, 742H1, 1155D9 and 465D1) integrations. However, two PERV-A loci (A1-6C and 1347C1) and one PERV-B locus (B3-7F) were absent in Duroc pigs. Moreover, two PERV-C loci (C2-6C and C4-2G) only existed in Korean native pigs and NIH miniature pigs. The results suggest that PERV insertional variations differ among pig breeds as well as among individuals within a breed. Also, the results presented here can be used for the selection of animals that do not have specific PERV integration for xenotransplantation research.

Defining the Sequence Elements and Candidate Genes for the Coloboma Mutation

The chicken coloboma mutation exhibits features similar to human congenital developmental malformations such as ocular coloboma, cleft-palate, dwarfism, and polydactyly. The coloboma-associated region and encoded genes were investigated using advanced genomic, genetic, and gene expression technologies. Initially, the mutation was linked to a 990 kb region encoding 11 genes; the application of the genetic and genomic tools led to a reduction of the linked region to 176 kb and the elimination of 7 genes. Furthermore, bioinformatics analyses of capture array-next generation sequence data identified genetic elements including SNPs, insertions, deletions, gaps, chromosomal rearrangements, and miRNA binding sites within the introgressed causative region relative to the reference genome sequence. Coloboma-specific variants within exons, UTRs, and splice sites were studied for their contribution to the mutant phenotype. Our compiled results suggest three genes for future studies. The three candidate genes, SLC30A5 (a zinc transporter), CENPH (a centromere protein), and CDK7 (a cyclin-dependent kinase), are differentially expressed (compared to normal embryos) at stages and in tissues affected by the coloboma mutation. Of these genes, two (SLC30A5 and CENPH) are considered high-priority candidate based upon studies in other vertebrate model systems.

Case study of sequence capture enrichment technology: identification of variation underpinning developmental syndromes in an amniote model

Chicken developmental mutants are valuable for discovering sequences and pathways controlling amniote development. Herein we applied the advanced technologies of targeted sequence genomic capture enrichment and next-generation sequencing to discover the causative element for three inherited mutations affecting craniofacial, limb and/or organ development. Since the mutations (coloboma, diplopodia-1 and wingless-2) were bred into a congenic line series and previously mapped to different chromosomes, each targeted mutant causative region could be compared to that of the other two congenic partners, thereby providing internal controls on a single array. Of the ~73 million 50-bp sequence reads, ~76% were specific to the enriched targeted regions with an average target coverage of 132-fold. Analysis of the three targeted regions (2.06 Mb combined) identified line-specific single nucleotide polymorphism (SNPs) and micro (1–3 nt) indels. Sequence content for regions indicated as gaps in the reference genome was generated, thus contributing to its refinement. Additionally, Mauve alignments were constructed and indicated putative chromosomal rearrangements. This is the first report of targeted capture array technology in an avian species, the chicken, an important vertebrate model; the work highlights the utility of employing advanced technologies in an organism with only a “draft stage” reference genome sequence.

The expression of preaxial polydactyly is influenced by modifying genetic elements and is not maintained by chromosomal inversion in an avian biomedical model

Polydactyly (Po) is a common mutation found in many vertebrates. The UCD-Po.003 congenic chicken line was previously characterized for Po inheritance (autosomal dominant) and the mutation was mapped to chromosome 2p. Here, we describe for the first time the range and variability of the phenotype in this congenic line. Further, we studied the hypothesis that a chromosomal inversion was responsible for the maintenance of a large (6.3 Mb) candidate gene region. Fluorescence in situ hybridization employing BACs encompassing a 10.7-Mb region of GGA2p showed that the Po chromosome was normal, i.e. exhibits the wild-type BAC order. Continued fine-mapping along with a change in breeding strategy reduced the size of the causative region to 1.43 Mb. Recent research indicates that the cause of preaxial Po resides within a 794-bp highly conserved zone of polarizing activity regulatory sequence (ZRS) element located in intron 5 of the LMBR1 gene; however, the ZRS polymorphism of interest is found in some but not all breeds of polydactylous chicken. Therefore, we sequenced the ZRS in 101 heterozygous and 30 unaffected (wild-type) individuals to establish the relevance of this region to the Po condition in the UCD-Po.003 congenic line. A single point mutation (C/A at coordinate GGA2p: 8,414,121) within the ZRS segregated with carrier status. The polydactylous UCD-Silkie line also maintains this SNP in addition to a single base deletion. An inheritance analysis of the phenotypic variation in UCD-Po.003 suggests recessive epistasis as the mode of inheritance for the additional modifying genetic elements, residing outside the ZRS, to impact the preaxial polydactyl phenotype. These results contribute to our understanding of the cause of Po in an important vertebrate model.