Genetic structure and origin of emu populations in Japanese farms inferred from large-scale SNP genotyping based on double-digest RAD-seq.

The emu is a novel poultry species in Japan. However, Japanese farmed emu populations have reduced genetic diversity owing to inbreeding. We have previously suggested that there are genetic resources in the Tohoku Safari Park (TSP) and Fuji/Kakegawa Kachoen Garden Park (FGP/KGP) to extend the genetic diversity of commercial emu farms based on microsatellite (SSR) and mitochondrial DNA. However, those markers provide relatively poor information. Thus, we investigated the genetic structure of farmed Japanese populations based on a large-scale genotyping system using RAD-seq and verified the usefulness of TSP and FGP/KGP as genetic resources for expanding genetic diversity. Admixture, phylogenetic, and principal component analyses based on 28,676 SNPs showed that TSP individuals were ancestors in the Okhotsk Emu Farm (OEF). FGP/KGP individuals showed a unique genetic component that differed from that of the others. We have previously reported that the mitochondrial haplotypes of FGP/KGP were shared with an isolated wild population in eastern Australia. These results suggest that FGP/KGP individuals originated from an eastern Australia isolated population different from other populations including ancestral of OEF/TSP. Our results would provide information for the development of Japanese emu farms and industry and for the conservation of genetic resources in the Australian wild emu.

Genetic relationships among emu populations in Japanese farms based on mitochondrial and microsatellite DNA polymorphisms.

Emus (Dromaius novaehollandiae) are expected to become a novel poultry species for producing eggs, meat, and oil. In our previous studies, Japanese emu populations were predicted to have reduced genetic diversity through inbreeding. For a sustainable emu industry in Japan, it is necessary to understand the current genetic structure and relationships in dispersed farms. In this study, we investigated the genetic structure and relationships of six Japanese emu farms based on mitochondrial DNA and microsatellite polymorphisms. We analyzed the DNA sequences of the mitochondrial D-loop region in 157 individuals and detected four haplotypes with four nucleotide substitution sites (Hap-a, Hap-b, Hap-c, and Hap-d). Analysis of molecular variance revealed that 43.6% of total variance was "among population," and the FST value was 0.436 with significant genetic differentiation (P < 0.001). In microsatellite analysis, the expected (HE ) and observed (HO ) heterozygosities were 0.53-0.64 and 0.44-0.59, respectively. Phylogenetic trees and STRUCTURE analysis revealed that the six Japanese farmed emu populations could be divided into four genetically differentiated groups. Therefore, we identified genetic resources that may be useful in extending the genetic diversity of Japanese farms and are predicted to contribute to the conservation and reconstruction of populations.

Gender-Difference in Hair Length as Revealed by Crispr-Based Production of Long-Haired Mice with Dysfunctional FGF5 Mutations.

Fibroblast growth factor 5 (FGF5) is an important molecule required for the transition from anagen to catagen phase of the mammalian hair cycle. We previously reported that Syrian hamsters harboring a 1-bp deletion in the Fgf5 gene exhibit excessive hair growth in males. Herein, we generated Fgf5 mutant mice using genome editing via oviductal nucleic acid delivery (GONAD)/improved GONAD (i-GONAD), an in vivo genome editing system used to target early embryos present in the oviductal lumen, to study gender differences in hair length in mutant mice. The two lines (Fgf5go-malc), one with a 2-bp deletion (c.552_553del) and the other with a 1-bp insertion (c.552_553insA) in exon 3 of Fgf5, were successfully established. Each mutation was predicted to disrupt a part of the FGF domain through frameshift mutation (p.Glu184ValfsX128 or p.Glu184ArgfsX128). Fgf5go-malc1 mice had heterogeneously distributed longer hairs than wild-type mice (C57BL/6J). Notably, this change was more evident in males than in females (p < 0.0001). Immunohistochemical analysis revealed the presence of FGF5 protein in the dermal papilla and outer root sheath of the hair follicles from C57BL/6J and Fgf5go-malc1 mice. Histological analysis revealed that the prolonged anagen phase might be the cause of accelerated hair growth in Fgf5go-malc1 mice.

Carcass traits and fat quality of breeding emu (Dromaius novaehollandiae) in Northern Japan.

Characterization of carcass traits and fat quality is important to effectively produce and genetically improve emus. We investigated carcass traits in 309 emus. The meat production of female emus showed a significantly higher value than that of males (P < 0.01). The fat weight of male (9.232 ± 3.156 kg) was larger than that of the female (7.772 ± 2.697 kg). The fat yield (fat weight per kg of body weight) was strongly correlated to body weight (r = 0.79 and r = 0.75 in male and female, respectively). The fat melting points of females and males were 19.19 ± 3.39°C and 19.39 ± 3.39°C, respectively, without significant difference. Since the fat melting point did not correlate to body and fat weights, we predicted that it was an independent trait from body growth and was highly influenced by genetic elements. Percentages of palmitic, stearic, oleic, linoleic, and α-linolenic acids were 22.27 ± 3.50%, 9.37 ± 1.90%, 54.11 ± 5.17%, 13.54 ± 7.80% and 0.71 ± 0.59%, respectively. Among them, linoleic acid contents showed a wide individual difference (range 0.3-19.9%). The oleic/stearic acid ratio showed a negative correlation to the fat melting point. These results suggest that the fat melting point is a good indicator of C18:1/C18:0 ratio in emu fat.

Development of 49 novel microsatellite markers from Next-generation sequencing data and a robust method for parentage tests in the emu (Dromaius novaehollandiae).

The emu is a useful and new breed of poultry, but their genetic improvement has not advanced yet due to their very recent domestication. Pedigree information is difficult to record because of their complex reproduction system (polyandry). To identify parent-offspring relationships in the emu, parentage test based on polymorphic DNA markers have to be developed. In this study, we isolated more than 25,000 microsatellite (simple sequence repeat, SSR) regions from Next-generation sequencing data via the QDD pipeline and developed 49 SSR markers with polymorphism in the Japanese farmed emu. The dinucleotide motifs, (AC)n, (AT)n and (AG)n, were the most frequently detected and were found on 10,167 (38.55%), 8,114 (30.76%) and 4,796 (18.18%) contigs, respectively. Forty-nine novel SSR markers were characterized in 20 individuals and showed NA ranged from 2 to 12, with an average of 4.2. HE/HO ranged from 0.389/0.071 to 0.702/1.000 with an average of 0.601/0.515. PIC value ranged from 0.059 to 0.886 with an average of 0.528, and 17 of 49 markers showed a higher polymorphism than 0.500. Thirty-four individuals were genotyped using 12 markers, and CERVUS simulations based on genotype showed that parents of all offspring were identified with 0.9995-1.0 probability. Thus, 49 novel SSR markers and a robust method for parentage test for the Japanese emu were developed.

Mitochondrial DNA variations in Japanese farmed emu populations.

The emu (Dromaius novaehollandiae) is a new poultry. In this study, we investigated the haplotype composition of mitochondrial DNA among emu populations farmed in Japan. We sequenced the D-loop region in 109 individuals, and detected four substitution sites and three haplotypes (Hap-a, -b, and -c). Hap-a was the most frequently observed haplotype in the Japanese populations. Although Hap-c was a rare haplotype in not only Japanese but also Australian populations, it was detected with high frequency in the Japanese farmed population. The AMOVA indicated that 9% of total variance was "among population". The FST value was 0.087 and genetic differentiation was significant (P<0.01). These results may contribute to conserving the genetic resources available for the Japanese emu industry.

c.753A>G genome editing of a Cdh23ahl allele delays age-related hearing loss and degeneration of cochlear hair cells in C57BL/6J mice.

