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.

Investigation on muscle fiber types and meat quality and estimation of their heritability and correlation coefficients with each other in four pig populations.

The aim of this study was to investigate the muscle fiber types and meat quality in four populations and estimate the heritability and correlation coefficients of those traits in Shanxia long black pig (SX). In this study, a total of 318 pigs were recorded for 16 traits of the muscle fiber types and meat quality in four populations, including 256 individuals from the new breed SX. The population had a significant effect on all recorded traits, and the meat quality of the Lulai black pig was better than the remaining populations. The heritability (h2 ) of meat quality traits was from 0.06 (pH at 24 h) to 0.47 (shearing force), and the muscle fiber types belonged to the traits with low to medium heritability. The density of total fiber had the highest h2 (0.40), while the percentage of type IIA had the lowest h2 (0.04). Most traits are phenotypically correlated with each other, but only a small proportion of traits are genetically correlated with each other. None fiber type genetically correlated with meat quality significantly, because the genetic correlation coefficients had large standard errors. These results provided some insights into genetic improvements for the meat quality in pig breeds and also indicated that the parameters of muscle fiber characteristics can explain parts of the variation in meat quality.

Extraction, structural characterization and immunoactivity of glucomannan type polysaccahrides from Lilium brownii var. viridulum Baker.

Homogeneous polysaccharide (LBP) was extracted and purified from the bulblets of Lilium brownii var. viridulum Baker with a molecular weight of 312 kDa. The monosaccharides are composed of mannose and glucose, and the corresponding molar ratios are 0.582 and 0.418, respectively. FT-IR, LC-MS, NMR, GC-MS and HPAEC were used to analyze the functional groups, glycosidic linkages and chemical structure of LBP, which was a 1-4-linked glucomannan and contained a dodecasaccharide repeating units of →4)-ß-D-Manp-(1 â†’ 4)-ß-D-Manp-(1 â†’ 4)-ß-D-Manp-(1 â†’ 4)-ß-D-Glcp-(1 â†’ 4)-ß-D-Manp-(1 â†’ 4)-ß-D-Manp-(1 â†’ 4)-ß-D-Glcp-(1 â†’ 4)-α-D-Glcp-(1 â†’ 4)-ß-D-Glcp-(1 â†’ 4)-ß-D-Glcp-(1 â†’ 4)-ß-D-Manp-(1 â†’ 4)-ß-D-Manp-(1 â†’ . In vitro experimental results showed that LBP had noble biocompatibility, and a low dose of 5 µg/mL LBP significantly up-regulated the mRNA expression of TNF-α, iNOS, IL-6, IL-1ß and Toll-like receptors family (TLRs) in RAW 264.7 cells. In conclusion, LBP played an important role in immunomodulation, and further studies on the specific immunomodulatory mechanisms of LBP on RAW 264.7 cells are still needed.

Association between MC1R gene and coat color segregation in Shanxia long black pig and Lulai black pig.

BACKGROUND: Coat color, as a distinct phenotypic characteristic of pigs, is often subject to preference and selection, such as in the breeding process of new breed. Shanxia long black pig was derived from an intercross between Berkshire boars and Licha black pig sows, and it was bred as a paternal strain with high-quality meat and black coat color. Although the coat color was black in the F1 generation of the intercross, it segregated in the subsequent generations. This study aims to decode the genetic basis of coat color segregation and develop a method to distinct black pigs from the spotted in Shanxia long black pig. RESULTS: Only a QTL was mapped at the proximal end of chromosome 6, and MC1R gene was picked out as functional candidate gene. A total of 11 polymorphic loci were identified in MC1R gene, and only the c.67_68insCC variant was co-segregating with coat color. This locus isn't recognized by any restriction endonuclease, so it can't be genotyped by PCR-RFLP. The c.370G > A polymorphic locus was also significantly associated with coat color, and has been in tightly linkage disequilibrium with the c.67_68insCC. Furthermore, it is recognized by BspHI. Therefore, a PCR-RFLP method was set up to genotype this locus. Besides the 175 sequenced individuals, another more 1,391 pigs were genotyped with PCR-RFLP, and all of pigs with GG (one band) were black. CONCLUSION: MC1R gene (c.67_68insCC) is the causative gene (mutation) for the coat color segregation, and the PCR-RFLP of c.370G > A could be used in the breeding program of Shanxia long black pig.

Genetic dissection and genomic prediction for pork cuts and carcass morphology traits in pig.

BACKGROUND: As pre-cut and pre-packaged chilled meat becomes increasingly popular, integrating the carcass-cutting process into the pig industry chain has become a trend. Identifying quantitative trait loci (QTLs) of pork cuts would facilitate the selection of pigs with a higher overall value. However, previous studies solely focused on evaluating the phenotypic and genetic parameters of pork cuts, neglecting the investigation of QTLs influencing these traits. This study involved 17 pork cuts and 12 morphology traits from 2,012 pigs across four populations genotyped using CC1 PorcineSNP50 BeadChips. Our aim was to identify QTLs and evaluate the accuracy of genomic estimated breed values (GEBVs) for pork cuts. RESULTS: We identified 14 QTLs and 112 QTLs for 17 pork cuts by GWAS using haplotype and imputation genotypes, respectively. Specifically, we found that HMGA1, VRTN and BMP2 were associated with body length and weight. Subsequent analysis revealed that HMGA1 primarily affects the size of fore leg bones, VRTN primarily affects the number of vertebrates, and BMP2 primarily affects the length of vertebrae and the size of hind leg bones. The prediction accuracy was defined as the correlation between the adjusted phenotype and GEBVs in the validation population, divided by the square root of the trait's heritability. The prediction accuracy of GEBVs for pork cuts varied from 0.342 to 0.693. Notably, ribs, boneless picnic shoulder, tenderloin, hind leg bones, and scapula bones exhibited prediction accuracies exceeding 0.600. Employing better models, increasing marker density through genotype imputation, and pre-selecting markers significantly improved the prediction accuracy of GEBVs. CONCLUSIONS: We performed the first study to dissect the genetic mechanism of pork cuts and identified a large number of significant QTLs and potential candidate genes. These findings carry significant implications for the breeding of pork cuts through marker-assisted and genomic selection. Additionally, we have constructed the first reference populations for genomic selection of pork cuts in pigs.

Dominance is common in mammals and is associated with trans-acting gene expression and alternative splicing.

BACKGROUND: Dominance and other non-additive genetic effects arise from the interaction between alleles, and historically these phenomena play a major role in quantitative genetics. However, most genome-wide association studies (GWAS) assume alleles act additively. RESULTS: We systematically investigate both dominance-here representing any non-additive within-locus interaction-and additivity across 574 physiological and gene expression traits in three mammalian stocks: F2 intercross pigs, rat heterogeneous stock, and mice heterogeneous stock. Dominance accounts for about one quarter of heritable variance across all physiological traits in all species. Hematological and immunological traits exhibit the highest dominance variance, possibly reflecting balancing selection in response to pathogens. Although most quantitative trait loci (QTLs) are detectable as additive QTLs, we identify 154, 64, and 62 novel dominance QTLs in pigs, rats, and mice respectively that are undetectable as additive QTLs. Similarly, even though most cis-acting expression QTLs are additive, gene expression exhibits a large fraction of dominance variance, and trans-acting eQTLs are enriched for dominance. Genes causal for dominance physiological QTLs are less likely to be physically linked to their QTLs but instead act via trans-acting dominance eQTLs. In addition, thousands of eQTLs are associated with alternatively spliced isoforms with complex additive and dominant architectures in heterogeneous stock rats, suggesting a possible mechanism for dominance. CONCLUSIONS: Although heritability is predominantly additive, many mammalian genetic effects are dominant and likely arise through distinct mechanisms. It is therefore advantageous to consider both additive and dominance effects in GWAS to improve power and uncover causality.

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.

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.

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.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity.

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by FST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

Inulin-type fructans obtained from Atractylodis Macrocephalae by water/alkali extraction and immunoregulatory evaluation.

Two homogenous polysaccharides extracted from Atractylodes macrocephala Koidz. were investigated by water extraction (AMP-FW) and alkali solution extraction (AMP-FA) after purification by anion exchange column and size exclusion chromatography. The molecular weight of AMP-FW and AMP-FA were 2874 Da and 3438 Da, respectively, estimated by high performance gel permeation chromatography (HPGPC). The monosaccharide compositions of AMP-FW and AMP-FA were glucose and fructose at a molar ratio of 0.11:0.89 determined by high performance anion exchange chromatography (HPAEC). The functional groups, glycosidic linkages and the chemical structure were characterized by FT-IR, GC-MS and NMR, which comprehensively indicated a similar inulin-type fructan structure of the two polysaccharides from A. macrocephala. However, the scanning electron microscopy (SEM) results showed different microstructures that irregular lamellar shape for the AMP-FW and spheroid shape for the AMP-FA. The further studies on immunomodulation showed that AMP-FW at 50 µg/mL could significantly (P < 0.05) stimulate RAW 264.7 cells by enhancing the mRNA expression of TNF-α and IL-1ß, which had a relative high immunomodulatory potential when compared to AMP-FA. Their activation on different toll-like receptors (TLR) also indicated their different roles in the immunoregulation. Overall, these findings reported here will serve as the basis for further structure-activity relationship studies.

Using phenotypic and genotypic big data to investigate the effect of muscle fiber characteristics on meat quality and eating quality traits in pigs.

Determining genetic correlations (GCs) between phenotypes that can be replicated across breeds or generations is important for animal breeding. A comprehensive and objective evaluation of this is dependent on enough variations in the studied phenotypes. To this end, we constructed a worldwide distributed eight-breeds crossbreed mosaic pig population and estimated the genetic and phenotypic correlations of muscle fiber characteristics (MFCs) with meat quality and eating quality traits using F6 and F7 generations (∼590 samples/generation) of this population. The GCs of the density of type IIA fibers and type IIB fibers with the loin-eye area, a*, color score, firmness score, and those of the proportions of the two fiber types with pH24h and b* were moderate to high (|rg| ≥ 0.3) in both populations. We also obtained moderate to high GCs of mean fiber density with five sensory quality traits. Our results provide an important reference for improving meat quality through the genetic regulation of MFCs.

Quantification and visualization of meat quality traits in pork using hyperspectral imaging.

Accurate and rapid determination of meat quality traits plays key roles in food industry and pig breeding. Currently, most of the spectroscopic instruments developed for meat quality determination can only obtain the spectral average value of the sample, so it is difficult to evaluate the spatial variation of meat quality traits. In this study, we evaluated the predictive potential of 14 meat quality traits based on large-scale VIS/NIR hyperspectral images collected by SpecimIQ. When predictions were based solely on hyperspectral data, the prediction accuracy (R2cv) for the majority of meat qualities ranged from 0.60 to 0.70. After adding texture information, the prediction accuracy of all traits is improved by different magnitudes (R2cv increases from 1.5% to 16.4%). Finally, the best model was utilized to visualize the spatial distribution of Fat (%) and Moisture (%) to assess their homogeneity. These results suggest that hyperspectral imaging has great potential for predicting and visualizing various meat qualities, as well as industrial applications for automated measurements.

Massive genome investigations reveal insights of prevalent introgression for environmental adaptation and triterpene biosynthesis in Ganoderma.

Genome introgression is one of the driving forces that can increase species and genetic diversity and facilitate the adaptive evolution of organisms and biodiversity conservation. However, the genomic introgression and its contribution to biodiversity of macrofungi are still unclear. The genus Ganoderma is a typical macrofungal group that plays crucial roles in forest ecosystem as saprophytic organisms and plant pathogens, and is also involved in human health as medicinal mushrooms. Most public Ganoderma genomes are fragmented, and reference genomes and whole-genome information of diverse germplasm resources for many Ganoderma species are lacking, thus hindering functional and evolutionary genomic investigations among Ganoderma species. In this study, we provide high-quality genomes of 10 Ganoderma species and whole-genome variants data of 224 individuals from various ecoregions, enabling us to infer the phylogeny of Ganoderma species and their historical population dynamics. Based on whole-genome variants, widespread and genome-wide introgression among Ganoderma species is revealed. Genes with significant introgression signals were related to stress response, digestive absorption, and secondary metabolite synthesis, factors that may contribute to environmental adaptation and important biocomponent metabolism. CYP512U6, an essential functional gene in the CYP450 family related to Ganoderma triterpene synthesis, was detected with significant introgression and selection signals combined with Ganoderma metabolomic analysis, indicating that both ancient gene exchange and recent domestication have contributed to the categories and content of secondary metabolites of Ganoderma. The reference genomes, whole-genome variants, and metabolite profiles could serve as abundant and valuable genetic resources for evolution, ecology, and conservation investigations of Ganoderma species and other macrofungi.

Whole-genome sequence-based association analyses on an eight-breed crossed heterogeneous stock of pigs reveal the genetic basis of skeletal muscle fiber characteristics.

Skeletal muscle fiber characteristics (MFCs) have been extensively studied due to their importance to human health and athletic ability, as well as to the quantity and quality of livestock meat production. Hence, we performed a genome-wide association study (GWAS) on nine muscle fiber traits by using whole genome sequence data in an eight-breed crossed heterogeneous stock pig population. This GWAS revealed 67 quantitative trait loci (QTLs) for these traits. The most significant GWAS signal was detected in the region of Sus scrofa chromosome 12 (SSC12) containing the MYH gene family. Notably, we identified a significant SNP rs322008693 (P = 7.52E-09) as the most likely causal mutation for the total number of muscle fibers (TNMF) QTL on SSC1. The results of EMSA and luciferase assays indicated that the rs322008693 SNP resided in a functional element. These findings provide valuable molecular markers for pig meat production selection as well as for deciphering the genetic mechanisms of the muscle fiber physiology.

Extensive evaluation of prediction performance for 15 pork quality traits using large scale VIS/NIRS data.

Visible and near-infrared spectroscopy (VIS/NIRS) has been extensively used in the livestock and food industries to quantify meat quality. Here, we collected VIS/NIRS data of 1206 pigs longissimus muscle, measured the corresponding 15 meat quality traits, and used seven models to predict these meat quality traits. The prediction performances of 7 models varied among predicted traits, with the Rcv2 of most traits above 0.9 in the best model. We have also established a new method, spectral-wide association analysis (SWAS), to select the feature wavelengths of measured traits. Results showed that the prediction performance is proportionate to the number of identified significant association wavelengths. We used the selected wavelengths to perform prediction again, and the prediction accuracy was similar to results with full wavelength using the best model, indicating effectiveness of feature wavelengths selection methods.

Pan-Genomes Provide Insights into the Genetic Basis of Auricularia heimuer Domestication.

In order to reveal the genetic variation signals of Auricularia heimuer that have occurred during their domestication and to find potential functional gene families, we constructed a monokaryotic pan-genome of A. heimuer representing four cultivated strains and four wild strains. The pan-genome contained 14,089 gene families, of which 67.56% were core gene families and 31.88% were dispensable gene families. We screened substrate utilization-related genes such as the chitinase gene ahchi1 of the glycoside hydrolase (GH) 18 family and a carbohydrate-binding module (CBM)-related gene from the dispensable families of cultivated populations. The genomic difference in the ahchi1 gene between the wild and cultivated genomes was caused by a 33 kb presence/absence variation (PAV). The detection rate of the ahchi1 gene was 93.75% in the cultivated population, significantly higher than that in the wild population (17.39%), indicating that it has been selected in cultivated strains. Principal component analysis (PCA) of the polymorphic markers in fragments near the ahchi1 gene was enriched in cultivated strains, and this was caused by multiple independent instances of artificial selection. We revealed for the first time the genetic basis of the ahchi1 gene in domestication, thereby providing a foundation for elucidating the potential function of the ahchi1 gene in the breeding of A. heimuer.

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.

Mapping short tandem repeats for liver gene expression traits helps prioritize potential causal variants for complex traits in pigs.

BACKGROUND: Short tandem repeats (STRs) were recently found to have significant impacts on gene expression and diseases in humans, but their roles on gene expression and complex traits in pigs remain unexplored. This study investigates the effects of STRs on gene expression in liver tissues based on the whole-genome sequences and RNA-Seq data of a discovery cohort of 260 F6 individuals and a validation population of 296 F7 individuals from a heterogeneous population generated from crosses among eight pig breeds. RESULTS: We identified 5203 and 5868 significantly expression STRs (eSTRs, FDR < 1%) in the F6 and F7 populations, respectively, most of which could be reciprocally validated (π1 = 0.92). The eSTRs explained 27.5% of the cis-heritability of gene expression traits on average. We further identified 235 and 298 fine-mapped STRs through the Bayesian fine-mapping approach in the F6 and F7 pigs, respectively, which were significantly enriched in intron, ATAC peak, compartment A and H3K4me3 regions. We identified 20 fine-mapped STRs located in 100 kb windows upstream and downstream of published complex trait-associated SNPs, which colocalized with epigenetic markers such as H3K27ac and ATAC peaks. These included eSTR of the CLPB, PGLS, PSMD6 and DHDH genes, which are linked with genome-wide association study (GWAS) SNPs for blood-related traits, leg conformation, growth-related traits, and meat quality traits, respectively. CONCLUSIONS: This study provides insights into the effects of STRs on gene expression traits. The identified eSTRs are valuable resources for prioritizing causal STRs for complex traits in pigs.