- Research article
- Open Access
De novo ZIC2 frameshift variant associated with frontonasal dysplasia in a Limousin calf
BMC Genomics volume 22, Article number: 1 (2021)
Bovine frontonasal dysplasias like arhinencephaly, synophthalmia, cyclopia and anophthalmia are sporadic congenital facial malformations. In this study, computed tomography, necropsy, histopathological examinations and whole genome sequencing on an Illumina NextSeq500 were performed to characterize a stillborn Limousin calf with frontonasal dysplasia. In order to identify private genetic and structural variants, we screened whole genome sequencing data of the affected calf and unaffected relatives including parents, a maternal and paternal halfsibling.
The stillborn calf exhibited severe craniofacial malformations. Nose and maxilla were absent, mandibles were upwardly curved and a median cleft palate was evident. Eyes, optic nerve and orbital cavities were not developed and the rudimentary orbita showed hypotelorism. A defect centrally in the front skull covered with a membrane extended into the intracranial cavity. Aprosencephaly affected telencephalic and diencephalic structures and cerebellum. In addition, a shortened tail was seen. Filtering whole genome sequencing data revealed a private frameshift variant within the candidate gene ZIC2 in the affected calf. This variant was heterozygous mutant in this case and homozygous wild type in parents, half-siblings and controls.
We found a novel ZIC2 frameshift mutation in an aprosencephalic Limousin calf. The origin of this variant is most likely due to a de novo mutation in the germline of one parent or during very early embryonic development. To the authors’ best knowledge, this is the first identified mutation in cattle associated with bovine frontonasal dysplasia.
Frontonasal dysplasias (FND) comprise a heterogeneous group of disorders with congenital polymalformations caused by abnormal median facial development. The absence of the olfactory tract (arhinencephaly), partial fusion (synophthalmia), complete fusion (cyclopia) or the total absence (anopthalmia) of ocular globes are phenotypes of FND found in cattle . FND can also include cerebral defects as holoprosencephaly (HPE) characterized by incomplete cleavage of the prosencephalon into two hemispheres . Additional organ defects like heart anomalies may occur .
In human and mouse, several genetic variants are associated with HPE and FND. Besides chromosomal defects, critical FND candidate genes were identified [4,5,6]. In human, the spectrum of Aristaless-Like Homeoboxprotein (ALX) related FNDs involves recessively inherited loss-of-function mutations in the ALX1, ALX3, and ALX4 genes . In Burmese cats, a 12-bp in-frame deletion in the ALX1 gene is responsible for severe craniofacial malformations and FND . In ruminants, sporadic cases of different types of cyclopia and FND occur in cattle [1, 9,10,11,12,13,14,15,16,17], buffalo [18,19,20], sheep [21,22,23] and goats  (Table 1). Whole genome scans for FND were not yet performed in ruminants.
The objective of this study was to give a full phenotypic characterization of a sporadic FND case and filter out private genetic and structural variants for this case using whole genome sequencing (WGS) data of the affected calf, his parents and two unaffected half-siblings.
The stillborn male Limousin calf had a weight of 52 kg and the crest-rump length was 103 cm (40.55 in.). The head showed severe craniofacial alterations such as abnormal development of the maxillary processes and face skull (Fig. 1). Nasal openings were missing and the region of the upper jaw and muzzle seemed to be constricted and covered with haired skin. Eyes, eyelids or eyelashes were not discernible. The tongue was protruded. A skin defect covered by a membrane was located centrally in the front of the skull. A very short tail was obvious.
Computed tomography (CT)
In the CT scan, the skull showed severe abnormalities compared to a healthy animal (Fig. 2). The maxilla was bilaterally completely absent. The mandibles were strongly curved upwards. The right part of the mandible was smaller than the left one. Four pairs of incisive teeth were completely developed, the molars and premolars were rudimentary and curved. The nasal bones, forming the nose in normal animals (os nasale, os incisivum, os lacrimale, os ethmoidale, os presphenoidale, choana (os palatinum, os spehnoidale) and vomer) were absent. The os frontale had open structures. The orbits were rudimentary and placed on the front side near to the midline of the face instead of the lateral sides of the head. The palate had a wide and complete median cleft.
The head showed severe deformations and absence of facial structures. The nose and the maxilla were not developed and the mandibles were curved upwards. Eyes, optic nerv and orbital cavities were absent. Centrally on the frontal skull a defect, measuring 7 × 5 × 0.5 cm in width, extended into the intracranial cavity. A thin brownish membrane covered this bone defect. The whole skull was heavily deformed and appeared shortened and squat. Only parts of the brain stem were developed. Aprosencephaly affected telencephalic (cerebral cortex, hippocampus and striatum) and diencephalic structures (globus pallidus, thalamus, hypothalamus and hypophysis) and cerebellum. The malformations of the forebrain included an absence of the prosencephalon (including telencephalon and diencephalon) and a hypoplastic mesencephalon and rhombencephalon. In the oral cavity, incisors and canini were developed, but premolars were missing. A median cleft palate of the soft and hard palate due to a complete median disclosure was evident. The vertebral column consisted of 7 cervical, 13 thoracic, 6 lumbar, 5 sacral, and only 4 coccygeal vertebral bodies.
The heart showed severe complex malformations with a persisting foramen ovale measuring 2 × 2 cm, a completely absent cardiac septum and only one single ventricle. Furthermore, a complete fetal atelectasis was found.
Histologically, the membrane covering the central defect on the skull was composed of collagen-rich connective tissue. Directly adherent to this membrane, parts of meninges and neuronal tissue were seen. Remnants of brain tissue consisted of mainly white matter with single scattered neurons.
The Limousin calf was the only one FND case in this herd with purebred Limousin (Fig. 3). The actual herd size when the FND case occurred was 60 females in reproductive age. The sire of the present case was a natural service purebred Limousin bull, which sired all heifers and cows of this herd for two consecutive years. An autosomal recessive mode of inheritance, where both, the dam and the sire, are heterozygous mutant, is very unlikely due to the sporadic occurrence of this case and 120 progeny within 2 years from matings of the same bull in this herd. A paternal dominant germline mutation also seems unlikely because of the large number of progeny and only one case among 120 births. A maternal germline mutation or a novel mutation in early embryonic development may be considered as possible hypotheses.
Whole genome sequencing
We screened 755 filtered genetic and structural variants (Additional file 1: Table S1) which were filtered out under the hypothesis of a recessive inheritance (heterozygous in both parents and homozygous mutant in the case) or a de novo dominant germline mutation (heterozygous in the case and homozygous wild type in both parents and half-sibs) (Fig. 3). There were no variants, which fulfilled the condition of a recessive inheritance. Then we searched for candidate genes (Additional file 2: Table S2) among the filtered variants. The sole candidate gene with a variant identified in the case was ZIC2. A 1-bp deletion (g.80722845TC>T; ARS-UCD1.2:g.76742066TC>T) in exon 4 of the candidate gene ZIC2 (Fig. 4) was heterozygous in the affected calf and homozygous wild type in both parents and all other Limousin herdmates and private controls. This variant was not yet deposited in dbSNP. Due to this heterozygous 1-bp deletion, a new stop codon (TAA) was predicted for position ZIC2:g.80723235 (ARS-UCD1.2:g.76742456). The wild type stop codon (TGA) at position ZIC2:g.80723063 (ARS-UCD1.2:g.76742284) was removed.
The protein effect of this variant was predicted to be probably damaging (0.991) using PolyPhen-2. The sequence of the mutant protein has a likely amino acid (aa) exchange of prolin to arginine at position 543 (p.Pro543fs; ARS-UCD1.2:p.Pro454fs) followed by an altered sequence of 127 aa (Additional file 3: Figure S1). The mutant ZIC2 protein was predicted to have additional 56 aa in comparison to the wild type protein and instead of 614 aa (ARS-UCD1.2:525aa) an increase to 670 aa (ARS-UCD1.2:581aa). The frameshift mutation was predicted to cause a truncation of the Zinc Finger Protein Zic and Gli (PTHR19818) domain at codon 543 using ORFfinder.
Congenital FND, arhinencephaly and severe defects of the central nervous system were the main and characteristic findings in the present case. Similar brain malformations like HPE and aprosencephaly were previously reported in a Hereford calf  and a lamb . Contrasting to the present case, bone structures of the head were completely absent in this lamb. A Brown Swiss crossbreed calf showed very similar malformations of facial structures with the present case . Muzzle and nasal bones were completely absent and the mandibles were strongly curved upwards as well as a membraned opening to the intracranial cavity were common findings among our case and this Brown Swiss crossbreed FND calf. Ophthalmological malformations characterized by hypotelorism without indications of cyclopia or synophthalmia in the present case were also often seen in previously reported FND cases in cattle [1, 12,13,14, 16]. A cleft palate was an accompanying sign in a Japanese Brown calf , Murrah Buffalo calf  as well in our case. We observed a shortened tail in the affected Limousin calf investigated here, but could not confirm a complete absence of the tail as previously described in a crossbreed calf with mild FND .
In human with arhinencephaly, synophthalmia and HPE also heart defects like inter-ventricular communication orifice were found similar to the complex heart malformations seen in the present case of our Limousin calf .
The WGS data analysis of the present case retrieved a heterozygous frameshift variant in ZIC2, which is a strong candidate gene for HPE in human and mouse. We validated the wild type genotypes in the parents, two half-siblings and controls of the affected Limousin calf but were not able to get Sanger sequences from a further tissue sample of the affected calf. The reason Sanger sequencing failed is that DNA degraded because sampling could only be done from the insufficiently cooled body about 2 weeks after birth of the affected calf. In cases of human HPE associated with ZIC2 mutations, only mild to moderate facial malformations were observed [25,26,27]. Despite the mild facial phenotypes, these patients often exhibit severe HPE and neurological impairments [25,26,27]. However, mouse models revealed severe facial malformations in embryos like as hypotelorism, synophthalmy and cyclopia . In human, patients with a 13q deletion syndrome including ZIC2, showed diverse phenotypes with congenital defects including craniofacial dysmorphy, HPE, aprosencephaly and heart malformations . The subfamily PTHR19818:SF18 is part of the ‘Zinc Finger Protein ZIC and GLI’ protein domain family. Due to RNA polymerase II regulatory region sequence and transcription factor activity, the molecular function of this protein domain consists in creating the zinc finger transcription factor that regulates mRNA transcription . Combined regulative function of GLI and ZIC genes responds to inductive signals and induces patterned neural cell differentiation . The zinc-finger domains of the Zic and Gli factors bind to the identical target sequences as the transcription factors, which mediate hedgehog signals . We suggest that the modified ZIC2 protein domain in this case altered mRNA transcription regulation and the neural cell differentiation in early embryonic development.
We excluded an autosomal recessive mode of inheritance for the present FND phenotype. Only the affected calf had a heterozygous mutant genotype for the ZIC2 candidate variant and all family members including parents were homozygous wild type. A paternal or maternal germline mutation may be likely the source of the ZIC2 frameshift or even a novel mutation in very early embryonic development may lead to a post-zygotic heterozygous mutant genotype . For the present case, the hypothesis of a paternal germline mutation seems less likely based on the pattern of the ZIC2 frameshift mutation in the Limousin family and the frequency of FND in about 120 progeny. In human, pathogenic ZIC2 variants are largely de novo and in addition, more frequently de novo than any other HPE-associated genes . In addition, ZIC2 is one of the most commonly heterozygous mutated genes in human HPE patients [35, 36]. A large spectrum of mutations were reported in human and most of them were located in zinc finger domains [34,35,36]. Mutations within the N-terminal region resulted more often in the most severe HPE cases (alobar form), in which the brain is not divided into hemispheres. Semilobar forms, characterized by incomplete forebrain division, were more frequently associated with mutations within the C-terminal region .
In the present study, we identified a de novo frameshift variant in ZIC2 for lethal frontonasal dysplasia in a Limousin calf. To the authors’ best knowledge, this is the first report of an FND associated mutation in cattle identified through WGS.
The examined male purebred Limousin calf was stillborn in July 2017 after an extended duration of pregnancy by 22 days. We got this stillborn calf from the private owner and herdbook breeder of this Limousin herd. According to the report of the owner, the calf was the first case of any skeletal malformations seen in this cattle herd. The sire of the affected calf was a natural service Limousin bull, registered in the Limousin herd book. This bull sired all heifers and cows of this herd in the calving period when the case was observed. The dam and the sire of the calf were in healthy condition and had normal appearance. Further, all half-siblings of the affected calf were in healthy condition and had no congenital abnormalities. We collected tissue samples from the affected calf and EDTA-blood samples from the vena jugularis of the sire and dam of the affected calf, and further of two maternal half-siblings, 9 paternal half-siblings and three unrelated control animals. Pedigree data from the herd with the affected calf were collected and edited for analysis.
The stillborn affected calf was submitted to CT and a patho-morphological examination. CT scans were acquired with a multislice helical CT scanner (Brilliance TM CT 16 BigBore, Philips Healthcare, The Netherlands) using conventional settings (120 kV and 315 mA) and 0.8 mm slice thickness. Necropsy and histopathological examinations were performed for the affected calf with particular reference on the malformations of the head and brain. Samples of the membrane covering the defect on the front skull were taken, fixed in 10% formalin and embedded in paraffin. Sections of 4 μm were stained with hematoxylin and eosin and were microscopically examined.
Whole genome sequencing
To identify potential causal variants, genomic DNA from the affected calf, its dam and sire, one female maternal and one male paternal unaffected half-siblings were isolated using a chloroform extraction. Quality control of the isolated DNA was performed and ensuing libraries from these samples were prepared according to the manufacturers protocols using the NEBNext Ultra DNA Library Prep Kit for Illumina (New England BioLabs, Ipswich, MA, USA).WGS was performed on an Illumina NextSeq500 (Illumina, San Diego, CA, USA) in a 2 × 150 bp paired-end mode. Quality control of the WGS data was done using fastqc 0.11.5  and the reads were trimmed using PRINSEQ (V 0.20.4) . The remaining data were mapped to the bovine reference genome UMD 3.1. (Ensembl) using BWA 0.7.13 . Mean coverage was 8.38X with a mean error rate (mismatches/bases mapped) of 1.84 × 10− 2 and a mean mapping rate of 96.62%. Bam-files were sorted, indexed and marked for duplicates using SAMtools 1.3.1  and Picard tools (http://broadinstitute.github.io/picard/, version 2.9.4). Variants were called with GATK, version 4.0 , using Base Quality Score Recalibration (BQSR), Haplotype Caller and Variant Recalibrator. Variant calling was done for the five Limousin samples and further 89 WGS data from cattle of the breeds Holstein, Fleckvieh, Braunvieh, Vorderwald, German Angus, Galloway, Limousin, Charolais, Hereford, Tyrolean Grey and Miniature Zebu. The resulting vcf-file was used for screening private variants in the FND case. Variants selected for analysis had a read depth of 2–999 and quality score values > 20. We first screened the bovine ALX1 gene (ENSBTAG00000014977, UMD3.1) on BTA 514991697–15,013,533 bp (ARS-UCD1.2: 14917352–14,939,191 bp) from the vcf-file using SAS, version 9.4 (Statistical Analysis System, Cary, NC, USA) to identify all coding and noncoding variants of the affected calf. All further variants filtered had high or moderate effects according to prediction toolbox SNPEff version 4.3 q (2017-08-30, SNPEff database UMD3.1.86) . We filtered for variants, which were homozygous mutant in the affected Limousin calf and heterozygous mutant in the dam and sire and heterozygous or homozygous wild type in both half-siblings using SAS, version 9.4. In addition, we screened for variants, which were heterozygous mutant in the affected Limousin calf but homozygous wild in both parents and both half-siblings. All filtered variants had to be homozygous wild type in all 89 private controls and not yet deposited in dbSNP. Filtered variants were screened for candidate genes associated with the terms FND, arhinencephaly, holosprosencephaly and cyclopia for all mammalian species retrieved from National Center for Biotechnology Information (NCBI, www.ncbi.nlm.nih.gov) and Online Mendelian Inheritance in Animals (OMIA, www.omia.org, date of access: 12 April 2020). In addition, we mapped the NGS data to ARS-UCD1.2 and visualized the bam-files at critical positions for possible causative variants using Integrated Genomics Viewer (IGV) . In order to check for large deletions/insertions and duplications as well as chromosomal aberrations, we used the breakpoint prediction framework LUMPY . We compared all WGS data from 89 private controls with data from the affected calf and its relatives. Structural variants were filtered out which were heterozygous mutant in the affected calf and homozygous wild type in all other animals or which were homozygous mutant in the affected calf and heterozygous in parents and not homozygous mutant in its half-sibs. These structural variants were investigated for their potential functional effects by comparison of their genomic positions with those regions harboring candidate genes for HPE in human and animals according to National Center for Biotechnology Information (NCBI, www.ncbi.nlm.nih.gov).
To verify the previously investigated protein effect of the ZIC2 variant, we applied the PolyPhen-2 (Polymorphism Phenotyping v2) tool. To visualize predicted changes in the ZIC2 protein, we used ORFfinder (NCBI, www.ncbi.nlm.nih.gov/orffinder). The resulting amino acid (aa) sequences of the wild type protein (NSBTAP00000035289, UMD3.1) and predicted mutant protein were compared and aligned using ClustalW2, version 2.1 . Predicted protein domains were identified using InterPro .
Availability of data and materials
Variants were submitted to European Variation Archive (https://www.ebi.ac.uk/eva/?Home) referred to as PRJEB36774 (ARS-UCD1.2:g.76742066TC>T). WGS data of the calves, the sire, the dam and the maternal and paternal half-siblings were deposited in NCBI Sequence Read Archive under the project number PRJNA526664 (SAMN11107014, SAMN11107015, SAMN11107016, SAMN11107017, SAMN11107018). Further WGS data were retrieved from Sequence Read Archive (SRA, NCBI).
- ALX :
Base quality score recalibration
Bos taurus autosome
Ethylenediamine tetraacetic acid
- Gli :
Integrated Genomics Viewer
- OMIA :
Online Mendelian Inheritance in Animals
Polymerase chain reaction
Statistical Analysis System
Whole genome sequencing
- ZIC2 :
Zinc finger family member 2
Schulze U, Distl O. Case report. Arhinia and cyclopia in a German Fleckvieh calf. DTW. 2006;113(6):236–9.
Cohen MM. Holoprosencephaly. Clinical, anatomic, and molecular dimensions. Birth Defects Res A Clin Mol Teratol. 2006;76(9):658–73.
Poaty H, Gentien D, Reyes C, Silou J. Macroscopic analysis of fetus having Arhinencephaly, Synophthalmia and Holoprosencephaly. Human Genet Embryol. 2016;6:135.
Roessler E, Muenke M. The molecular genetics of holoprosencephaly. Am J Med Genet C: Semin Med Genet. 2010;1:52–61.
Dubourg C, Bendavid C, Pasquier L, Henry C, Odent S, David V. Holoprosencephaly. Orphanet J Rare Dis. 2007;2(1):8.
Farlie PG, Baker NL, Yap P, Tan TY. Frontonasal dysplasia: towards an understanding of molecular and developmental Aetiology. Mol Syndromol. 2016;7(6):312–21.
Uz E, Alanay Y, Aktas D, Vargel I, Gucer S, Tuncbilek G, et al. Disruption of ALX1 causes extreme microphthalmia and severe facial clefting: expanding the spectrum of autosomal-recessive ALX-related frontonasal dysplasia. Am J Hum Genet. 2010;86(5):789–96.
Lyons LA, Erdman CA, Grahn RA, Hamilton MJ, Carter MJ, Helps CR, et al. Aristaless-like Homeobox protein 1 (ALX1) variant associated with craniofacial structure and frontonasal dysplasia in Burmese cats. Dev Biol. 2016;409(2):451–8.
Cho D-Y, Zeman D, Miller J. Holoprosencephaly in a bovine calf. Acta Neuropathol. 1985;67(3):322–5.
Hammoda A, Abdoud M. Cyclopia in a Friesian calf. Egypt J Phytopathol. 1989;17(1):76–82.
Moritomo Y, Tsuda T, Miyamoto H. Craniofacial skeletal abnormalities in anomalous calves with clefts of the face. J Vet Med Sci. 1999;61(10):1147–52.
Kim B-S, Park H-J, Cho J-W, Berzina D, Zeeshan M, Ejaz S, et al. Cyclopia in a calf. Korean J Vet Serv. 2006;29(4):503–6.
ÖZcan K, Gürbulak K, Takci I, ÖZen H, Kacar C. Pancarci̇ M. atypical cyclopia in a brown swiss cross calf: a case report. Anat Histol Embryol. 2006;35(3):152–4.
Malik AA, Sheikh GG, Lone FA, Islam R, Khatun A. Delivery of Cyclopia monster fetus from a crossbred heifer suffering from dystocia. J Adv Vet Anim Res. 2013;3(4):167–8.
Osman HE, Shuaib Y, Ibrahim Z. Congenital deformities in the Sudan 5. Craniofacial deformities and coccygeal agenesis in anomalous bovine calf. Assiut vet. Med. J; 2013.
Nourani H, Karimi I, Vardanjani HR. Synophthalmia in a Holstein cross calf, vol. 333. Urmia: Veterinary Research Forum Faculty of Veterinary Medicine, Urmia University; 2014.
Agerholm JS, McEvoy FJ, Heegaard S, Charlier C, Jagannathan V, Drögemüller C. A de novo missense mutation of FGFR2 causes facial dysplasia syndrome in Holstein cattle. BMC Genet. 2017;18(1):74.
Pandey A, Saini G, Chander S, Chaudhary R, Jakhar P, Singh M, et al. Dystocia due to abnormal calf in a buffalo: a case report. Buffalo Bull. 2010;29(4):315–7.
Singh H, Gupta G, Jan M, Nabi S, Singh J, Dey S. Atypical cyclopia in a buffalo calf. Buff Bull. 2013;32:1.
Albarella S, Ciotola F, D’Anza E, Coletta A, Zicarelli L, Peretti V. Congenital malformations in river Buffalo (Bubalus bubalis). Animals. 2017;7(2):9.
Brachthäuser L, Klumpp S, Hecht W, Kuchelmeister K, Reinacher M, Ebbert W, et al. Aprosencephaly with otocephaly in a lamb (Ovis aries). Vet Pathol. 2012;49(6):1043–8.
Dantas F, Medeiros G, Figueiredo A, Thompson K, Riet-Correa F. Skeletal dysplasia with craniofacial deformity and disproportionate dwarfism in hair sheep of northeastern Brazil. J Comp Pathol. 2014;150(2):245–52.
Welch K, Panter K, Lee S, Gardner D, Stegelmeier B, Cook D. Cyclopamine induced synophthalmia in sheep: defining a critical window and toxicokinetic evaluation. J Appl Toxicol. 2009;29(5):414–21.
Rashed R, Al-kafafy M, Abdellah B, Sharshar A, Shoghy K, Erasha A, et al. Cyclopia of goat: micro and macroscopic, radiographic and computed tomographic studies. Alexandria J Vet Sci. 2014;42:1.
Brown LY, Odent S, David V, Blayau M, Dubourg C, Apacik C, et al. Holoprosencephaly due to mutations in ZIC2: alanine tract expansion mutations may be caused by parental somatic recombination. Hum Mol Genet. 2001;10(8):791–6.
Dubourg C, Lazaro L, Pasquier L, Bendavid C, Blayau M, Duff FL, et al. Molecular screening of SHH, ZIC2, SIX3, and TGIF genes in patients with features of holoprosencephaly spectrum: mutation review and genotype–phenotype correlations. Hum Mutat. 2004;24(1):43–51.
Xiong J, Xiang B, Chen X, Cai T. Case report: a novel mutation in ZIC2 in an infant with microcephaly, holoprosencephaly, and arachnoid cyst. Medicine. 2019;98:10.
Warr N, Powles-Glover N, Chappell A, Robson J, Norris D, Arkell RM. Zic2-associated holoprosencephaly is caused by a transient defect in the organizer region during gastrulation. Hum Mol Genet. 2008;17(19):2986–96.
Quélin C, Bendavid C, Dubourg C, De La Rochebrochard C, Lucas J, Henry C, et al. Twelve new patients with 13q deletion syndrome: genotype–phenotype analyses in progress. Eur J Med Genet. 2009;52(1):41–6.
Lawson MJ, Zhang LJN. Sexy gene conversions: locating gene conversions on the X-chromosome. Nucleic Acids Res. 2009;37(14):4570–9.
Brewster R, Lee J, Altaba ARJN. Gli/Zic factors pattern the neural plate by defining domains of cell differentiation. Nature. 1998;393(6685):579.
Mizugishi K, Aruga J, Nakata K, Mikoshiba K. Molecular properties of Zic proteins as transcriptional regulators and their relationship to GLI proteins. J Biol Chem. 2001;276(3):2180–8.
Acuna-Hidalgo R, Bo T, Kwint MP, Van De Vorst M, Pinelli M, Veltman JA, et al. Post-zygotic point mutations are an underrecognized source of de novo genomic variation. Am J Hum Genet. 2015;97(1):67–74.
Solomon BD, Lacbawan F, Mercier S, Clegg NJ, Delgado MR, Rosenbaum K, et al. Mutations in ZIC2 in human holoprosencephaly: description of a novel ZIC2 specific phenotype and comprehensive analysis of 157 individuals. J Med Genet. 2010;47(8):513–24.
Roessler E, Hu P, Marino J, Hong S, Hart R, Berger S, et al. Common genetic causes of holoprosencephaly are limited to a small set of evolutionarily conserved driver genes of midline development coordinated by TGF-β, hedgehog, and FGF signaling. Hum Mutat. 2018;39(10):1416–27.
Roessler E, Lacbawan F, Dubourg C, Paulussen A, Herbergs J, Hehr U, et al. The full spectrum of holoprosencephaly-associated mutations within the ZIC2 gene in humans predicts loss-of-function as the predominant disease mechanism. Hum Mutat. 2009;30(4):E541–54.
Andrews S. A quality control tool for high throughput sequence data. In: Reference Source. http://www.bioinformatics.babraham.ac.uk/projects/fastqc; 2010.
Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4.
Li H, Durbin R. Fast and accurate short read alignment with burrows–wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6(2):80–92.
Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24.
Layer RM, Chiang C, Quinlan AR, Hall IM. LUMPY: a probabilistic framework for structural variant discovery. Genome Biol. 2014;15(6):R84.
Larkin MA, Blackshields G, Brown N, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.
Mitchell AL, Attwood TK, Babbitt PC, Blum M, Bork P, Bridge A, et al. InterPro in 2019: improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 2019;47(D1):D351–60.
The authors are grateful to the Limousin breeder donating samples and providing data for this research project. We thank J. Wrede for his help in bioinformatics and M. Drabert, H. Klippert-Hasberg and N. Wagner for support in DNA preparation. We thank H. Klippert-Hasberg for expert technical assistance in library preparation and next generation sequencing. We gratefully acknowledge The North-German Supercomputing Alliance (HLRN) for providing high performance computing infrastructure and resources.
Not applicable. Open Access funding enabled and organized by Projekt DEAL.
Ethics approval and consent to participate
All animal work has been conducted according to the national and international guidelines for animal welfare. The present study had been specifically approved by the Institutional Animal Care and Use Committee of Lower Saxony, the State Veterinary Office Niedersächsisches Landesamt für Verbraucherschutz und Lebensmittelsicherheit, Oldenburg, Germany (registration number 33.9–42502-05-04A247). We obtained written informed approval from the owner of the Limousin cattle to use samples and photographs from his animals for current research, publication and further investigations. Samples of controls from other breeds were taken from our biobank and all the owners of these animals gave written informed consent to use these samples.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Filtering result of whole genome sequencing data revealed 755 variants. Only one variant (printed in bold) on BTA 12 within the critical candidate gene ZIC2, which was associated with holoprosencephaly in mammalian animals, exclusively occurred in the affected calf. The dam, the sire, one male paternal (a) and one female maternal half-sibling (b) were homozygous wild type, as well as further 89 controls of the breeds Holstein, Fleckvieh, Braunvieh, Vorderwald, German Angus, Galloway, Limousin, Charolais, Hereford, Tyrolean Grey and Miniature Zebu.
Candidate genes for frontonasal dysplasia (FND), arhinencephaly, holoprosencephaly and cyclopia in mammalian animals according to NCBI. The bovine orthologues gene were presented with the chromosomal position according to UMD3.1. The candidate gene ZIC2 which was filtered out of whole genome sequencing data is in bold.
Protein sequences of the wild type and mutant ZIC2 protein (NSBTAP00000035289) are compared. The amino acid frameshift of prolin to arginine at position 543 is red framed. The extension of the protein sequence over 127 amino acids for the mutant protein is given. Asterisks mark identical amino acids between the wild type and mutant sequences.
About this article
Cite this article
Braun, M., Lehmbecker, A., Eikelberg, D. et al. De novo ZIC2 frameshift variant associated with frontonasal dysplasia in a Limousin calf. BMC Genomics 22, 1 (2021). https://doi.org/10.1186/s12864-020-07350-y
- Stillborn: arhinencephaly
- de novo mutation