C57BL/6J mice have long been studied as a model of age-related hearing loss (ARHL). In C57BL/6J mice, ARHL begins in the high-frequency range at 3 months of age and spreads toward low frequencies by 10 months of age. We previously confirmed that c.753A>G genome editing of an ahl allele (c.753A) in the cadherin 23 gene (Cdh23) suppressed the onset of ARHL until 12 months of age. We further investigated the hearing phenotypes of the original and genome-edited C57BL/6J-Cdh23+/+ (c.753G/G) mice until 24 months of age. The hearing tests revealed that most of the C57BL/6J mice maintained good hearing levels until 14 months of age following genome editing of a Cdh23ahl allele. However, the hearing levels of the C57BL/6J-Cdh23+/+ mice gradually declined, and severe ARHL developed with increasing age. ARHL in the C57BL/6J mice was correlated with degeneration of the stereocilia in cochlear hair cells. The stereocilia degeneration was rescued in the C57BL/6J-Cdh23+/+ mice at 12 months of age, but the stereocilia bundles exhibited abnormal phenotypes similar to those of the original C57BL/6J mice at more advanced ages. Therefore, genome editing of Cdh23ahl did not completely suppress ARHL in C57BL/6J mice. We also compared the hearing levels of C57BL/6J-Cdh23+/+ mice with those of C3H/HeN and MSM/Ms mice, which carry the Cdh23+ allele. The severity and onset patterns of ARHL in the C57BL/6J-Cdh23+/+ mice differed from those observed in other Cdh23+/+ mice. Therefore, we hypothesize that other susceptible and/or resistant alleles of ARHL exist in the genetic backgrounds of these mice.

Development and characterization of ten novel microsatellite loci for the emu (Dromaius novaehollandiae) and genetic diversity of Japanese farm populations.

The emu (Dromaius novaehollandiae) is a useful poultry animal farmed for fat, meat, and eggs. Genetic structure and relationships among farmed emu populations in Japan are unknown and the number of microsatellite markers for genetic analysis of the emu is insufficient. In this study, we isolated 16 microsatellites from the emu genome and developed ten new microsatellite markers. These microsatellite markers were used to characterize three farm emu populations in Japan. The number of alleles ranged from 3 to 13 and the expected (HE) and observed heterozygosity (HO) of these microsatellite loci was 0.187-0.802 and 0.179-0.647, respectively. The polymorphic information content ranged from 0.176 to 0.786. Positive inbreeding coefficient (FIS) values were detected in all tested populations, and they ranged from 0.027 to 0.540. These results suggest that farm populations of the emu in Japan resulted from inbreeding. The fixation index (FST) values ranged from 0.026 to 0.061, and phylogenetic trees and population structure analysis confirmed no definitive genetic differentiation among the three populations. Therefore, these populations are at a relatively low level of genetic differentiation at present. The microsatellite markers developed in our study can be utilized for genetic analysis and preservation of genetic resources in the emu.

Genetic modifiers of rodent animal models: the role in cataractogenesis.

Visual impairment leads to a decrease in quality of life. Cataract is the most commonly observed ocular disease in humans that causes vision disorders. The risk factors associated with cataract development include aging, infections, eye injuries, environmental causes, such as radiation and exposure to ultraviolet rays in sunlight, and genetic mutations. Additionally, several cataract patients display phenotypic heterogeneity, suggesting the role of genetic modifiers in the modulation of severity and onset time of cataractogenesis. However, the genetic modifiers associated with cataract have not been identified in humans yet. In contrast, the identification and mapping of genetic modifiers have been successfully carried out in mice and rats. In this review, we focus on the genetic modifiers of cataract in the rodent models.

OHC-TRECK: A Novel System Using a Mouse Model for Investigation of the Molecular Mechanisms Associated with Outer Hair Cell Death in the Inner Ear.

Outer hair cells (OHCs) are responsible for the amplification of sound, and the death of these cells leads to hearing loss. Although the mechanisms for sound amplification and OHC death have been well investigated, the effects on the cochlea after OHC death are poorly understood. To study the consequences of OHC death, we established an OHC knockout system using a novel mouse model, Prestin-hDTR, which uses the prestin promoter to express the human diphtheria toxin (DT) receptor gene (hDTR). Administration of DT to adult Prestin-hDTR mice results in the depletion of almost all OHCs without significant damage to other cochlear and vestibular cells, suggesting that this system is an effective tool for the analysis of how other cells in the cochlea and vestibula are affected after OHC death. To evaluate the changes in the cochlea after OHC death, we performed differential gene expression analysis between the untreated and DT-treated groups of wild-type and Prestin-hDTR mice. This analysis revealed that genes associated with inflammatory/immune responses were significantly upregulated. Moreover, we found that several genes linked to hearing loss were strongly downregulated by OHC death. Together, these results suggest that this OHC knockout system is a useful tool to identify biomarkers associated with OHC death.

i-GONAD: a robust method for in situ germline genome engineering using CRISPR nucleases.

We present a robust method called improved-Genome editing via Oviductal Nucleic Acids Delivery (i-GONAD) that delivers CRISPR ribonucleoproteins to E0.7 embryos via in situ electroporation. The method generates mouse models containing single-base changes, kilobase-sized deletions, and knock-ins. The efficiency of i-GONAD is comparable to that of traditional microinjection methods, which rely on ex vivo handling of zygotes and require recipient animals for embryo transfer. In contrast, i-GONAD avoids these technically difficult steps, and it can be performed at any laboratory with simple equipment and technical expertise. Further, i-GONAD-treated females retain reproductive function, suggesting future use of the method for germline gene therapy.

Pde6brd1 mutation modifies cataractogenesis in Foxe3rct mice.

The Foxe3rct mutation, which causes early-onset cataracts, is a recessive mutation found in SJL/J mice. A previous study reported that cataract phenotypes are modified by the genetic background of mouse inbred strains and that the Pde6brd1 mutation, which induced degeneration of the photoreceptor cells, is a strong candidate genetic modifier to accelerate the severity of cataractogenesis of Foxe3rct mice. We created congenic mice by transferring a genomic region including the Foxe3rct mutation to the B6 genetic background, which does not carry the Pde6brd1 mutation. In the congenic mice, the cataract phenotypes became remarkably mild, and the development of cataracts was suppressed for a long time. Moreover, we created transgenic mice by injecting BAC clones including the wild-type Pde6b gene into the eggs of SJL-Foxe3rct mice. Although the resistant effect for cataract phenotypes in transgenic mice was less than that in congenic mice, the severity and onset time of cataract phenotypes were clearly improved and delayed, respectively, compared with the phenotypes of the original SJL-Foxe3rct mice. These results clearly show that the development of early-onset cataracts requires at least two mutant alleles of Foxe3rct and Pde6brd1, and another modifier associated with the severity of cataract phenotypes in Foxe3rct mice underlies the genetic backgrounds in mice.

A novel splice site mutation of myosin VI in mice leads to stereociliary fusion caused by disruption of actin networks in the apical region of inner ear hair cells.

An unconventional myosin encoded by the myosin VI gene (MYO6) contributes to hearing loss in humans. Homozygous mutations of MYO6 result in nonsyndromic profound congenital hearing loss, DFNB37. Kumamoto shaker/waltzer (ksv) mice harbor spontaneous mutations, and homozygous mutants exhibit congenital defects in balance and hearing caused by fusion of the stereocilia. We identified a Myo6c.1381G>A mutation that was found to be a p.E461K mutation leading to alternative splicing errors in Myo6 mRNA in ksv mutants. An analysis of the mRNA and protein expression in animals harboring this mutation suggested that most of the abnormal alternatively spliced isoforms of MYO6 are degraded in ksv mice. In the hair cells of ksv/ksv homozygotes, the MYO6 protein levels were significantly decreased in the cytoplasm, including in the cuticular plates. MYO6 and stereociliary taper-specific proteins were mislocalized along the entire length of the stereocilia of ksv/ksv mice, thus suggesting that MYO6 attached to taper-specific proteins at the stereociliary base. Histological analysis of the cochlear hair cells showed that the stereociliary fusion in the ksv/ksv mutants, developed through fusion between stereociliary bundles, raised cuticular plate membranes in the cochlear hair cells and resulted in incorporation of the bundles into the sheaths of the cuticular plates. Interestingly, the expression of the stereociliary rootlet-specific TRIO and F-actin binding protein (TRIOBP) was altered in ksv/ksv mice. The abnormal expression of TRIOBP suggested that the rootlets in the hair cells of ksv/ksv mice had excessive growth. Hence, these data indicated that decreased MYO6 levels in ksv/ksv mutants disrupt actin networks in the apical region of hair cells, thereby maintaining the normal structure of the cuticular plates and rootlets, and additionally provided a cellular basis for stereociliary fusion in Myo6 mutants.

