Three amphioxus reference genomes reveal gene and chromosome evolution of chordates.

The slow-evolving invertebrate amphioxus has an irreplaceable role in advancing our understanding of the vertebrate origin and innovations. Here we resolve the nearly complete chromosomal genomes of three amphioxus species, one of which best recapitulates the 17 chordate ancestor linkage groups. We reconstruct the fusions, retention, or rearrangements between descendants of whole-genome duplications, which gave rise to the extant microchromosomes likely existed in the vertebrate ancestor. Similar to vertebrates, the amphioxus genome gradually establishes its three-dimensional chromatin architecture at the onset of zygotic activation and forms two topologically associated domains at the Hox gene cluster. We find that all three amphioxus species have ZW sex chromosomes with little sequence differentiation, and their putative sex-determining regions are nonhomologous to each other. Our results illuminate the unappreciated interspecific diversity and developmental dynamics of amphioxus genomes and provide high-quality references for understanding the mechanisms of chordate functional genome evolution.

Creating a functional single-chromosome yeast.

Eukaryotic genomes are generally organized in multiple chromosomes. Here we have created a functional single-chromosome yeast from a Saccharomyces cerevisiae haploid cell containing sixteen linear chromosomes, by successive end-to-end chromosome fusions and centromere deletions. The fusion of sixteen native linear chromosomes into a single chromosome results in marked changes to the global three-dimensional structure of the chromosome due to the loss of all centromere-associated inter-chromosomal interactions, most telomere-associated inter-chromosomal interactions and 67.4% of intra-chromosomal interactions. However, the single-chromosome and wild-type yeast cells have nearly identical transcriptome and similar phenome profiles. The giant single chromosome can support cell life, although this strain shows reduced growth across environments, competitiveness, gamete production and viability. This synthetic biology study demonstrates an approach to exploration of eukaryote evolution with respect to chromosome structure and function.

Phenotypic correlations of carpal gland diverticular number with production traits and its genome-wide association analysis in multiple pig populations.

Pig carpal glands play crucial roles in territorial recognition, reproductive behavior, and information exchange; however, their effects on production traits and underlying genetic mechanisms remain unclear. In this study, 1028 pigs from six populations were counted for the carpal gland diverticular numbers (CGDNs) on the left (CGDNL) and right (CGDNR) legs, and their carcass and meat quality traits were assessed. The CGDNs were significantly different among the populations, and Licha Black pigs had a lower CGDN than the Bama Xiang breed. It was also significantly different between sexes, with males having more diverticula than females (p ≤ 0.0391). Moreover, the number was asymmetric, with CGDNR being significantly higher than CGDNL. Notably, CGDNs was significantly correlated with each other in phenotype and genetics and with 24-h pH, 24-h meat color score, 24-h marbling score, fat content, moisture content, sodium salt content, and saturated fatty acid content in phenotype. Furthermore, genome-wide association analyses identified seven SNPs in association with CGDNs at a 5% genome-wide significance level, all of which were located in a 1.78-Mb (35.347-37.129 Mb) region on chromosome 1. CNC10010837 and CNC10010840 were the top SNPs: both had an additive effect of 0.789 ± 0.120 on CGDNR with p = 8.31E-10. These findings provide important insights into the functions and underlying genetic mechanisms of swine carpal glands.

Genomic prediction in pigs using data from a commercial crossbred population: insights from the Duroc x (Landrace x Yorkshire) three-way crossbreeding system.

