RESUMO
We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or â¼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.
Assuntos
Genoma , Primatas , Animais , Humanos , Sequência de Bases , Primatas/classificação , Primatas/genética , Evolução Biológica , Análise de Sequência de DNA , Variação Estrutural do GenomaRESUMO
Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size1. As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions2,3. Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome4,5. We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA.
Assuntos
Centrômero , Evolução Molecular , Variação Genética , Animais , Humanos , Centrômero/genética , Centrômero/metabolismo , Proteína Centromérica A/metabolismo , Metilação de DNA/genética , DNA Satélite/genética , Cinetocoros/metabolismo , Macaca/genética , Pan troglodytes/genética , Polimorfismo de Nucleotídeo Único/genética , Pongo/genética , Masculino , Feminino , Padrões de Referência , Imunoprecipitação da Cromatina , Haplótipos , Mutação , Amplificação de Genes , Alinhamento de Sequência , Cromatina/genética , Cromatina/metabolismo , Especificidade da EspécieRESUMO
Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.
Assuntos
Hominidae , Cromossomo X , Cromossomo Y , Animais , Feminino , Masculino , Gorilla gorilla/genética , Hominidae/genética , Hominidae/classificação , Hylobatidae/genética , Pan paniscus/genética , Pan troglodytes/genética , Filogenia , Pongo abelii/genética , Pongo pygmaeus/genética , Telômero/genética , Cromossomo X/genética , Cromossomo Y/genética , Evolução Molecular , Variações do Número de Cópias de DNA/genética , Humanos , Espécies em Perigo de Extinção , Padrões de ReferênciaRESUMO
TBC1D3 is a primate-specific gene family that has expanded in the human lineage and has been implicated in neuronal progenitor proliferation and expansion of the frontal cortex. The gene family and its expression have been challenging to investigate because it is embedded in high-identity and highly variable segmental duplications. We sequenced and assembled the gene family using long-read sequencing data from 34 humans and 11 non-human primate species. Our analysis shows that this particular gene family has independently duplicated in at least five primate lineages, and the duplicated loci are enriched at sites of large-scale chromosomal rearrangements on Chromosome 17. We find that all human copy number variation maps to two distinct clusters located at Chromosome 17q12 and that humans are highly structurally variable at this locus, differing by as many as 20 copies and ~1 Mbp in length depending on haplotypes. We also show evidence of positive selection, as well as a significant change in the predicted human TBC1D3 protein sequence. Lastly, we find that, despite multiple duplications, human TBC1D3 expression is limited to a subset of copies and, most notably, from a single paralog group: TBC1D3-CDKL These observations may help explain why a gene potentially important in cortical development can be so variable in the human population.
RESUMO
The deleted in azoospermia like (DAZL) is required for germ cells development and maintenance. In chickens, the mRNA and protein of DAZL, a representative maternally inherited germ plasm factor, are detected in the germ plasm of oocyte, zygote, and all stages of the intrauterine embryos. However, it is still insufficient to explain the origin and specification process of chicken germ cells, because the stage at which the zygotic transcription of DAZL occurs and the stage at which the maternal DAZL RNA/protein clears have not yet been fully identified. Moreover, a comprehensive understanding of the expression of DAZL interacting genes during the germ cells specification and development and zygotic genome activation (ZGA) is lacking in chickens. In this study, we identified a set of DAZL interacting genes in chickens using in silico prediction method. Then, we analyzed the whole-transcriptome sequencing (WTS)-based expression of DAZL and its interacting genes in the chicken oocyte, zygote, and Eyal-Giladi and Kochav (EGK) stage embryos (EGK.I to EGK.X). In the results, DAZL transcripts are increased in the zygote (onset of transcription), maintained the increased level until EGK.VI, and decreased from EGK.VIII (possible clearance of maternal RNAs). Among the DAZL interacting genes, most of them are increased either at 1st ZGA or 2nd ZGA, indicating their involvement in germ cells specification and development.
Assuntos
Galinhas/genética , Proteínas de Ligação a RNA/genética , Animais , Diferenciação Celular/fisiologia , Embrião de Galinha , Galinhas/crescimento & desenvolvimento , Galinhas/metabolismo , Epistasia Genética , Feminino , Regulação da Expressão Gênica no Desenvolvimento , Células Germinativas/crescimento & desenvolvimento , Células Germinativas/metabolismo , Proteínas de Ligação a RNA/metabolismo , Análise de Sequência de RNA/métodos , Análise de Sequência de RNA/veterinária , ZigotoRESUMO
Tandem repeats (TRs) are widespread in the genomes of all living organisms. In eukaryotes, they are found in both coding and noncoding regions and have potential roles in the regulation of cellular processes such as transcription, translation and in the modification of protein structure. Recent studies have highlighted TRs as a key regulator of gene expression and a potential contributor to human evolution. Thus, TRs are emerging as an important source of variation that can result in differential gene expression at intra- and inter-species levels. In this study, we performed a genome-wide survey to identify TRs that have emerged in the human lineage. We further examined these loci to explore their potential functional significance for human evolution. We identified 152 human-specific TR (HSTR) loci containing a repeat unit of more than ten bases, with most of them showing a repeat count of two. Gene set enrichment analysis showed that HSTR-associated genes were associated with biological functions in brain development and synapse function. In addition, we compared gene expression of human HSTR loci with orthologues from non-human primates (NHP) in seven different tissues. Strikingly, the expression level of HSTR-associated genes in brain tissues was significantly higher in human than in NHP. These results suggest the possibility that de novo emergence of TRs could have resulted in altered gene expression in humans within a short-time frame and contributed to the rapid evolution of human brain function.
Assuntos
Encéfalo/metabolismo , Regulação da Expressão Gênica , Especificidade de Órgãos/genética , Sequências de Repetição em Tandem/genética , Animais , Sequência de Bases , Evolução Molecular , Genoma Humano/genética , Humanos , Taxa de Mutação , Primatas/genética , Homologia de Sequência do Ácido NucleicoRESUMO
BACKGROUND: Indigenous cattle in Africa have adapted to various local environments to acquire superior phenotypes that enhance their survival under harsh conditions. While many studies investigated the adaptation of overall African cattle, genetic characteristics of each breed have been poorly studied. RESULTS: We performed the comparative genome-wide analysis to assess evidence for subspeciation within species at the genetic level in trypanotolerant N'Dama cattle. We analysed genetic variation patterns in N'Dama from the genomes of 101 cattle breeds including 48 samples of five indigenous African cattle breeds and 53 samples of various commercial breeds. Analysis of SNP variances between cattle breeds using wMI, XP-CLR, and XP-EHH detected genes containing N'Dama-specific genetic variants and their potential associations. Functional annotation analysis revealed that these genes are associated with ossification, neurological and immune system. Particularly, the genes involved in bone formation indicate that local adaptation of N'Dama may engage in skeletal growth as well as immune systems. CONCLUSIONS: Our results imply that N'Dama might have acquired distinct genotypes associated with growth and regulation of regional diseases including trypanosomiasis. Moreover, this study offers significant insights into identifying genetic signatures for natural and artificial selection of diverse African cattle breeds.
Assuntos
Bovinos/genética , Bovinos/parasitologia , Genômica , Polimorfismo de Nucleotídeo Único , Trypanosoma/fisiologia , Animais , Doenças dos Bovinos/imunologia , Doenças dos Bovinos/parasitologia , Códon sem Sentido , Resistência à Doença/genética , Mutação de Sentido Incorreto , Especificidade da EspécieRESUMO
Segmental duplications (SDs) contribute significantly to human disease, evolution, and diversity yet have been difficult to resolve at the sequence level. We present a population genetics survey of SDs by analyzing 170 human genome assemblies where the majority of SDs are fully resolved using long-read sequence assembly. Excluding the acrocentric short arms, we identify 173.2 Mbp of duplicated sequence (47.4 Mbp not present in the telomere-to-telomere reference) distinguishing fixed from structurally polymorphic events. We find that intrachromosomal SDs are among the most variable with rare events mapping near their progenitor sequences. African genomes harbor significantly more intrachromosomal SDs and are more likely to have recently duplicated gene families with higher copy number when compared to non-African samples. A comparison to a resource of 563 million full-length Iso-Seq reads identifies 201 novel, potentially protein-coding genes corresponding to these copy number polymorphic SDs.
RESUMO
TBC1D3 is a primate-specific gene family that has expanded in the human lineage and has been implicated in neuronal progenitor proliferation and expansion of the frontal cortex. The gene family and its expression have been challenging to investigate because it is embedded in high-identity and highly variable segmental duplications. We sequenced and assembled the gene family using long-read sequencing data from 34 humans and 11 nonhuman primate species. Our analysis shows that this particular gene family has independently duplicated in at least five primate lineages, and the duplicated loci are enriched at sites of large-scale chromosomal rearrangements on chromosome 17. We find that most humans vary along two TBC1D3 clusters where human haplotypes are highly variable in copy number, differing by as many as 20 copies, and structure (structural heterozygosity 90%). We also show evidence of positive selection, as well as a significant change in the predicted human TBC1D3 protein sequence. Lastly, we find that, despite multiple duplications, human TBC1D3 expression is limited to a subset of copies and, most notably, from a single paralog group: TBC1D3-CDKL. These observations may help explain why a gene potentially important in cortical development can be so variable in the human population.
RESUMO
We present haplotype-resolved reference genomes and comparative analyses of six ape species, namely: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. We achieve chromosome-level contiguity with unparalleled sequence accuracy (<1 error in 500,000 base pairs), completely sequencing 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, providing more in-depth evolutionary insights. Comparative analyses, including human, allow us to investigate the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference. This includes newly minted gene families within lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes, and subterminal heterochromatin. This resource should serve as a definitive baseline for all future evolutionary studies of humans and our closest living ape relatives.
RESUMO
The biological role of the repetitive DNA sequences in the human genome remains an outstanding question. Recent long-read human genome assemblies have allowed us to identify a function for one of these repetitive regions. We have uncovered a tandem array of conserved primate-specific retrogenes encoding the protein Elongin A3 (ELOA3), a homolog of the RNA polymerase II (RNAPII) elongation factor Elongin A (ELOA). Our genomic analysis shows that the ELOA3 gene cluster is conserved among primates and the number of ELOA3 gene repeats is variable in the human population and across primate species. Moreover, the gene cluster has undergone concerted evolution and homogenization within primates. Our biochemical studies show that ELOA3 functions as a promoter-associated RNAPII pause-release elongation factor with distinct biochemical and functional features from its ancestral homolog, ELOA. We propose that the ELOA3 gene cluster has evolved to fulfil a transcriptional regulatory function unique to the primate lineage that can be targeted to regulate cellular hyperproliferation.
Assuntos
Fatores de Alongamento de Peptídeos , RNA Polimerase II , Animais , Humanos , RNA Polimerase II/genética , RNA Polimerase II/metabolismo , Fatores de Alongamento de Peptídeos/genética , Primatas/genética , Elonguina/genética , Família Multigênica , Sequências de Repetição em Tandem/genéticaRESUMO
Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.
RESUMO
BACKGROUND: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic began in 2019 but it remains as a serious threat today. To reduce and prevent spread of the virus, multiple vaccines have been developed. Despite the efforts in developing vaccines, Omicron strain of the virus has recently been designated as a variant of concern (VOC) by the World Health Organization (WHO). OBJECTIVE: To develop a vaccine candidate against Omicron strain (B.1.1.529, BA.1) of the SARS-CoV-19. METHODS: We applied reverse vaccinology methods for BA.1 and BA.2 as the vaccine target and a control, respectively. First, we predicted MHC I, MHC II and B cell epitopes based on their viral genome sequences. Second, after estimation of antigenicity, allergenicity and toxicity, a vaccine construct was assembled and tested for physicochemical properties and solubility. Third, AlphaFold2, RaptorX and RoseTTAfold servers were used to predict secondary structures and 3D structures of the vaccine construct. Fourth, molecular docking analysis was performed to test binding of our construct with angiotensin converting enzyme 2 (ACE2). Lastly, we compared mutation profiles on the epitopes between BA.1, BA.2, and wild type to estimate the efficacy of the vaccine. RESULTS: We collected a total of 10 MHC I, 9 MHC II and 5 B cell epitopes for the final vaccine construct for Omicron strain. All epitopes were predicted to be antigenic, non-allergenic and non-toxic. The construct was estimated to have proper stability and solubility. The best modelled tertiary structures were selected for molecular docking analysis with ACE2 receptor. CONCLUSIONS: These results suggest the potential efficacy of our newly developed vaccine construct as a novel vaccine candidate against Omicron strain of the coronavirus.
Assuntos
COVID-19 , Vacinas Virais , Enzima de Conversão de Angiotensina 2 , COVID-19/prevenção & controle , Vacinas contra COVID-19 , Epitopos de Linfócito B/química , Epitopos de Linfócito B/genética , Epitopos de Linfócito T/química , Epitopos de Linfócito T/genética , Humanos , Simulação de Acoplamento Molecular , SARS-CoV-2/genética , Desenvolvimento de Vacinas , Vacinologia/métodos , Vacinas Virais/química , Vacinas Virais/genéticaRESUMO
The little skate Leucoraja erinacea, a cartilaginous fish, displays pelvic fin driven walking-like behavior using genetic programs and neuronal subtypes similar to those of land vertebrates. However, mechanistic studies on little skate motor circuit development have been limited, due to a lack of high-quality reference genome. Here, we generated an assembly of the little skate genome, with precise gene annotation and structures, which allowed post-genome analysis of spinal motor neurons (MNs) essential for locomotion. Through interspecies comparison of mouse, skate and chicken MN transcriptomes, shared and divergent gene expression profiles were identified. Comparison of accessible chromatin regions between mouse and skate MNs predicted shared transcription factor (TF) motifs with divergent ones, which could be used for achieving differential regulation of MN-expressed genes. A greater number of TF motif predictions were observed in MN-expressed genes in mouse than in little skate. These findings suggest conserved and divergent molecular mechanisms controlling MN development of vertebrates during evolution, which might contribute to intricate gene regulatory networks in the emergence of a more sophisticated motor system in tetrapods.
Assuntos
Rajidae , Animais , Camundongos , Cromatina/metabolismo , Neurônios Motores , Rajidae/genética , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismo , Caminhada , GenomaRESUMO
Water-to-land transition has been independently evolved in multiple vertebrate lineages including the most recent common ancestor of tetrapod and multiple fish clades, and among them, mudskippers uniquely adapted to the mudflat. Even though physiological and morphological adaptation of mudskippers is thought to resemble that of the ancestral tetrapod, it is unclear if they share genome-wide evolutionary signatures. To detect potential signatures of positive selection in mudskipper and tetrapods, we analyzed 4118 singleton orthologues of terrestrial tetrapods, coelacanth, mudskipper, and fully aquatic fishes. Among positively selected genes identified in mudskipper and tetrapod lineages, genes involved in immune responses, mitochondrial oxidative phosphorylation, and kidney development were detected. On the other hand, tetrapod-specific and mudskipper-specific positively selected genes were functionally enriched for DNA repair processes, which could be associated with higher exposure to UV light. We also performed gene family analysis and discovered convergent contraction of eight gene families, including ßγ-crystallin coding genes in both tetrapod and mudskipper lineages. Findings of this study suggest the similar genetic adaptation against environmental constraints between the ancient tetrapod and mudskippers for their land adaptation.
RESUMO
Stress resistance is an important trait expected of lactic acid bacteria used in food manufacturing. Among the various sources of stress, high temperature is a key factor that interrupts bacterial growth. In this regards, constant efforts are made for the development of heat-resistant strains, but few studies were done accompanying genomic analysis to identify the causal factors of the resistance mechanisms. Furthermore, it is also thought that tolerance to multiple stresses are equally important. Herein, we isolated one Enterococcus faecium strain named BIOPOP-3 and completed a full-length genome sequence. Using this strain, a two-step adaptive laboratory evolution (ALE) method was applied to obtain a heat-resistant strain, BIOPOP-3 ALE. After sequencing the whole genome, we compared the two full-length sequences and identified one non-synonymous variant and four indel variants that could potentially confer heat resistance, which were technically validated by resequencing. We experimentally verified that the evolved strain was significantly enhanced in not only heat resistance but also acid and bile resistance. We demonstrated that the developed heat-resistant strain can be applied in animal feed manufacturing processes. The multi-stress-resistant BIOPOP-3 ALE strain developed in this study and the two-step ALE method are expected to be widely applied in industrial and academic fields. In addition, we expect that the identified variants which occurred specifically in heat-resistant strain will enhance molecular biological understanding and be broadly applied to the biological engineering field.
RESUMO
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
RESUMO
Cattle pastoralism plays a central role in human livelihood in Africa. However, the genetic history of its success remains unknown. Here, through whole-genome sequence analysis of 172 indigenous African cattle from 16 breeds representative of the main cattle groups, we identify a major taurine × indicine cattle admixture event dated to circa 750-1,050 yr ago, which has shaped the genome of today's cattle in the Horn of Africa. We identify 16 loci linked to African environmental adaptations across crossbred animals showing an excess of taurine or indicine ancestry. These include immune-, heat-tolerance- and reproduction-related genes. Moreover, we identify one highly divergent locus in African taurine cattle, which is putatively linked to trypanotolerance and present in crossbred cattle living in trypanosomosis-infested areas. Our findings indicate that a combination of past taurine and recent indicine admixture-derived genetic resources is at the root of the present success of African pastoralism.
Assuntos
Adaptação Fisiológica/genética , Cruzamento , Bovinos , Genoma , Sequenciamento Completo do Genoma , África , Alelos , Animais , Bovinos/genética , Genótipo , Temperatura Alta/efeitos adversos , Mosaicismo , Polimorfismo de Nucleotídeo Único , Reprodução/genética , Sequenciamento Completo do Genoma/veterináriaRESUMO
Diseases prediction has been performed by machine learning approaches with various biological data. One of the representative data is the gut microbial community, which interacts with the host's immune system. The abundance of a few microorganisms has been used as markers to predict diverse diseases. In this study, we hypothesized that multi-classification using machine learning approach could distinguish the gut microbiome from following six diseases: multiple sclerosis, juvenile idiopathic arthritis, myalgic encephalomyelitis/chronic fatigue syndrome, acquired immune deficiency syndrome, stroke and colorectal cancer. We used the abundance of microorganisms at five taxonomy levels as features in 696 samples collected from different studies to establish the best prediction model. We built classification models based on four multi-class classifiers and two feature selection methods including a forward selection and a backward elimination. As a result, we found that the performance of classification is improved as we use the lower taxonomy levels of features; the highest performance was observed at the genus level. Among four classifiers, LogitBoost-based prediction model outperformed other classifiers. Also, we suggested the optimal feature subsets at the genus-level obtained by backward elimination. We believe the selected feature subsets could be used as markers to distinguish various diseases simultaneously. The finding in this study suggests the potential use of selected features for the diagnosis of several diseases.
Assuntos
Doença/classificação , Microbioma Gastrointestinal/genética , Metagenômica/métodos , Algoritmos , Doença/genética , Previsões/métodos , Humanos , Aprendizado de Máquina , Máquina de Vetores de SuporteRESUMO
Red sea bream, a popular fish resource in Korea and Japan, is being bred in fish farms of the two countries. It is hypothesized that the genomes of red sea bream are influenced by decades of artificial selection. This study investigates the impact of artificial selection on genomes of red sea bream. Whole genome sequencing was conducted for 40 samples of red sea bream either from Ehime, Nagasaki and Tongyeong fish farms or from the wild. Population stratification based on whole genome data was investigated and the genomic regions of fish farm populations under selection were identified using XP-EHH and relative nucleotide diversity. Gene ontology analysis revealed that different functions were enriched in different fish farms. In conclusion, this study highlights the difference between independently cultured red sea bream populations by showing that influence of artificial selection acted upon completely different genes related to different functions including metabolic and developmental processes.