A novel missense mutation of Mip causes semi-dominant cataracts in the Nat mouse.

Major intrinsic protein of lens fiber (MIP) is one of the proteins essential for maintaining lens transparency while also contributing to dominant cataracts in humans. The Nodai cataract (Nat) mice harbor a spontaneous mutation in Mip and develop early-onset nuclear cataracts. The Nat mutation is a c.631G>A mutation (MipNat), resulting in a glycine-to-arginine substitution (p.Gly211Arg) in the sixth transmembrane domain. The MipNat/Nat homozygotes exhibit congenital cataracts caused by the degeneration of lens fiber cells. MIP normally localizes to the lens fiber cell membranes. However, the MipNat/Nat mice were found to lack an organelle-free zone, and the MIP was mislocalized to the nuclear membrane and perinuclear region. Furthermore, the MipNat/+ mice exhibited milder cataracts than MipNat/Nat mice due to the slight degeneration of the lens fiber cells. Although there were no differences in the localization of MIP to the membranes of lens fiber cells in MipNat/+ mice compared to that in wild-type mice, the protein levels of MIP were significantly reduced in the eyes. These findings suggest that cataractogenesis in MipNat mutants are caused by defects in MIP expression. Overall, the MipNat mice offer a novel model to better understand the phenotypes and mechanisms for the development of cataracts in patients that carry missense mutations in MIP.

GONAD: A Novel CRISPR/Cas9 Genome Editing Method that Does Not Require Ex Vivo Handling of Embryos.

Transgenic technologies used for creating a desired genomic change in animals involve three critical steps: isolation of fertilized eggs, microinjection of transgenic DNA into them and their subsequent transfer to recipient females. These ex vivo steps have been widely used for over 3 decades and they were also readily adapted for the latest genome editing technologies such as ZFNs, TALENs, and CRISPR/Cas9 systems. We recently developed a method called GONAD (Genome editing via Oviductal Nucleic Acids Delivery) that does not require all the three critical steps of transgenesis and therefore relieves the bottlenecks of widely used animal transgenic technologies. Here we provide protocols for the GONAD system.

A 1-bp deletion in Fgf5 causes male-dominant long hair in the Syrian hamster.

Hair length in mammals is generally regulated by the hair cycle, and its disruption leads to abnormal hair morphogenesis in several species. FGF5, one of the hair cycle regulators, has a role in inducing catagen, and that mutation causes abnormal hair length in both sexes in humans, mice, dogs, and cats. Male-dominant long-haired coat (MALC) is an inbred strain of Syrian hamster exhibiting spontaneous long hair in males. After castration, MALC exhibited significantly shorter hair than the control individuals, but testosterone administration to castrated MALC showed reversion to the original phenotype. Moreover, flutamide administration led to MALC phenotype repression. Histological analysis revealed that hair follicle regression was shown in the wild-type 4 weeks after depilation, but that of MALC remained in the anagen phase. We detected a c.546delG of Fgf5 in MALC (Fgf5malc) that might lead to truncation resulting from a frame shift in FGF5 (p.Arg184GlyfsX6). Additionally, homozygous Fgf5malc was only detected in long-haired (Slc:Syrian×MALC)F2 and (J-2-Nn×MALC)F2 progenies, and all homozygous wild and heterozygous Fgf5malc individuals showed normal hair length. Thus, Fgf5malc leads to male-dominant long hair via a prolonged anagen phase which is affected by testosterone in hamsters. To our knowledge, this report is the first to present the sexual dimorphism of hair length caused by the Fgf5 mutation.

GONAD: Genome-editing via Oviductal Nucleic Acids Delivery system: a novel microinjection independent genome engineering method in mice.

Microinjection is considered the gold standard technique for delivery of nucleic acids (NAs; transgenes or genome editing tools such as CRISPR/Cas9 systems) into embryos, for creating genetically modified organisms. It requires sophisticated equipment as well as well-trained and highly skilled personnel to perform the micro-injection technique. Here, we describe a novel and simple microinjection-independent technique, called Genome-editing via Oviductal Nucleic Acids Delivery (GONAD). Using GONAD, we show that NAs (e.g., eGFP mRNA or Cas9 mRNA/sgRNAs) can be effectively delivered to pre-implantation embryos within the intact mouse oviduct by a simple electroporation method, and result in the desired genetic modification in the embryos. Thus GONAD can bypass many complex steps in transgenic technology such as isolation of zygotes, microinjection of NAs into them, and their subsequent transfer to pseudo-pregnant animals. Furthermore, this method can potentially be used for genome editing in species other than mice.