BACKGROUND: Genomic selection is widely applied for genetic improvement in livestock crossbreeding systems to select excellent nucleus purebred (PB) animals and to improve the performance of commercial crossbred (CB) animals. Most current predictions are based solely on PB performance. Our objective was to explore the potential application of genomic selection of PB animals using genotypes of CB animals with extreme phenotypes in a three-way crossbreeding system as the reference population. Using real genotyped PB as ancestors, we simulated the production of 100,000 pigs for a Duroc x (Landrace x Yorkshire) DLY crossbreeding system. The predictive performance of breeding values of PB animals for CB performance using genotypes and phenotypes of (1) PB animals, (2) DLY animals with extreme phenotypes, and (3) random DLY animals for traits of different heritabilities ([Formula: see text] = 0.1, 0.3, and 0.5) was compared across different reference population sizes (500 to 6500) and prediction models (genomic best linear unbiased prediction (GBLUP) and Bayesian sparse linear mixed model (BSLMM)). RESULTS: Using a reference population consisting of CB animals with extreme phenotypes showed a definite predictive advantage for medium- and low-heritability traits and, in combination with the BSLMM model, significantly improved selection response for CB performance. For high-heritability traits, the predictive performance of a reference population of extreme CB phenotypes was comparable to that of PB phenotypes when the effect of the genetic correlation between PB and CB performance ([Formula: see text]) on the accuracy obtained with a PB reference population was considered, and the former could exceed the latter if the reference size was large enough. For the selection of the first and terminal sires in a three-way crossbreeding system, prediction using extreme CB phenotypes outperformed the use of PB phenotypes, while the optimal design of the reference group for the first dam depended on the percentage of individuals from the corresponding breed that the PB reference data comprised and on the heritability of the target trait. CONCLUSIONS: A commercial crossbred population is promising for the design of the reference population for genomic prediction, and selective genotyping of CB animals with extreme phenotypes has the potential for maximizing genetic improvement for CB performance in the pig industry.

Genetic dissection of 26 meat cut, meat quality and carcass traits in four pig populations.

BACKGROUND: Currently, meat cut traits are integrated in pig breeding objectives to gain extra profit. However, little is known about the heritability of meat cut proportions (MCP) and their correlations with other traits. The aims of this study were to assess the heritability and genetic correlation of MCP with carcass and meat quality traits using single nucleotide polymorphism chips and conduct a genome-wide association study (GWAS) to identify candidate genes for MCP. RESULTS: Seventeen MCP, 12 carcass, and seven meat quality traits were measured in 2012 pigs from four populations (Landrace; Yorkshire; Landrace and Yorkshire hybrid pigs; Duroc, and Landrace and Yorkshire hybrid pigs). Estimates of the heritability for MCP ranged from 0.10 to 0.55, with most estimates being moderate to high and highly consistent across populations. In the combined population, the heritability estimates for the proportions of scapula bone, loin, back fat, leg bones, and boneless picnic shoulder were 0.44 ± 0.04, 0.36 ± 0.04, 0.44 ± 0.04, 0.38 ± 0.04, and 0.39 ± 0.04, respectively. Proportion of middle cuts was genetically significantly positively correlated with intramuscular fat content and backfat depth. Proportion of ribs was genetically positively correlated with carcass oblique length and straight length (0.35 ± 0.08 to 0.45 ± 0.07) and negatively correlated with backfat depth (- 0.26 ± 0.10 to - 0.45 ± 0.10). However, weak or nonsignificant genetic correlations were observed between most MCP, indicating their independence. Twenty-eight quantitative trait loci (QTL) for MCP were detected by GWAS, and 24 new candidate genes related to MCP were identified, which are involved with growth, height, and skeletal development. Most importantly, we found that the development of the bones in different parts of the body may be regulated by different genes, among which HMGA1 may be the strongest candidate gene affecting forelimb bone development. Moreover, as previously shown, VRTN is a causal gene affecting vertebra number, and BMP2 may be the strongest candidate gene affecting hindlimb bone development. CONCLUSIONS: Our results indicate that breeding programs for MCP have the potential to enhance carcass composition by increasing the proportion of expensive cuts and decreasing the proportion of inexpensive cuts. Since MCP are post-slaughter traits, the QTL and candidate genes related to these traits can be used for marker-assisted and genomic selection.

Chromosomal genome and population genetic analyses to reveal genetic architecture, breeding history and genes related to cadmium accumulation in Lentinula edodes.

BACKGROUND: Lentinula edodes (Berk.) is the second most productive mushroom in the world. It contains compounds effective for antiviral, antitumor, antioxidant and immune regulation. Although genomes have previously been reported for this species, a high-quality chromosome-level reference for L. edodes is unavailable. This hinders detailed investigation of population genetics, breeding history of strains and genes related to environmental stress responses. RESULTS: A high-quality chromosome-level genome was constructed. We separated a monokaryon from protoplasts of the commercial L. edodes strain L808 and assembled the genome of L. edodes using PacBio long-read and Illumina short-read sequencing, along with the high-throughput chromatin conformation capture (Hi-C) technique. We assembled a 45.87 Mb genome, and 99% of the sequences were anchored onto 10 chromosomes. The contig and scaffold N50 length were 2.17 and 4.94 Mb, respectively. Over 96% of the complete Benchmarking Universal Single-Copy Orthologs (BUSCO) were identified, and 9853 protein-coding genes were predicted. We performed population genome resequencing using 34 wild strains and 65 commercial cultivars of L. edodes originating from China, Japan, the United States and Australia. Based on whole-genome variants, we showed substantial differences in the Chinese wild population, which divided into different branches according to the main areas of their geographical distribution. We also determined the breeding history of L. edodes at the molecular level, and demonstrated that the cultivated strains in China mainly originated from wild strains from China and Northeast Asia. Phenotypic analysis showed that 99 strains exhibited differences on the Cd accumulation. Three significant loci in the of L. edodes genome were identified using the genome-wide association study (GWAS) of Cd accumulation traits. Functional genes associated with Cd accumulation traits were related to DNA ligase and aminoacyl tRNA synthetase, indicating that DNA damage repair and in vivo protein translation may be responses to Cd stress. CONCLUSIONS: A high-quality chromosome-level genome and population genetic data of L. edodes provide genetic resources for functional genomic, evolutionary and artificial breeding studies for L. edodes.

A further look at quantitative trait loci for growth and fatness traits in a White Duroc × Erhualian F3 intercross population.

Genetic analysis of porcine growth and fatness traits is beneficial to the swine industry and provides a reference to understand human obesity. Here, we obtained 29 growth and fatness traits for 473 individuals from a White Duroc × Erhualian F3 intercross population. Basic statistical analyses showed that: (1) Positive correlations between different-stage body weights were detected, the shorter the time interval the stronger the correlation. (2) Strong correlations existed in the paired fatness traits. (3) With the growth of age, the correlation between fatness and body weight was increasing. All pigs were genotyped by Illumina 50 K SNP chips and their whole-genome genotypes were imputed referred to 109 re-sequencing data. We performed common and imputation-based GWASs for these traits. Two genome-wide significant loci on swine chromosome (SSC) 4 and 7 were repeatedly detected. The strongest association (P = 3.24 × 10-19) was detected at 31.96 Mb on SSC7 for leaf fat weight. On this locus, seven major haplotypes were identified, of which two were novel and had an increasing-fatness effect. In the imputation-based GWAS, three new loci were identified. Our findings provide further insights into and enhance our understanding of genetic mechanism of porcine growth and fat deposition.

Identification and validation of a regulatory mutation upstream of the BMP2 gene associated with carcass length in pigs.

BACKGROUND: Carcass length is very important for body size and meat production for swine, thus understanding the genetic mechanisms that underly this trait is of great significance in genetic improvement programs for pigs. Although many quantitative trait loci (QTL) have been detected in pigs, very few have been fine-mapped to the level of the causal mutations. The aim of this study was to identify potential causal single nucleotide polymorphisms (SNPs) for carcass length by integrating a genome-wide association study (GWAS) and functional assays. RESULTS: Here, we present a GWAS in a commercial Duroc × (Landrace × Yorkshire) (DLY) population that reveals a prominent association signal (P = 4.49E-07) on pig chromosome 17 for carcass length, which was further validated in two other DLY populations. Within the detected 1 Mb region, the BMP2 gene stood out as the most likely causal candidate because of its functions in bone growth and development. Whole-genome gene expression studies showed that the BMP2 gene was differentially expressed in the cartilage tissues of pigs with extreme carcass length. Then, we genotyped an additional 267 SNPs in 500 selected DLY pigs, followed by further whole-genome SNP imputation, combined with deep genome resequencing data on multiple pig breeds. Reassociation analyses using genotyped and imputed SNP data revealed that the rs320706814 SNP, located approximately 123 kb upstream of the BMP2 gene, was the strongest candidate causal mutation, with a large association with carcass length, with a ~ 4.2 cm difference in length across all three DLY populations (N = 1501; P = 3.66E-29). This SNP segregated in all parental lines of the DLY (Duroc, Large White and Landrace) and was also associated with a significant effect on body length in 299 pure Yorkshire pigs (P = 9.2E-4), which indicates that it has a major value for commercial breeding. Functional assays showed that this SNP is likely located within an enhancer and may affect the binding affinity of transcription factors, thereby regulating BMP2 gene expression. CONCLUSIONS: Taken together, these results suggest that the rs320706814 SNP on pig chromosome 17 is a putative causal mutation for carcass length in the widely used DLY pigs and has great value in breeding for body size in pigs.

