Next Article in Journal
On the Population Dynamics of Junk: A Review on the Population Genomics of Transposable Elements
Previous Article in Journal
Genome Analysis of Hypomyces perniciosus, the Causal Agent of Wet Bubble Disease of Button Mushroom (Agaricus bisporus)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 6
10.3390/genes10060418
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Loci Controlling the Dwarfism Trait in the White Sailfin Molly (Poecilia latipinna) Using Genome-Wide Association Studies Based on Genotyping-By-Sequencing

by
Fan Shao
†,
Jing Liu
†,
Mengyuan Ren
†,
Junying Li
,
Haigang Bao
* and
Changxin Wu
National Engineering Laboratory for Animal Breeding, Beijing Key Laboratory of Animal Genetic Improvement, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
*
Author to whom correspondence should be addressed.
†
These authors made equal contributions to this work.
Genes 2019, 10(6), 418; https://doi.org/10.3390/genes10060418
Submission received: 22 March 2019 / Revised: 17 May 2019 / Accepted: 27 May 2019 / Published: 30 May 2019
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Dwarfism is a condition defined by low harvest weight in fish, but also results in strange body figures which may have potential for the selective breeding of new ornamental fish strains. The objectives of this study are to reveal the physiological causes of dwarfism and identify the genetic loci controlling this trait in the white sailfin molly. Skeletons of dwarf and normal sailfin mollies were observed by X-ray radioscopy and skeletal staining. Genome-wide association studies based on genotyping-by-sequencing (n = 184) were used to map candidate genomic regions associated with the dwarfism trait. Quantitative real-time PCR was performed to determine the expression level of candidate genes in normal (n = 8) and dwarf (n = 8) sailfin mollies. We found that the dwarf sailfin molly has a short and dysplastic spine in comparison to the normal fish. Two regions, located at NW_015112742.1 and NW_015113621.1, were significantly associated with the dwarfism trait. The expression level of three candidate genes, ADAMTS like 1, Larp7 and PPP3CA, were significantly different between the dwarf and normal sailfin mollies in the hepatopancreas, with PPP3CA also showing significant differences in the vertebrae and Larp7 showing significant differences in the muscle. This study identified genomic regions and candidate genes associated with the dwarfism trait in the white sailfin molly and would provide a reference to determine dwarf-causing variations.

1. Introduction

Dwarf ornamental fish, such as the balloon molly and dwarf swordtail, are popular in China for their lovely body shapes and graceful activities. By analyzing the skeletal morphology of the dwarf and normal swordtail fishes, Li et al. found the spine length of the dwarf swordtail to be shorter than that of the normal fish, and speculated that dwarfism in the swordtail fish was caused by a developmental defect in the spine [1].
Spinal deformities are a frequent phenomenon in many different fish species and can be caused by many factors during fish development, such as bacterial and parasitic infections, water environment, elevated egg incubation temperatures, pollution, nutrients, inappropriate light regimes and vaccination, etc. [2,3,4,5,6,7]. Spinal deformities often result in increased metabolic costs in the deformed fish, leading to reduced growth and low harvest weight [8]. The more deformed the spine is, the smaller the size and weight of the fish [9]. Genetic variations may play an important role in the growth and development of animal spines [10,11]. McMenamin et al. revealed a premature stop codon in the gh1 gene of zebrafish which retarded growth [12]. Inohaya et al. reported that a DNA fragment insertion in the promoter region of the wnt4b gene caused fused vertebrae and shortened body length in medaka [13]. Li et al. suggested that the dwarf trait of the swordtail fish could be inherited and, in accordance with the basic characteristics of quantitative trait inheritance, might be controlled by multiple quantitative trait genes [1].
Genome-wide association studies (GWAS) present a popular and powerful method to map candidate loci associated with target traits [14]. High-throughput single nucleotide polymorphism (SNP) chips with thousands of SNPs uniformly distributed along the genome are usually needed to genotype SNPs and improve the effect of GWAS [15]. However, for species without economical high-density SNP chips, a simple, low-cost and efficient method, such as genotyping-by-sequencing (GBS), might be the best technology, based on second-generation sequencing technology, to obtain enough SNPs to uniformly cover the entire genome [16].
The sailfin molly (Poecilia latipinna) is a type of ornamental fish that originated from America or Mexico [17]. Generally, two types of figures are observed in the sailfin molly, one is dwarf and the other is normal. The dwarf sailfin molly, also called the balloon belly molly, is more popular than the normal sailfin molly, but the mechanism of how the balloon phenotype comes into being is still unknown. In the present study, the skeleton shape of the dwarf and normal sailfin mollies was observed and GWAS based on GBS was performed on 184 white sailfin mollies to map the genomic regions associated with the dwarfism trait. We also determined the expression levels of candidate genes within important genomic regions to screen for candidate genes which could be important in controlling the dwarfism trait.

