Regulatory variants at 2q33.1 confer schizophrenia risk by modulating distal gene TYW5 expression.

Genome-wide association studies have shown that genetic variants at 2q33.1 are strongly associated with schizophrenia. However, potential causal variants in this locus and their roles in schizophrenia remain unknown. Here, we identified two functional variants (rs796364 and rs281759) that disrupt CTCF, RAD21 and FOXP2 binding at 2q33.1. We systematically investigated the regulatory mechanisms of these two variants with serial experiments, including reporter gene assays and electrophoretic mobility shift assay. Intriguingly, these two single nucleotide polymorphisms physically interacted with TYW5 and showed the most significant associations with TYW5 expression in human brain. Consistently, CRISPR-Cas9-mediated genome editing confirmed the regulatory effect of the two single nucleotide polymorphisms on TYW5 expression. Additionally, expression analysis indicated that TYW5 was significantly upregulated in brains of schizophrenia cases compared with controls, suggesting that rs796364 and rs281759 might confer schizophrenia risk by modulating TYW5 expression. We over-expressed TYW5 in mouse neural stem cells and rat primary neurons to mimic its upregulation in schizophrenia and found significant alterations in the proliferation and differentiation of neural stem cells, as well as dendritic spine density following TYW5 overexpression, indicating its important roles in neurodevelopment and spine morphogenesis. Furthermore, we independently confirmed the association between rs796364 and schizophrenia in a Chinese cohort of 8202 subjects. Finally, transcriptome analysis revealed that TYW5 affected schizophrenia-associated pathways. These lines of evidence consistently revealed that rs796364 and rs281759 might contribute to schizophrenia risk by regulating the expression of TYW5, a gene whose expression dysregulation affects two important schizophrenia pathophysiological processes (i.e. neurodevelopment and dendritic spine formation).

Functional genomics elucidates regulatory mechanisms of Parkinson's disease-associated variants.

BACKGROUND: Genome-wide association studies (GWASs) have identified multiple risk loci for Parkinson's disease (PD). However, identifying the functional (or potential causal) variants in the reported risk loci and elucidating their roles in PD pathogenesis remain major challenges. To identify the potential causal (or functional) variants in the reported PD risk loci and to elucidate their regulatory mechanisms, we report a functional genomics study of PD. METHODS: We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) (from neuronal cells and human brain tissues) data and GWAS-identified single-nucleotide polymorphisms (SNPs) in PD risk loci. We then conducted a series of experiments and analyses to validate the regulatory effects of these (i.e., functional) SNPs, including reporter gene assays, allele-specific expression (ASE), transcription factor (TF) knockdown, CRISPR-Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis. RESULTS: We identified 44 SNPs (from 11 risk loci) affecting the binding of 12 TFs and we validated the regulatory effects of 15 TF binding-disrupting SNPs. In addition, we also identified the potential target genes regulated by these TF binding-disrupting SNPs through eQTL analysis. Finally, we showed that 4 eQTL genes of these TF binding-disrupting SNPs were dysregulated in PD cases compared with controls. CONCLUSION: Our study systematically reveals the gene regulatory mechanisms of PD risk variants (including widespread disruption of CTCF binding), generates the landscape of potential PD causal variants, and pinpoints promising candidate genes for further functional characterization and drug development.

Regulatory mechanisms of major depressive disorder risk variants.

Major depressive disorder (MDD) is one of the most prevalent psychiatric disorders and a leading cause of disability worldwide. Though recent genome-wide association studies (GWAS) have identified multiple risk variants for MDD, how these variants confer MDD risk remains largely unknown. Here we systematically characterize the regulatory mechanism of MDD risk variants using a functional genomics approach. By integrating chromatin immunoprecipitation sequencing (ChIP-Seq) (from human brain tissues or neuronal cells) and position weight matrix (PWM) data, we identified 34 MDD risk SNPs that disrupt the binding of 15 transcription factors (TFs). We verified the regulatory effect of the TF binding-disrupting SNPs with reporter gene assays, allelic-specific expression analysis, and CRISPR-Cas9-mediated genome editing. Expression quantitative trait loci (eQTL) analysis identified the target genes that might be regulated by these regulatory risk SNPs. Finally, we found that NEGR1 (regulated by the TF binding-disrupting MDD risk SNP rs3101339) was dysregulated in the brains of MDD cases compared with controls, implying that rs3101339 may confer MDD risk by affecting NEGR1 expression. Our findings reveal how genetic variants contribute to MDD risk by affecting TF binding and gene regulation. More importantly, our study identifies the potential MDD causal variants and their target genes, thus providing pivotal candidates for future mechanistic study and drug development.

Functional variants fine-mapping and gene function characterization provide insights into the role of ZNF323 in schizophrenia pathogenesis.

Schizophrenia is a severe mental disease characterized with positive symptoms, negative symptoms, and cognitive impairments. Although recent genome-wide association studies (GWASs) have identified over 145 risk loci for schizophrenia, pinpointing the causal variants and genes at the reported loci and elucidating their roles in schizophrenia remain major challenges. Here we identify a functional single-nucleotide polymorphism (SNP; rs213237) in ZNF323 promoter by using functional fine-mapping. We found that allelic differences at rs213237 affected the ZNF323 promoter activity significantly. Consistently, expression quantitative trait loci (eQTL) analysis showed that rs213237 was significantly associated with ZNF323 expression in diverse human brain tissues, suggesting that rs213237 may contribute to schizophrenia risk through regulating ZNF323 expression. Interestingly, we found that ZNF323 protein was localized in the nucleus and knockdown of ZNF323 in macaque neural stem cells (mNSCs) significantly impaired proliferation and survival of mNSCs. We further showed that stable knockdown of ZNF323 in SH-SY5Y cells resulted in significant decrease of the tyrosine hydroxylase (TH) protein expression. Finally, transcriptome analysis revealed that ZNF323 may regulate pivotal schizophrenia risk genes (including VIPR2 and NPY) and schizophrenia-associated pathways (including PI3K-AKT and NOTCH signaling pathways), suggesting that ZNF323 may be a major regulator of schizophrenia risk genes. Our study reveals how a genetic variant in ZNF323 promoter contributes to schizophrenia risk through regulating ZNF323 expression. More importantly, our findings demonstrate that ZNF323 may have a pivotal role in schizophrenia pathogenesis through regulating schizophrenia risk genes and schizophrenia-associated biological processes (including neurodevelopment, PI3K-AKT, and NOTCH signaling pathways).

Correction to: ZIKV infection effects changes in gene splicing, isoform composition and lncRNA expression in human neural progenitor cells.

In the original publication of this article [1], two grants from National Science Foundation of China and Yunnan Provincial Government (U1602226) and by National Science Foundation of China (2016YFC1200404) were omitted in the 'Funding' section. The correct 'Funding' section is below.

ZIKV infection effects changes in gene splicing, isoform composition and lncRNA expression in human neural progenitor cells.

BACKGROUND: The Zika virus (ZIKV) is a mosquito-borne flavivirus that causes microcephaly and Guillain-Barré syndrome in infected individuals. To obtain insights into the mechanism of ZIKV infection and pathogenesis, we analyzed the transcriptome of ZIKV infected human neural progenitor cells (hNPCs) for changes in alternative splicing (AS), gene isoform (ISO) composition and long noncoding RNAs (lncRNAs) expression. METHODS: We analyzed differentially expressed lncRNAs, AS, ISO from RNA-seq data in ZIKV infected hNPCs. RESULTS: We obtained 149 differentially expressed lncRNAs, including potential viral targets to modulate cellular processes such as cell cycle, apoptosis and immune response. The infection induced 262 cases of AS occurring in 229 genes, which were enriched in cell death, RNA processing, transport, and neuron development. Among 691 differentially expressed ISOs, upregulated ISOs were enriched in signaling, regulation of transcription, and amino acid biosynthesis, while downregulated ISOs were mostly enriched in cell cycle. Importantly, these analyses revealed specific links between ZIKV induced changes in cellular pathways and the type of changes in the host transcriptome, suggesting important regulatory mechanisms. CONCLUSIONS: Our analyses revealed candidate lncRNAs, AS events and ISOs which may function in ZIKV infection induced cell cycle disruption, apoptosis and attenuation of neurogenesis, and shed light on the roles of lncRNAs, AS and ISOs in virus-host interactions, and would facilitate future studies of ZIKV infection and pathogenesis.

Adaptive evolution of interleukin-3 (IL3), a gene associated with brain volume variation in general human populations.

Greatly expanded brain volume is one of the most characteristic traits that distinguish humans from other primates. Recent studies have revealed genes responsible for the dramatically enlarged human brain size (i.e., the microcephaly genes), and it has been well documented that many microcephaly genes have undergone accelerated evolution along the human lineage. In addition to being far larger than other primates, human brain volume is also highly variable in general populations. However, the genetic basis underlying human brain volume variation remains elusive and it is not known whether genes regulating human brain volume variation also have experienced positive selection. We have previously shown that genetic variants (near the IL3 gene) on 5q33 were significantly associated with brain volume in Chinese population. Here, we provide further evidence that support the significant association of genetic variants on 5q33 with brain volume. Bioinformatic analyses suggested that rs31480 is likely to be the causal variant among the studied SNPs. Molecular evolutionary analyses suggested that IL3 might have undergone positive selection in primates and humans. Neutrality tests further revealed signatures of positive selection of IL3 in Han Chinese and Europeans. Finally, extended haplotype homozygosity (EHH) and relative EHH analyses showed that the C allele of SNP rs31480 might have experienced recent positive selection in Han Chinese. Our results suggest that IL3 is an important genetic regulator for human brain volume variation and implied that IL3 might have experienced weak or modest positive selection in the evolutionary history of humans, which may be due to its contribution to human brain volume.

Functional prediction of differentially expressed lncRNAs in HSV-1 infected human foreskin fibroblasts.

BACKGROUND: One of the most important functions of long noncoding RNAs (lncRNAs) is to control protein coding gene transcription by acting locally in cis, or remotely in trans. Herpes Simplex Virus type I (HSV-1) latently infects over 80 % of the population, its reactivation from latency usually results in productive infections in human epithelial cells, and is responsible for the common cold sores and genital Herpes. HSV-1 productive infection leads to profound changes in the host cells, including the host transcriptome. However, how genome wide lncRNAs expressions are affected by the infection and how lncRNAs expression relates to protein coding gene expression have not been analyzed. METHODS: We analyzed differentially expressed lncRNAs and their potential targets from RNA-seq data in HSV-1 infected human foreskin fibroblast (HFF) cells. Based on correlations of expression patterns of differentially expressed protein-coding genes and lncRNAs, we predicted that these lncRNAs may regulate, either in cis or in trans, the expression of many cellular protein-coding genes. RESULTS: Here we analyzed HSV-1 infection induced, differentially expressed lncRNAs and predicted their target genes. We detected 208 annotated and 206 novel differentially expressed lncRNAs. Gene Ontology and Pathway enrichment analyses revealed potential lncRNA targets, including genes in chromatin assembly, genes in neuronal development and neurodegenerative diseases and genes in the immune response, such as Toll-like receptor signaling and RIG-I-like receptor signaling pathways. CONCLUSIONS: We found that differentially expressed lncRNAs may regulate the expression of many cellular protein-coding genes involved in pathways from native immunity to neuronal development, thus revealing important roles of lncRNAs in the regulation of host transcriptional programs in HSV-1 infected human cells.

Functional genomic analysis delineates regulatory mechanisms of GWAS-identified bipolar disorder risk variants.

BACKGROUND: Genome-wide association studies (GWASs) have identified multiple risk loci for bipolar disorder (BD). However, pinpointing functional (or causal) variants in the reported risk loci and elucidating their regulatory mechanisms remain challenging. METHODS: We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) data from human brain tissues (or neuronal cell lines) and position weight matrix (PWM) data to identify functional single-nucleotide polymorphisms (SNPs). Then, we verified the regulatory effects of these transcription factor (TF) binding-disrupting SNPs (hereafter referred to as "functional SNPs") through a series of experiments, including reporter gene assays, allele-specific expression (ASE) analysis, TF knockdown, CRISPR/Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis. Finally, we overexpressed PACS1 (whose expression was most significantly associated with the identified functional SNPs rs10896081 and rs3862386) in mouse primary cortical neurons to investigate if PACS1 affects dendritic spine density. RESULTS: We identified 16 functional SNPs (in 9 risk loci); these functional SNPs disrupted the binding of 7 TFs, for example, CTCF and REST binding was frequently disrupted. We then identified the potential target genes whose expression in the human brain was regulated by these functional SNPs through eQTL analysis. Of note, we showed dysregulation of some target genes of the identified TF binding-disrupting SNPs in BD patients compared with controls, and overexpression of PACS1 reduced the density of dendritic spines, revealing the possible biological mechanisms of these functional SNPs in BD. CONCLUSIONS: Our study identifies functional SNPs in some reported risk loci and sheds light on the regulatory mechanisms of BD risk variants. Further functional characterization and mechanistic studies of these functional SNPs and candidate genes will help to elucidate BD pathogenesis and develop new therapeutic approaches and drugs.

The schizophrenia-associated missense variant rs13107325 regulates dendritic spine density.

The missense variant rs13107325 (C/T, p.Ala391Thr) in SLC39A8 consistently showed robust association with schizophrenia in recent genome-wide association studies (GWASs), suggesting the potential pathogenicity of this non-synonymous risk variant. Nevertheless, how this missense variant confers schizophrenia risk remains unknown. Here we constructed a knock-in mouse model (by introducing a threonine at the 393th amino acid of mouse SLC39A8 (SLC39A8-p.393T), which corresponds to rs13107325 (p.Ala391Thr) of human SLC39A8) to explore the potential roles and biological effects of this missense variant in schizophrenia pathogenesis. We assessed multiple phenotypes and traits (associated with rs13107325) of the knock-in mice, including body and brain weight, concentrations of metal ions (including cadmium, zinc, manganese, and iron) transported by SLC39A8, blood lipids, proliferation and migration of neural stem cells (NSCs), cortical development, behaviors and cognition, transcriptome, dendritic spine density, and synaptic transmission. Many of the tested phenotypes did not show differences in SLC39A8-p.393T knock-in and wild-type mice. However, we found that zinc concentration in brain and blood of SLC39A8-p.393T knock-in mice was dysregulated compared with wild-types, validating the functionality of rs13107325. Further analysis indicated that cortical dendritic spine density of the SLC39A8-p.393T knock-in mice was significantly decreased compared with wild-types, indicating the important role of SLC39A8-p.393T in dendritic spine morphogenesis. These results indicated that SLC39A8-p.393T knock-in resulted in decreased dendritic spine density, thus mimicking the dendritic spine pathology observed in schizophrenia. Our study indicates that rs13107325 might confer schizophrenia risk by regulating zinc concentration and dendritic spine density, a featured characteristic that was frequently reported to be decreased in schizophrenia.

A functional missense variant in ITIH3 affects protein expression and neurodevelopment and confers schizophrenia risk in the Han Chinese population.

The Psychiatric Genomics Consortium (PGC) has recently identified 10 potential functional coding variants for schizophrenia. However, how these coding variants confer schizophrenia risk remains largely unknown. Here, we investigate the associations between eight potential functional coding variants identified by PGC and schizophrenia in a large Han Chinese sample (n = 4022 cases and 9270 controls). Among the eight tested single nucelotide polymorphisms (SNPs), rs3617 (a missense variant, p.K315Q in the ITIH3 gene) showed genome-wide significant association with schizophrenia in the Han Chinese population (P = 8.36 × 10-16), with the same risk allele as in PGC. Interestingly, rs3617 is located in a genomic region that is highly evolutionarily conserved, and its schizophrenia risk allele (C allele) was associated with lower ITIH3 mRNA and protein expression. Intriguingly, mouse neural stem cells stably overexpressing ITIH3 with different alleles of rs3617 exhibited significant differences in proliferation, migration, and differentiation, suggesting the impact of rs3617 on neurodevelopment. Subsequent transcriptome analysis found that the differentially expressed genes in neural stem cells stably overexpressing different alleles of rs3617 were significantly enriched in schizophrenia-related pathways, including cell adhesion, synapse assembly, MAPK and PI3K-AKT pathways. Our study provides convergent lines of evidence suggesting that rs3617 in ITIH3 likely affects protein function and neurodevelopment and thereby confers risk of schizophrenia.

Functional genomics reveal gene regulatory mechanisms underlying schizophrenia risk.