A worldwide map of swine short tandem repeats and their associations with evolutionary and environmental adaptations.

BACKGROUND: Short tandem repeats (STRs) are genetic markers with a greater mutation rate than single nucleotide polymorphisms (SNPs) and are widely used in genetic studies and forensics. However, most studies in pigs have focused only on SNPs or on a limited number of STRs. RESULTS: This study screened 394 deep-sequenced genomes from 22 domesticated pig breeds/populations worldwide, wild boars from both Europe and Asia, and numerous outgroup Suidaes, and identified a set of 878,967 polymorphic STRs (pSTRs), which represents the largest repository of pSTRs in pigs to date. We found multiple lines of evidence that pSTRs in coding regions were affected by purifying selection. The enrichment of trinucleotide pSTRs in coding sequences (CDS), 5'UTR and H3K4me3 regions suggests that trinucleotide STRs serve as important components in the exons and promoters of the corresponding genes. We demonstrated that, compared to SNPs, pSTRs provide comparable or even greater accuracy in determining the breed identity of individuals. We identified pSTRs that showed significant population differentiation between domestic pigs and wild boars in Asia and Europe. We also observed that some pSTRs were significantly associated with environmental variables, such as average annual temperature or altitude of the originating sites of Chinese indigenous breeds, among which we identified loss-of-function and/or expanded STRs overlapping with genes such as AHR, LAS1L and PDK1. Finally, our results revealed that several pSTRs show stronger signals in domestic pig-wild boar differentiation or association with the analysed environmental variables than the flanking SNPs within a 100-kb window. CONCLUSIONS: This study provides a genome-wide high-density map of pSTRs in diverse pig populations based on genome sequencing data, enabling a more comprehensive characterization of their roles in evolutionary and environmental adaptation.

A whole genome sequence association study of muscle fiber traits in a White Duroc×Erhualian F2 resource population.

OBJECTIVE: Muscle fiber types, numbers and area are crucial aspects associated with meat production and quality. However, there are few studies of pig muscle fibre traits in terms of the detection power, false discovery rate and confidence interval precision of whole-genome quantitative trait loci (QTL). We had previously performed genome scanning for muscle fibre traits using 183 microsatellites and detected 8 significant QTLs in a White Duroc× Erhualian F2 population. The confidence intervals of these QTLs ranged between 11 and 127 centimorgan (cM), which contained hundreds of genes and hampered the identification of QTLs. A whole-genome sequence imputation of the population was used for fine mapping in this study. METHODS: A whole-genome sequences association study was performed in the F2 population. Genotyping was performed for 1,020 individuals (19 F0, 68 F1, and 933 F2). The whole-genome variants were imputed and 21,624,800 single nucleotide polymorphisms (SNPs) were identified and examined for associations to 11 longissimus dorsi muscle fiber traits. RESULTS: A total of 3,201 significant SNPs comprising 7 novel QTLs showing associations with the relative area of fiber type I (I_RA), the fiber number per square centimeter (FN) and the total fiber number (TFN). Moreover, one QTL on pig chromosome 14 was found to affect both FN and TFN. Furthermore, four plausible candidate genes associated with FN (kinase non-catalytic C-lobe domain containing [KNDC1]), TFN (KNDC1), and I_RA (solute carrier family 36 member 4, contactin associated protein like 5, and glutamate metabotropic receptor 8) were identified. CONCLUSION: An efficient and powerful imputation-based association approach was utilized to identify genes potentially associated with muscle fiber traits. These identified genes and SNPs could be explored to improve meat production and quality via marker-assisted selection in pigs.

Muscle glycogen level and occurrence of acid meat in commercial hybrid pigs are regulated by two low-frequency causal variants with large effects and multiple common variants with small effects.

BACKGROUND: Meat production from the commercial crossbred Duroc × (Landrace × Yorkshire) (DLY) pig is predominant in the pork industry, but its meat quality is often impaired by low ultimate pH (pHu). Muscle glycogen level at slaughter is closely associated with pHu and meat technological quality, but its genetic basis remains elusive. The aim of this study was to identify genes and/or causative mutations associated with muscle glycogen level and other meat quality traits by performing a genome-wide association study (GWAS) and additional analyses in a population of 610 DLY pigs. RESULTS: Our initial GWAS identified a genome-wide significant (P = 2.54e-11) quantitative trait locus (QTL) on SSC15 (SSC for Sus scrofa chromosome) for the level of residual glycogen and glucose (RG) in the longissimus muscle at 45 min post-mortem. Then, we demonstrated that a low-frequency (minor allele frequency = 0.014) R200Q missense mutation in the PRKAG3 (RN) gene caused this major QTL effect on RG. Moreover, we showed that the 200Q (RN-) allele was introgressed from the Hampshire breed into more than one of the parental breeds of the DLY pigs. After conditioning on R200Q, re-association analysis revealed three additional QTL for RG on SSC3 and 4, and on an unmapped scaffold (AEMK02000452.1). The SSC3 QTL was most likely caused by a splice mutation (g.8283C>A) in the PHKG1 gene that we had previously identified. Based on functional annotation, the genes TMCO1 on SSC4 and CKB on the scaffold represent promising candidate genes for the other two QTL. There were significant interaction effects of the GWAS tag SNPs at those two loci with PRKAG3 R200Q on RG. In addition, a number of common variants with potentially smaller effects on RG (P < 10-4) were uncovered by a second conditional GWAS after adjusting for the two causal SNPs, R200Q and g.8283C>A. CONCLUSIONS: We found that the RN- allele segregates in the parental lines of our DLY population and strongly influences its meat quality. Our findings also indicate that the genetic basis of RG in DLY can be mainly attributed to two major genes (PRKAG3 and PHKG1), along with many minor genes.

An association study using imputed whole-genome sequence data identifies novel significant loci for growth-related traits in a Duroc × Erhualian F2 population.

The average daily gain (ADG) and body weight (BW) are very important traits for breeding programs and for the meat production industry, which have attracted many researchers to delineate the genetic architecture behind these traits. In the present study, single- and multi-trait genome-wide association studies (GWAS) were performed between imputed whole-genome sequence data and the traits of the ADG and BW at different stages in a large-scale White Duroc × Erhualian F2 population. A bioinformatics annotation analysis was used to assist in the identification of candidate genes that are associated with these traits. Five and seven genome-wide significant quantitative trait loci (QTLs) were identified by single- and multi-trait GWAS, respectively. Furthermore, more than 40 genome-wide suggestive loci were detected. On the basis of the whole-genome sequence association study and the bioinformatics analysis, NDUFAF6, TNS1 and HMGA1 stood out as the strongest candidate genes. The presented single- and multi-trait GWAS analysis using imputed whole-genome sequence data identified several novel QTLs for pig growth-related traits. Integrating the GWAS with bioinformatics analysis can facilitate the more accurate identification of candidate genes. Higher imputation accuracy, time-saving algorithms, improved models and comprehensive databases will accelerate the identification of causal genes or mutations, which will contribute to genomic selection and pig breeding in the future.

Unravelling the genetic loci for growth and carcass traits in Chinese Bamaxiang pigs based on a 1.4 million SNP array.

Bamaxiang pig is from Guangxi province in China, characterized by its small body size and two-end black coat colour. It is an important indigenous breed for local pork market and excellent animal model for biomedical research. In this study, we performed genomewide association studies (GWAS) on 43 growth and carcass traits in 315 purebred Bamaxiang pigs based on a 1.4 million SNP array. We observed considerable phenotypic variability in the growth and carcass traits in the Bamaxiang pigs. The corresponding SNP based heritability varied greatly across the 43 traits and ranged from 9.0% to 88%. Through a conditional GWAS, we identified 53 significant associations for 35 traits at p value threshold of 10-6 . Among which, 26 associations on chromosome 3, 7, 14 and X passed a genomewide significance threshold of 5 × 10-8 . The most remarkable loci were at around 30.6 Mb on chromosome 7, which had growth stage-dependent effects on body lengths and cannon circumferences and showed large effects on multiple carcass traits. We discussed HMGA1 NUDT3, EIF2AK1, TMEM132C and AFF2 that near the lead SNP of significant loci as plausible candidate genes for corresponding traits. We also showed that including phenotypic covariate in GWAS can help to reveal additional significant loci for the target traits. The results provide insight into the genetic architecture of growth and carcass traits in Bamaxiang pigs.

Comparative analysis of the GBLUP, emBayesB, and GWAS algorithms to predict genetic values in large yellow croaker (Larimichthys crocea).

BACKGROUND: The advances of sequencing technology accelerate the development of theory of molecular quantitative genetics such as QTL mapping, genome-wide association study and genomic selection. This paper was designed to study genomic selection in large yellow croaker breeding. The aims of this study were: (i) to estimate heritability values of traits in large yellow croaker; (ii) to assess feasibility of genomic selection in the traits of growth rate and meat quality; (iii) to compare predictive accuracies affected by different algorithms and training sizes, and to find what training sizes could reach ideal accuracies; (iv) to compare results of GWAS with genomic prediction, and to assess feasibility of pre-selection of significant SNPs in genomic selection. 500 individuals were tested in the trait of body weight and body length, while 176 were tested in the percentage of n-3 highly unsaturated fatty acids (n-3HUFA) in muscle. GBLUP and emBayesB were used to perform genomic prediction. RESULTS: Genotyping-By-Sequencing method was used to construct the libraries for the NGS sequencing and find ~30,000 SNPs. Heritability estimates were 0.604, 0.586 and 0.438 for trait of body weight, body length and n-3HUFA, respectively. The predictive abilities estimated by GBLUP showed higher than that by emBayesB in traits of body weight and body length. However, the result was just the opposite in n-3HUFA. According to fit the curve of predictive accuracy, we estimated that at least 1000 individuals in training set could reach an accuracy of 0.8 in body weight and body length. GBLUP, emBayesB and GWAS could not always find significant SNPs associated with phenotypes consistently. Significant SNPs were selected by emBayesB could obtain the largest proportions to explain total additive genetic variances. CONCLUSIONS: This research showed that genomic selection was feasible in large yellow croaker breeding. We suggest doing a test before deciding to use which algorithm in specific trait in genomic prediction. We estimated required training sizes to reach ideal predictive accuracies and assessed feasibility of pre-selection of SNPs successfully. Because of high mortality rate of fish and high cost in genomic sequencing, genomic selection may be more suitable for applying on some traits which cannot be measured on candidates directly.

Genetic architecture of fatty acid composition in the longissimus dorsi muscle revealed by genome-wide association studies on diverse pig populations.

BACKGROUND: Fatty acid composition in muscle is an important factor that affects the nutritive value and taste of pork. To investigate the genetic architecture of fatty acid composition of pork, we measured fatty acid contents in longissimus dorsi muscle of 1244 pigs from three divergent populations and conducted genome-wide association studies (GWAS) for fatty acid contents. RESULTS: We detected 26 genome-wide significant quantitative trait loci (QTL) on eight chromosomes (SSC for Sus scrofa) for eight fatty acids. These loci not only replicated previously reported QTL for C18:0 on SSC14 and C20:0 on SSC16, but also included several novel QTL such as those for C20:1 on SSC7, C14:0 on SSC9, and C14:0, C16:0 and C16:1 on SSC12. Furthermore, we performed a meta-analysis of GWAS on five populations, including the three populations that were investigated in this study and two additional populations that we had previously examined. This enhanced the strength of the associations detected between fatty acid composition and several marker loci, especially for those for C18:0 on SSC14 and C20:0 on SSC16. The genes ELOVL5, ELOVL6, ELOVL7, FASN, SCD and THRSP, which have functions that are directly relevant to fatty acid metabolism, are proximal to the top associated markers at six significant QTL. CONCLUSIONS: The findings improve our understanding of the genetic architecture of fatty acid composition in pork and contribute to further fine-map and characterize genes that influence fatty acid composition.

A genome-wide association study of copy number variations with umbilical hernia in swine.

Umbilical hernia (UH) is one of the most common congenital defects in pigs, leading to considerable economic loss and serious animal welfare problems. To test whether copy number variations (CNVs) contribute to pig UH, we performed a case-control genome-wide CNV association study on 905 pigs from the Duroc, Landrace and Yorkshire breeds using the Porcine SNP60 BeadChip and penncnv algorithm. We first constructed a genomic map comprising 6193 CNVs that pertain to 737 CNV regions. Then, we identified eight CNVs significantly associated with the risk for UH in the three pig breeds. Six of seven significantly associated CNVs were validated using quantitative real-time PCR. Notably, a rare CNV (CNV14:13030843-13059455) encompassing the NUGGC gene was strongly associated with UH (permutation-corrected P = 0.0015) in Duroc pigs. This CNV occurred exclusively in seven Duroc UH-affected individuals. SNPs surrounding the CNV did not show association signals, indicating that rare CNVs may play an important role in complex pig diseases such as UH. The NUGGC gene has been implicated in human omphalocele and inguinal hernia. Our finding supports that CNVs, including the NUGGC CNV, contribute to the pathogenesis of pig UH.

Production of Transgenic Pigs with an Introduced Missense Mutation of the Bone Morphogenetic Protein Receptor Type IB Gene Related to Prolificacy.

In the last few decades, transgenic animal technology has witnessed an increasingly wide application in animal breeding. Reproductive traits are economically important to the pig industry. It has been shown that the bone morphogenetic protein receptor type IB (BMPR1B) A746G polymorphism is responsible for the fertility in sheep. However, this causal mutation exits exclusively in sheep and goat. In this study, we attempted to create transgenic pigs by introducing this mutation with the aim to improve reproductive traits in pigs. We successfully constructed a vector containing porcine BMPR1B coding sequence (CDS) with the mutant G allele of A746G mutation. In total, we obtained 24 cloned male piglets using handmade cloning (HMC) technique, and 12 individuals survived till maturation. A set of polymerase chain reactions indicated that 11 of 12 matured boars were transgene-positive individuals, and that the transgenic vector was most likely disrupted during cloning. Of 11 positive pigs, one (No. 11) lost a part of the terminator region but had the intact promoter and the CDS regions. cDNA sequencing showed that the introduced allele (746G) was expressed in multiple tissues of transgene-positive offspring of No.11. Western blot analysis revealed that BMPR1B protein expression in multiple tissues of transgene-positive F1 piglets was 0.5 to 2-fold higher than that in the transgene-negative siblings. The No. 11 boar showed normal litter size performance as normal pigs from the same breed. Transgene-positive F1 boars produced by No. 11 had higher semen volume, sperm concentration and total sperm per ejaculate than the negative siblings, although the differences did not reached statistical significance. Transgene-positive F1 sows had similar litter size performance to the negative siblings, and more data are needed to adequately assess the litter size performance. In conclusion, we obtained 24 cloned transgenic pigs with the modified porcine BMPR1B CDS using HMC. cDNA sequencing and western blot indicated that the exogenous BMPR1B CDS was successfully expressed in host pigs. The transgenic pigs showed normal litter size performance. However, no significant differences in litter size were found between transgene-positive and negative sows. Our study provides new insight into producing cloned transgenic livestock related to reproductive traits.

Rapid construction of genome map for large yellow croaker (Larimichthys crocea) by the whole-genome mapping in BioNano Genomics Irys system.

BACKGROUND: Large yellow croaker (Larimichthys crocea) is an important commercial fish in China and East-Asia. The annual product of the species from the aqua-farming industry is about 90 thousand tons. In spite of its economic importance, genetic studies of economic traits and genomic selections of the species are hindered by the lack of genomic resources. Specifically, a whole-genome physical map of large yellow croaker is still missing. The traditional BAC-based fingerprint method is extremely time- and labour-consuming. Here we report the first genome map construction using the high-throughput whole-genome mapping technique by nanochannel arrays in BioNano Genomics Irys system. RESULTS: For an optimal marker density of ~10 per 100 kb, the nicking endonuclease Nt.BspQ1 was chosen for the genome map generation. 645,305 DNA molecules with a total length of ~112 Gb were labelled and detected, covering more than 160X of the large yellow croaker genome. Employing IrysView package and signature patterns in raw DNA molecules, a whole-genome map of large yellow croaker was assembled into 686 maps with a total length of 727 Mb, which was consistent with the estimated genome size. The N50 length of the whole-genome map, including 126 maps, was up to 1.7 Mb. The excellent hybrid alignment with large yellow croaker draft genome validated the consensus genome map assembly and highlighted a promising application of whole-genome mapping on draft genome sequence super-scaffolding. The genome map data of large yellow croaker are accessible on lycgenomics.jmu.edu.cn/pm. CONCLUSION: Using the state-of-the-art whole-genome mapping technique in Irys system, the first whole-genome map for large yellow croaker has been constructed and thus highly facilitates the ongoing genomic and evolutionary studies for the species. To our knowledge, this is the first public report on genome map construction by the whole-genome mapping for aquatic-organisms. Our study demonstrates a promising application of the whole-genome mapping on genome maps construction for other non-model organisms in a fast and reliable manner.

Genome-wide association analysis reveals genetic loci and candidate genes for meat quality traits in Chinese Laiwu pigs.

Meat quality traits have economically significant impacts on the pig industry, and can be improved using molecular approaches in pig breeding. Since 1994 when the first genome-wide scan for quantitative trait loci (QTLs) in pig was reported, over the past two decades, numerous QTLs have been identified for meat quality traits by family based linkage analyses. However, little is known about the genetic variants for meat quality traits in Chinese purebred or outbred populations. To unveil it, we performed a genome-wide association study for 10 meat quality traits in Chinese purebred Laiwu pigs. In total, 75 significant SNPs (P < 1.01 × 10(-6)) and 33 suggestive SNPs (P < 2.03 × 10(-5)) were identified. On SSC12, a region between 56.22 and 61.49 Mb harbored a cluster of SNPs that were associated with meat color parameters (L*, lightness; a*, redness; b*, yellowness) and moisture content of longissimus muscle (LM) and semimembranosus muscle at the genome-wide significance level. A region on SSC4 also has pleiotropic effects on moisture content and drip loss of LM. In addition, this study revealed at least five novel QTLs and several candidate genes including 4-linked MYH genes (MYH1, MYH2, MYH3, and MYH13), MAL2, LPAR1, and PRKAG3 at four significant loci. Except for the SSC12 QTL, other QTLs are likely tissue-specific. These results provide new insights into the genetic basis of meat quality traits in Chinese Laiwu pigs and some significant SNPs reported here could be incorporated into the selection programs involving this breed.

Molecular cloning and functional characterization of a QM protein in large yellow croaker (Larimichthys crocea).

Since it was proposed to be a tumor suppressor in 1991, QM protein has attracted intensive and wide attention in plants, animals and fungi research fields. Up to date, however, the function of QM protein in fish immunity remains unknown. In this investigation, a QM gene (named as LycQM gene) was cloned from large yellow croaker (Larimichthys crocea), and LycQM protein was expressed in Escherichia coli and purified. The LycQM gene was ubiquitously transcribed in multi-tissues, including spleen, muscle, heart, liver, intestine, blood and head kidney. By quantitative real-time RT-PCR analysis, we found the highest and the lowest expression level of LycQM gene in head kidney and in heart, respectively. Time course analysis showed that LycQM expression was obviously up-regulated in blood and head kidney after immunization with polyinosinic polycytidynic acid (poly I:C), formalin-inactive Gram-negative bacterium Vibrio parahaemolyticus and bacterial lipopolysaccharides (LPS). Moreover, as demonstrated by RNAi assays, LycQM protein could regulate the activity of phenoloxidase, a key enzyme in the proPO activation system of immunity. These results suggested that LycQM protein might play an important role in the immune response against microorganisms in large yellow croaker.