Expression of truncated PITX3 in the developing lens leads to microphthalmia and aphakia in mice.

Microphthalmia is a severe ocular disorder, and this condition is typically caused by mutations in transcription factors that are involved in eye development. Mice carrying mutations in these transcription factors would be useful tools for defining the mechanisms underlying developmental eye disorders. We discovered a new spontaneous recessive microphthalmos mouse mutant in the Japanese wild-derived inbred strain KOR1/Stm. The homozygous mutant mice were histologically characterized as microphthalmic by the absence of crystallin in the lens, a condition referred to as aphakia. By positional cloning, we identified the nonsense mutation c.444C>A outside the genomic region that encodes the homeodomain of the paired-like homeodomain transcription factor 3 gene (Pitx3) as the mutation responsible for the microphthalmia and aphakia. We examined Pitx3 mRNA expression of mutant mice during embryonic stages using RT-PCR and found that the expression levels are higher than in wild-type mice. Pitx3 over-expression in the lens during developmental stages was also confirmed at the protein level in the microphthalmos mutants via immunohistochemical analyses. Although lens fiber differentiation was not observed in the mutants, strong PITX3 protein signals were observed in the lens vesicles of the mutant lens. Thus, we speculated that abnormal PITX3, which lacks the C-terminus (including the OAR domain) as a result of the nonsense mutation, is expressed in mutant lenses. We showed that the expression of the downstream genes Foxe3, Prox1, and Mip was altered because of the Pitx3 mutation, with large reductions in the lens vesicles in the mutants. Similar profiles were observed by immunohistochemical analysis of these proteins. The expression profiles of crystallins were also altered in the mutants. Therefore, we speculated that the microphthalmos/aphakia in this mutant is caused by the expression of truncated PITX3, resulting in the abnormal expression of downstream targets and lens fiber proteins.

A 5-bp insertion in Mip causes recessive congenital cataract in KFRS4/Kyo rats.

We discovered a new cataract mutation, kfrs4, in the Kyoto Fancy Rat Stock (KFRS) background. Within 1 month of birth, all kfrs4/kfrs4 homozygotes developed cataracts, with severe opacity in the nuclei of the lens. In contrast, no opacity was observed in the kfrs4/+ heterozygotes. We continued to observe these rats until they reached 1 year of age and found that cataractogenesis did not occur in kfrs4/+ rats. To define the histological defects in the lenses of kfrs4 rats, sections of the eyes of these rats were prepared. Although the lenses of kfrs4/kfrs4 homozygotes showed severely disorganised fibres and vacuolation, the lenses of kfrs4/+ heterozygotes appeared normal and similar to those of wild-type rats. We used positional cloning to identify the kfrs4 mutation. The mutation was mapped to an approximately 9.7-Mb region on chromosome 7, which contains the Mip gene. This gene is responsible for a dominant form of cataract in humans and mice. Sequence analysis of the mutant-derived Mip gene identified a 5-bp insertion. This insertion is predicted to inactivate the MIP protein, as it produces a frameshift that results in the synthesis of 6 novel amino acid residues and a truncated protein that lacks 136 amino acids in the C-terminal region, and no MIP immunoreactivity was observed in the lens fibre cells of kfrs4/kfrs4 homozygous rats using an antibody that recognises the C- and N-terminus of MIP. In addition, the kfrs4/+ heterozygotes showed reduced expression of Mip mRNA and MIP protein and the kfrs4/kfrs4 homozygotes showed no expression in the lens. These results indicate that the kfrs4 mutation conveys a loss-of-function, which leads to functional inactivation though the degradation of Mip mRNA by an mRNA decay mechanism. Therefore, the kfrs4 rat represents the first characterised rat model with a recessive mutation in the Mip gene.