2. Materials and Methods

2.1. Ethics Statement

All experimental procedures and animal use were approved by the Ethics Review Committee for Laboratory Animal Welfare and Animal Experiment of China Agricultural University (approval number: CAU20170502-3).

2.2. Animals and Phenotyping

We purchased 184 white sailfin mollies at five to six months of age from Shanghai Yuyi Tropical Ornamental Fish Farm Corporation. Classification of body shape was based on the general market classification criteria. The torso length (TL) and body height (H) of each white sailfin molly were measured using a Vernier caliper (Figure 1) and the ratio of TL and H (TL/H) was calculated and used as an index for dwarfism and for subsequent analysis. The least squares algorithm of the SAS software package (SAS 9.2, SAS Institute Inc., Cary, NC, USA) was used to analyze the phenotypic values, which were then expressed as mean ± SE.

2.3. Comparison of Skeletal Structures between the Dwarf Molly and the Normal Molly

Several specimens of each body type were selected at random and their skeletal features were observed by an X-ray radioscopy inspection system (two specimens of each body type; FX Pro/FX, Carestream, USA) or skeletal staining [18] (three and five specimens of the normal and dwarf mollies, respectively).

2.4. Genotyping

The tail fin of each fish was collected to extract DNA. DNA extraction was performed using the TIANamp Genomic DNA Kit (cat. #DP304-02, TIANGEN, Beijing, China) according to the protocol supplied. All samples were genotyped using the genotyping-by-sequencing method [19] with some revisions. Briefly, genomic DNA was incubated at 37 °C with HaeIII, EcoRI, MseI, T4 ligase, ATP and a MseI Y adapter. The oligo nucleotide sequences containing barcodes are as follows: 5′-ACACTCTTTCCCTACACGACGCTCTTCCGATCTxxxx and 5′-TAAyyyyAGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGT. The common adapters are as follows: 5′-TAAAGATCGGAAGAGCGGTTCAGCAGGAATGCCGAG and 5′-CTCGGCATTCCTGCTGAACCGCTCTTCCGATCT. After the restriction-ligation reactions, samples were further digested with the restriction enzymes HaeIII and EcoRI at 37 °C and then purified with an Agencourt AMPure XP (Beckman, Brea, CA, USA). PCR reactions were performed on the purified samples using a Phusion Master Mix universal primer and index primer. The polymerase chain reaction (PCR) products were purified again, and all PCR products were pooled to generate one library that included all 184 individuals. Fragments of 375–400 bp (with indexes and adapters) were isolated using a Gel Extraction Kit (Qiagen, Hilden, Germany). After being purified again, the fragments were diluted for sequencing and pair-end sequencing (PE, 150 bp) was performed using the Illumina HiSeq 4000 sequencing platform at Novogene Bioinformatics Technology Co. Ltd. (Beijing, China). After quality control, Illumina sequencing short reads were processed with Burrows-Wheeler Aligner’s maximal exact match (BWA-mem) [20], smatools-bcftools pipeline [21], vcftools [22] and several in-house scripts to call and filter SNPs and genotype all samples.

2.5. Population Structural and Linkage Disequilibrium Analyses

The population structure was estimated by principal component analyses (PCA) using the software package GCTA [23] and linkage disequilibrium (LD) analysis was performed with PopLDdecay [24].

2.6. Genome-Wide Association Studies

We used GEMMA software [25] to perform genome-wide association studies for the ratio of TL and H using a generalized linear mixed model. The population genetic structure was considered to be a fixed effect and individual kinship was regarded as a random effect in the mixed linear model. The linear mixed model was constructed as follows:
y = Xα + Zβ +ξ+ e
where y is the vector of the ratio of TL/H, X is the covariance matrix of the fixed effect, α is the estimated vector of the fixed effect, Z is the indicator matrix of the SNP, β is the effect of the SNP, ξ ~ N(0, Kφ2) is the polygenic effect, K is the genomic kinship matrix to adjust for kinship among samples, φ2 is the additive genetic variance and e ~ (0, δe2) is the random residual variance. The significance threshold was set to −log(1/n) [26].

2.7. RNA Extraction and Quantitative Real-Time PCR

We selected 16 white sailfin mollies of about six months in age, including eight dwarf and eight normal sailfin mollies, at random from an ornamental fish market in Beijing. All fish samples were collected and raised in the same tank with fresh water at a temperature of 25 °C. After about two weeks the mollies were killed by rapidly cutting off the head, and the hepatopancreas and muscle were separated out immediately and stored at −80 °C for RNA extraction. Vertebrae from the end of the head to the anus were also cut off and dipped into trizol for about 25 min to remove the remaining muscle. After removing all muscle, the remaining vertebrae were stored at −80 °C for RNA extraction, or dipped into a new trizol solution for homogenation. After homogenation, the total RNA was extracted using the RNAprep Pure Tissue Kit (cat. #DP431, TIANGEN) according to the protocol supplied. The cDNA was obtained using a FastQuant RT Kit (cat. #KR106, TIANGEN) and quantitative real time PCR (qRT-PCR) was performed with 20 µL of the reaction system using SuperReal PreMix Plus (SYBR Green; cat. #FP205, TIANGEN) according to the protocols supplied. The thermal cycling procedure of qRT-PCR was as follows: 95 ℃ for 15 min followed by 40 cycles of 95 °C for 20 s, 57 °C for 20 s and 72 °C for 20 s, followed by a melting/dissociation curve stage for 6 s. Each sample had three technical duplications, and the expression levels of the target genes were adjusted by β-actin RNA expression. The primer sequences used in the present study are listed in Table S1. The 2−∆∆Ct method was used to calculate the relative expression of each target gene [27]. Differences in expression were tested by t-tests using GraphPad Prism (Version 6; GraphPad Software, La Jolla, CA, USA). The significance thresholds were set at p < 0.05 and p < 0.01.

3. Results and Discussion

3.1. Phenotypic Data of the White Molly

After discarding a few ambiguous samples, 184 white mollies, including 92 normal white sailfin mollies (NWMs, 46 females and 46 males) and 92 dwarf white sailfin mollies (DWMs, 46 females and 46 males), were collected. The TL, H and TL/H were measured and calculated as shown in Table 1. From Table 1, we can see that the average TL/H of NWMs was 2.2256, which is significantly higher than that of DWMs (p < 0.01).

3.2. Comparison of Skeletal Structures between the Dwarf Molly and the Normal Molly

In appearance, NWMs are slender and flat while DWMs have short bodies and convexly arched backs. The X-ray radioscopy inspection system and skeletal staining [18] were used to observe the difference in skeletal structure between the NWMs and DWMs. As shown in Figure 2, there were differences in the skeletal structures between the NWMs and DWMs, and the most obvious difference was in the spine. DWMs were short and their vertebrae were dysplastic in comparison with NWMs. The skeletal stainings were clear enough to count the vertebrae, as shown in Figure 2B,C. The white sailfin molly has 27–29 vertebrae and no significant difference (p = 0.16) was found in vertebrae number between the NWMs (n = 3) and DWMs (n = 5) in the present study.
Dwarfism is a manifestation of abnormal skeletal development [2,28]. Skeletal dysplasia of tibia, cartilage and the spine can cause animal dwarfism [10,29,30]. Spinal deformities may play an important role in fish dwarfism [1]. Body deformity and dwarfism were often found to be associated with scoliosis, as demonstrated by Afonso et al. [31]. However, in the present study we did not find curved spines in the DWMs. Since no other reasonable factor was found to explain the dwarfism, we hypothesize that the dysplastic vertebrae seen in the DWMs are the main cause of dwarfism in the white sailfin molly, as has been suggested previously in a study of dwarfism in swordtail fish [1].

3.3. Sequencing Data Statistics and Quality Assessment

A summary of the next-generation sequencing is shown in Table S2. In total, 83.615 Gb of raw data and 83.612 Gb of clean data were obtained from 184 samples of white sailfin mollies, averaging 0.454 Gb per sample. For each sample, the GC distribution of the sequencing data was normal, and the value of Q30 was larger than 90.08% (Table S2). The results show that the sequencing quality was high and the data could be used for subsequent analysis.

3.4. Single Nucleotide Polymorphism Filtering and Population Component Analysis

In total, 3,021,852 SNPs were obtained from NWMs and DWMs using GBS technology. After removing the SNPs which fell outside of the following thresholds: sequencing depth below 4×, detection rate less than 20%, minimum allele frequency less than 0.05 and Hardy–Weiberg equilibrium (HWE) less than 0.00001, we obtained 113,814 high-quality SNPs for subsequent analysis with an SNP density of about 1.7 SNP per 10 kb in the whole genome. LD analysis was also performed with PopLDdecay [24] and the results are shown in Figure S1. With LD reduced to 50% of the highest original value (r2 > 0.30) set as the threshold, the LD decay distance of the population is about 100 kb (Figure S1). Principle component analysis (PCA) is a useful method to reveal population stratification, which can cause false positives in GWAS [32]. In the present study, PCA was performed with GCTA [23] and the results are shown in Figure 3. From Figure 3, we can see significant population stratification, and so a univariate linear mixed model in GEMMA was performed to adjust for population stratification in the subsequent GWAS.

3.5. Genome-Wide Association Studies

In order to correct for the influence of individual size on phenotypic judgement, TL/H was calculated and used as the index of dwarfism in this study, and GEMMA [25] was used to test the association between the 113,814 SNPs and TL/H. Two regions on NW_015112742.1 and NW_015113621.1 were significantly associated with the trait (Figure 4, Table S3). The percent variance explained by the two loci was also calculated with GCTA [23] and the result is shown in Table S4. From Table S4, we can see that 99.3% of the phenotype variance can be explained by these two loci, which increases the credibility of our mapping results.

3.6. Prediction of Candidate Genes

To determine the candidate TL/H-related genes, significant SNPs were annotated with ANNOVAR [33] based on P_latipinna-1.0 (GCA_001443285.1; ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/001/443/285/GCF_001443285.1_P_latipinna-1.0/). Five genes, namely ADAMTS like 1, NLRP 12 like, Larp 7, PPP3CA and Bank 1, were annotated by the ranges from 20 kb upstream and downstream of the physical locations of those significant SNPs (Table S3).
Trait-causing genes often induce changes in gene expression levels, and so the gene expression levels of the five genes, except Bank 1, in the candidate region were determined by qRT-PCR (Figure 5). Bank 1 (B-cell scaffold protein with ankyrin repeats 1) encodes a B-cell-specific scaffold protein, a novel substrate of tyrosine kinases, and its overexpression in B-cells can enhance calcium mobilization induced by B-cell antigen receptor (BCR) stimulation [34]. It has been reported that Bank 1 is associated with systemic lupus erythematosus (SLE) [35,36]. The expression levels of the Bank 1 gene were too low to be accurately determined in the three tissues examined in the present study, so no expression data of the Bank 1 gene is shown here. There was no difference in the expression level of the NLRP 12-like gene between DWMs and NWMs in the vertebrae, hepatopancreas or muscle (Figure 5A). ADAMTS-like proteins may be enhancers of the ADAMTS proteases [37]. The ADAMTS-like 1 gene, also known as C9orf94, PUNCTIN, ADAMTSR1 or ADAMTSL-1 in humans, encodes a novel ADAMTS-like molecule which lacks the typical domains of the ADAMTS family and is composed of only ADAMTS ancillary domains in the extracellular matrix [37]. A significant difference (p < 0.05) between the expression levels of ADAMTS-like 1 in DWMs and NWMs was observed in the hepatopancreas, with the value observed for NWMs being greater than that observed for DWMs, but no difference was found in the vertebrae or muscle (Figure 5B). The Larp7 gene encodes a member of the ribonucleoprotein domain family 7, which is a unit of the ribonucleoprotein complex 7SK snRNP. Some studies have suggested that Larp7 mutations cause dwarf phenotypes [38,39,40]. Alamazi et al. [38] found that a seven-base homozygous repeat (c.1024_1030dupAAGGATA, p.T344K) on exon 8 of the Larp7 gene, which causes a frameshift mutation and prematurely terminates translation, is an underlying mechanism of primordial dwarfism (PD), which is characterized by embryonic growth degeneration, facial deformity and mental impairment. Recently, two more mutations (c.1091_1094delCGGT and c.1045_1051dupAAGGATA) on the Larp7 gene have been associated with PD [39]. In zebrafish, Larp7 has also been found to play a crucial role in embryo development [40]. In the present study, the expression of the Larp7 gene in DWMs was significantly higher (p < 0.05) than that in NWMs in the hepatopancreas and muscle, but not in the vertebrae (Figure 5C). PPP3CA encodes the alpha isoform of a subunit of calcineurin and is involved in neurodevelopmental diseases and systemic lupus erythematosus [35,36]. The expression of the PPP3CA gene in the hepatopancreas and vertebrae of DWMs was significantly (p < 0.05) or extremely significantly (p < 0.01) higher than that in NWMs (Figure 5D). The data also showed that DWMs possessed a higher expression level of the PPP3CA gene in muscle compared with NWMs, despite the fact that the difference was not significant (Figure 5D; p = 0.16).
Some factors, such as fish developmental stage and potential differences in food intake between the two phenotypes, may affect gene expression levels. Although there is no significant difference between the two phenotypes of mollies in the expression of the NLRP 12-like gene at an age of six months, there may be differences at other stages, so we cannot exclude the role of NLRP 12-like in phenotypic formation based solely on this study. In comparison with NWMs, dwarf individuals may not be able to easily get food due to their deformed mouths and inflexible bodies, which can have impacts on gene expression. In order to reduce this effect, for the two-week duration of this study we fed all fish two times a day and evenly distributed enough food in the pool to ensure that all individuals could get enough food.

4. Conclusions

We found that the white sailfin molly has 27–29 vertebrae in our study. The dwarf sailfin molly has a short and dysplastic spine in comparison with the normal sailfin molly, which indicated that the arrested longitudinal development of the vertebrae could be the main cause of the dwarf trait. Six SNPs significantly associated with pygmy traits were screened in the white sailfin molly, located at NW_015112742.1 and NW_015113621.1. According to the annotation of gene function and relative expression levels, we suggest PPP3CA and Larp7 as the most important candidate dwarf-causing genes to be further investigated.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/10/6/418/s1, Table S1: The primer sequences for quantitative real-time PCR. Table S2: Summary of high-throughput sequencing data and quality control. Table S3: Significant SNPs associated with TL/H and candidate genes. Table S4: Estimated variance explained by the significant SNPs. Figure S1: Linkage disequilibrium (LD) decay distance of the population.

Author Contributions

Conceptualization, C.W. and H.B.; Methodology, H.B., F.S., M.R., J.L. (Jing Liu); Validation, J.L. (Jing Liu), M.R., H.B.; Formal Analysis, F.S., H.B., M.R., J.L. (Jing Liu); Investigation, H.B., F.S., M.R., J.L. (Jing Liu); Resources, J.L. (Junying Li), C.W.; Data Curation, H.B., F.S., M.R., J.L. (Jing Liu); Writing—Original Draft Preparation, J.L. (Jing Liu), H.B., F.S., M.R.; Writing—Review and Editing, H.B., J.L. (Jing Liu), F.S., M.R., J.L. (Junying Li), C.W.; Visualization, F.S., M.R. J.L. (Jing Liu); Supervision, H.B.; Project Administration, H.B.; Funding Acquisition, C.W.

Funding

This research received no external funding.

Acknowledgments

We would like to thank Pengfei Ji for his help in data analyses. This study was funded by the Science and Technology Foundation of China Agricultural University (Grand No. 88057007).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, K.B.; Liu, C.; Chang, O.Q.; Wang, F.; Ma, B.Y.; Wu, S.Q. The skeletal defects and inheritance in dwarf mutants of Xiphophorus helleri. J. Fish. China 2011, 35, 809–815. [Google Scholar] [CrossRef]
  2. Witten, P.E.; Gil-Martens, L.; Huysseune, A.; Takle, H.; Hjelde, K. Towards a classification and an understanding of developmental relationships of vertebral body malformations in Atlantic salmon (Salmo salar L.). Aquaculture 2009, 295, 6–14. [Google Scholar] [CrossRef]
  3. Haller, G.; McCall, K.; Jenkitkasemwong, S.; Sadler, B.; Antunes, L.; Nikolov, M.; Whittle, J.; Upshaw, Z.; Shin, J.; Baschal, E.; et al. A missense variant in SLC39A8 is associated with severe idiopathic scoliosis. Nat. Commun. 2018, 9, 1. [Google Scholar] [CrossRef] [PubMed]
  4. Polak-Juszczak, L. Impact of strontium on skeletal deformities in Baltic cod (Gadus morhua callaris L.). Chemosphere 2011, 83, 486–491. [Google Scholar] [CrossRef] [PubMed]
  5. Munday, J.S.; Perrott, M.R.; Symonds, J.E.; Walker, S.P.; Preece, M.A.; Davie, P.S. Prevalence of spinal abnormalities in Chinook salmon smolt and influence of early rearing temperature and growth rates. J. Fish Dis. 2018, 41, 1111–1116. [Google Scholar] [CrossRef]
  6. Kupsco, A.; Schlenk, D. Molecular mechanisms of selenium-Induced spinal deformities in fish. Aquat. Toxicol. 2016, 179, 143–150. [Google Scholar] [CrossRef]
  7. Torres-Dowdall, J.; Karagic, N.; Plath, M.; Riesch, R. Evolution in caves: Selection from darkness causes spinal deformities in teleost fishes. Biol. Lett. 2018, 14, 20180197. [Google Scholar] [CrossRef]
  8. Chin, H.N.; Loh, R.; Hong, Y.C.; Gibson-Kueh, S. Case studies of spinal deformities in ornamental koi, Cyprinus carpio L. J. Fish Dis. 2017, 40, 65–71. [Google Scholar] [CrossRef]
  9. Tong, X.H.; Liu, Q.H.; Xu, S.H.; Ma, D.Y.; Xiao, Z.Z.; Xiao, Y.S.; Li, J. Skeletal development and abnormalities of the vertebral column and of the fins in hatchery-reared turbot Scophthalmus maximus. J. Fish. Biol. 2012, 80, 486–502. [Google Scholar] [CrossRef]
  10. Gedeon, A.K.; Tiller, G.E.; Le Merrer, M.; Heuertz, S.; Tranebjaerg, L.; Chitayat, D.; Robertson, S.; Glass, I.A.; Savarirayan, R.; Cole, W.G.; et al. The molecular basis of X-linked spondyloepiphyseal dysplasia tarda. Am. J. Hum. Genet. 2001, 68, 1386–1397. [Google Scholar] [CrossRef]
  11. Xia, X.Y.; Cui, Y.X.; Zhou, Y.C.; Zhou, X.; Shi, Y.C.; Wei, L.; Li, X.J.; Huang, Y.F.; Huang, T.T. A novel insertion mutation in the SEDL gene results in X-linked spondyloepiphyseal dysplasia tarda in a large Chinese pedigree. Clin. Chim. Acta 2009, 410, 39–42. [Google Scholar] [CrossRef] [PubMed]
  12. McMenamin, S.K.; Minchin, J.E.N.; Gordon, T.N.; Rawls, J.F.; Parichy, D.M. Dwarfism and increased adiposity in the gh1 mutant Zebrafish vizzini. Endocrinology 2013, 154, 1476–1487. [Google Scholar] [CrossRef] [PubMed]
  13. Inohaya, K.; Takano, Y.; Kudo, A. Production of Wnt4b by floor plate cells is essential for the segmental patterning of the vertebral column in medaka. Development 2010, 137, 1807–1813. [Google Scholar] [CrossRef] [PubMed]
  14. Tenghe, A.M.M.; Bouwman, A.C.; Berglund, B.; Strandberg, E.; de Koning, D.J.; Veerkamp, R.F. Genome-wide association study for endocrine fertility traits using single nucleotide polymorphism arrays and sequence variants in dairy cattle. J. Dairy Sci. 2016, 99, 5470–5485. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Edea, Z.; Hong, J.K.; Jung, J.H.; Jung, J.H.; Kim, D.W.; Kim, Y.M.; Kim, E.S.; Shin, S.S.; Jung, Y.C.; Kim, K.S. Detecting selection signatures between Duroc and Duroc synthetic pig populations using high-density SNP chip. Anim. Genet. 2017, 48, 473–477. [Google Scholar] [CrossRef] [PubMed]
  16. Siadjeu, C.; Mayland-Quellhorst, E.; Albach, D.C. Genetic diversity and population structure of trifoliate yam (Dioscorea dumetorum Kunth) in Cameroon revealed by genotyping-by-sequencing (GBS). BMC Plant Biol. 2018, 18, 359–372. [Google Scholar] [CrossRef] [PubMed]
  17. Apodaca, J.J.; Trexler, J.C.; Jue, N.K.; Schrader, M.; Travis, J. Large-Scale Natural Disturbance Alters Genetic Population Structure of the Sailfin Molly, Poecilia latipinna. Am. Nat. 2013, 181, 254–263. [Google Scholar] [CrossRef]
  18. van Eeden, F.J.; Granato, M.; Schach, U.; Brand, M.; Furutani-Seiki, M.; Haffter, P.; Hammerschmidt, M.; Heisenberg, C.P.; Jiang, Y.J.; Kane, D.A.; et al. Mutations affecting somite formation and patterning in the zebrafish, Danio rerio. Development 1996, 123, 153–164. [Google Scholar]
  19. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef]
  20. Li, H.; Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 2010, 26, 589–595. [Google Scholar] [CrossRef]
  21. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; 1000 Genome Project Data Processing Subgroup. The Sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  22. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef] [PubMed]
  23. Yang, J.; Lee, S.H.; Goddard, M.E.; Visscher, P.M. GCTA: A Tool for Genome-wide Complex Trait Analysis. Am. J. Hum. Genet. 2011, 88, 76–82. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhang, C.; Dong, S.S.; Xu, J.Y.; He, W.M.; Yang, T.L. PopLDdecay: A fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 2019, 35, 1786–1788. [Google Scholar] [CrossRef] [PubMed]
  25. Zhou, X.; Stephens, M. Genome-wide efficient mixed-model analysis for association studies. Nat. Genet. 2012, 44, 821–824. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sonah, H.; O’Donoughue, L.; Cober, E.; Rajcan, I.; Belzile, F. Identification of loci governing eight agronomic traits using a GBS-GWAS approach and validation by QTL mapping in soya bean. Plant Biotechnol. J. 2015, 13, 211–221. [Google Scholar] [CrossRef] [PubMed]
  27. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  28. Xie, F.; Zhu, M.; Cao, Y.L.; Xiong, F. Correlation of SHOX gene deficiency and mutation with skeleton X-ray characteristics of idiopathic short stature patients. J. Third Mil. Med. Univ. 2011, 06, 604–607. [Google Scholar]
  29. Zhou, Z.H.; Shu, D.M.; Guo, K.Y.; Yang, C.F.; Xiong, Y. Study on the Inheritance Law of Yellow Feather Dwarf Chicken. Chin. J. Anim. Sci. 1998, 03, 26–27. [Google Scholar]
  30. Yun, S.F.; Yan, Z.L.; Li, C.P. The present situation of pathogenesis and genetics of short stature. Chin. J. Birth Health Hered. 2012, 03, 137–139. [Google Scholar]
  31. Afonso, J.M.; Montero, D.; Robaina, L.; Astorga, N.; Izquierdo, M.S.; Ginés, R. Association of a lordosis-scoliosis-kyphosis deformity in gilthead seabream (Sparus aurata) with family structure. Fish Physiol. Biochem. 2000, 22, 159–163. [Google Scholar] [CrossRef]
  32. Reich, D.; Price, A.L.; Patterson, N. Principal component analysis of genetic data. Nat. Genet. 2008, 40, 491–492. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
  34. Yokoyama, K.; Su, I.I.; Tezuka, T.; Yasuda, T.; Mikoshiba, K.; Tarakhovsky, A.; Yamamoto, T. BANK regulates BCR-induced calcium mobilization by promoting tyrosine phosphorylation of IP(3) receptor. EMBO J. 2002, 21, 83–92. [Google Scholar] [CrossRef] [PubMed]
  35. Wilbe, M.; Jokinen, P.; Truvé, K.; Seppala, E.H.; Karlsson, E.K.; Biagi, T.; Hughes, A.; Bannasch, D.; Andersson, G.; Hansson-Hamlin, H.; et al. Genome-wide association mapping identifies multiple loci for a canine SLE-related disease complex. Nat. Genet. 2010, 42, 250–254. [Google Scholar] [CrossRef] [PubMed]
  36. Chang, Y.K.; Yang, W.; Zhao, M.; Mok, C.C.; Chan, T.M.; Wong, R.W.; Lee, K.W.; Mok, M.Y.; Wong, S.N.; Ng, I.O.; et al. Association of BANK1 and TNFSF4 with systemic lupus erythematosus in Hong Kong Chinese. Genes Immun. 2009, 10, 414–420. [Google Scholar] [CrossRef] [PubMed]
  37. Hirohata, S.; Wang, L.W.; Miyagi, M.; Yan, L.; Seldin, M.F.; Keene, D.R.; Crabb, J.W.; Apte, S.S. Punctin, a novel ADAMTS-like molecule, ADAMTSL-1, in extracellular matrix. J. Biol. Chem. 2002, 277, 12182–12189. [Google Scholar] [CrossRef] [PubMed]
  38. Alazami, A.M.; Al-Owain, M.; Alzahrani, F.; Shuaib, T.; Al-Shamrani, H.; Al-Falki, Y.H.; Al-Qahtani, S.M.; Alsheddi, T.; Colak, D.; Alkuraya, F.S. Loss of function mutation in LARP7, chaperone of 7SK ncRNA, causes a syndrome of facial dysmorphism, intellectual disability, and primordial dwarfism. Hum. Mutat. 2012, 33, 1429–1434. [Google Scholar] [CrossRef]
  39. Hollink, I.H.; Alfadhel, M.; Al-Wakeel, A.S.; Ababneh, F.; Pfundt, R.; de Man, S.A.; Jamra, R.A.; Rolfs, A.; Bertoli-Avella, A.M.; van de Laar, I.M. Broadening the phenotypic spectrum of pathogenic LARP7 variants: Two cases with intellectual disability, variable growth retardation and distinct facial features. J. Hum. Genet. 2016, 61, 229–233. [Google Scholar] [CrossRef]
  40. Barboric, M.; Lenasi, T.; Chen, H.; Johansen, E.B.; Guo, S.; Peterlin, B.M. 7SK snRNP/P-TEFb couples transcription elongation with alternative splicing and is essential for vertebrate development. Proc. Natl. Acad. Sci. USA 2009, 106, 7798–7803. [Google Scholar] [CrossRef] [Green Version]