Genome-wide association studies (GWASs) have identified over 180 independent schizophrenia risk loci. Nevertheless, how the risk variants in the reported loci confer schizophrenia susceptibility remains largely unknown. Here we systematically investigate the gene regulatory mechanisms underpinning schizophrenia risk through integrating data from functional genomics (including 30 ChIP-Seq experiments) and position weight matrix (PWM). We identify 132 risk single nucleotide polymorphisms (SNPs) that disrupt transcription factor binding and we find that 97 of the 132 TF binding-disrupting SNPs are associated with gene expression in human brain tissues. We validate the regulatory effect of some TF binding-disrupting SNPs with reporter gene assays (9 SNPs) and allele-specific expression analysis (10 SNPs). Our study reveals gene regulatory mechanisms affected by schizophrenia risk SNPs (including widespread disruption of POLR2A and CTCF binding) and identifies target genes for mechanistic studies and drug development. Our results can be accessed and visualized at SZDB database ( http://www.szdb.org/ ).

The integrated landscape of causal genes and pathways in schizophrenia.

Genome-wide association studies (GWAS) have identified more than 100 loci that show robust association with schizophrenia risk. However, due to the complexity of linkage disequilibrium and gene regulatory, it is challenging to pinpoint the causal genes at the risk loci and translate the genetic findings from GWAS into disease mechanism and clinical treatment. Here we systematically predicted the plausible candidate causal genes for schizophrenia at genome-wide level. We utilized different approaches and strategies to predict causal genes for schizophrenia, including Sherlock, SMR, DAPPLE, Prix Fixe, NetWAS, and DEPICT. By integrating the results from different prediction approaches, we identified six top candidates that represent promising causal genes for schizophrenia, including CNTN4, GATAD2A, GPM6A, MMP16, PSMA4, and TCF4. Besides, we also identified 35 additional high-confidence causal genes for schizophrenia. The identified causal genes showed distinct spatio-temporal expression patterns in developing and adult human brain. Cell-type-specific expression analysis indicated that the expression level of the predicted causal genes was significantly higher in neurons compared with oligodendrocytes and microglia (P < 0.05). We found that synaptic transmission-related genes were significantly enriched among the identified causal genes (P < 0.05), providing further support for the dysregulation of synaptic transmission in schizophrenia. Finally, we showed that the top six causal genes are dysregulated in schizophrenia cases compared with controls and knockdown of these genes impaired the proliferation of neuronal cells. Our study depicts the landscape of plausible schizophrenia causal genes for the first time. Further genetic and functional validation of these genes will provide mechanistic insights into schizophrenia pathogenesis and may facilitate to provide potential targets for future therapeutics and diagnostics.

Comprehensive integrative analyses identify GLT8D1 and CSNK2B as schizophrenia risk genes.

Recent genome-wide association studies (GWAS) have identified multiple risk loci that show strong associations with schizophrenia. However, pinpointing the potential causal genes at the reported loci remains a major challenge. Here we identify candidate causal genes for schizophrenia using an integrative genomic approach. Sherlock integrative analysis shows that ALMS1, GLT8D1, and CSNK2B are schizophrenia risk genes, which are validated using independent brain expression quantitative trait loci (eQTL) data and integrative analysis method (SMR). Consistently, gene expression analysis in schizophrenia cases and controls further supports the potential role of these three genes in the pathogenesis of schizophrenia. Finally, we show that GLT8D1 and CSNK2B knockdown promote the proliferation and inhibit the differentiation abilities of neural stem cells, and alter morphology and synaptic transmission of neurons. These convergent lines of evidence suggest that the ALMS1, CSNK2B, and GLT8D1 genes may be involved in pathophysiology of schizophrenia.

Author Correction: Comprehensive integrative analyses identify GLT8D1 and CSNK2B as schizophrenia risk genes.

In the original version of this Article, the affiliation details for Qiushuo Shen incorrectly omitted 'Kunming College of Life Science, University of Chinese Academy of Sciences, Kunming, 650204, China'. This has now been corrected in both the PDF and HTML versions of the Article.

Cellular responses to HSV-1 infection are linked to specific types of alterations in the host transcriptome.

Pathogen invasion triggers a number of cellular responses and alters the host transcriptome. Here we report that the type of changes to cellular transcriptome is related to the type of cellular functions affected by lytic infection of Herpes Simplex Virus type I in Human primary fibroblasts. Specifically, genes involved in stress responses and nuclear transport exhibited mostly changes in alternative polyadenylation (APA), cell cycle genes showed mostly alternative splicing (AS) changes, while genes in neurogenesis, rarely underwent these changes. Transcriptome wide, the infection resulted in 1,032 cases of AS, 161 incidences of APA, 1,827 events of isoform changes, and up regulation of 596 genes and down regulations of 61 genes compared to uninfected cells. Thus, these findings provided important and specific links between cellular responses to HSV-1 infection and the type of alterations to the host transcriptome, highlighting important roles of RNA processing in virus-host interactions.

Identification of SLC25A37 as a major depressive disorder risk gene.

Major depressive disorder (MDD) is one of the most prevalent and disabling mental disorders, but the genetic etiology remains largely unknown. We performed a meta-analysis (14,543 MDD cases and 14,856 controls) through combining the GWAS data from the Major Depressive Disorder Working Group of the Psychiatric GWAS Consortium and the CONVERGE consortium and identified seven SNPs (four of them located in the downstream of SCL25A37) that showed suggestive associations (P < 5.0 × 10-7) with MDD. Systematic integration (Sherlock integrative analysis) of brain eQTL and GWAS meta-analysis identified SCL25A37 as a novel MDD risk gene (P = 2.22 × 10-6). A cis SNP (rs6983724, ∼28 kb downstream of SCL25A37) showed significant association with SCL25A37 expression (P = 1.19 × 10-9) and suggestive association with MDD (P = 1.65 × 10-7). We validated the significant association between rs6983724 and SCL25A37 expression in independent expression datasets. Finally, we found that SCL25A37 is significantly down-regulated in hippocampus and blood of MDD patients (P = 3.49 × 10-3 and P = 2.66 × 10-13, respectively). Our findings implicate that the SCL25A37 is a MDD susceptibility gene whose expression may influence MDD risk. The consistent down-regulation of SCL25A37 in MDD patients in three independent samples suggest that SCL25A37 may be used as a potential biomarker for MDD diagnosis. Further functional characterization of SCL25A37 may provide a potential target for future therapeutics and diagnostics.

Recent Positive Selection Drives the Expansion of a Schizophrenia Risk Nonsynonymous Variant at SLC39A8 in Europeans.

Natural selection has played important roles in optimizing complex human adaptations. However, schizophrenia poses an evolutionary paradox during human evolution, as the illness has strongly negative effects on fitness, but persists with a prevalence of ~0.5% across global populations. Recent studies have identified numerous risk variations in diverse populations, which might be able to explain the stable and high rate of schizophrenia morbidity in different cultures and regions, but the questions about why the risk alleles derived and maintained in human gene pool still remain unsolved. Here, we studied the evolutionary pattern of a schizophrenia risk variant rs13107325 (P < 5.0 × 10(-8) in Europeans) in the SLC39A8 gene. We found the SNP is monomorphic in Asians and Africans with risk (derived) T-allele totally absent, and further evolutionary analyses showed the T-allele has experienced recent positive selection in Europeans. Subsequent exploratory analyses implicated that the colder environment in Europe was the likely selective pressures, ie, when modern humans migrated "out of Africa" and moved to Europe mainland (a colder and cooler continent than Africa), new alleles derived due to positive selection and protected humans from risk of hypertension and also helped them adapt to the cold environment. The hypothesis was supported by our pleiotropic analyses with hypertension and energy intake as well as obesity in Europeans. Our data thus provides an intriguing example to illustrate a possible mechanism for maintaining schizophrenia risk alleles in the human gene pool, and further supported that schizophrenia is likely a product caused by pleiotropic effect during human evolution.

Comparative population genomics reveals genetic basis underlying body size of domestic chickens.

Body size is the most important economic trait for animal production and breeding. Several hundreds of loci have been reported to be associated with growth trait and body weight in chickens. The loci are mapped to large genomic regions due to the low density and limited number of genetic markers in previous studies. Herein, we employed comparative population genomics to identify genetic basis underlying the small body size of Yuanbao chicken (a famous ornamental chicken) based on 89 whole genomes. The most significant signal was mapped to the BMP10 gene, whose expression was upregulated in the Yuanbao chicken. Overexpression of BMP10 induced a significant decrease in body length by inhibiting angiogenic vessel development in zebrafish. In addition, three other loci on chromosomes 1, 2, and 24 were also identified to be potentially involved in the development of body size. Our results provide a paradigm shift in identification of novel loci controlling body size variation, availing a fast and efficient strategy. These loci, particularly BMP10, add insights into ongoing research of the evolution of body size under artificial selection and have important implications for future chicken breeding.