Genes 10 00418 g001 550
Figure 1. Diagram of the measured body traits.
Figure 1. Diagram of the measured body traits.
Genes 10 00418 g001
Genes 10 00418 g002 550
Figure 2. Skeletons of the normal molly and the dwarf molly. (A) Skeletons of normal and dwarf mollies under X-ray radioscopy; (B) skeleton of a normal molly using skeletal staining; (C) skeleton of a dwarf molly using skeletal staining.
Figure 2. Skeletons of the normal molly and the dwarf molly. (A) Skeletons of normal and dwarf mollies under X-ray radioscopy; (B) skeleton of a normal molly using skeletal staining; (C) skeleton of a dwarf molly using skeletal staining.
Genes 10 00418 g002
Genes 10 00418 g003 550
Figure 3. Principal component analysis performed with the GCTA package. Blue triangles represent NWMs and red triangles represent DWMs. Samples with the same body shape tend to come together, which implies that population stratification should be considered in the subsequent genome-wide association studies (GWAS). PC1: Principle component 1; PC2: Principle component 2.
Figure 3. Principal component analysis performed with the GCTA package. Blue triangles represent NWMs and red triangles represent DWMs. Samples with the same body shape tend to come together, which implies that population stratification should be considered in the subsequent genome-wide association studies (GWAS). PC1: Principle component 1; PC2: Principle component 2.
Genes 10 00418 g003
Genes 10 00418 g004 550
Figure 4. Manhattan plot (A) and quantile–quantile (Q–Q) plot (B) for the association of all single nucleotide polymorphisms (SNPs) with the ratio of tail fin length to height (TL/H). Regional plots (C,D) of the genome-wide significant SNPs. Six SNPs reached genome-wide significant level. In plots A, C and D the red line shows the p-value threshold of 8.79 × 10−6 (−log10(p-value) > 5.06).
Figure 4. Manhattan plot (A) and quantile–quantile (Q–Q) plot (B) for the association of all single nucleotide polymorphisms (SNPs) with the ratio of tail fin length to height (TL/H). Regional plots (C,D) of the genome-wide significant SNPs. Six SNPs reached genome-wide significant level. In plots A, C and D the red line shows the p-value threshold of 8.79 × 10−6 (−log10(p-value) > 5.06).
Genes 10 00418 g004
Genes 10 00418 g005 550
Figure 5. Comparison of the relative expression levels of genes in the candidate region in different tissues of the white normal sailfin molly (black bars) and the white dwarf sailfin molly (gray bars). (A) Relative expression level of Larp7; (B) relative expression level of ADAMTS like 1; (C) relative expression level of NLRP 12-like; (D) relative expression level of PPP3CA. Each bar represents the mean ± SE for each molly group (n = 8 for each group). * The values of gene expression levels are significantly different from each other (p < 0.05). ** The values of gene expression levels are significantly different from each other (p < 0.01).
Figure 5. Comparison of the relative expression levels of genes in the candidate region in different tissues of the white normal sailfin molly (black bars) and the white dwarf sailfin molly (gray bars). (A) Relative expression level of Larp7; (B) relative expression level of ADAMTS like 1; (C) relative expression level of NLRP 12-like; (D) relative expression level of PPP3CA. Each bar represents the mean ± SE for each molly group (n = 8 for each group). * The values of gene expression levels are significantly different from each other (p < 0.05). ** The values of gene expression levels are significantly different from each other (p < 0.01).
Genes 10 00418 g005
Table
Table 1. Phenotype differences between normal and dwarf white sailfin mollies.
Table 1. Phenotype differences between normal and dwarf white sailfin mollies.
CharacterNWMDWM
Torso Lenght (cm)3.6137 ± 0.0451 **2.1828 ± 0.0451
Height (cm)1.6264 ± 0.0211 *1.5666 ± 0.0211
TL/H2.2256 ± 0.0244 **1.4097 ± 0.0244
Note: NWM is the abbreviation of normal white sailfin molly and DWM means dwarf white sailfin molly. * Indicates a significant difference (p < 0.05) between NWMs and DWMs; ** Indicates an extremely significant difference (p < 0.01) between NWMs and DWMs.

Share and Cite

MDPI and ACS Style

Shao, F.; Liu, J.; Ren, M.; Li, J.; Bao, H.; Wu, C. Identification of Loci Controlling the Dwarfism Trait in the White Sailfin Molly (Poecilia latipinna) Using Genome-Wide Association Studies Based on Genotyping-By-Sequencing. Genes 2019, 10, 418. https://doi.org/10.3390/genes10060418

AMA Style

Shao F, Liu J, Ren M, Li J, Bao H, Wu C. Identification of Loci Controlling the Dwarfism Trait in the White Sailfin Molly (Poecilia latipinna) Using Genome-Wide Association Studies Based on Genotyping-By-Sequencing. Genes. 2019; 10(6):418. https://doi.org/10.3390/genes10060418

Chicago/Turabian Style

Shao, Fan, Jing Liu, Mengyuan Ren, Junying Li, Haigang Bao, and Changxin Wu. 2019. "Identification of Loci Controlling the Dwarfism Trait in the White Sailfin Molly (Poecilia latipinna) Using Genome-Wide Association Studies Based on Genotyping-By-Sequencing" Genes 10, no. 6: 418. https://doi.org/10.3390/genes10060418

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop