Integrated GWAS and transcriptome analysis reveals key genes associated with muscle fibre and fat traits in Gushi chicken.

1. In the following experiment meat quality traits of a Gushi-Anka F2 resource population were measured, and their heritability estimated. Intramuscular fat (IMF) had medium heritability (0.35) but leg muscle fibre density (LMD), leg muscle fibre diameter (LMF), breast muscle fibre density (BMD), fresh fat content (FFA), and absolute dry fat content (AFC) had low heritability (0-0.2). The IMF presented the most important genetic additive effect among the poultry meat quality-related traits studied.2. The phenotypic data of meat quality traits in the Gushi-Anka F2 resource population were combined with genotyping by sequencing (GBS) data to obtain genotype data. Six meat quality traits in 734 birds were analysed by GWAS. Based on these variants, 83 significant (-log10(p) > 4.42) single nucleotide polymorphisms and four quantitative trait loci (QTL) regions corresponding to 175 genes were identified. Further linkage disequilibrium (LD) analysis was conducted on chromosome 13 (Chr13) and chromosome 27 (Chr27) QTL regions.3. Based on the transcriptome data and GWAS results, 12 shared genes - ITGB3, DNAJC27, ETV4, C7orf50, FKBP1B, G3BP1, IGF2BP1, KCNH6, LOC416263, SCARA5, SMIM5 and TBL1XR1 were identified as candidate genes influencing muscle fibre and fat traits.

Deciphering the genetic basis of behavioral traits in dogs: Observed-trait GWAS and latent-trait GWAS analysis reveal key genes and variants.

Dogs exhibit remarkable phenotypic diversity, particularly in behavioral traits, making them an excellent model for studying the genetic basis of complex behaviors. Behavioral traits such as aggression and fear are highly heritable among different dog breeds, but their genetic basis is largely unknown. We used the genome-wide association study (GWAS) to identify candidate genes associated with nine behavioral traits including; stranger-directed aggression (SDA), owner-directed aggression (ODA), dog-directed aggression (DDA), stranger-directed fear (SDF), nonsocial fear (NF), dog-directed fear (DDF), touch sensitivity (TS), separation-related behavior (SRB) and attachment attention-seeking (AAS). The observed behavioral traits were collected from 38,714 to 40,460 individuals across 108 modern dog breeds. We performed a GWAS based on a latent trait extracted using the confirmatory factor analysis (CFA) method with nine observable behavioral traits and compared the results with those from the GWAS of the observed traits. Using both observed-trait and latent-trait GWAS, we identified 41 significant SNPs that were common between both GWAS methods, of which 26 were pleiotropic, as well as 10 SNPs unique to the latent-trait GWAS, and 5 SNPs unique to the observed-trait GWAS discovered. These SNPs were associated with 21 genes in latent-trait GWAS and 22 genes in the observed-trait GWAS, with 19 genes shared by both. According to previous studies, some of the genes from this study have been reported to be related to behavioral and neurological functions in dogs. In the human population, these identified genes play a role in either the formation of the nervous system or are linked to various mental health conditions. Taken together, our findings suggest that latent-trait GWAS for behavioral traits in dogs identifies significant latent genes that are neurologically prioritized.

Comparative transcriptome analysis and meta-QTLs mapping reveal the regulatory mechanism of cold tolerance in rice at the budding stage.

Rice (Oryza sativa L.) is one of the most extensively farmed food crops, but its development and productivity are significantly impacted by cold stress during the budding period. In this study, transcriptome sequencing was conducted on two types of rice: the cold-sensitive indica rice A117 and the substantially cold-tolerant japonica rice B106 under control and cold treatments. Differentially expressed genes between the two materials under cold conditions were analyzed using GO and KEGG enrichment analyses. The results revealed that processes such as the TCA cycle, glycolysis/glycogenesis, oxidative phosphorylation, and glutathione metabolism contribute to B106's cold tolerance. Additionally, an enrichment analysis of cold-induced genes in each material and shared genes identified significant enrichment in pathways such as glutathione metabolism, phenylpropanoid biosynthesis, and photosynthesis-antenna proteins. Initial cold tolerance QTLs at the rice bud stage were collected from published literature, and meta-QTL mapping identified 9 MQTLs. Gene expression profiling led to the identification of 75 potential DEGs within the 9 MQTLs region, from which four candidate genes (Os02g0194100, Os03g0802500, Os05g0129000, and Os07g0462000) were selected using qRT-PCR and gene annotation. These findings provide genetic resources for further research on the molecular mechanisms underlying rice's response to cold stress during the bud stage.

Genetic Foundations of Nellore Traits: A Gene Prioritization and Functional Analyses of Genome-Wide Association Study Results.

The main goal of this study was to pinpoint functional candidate genes associated with multiple economically important traits in Nellore cattle. After quality control, 1830 genomic regions sourced from 52 scientific peer-reviewed publications were used in this study. From these, a total of 8569 positional candidate genes were annotated for reproduction, 11,195 for carcass, 5239 for growth, and 3483 for morphological traits, and used in an over-representation analysis. The significant genes (adjusted p-values < 0.05) identified in the over-representation analysis underwent prioritization analyses, and enrichment analysis of the prioritized over-represented candidate genes was performed. The prioritized candidate genes were GFRA4, RFWD3, SERTAD2, KIZ, REM2, and ANKRD34B for reproduction; RFWD3, TMEM120A, MIEF2, FOXRED2, DUSP29, CARHSP1, OBI1, JOSD1, NOP58, and LOXL1-AS1 for the carcass; ANKRD34B and JOSD1 for growth traits; and no genes were prioritized for morphological traits. The functional analysis pinpointed the following genes: KIZ (plays a crucial role in spindle organization, which is essential in forming a robust mitotic centrosome), DUSP29 (involved in muscle cell differentiation), and JOSD1 (involved in protein deubiquitination, thereby improving growth). The enrichment of the functional candidate genes identified in this study highlights that these genes play an important role in the expression of reproduction, carcass, and growth traits in Nellore cattle.

Quantitative Trait Loci Mapping of Heading Date in Wheat under Phosphorus Stress Conditions.

Wheat (Triticum aestivum L.) is a crucial cereal crop, contributing around 20% of global caloric intake. However, challenges such as diminishing arable land, water shortages, and climate change threaten wheat production, making yield enhancement crucial for global food security. The heading date (HD) is a critical factor influencing wheat's growth cycle, harvest timing, climate adaptability, and yield. Understanding the genetic determinants of HD is essential for developing high-yield and stable wheat varieties. This study used a doubled haploid (DH) population from a cross between Jinmai 47 and Jinmai 84. QTL analysis of HD was performed under three phosphorus (P) treatments (low, medium, and normal) across six environments, using Wheat15K high-density SNP technology. The study identified 39 QTLs for HD, distributed across ten chromosomes, accounting for 2.39% to 29.52% of the phenotypic variance. Notably, five stable and major QTLs (Qhd.saw-3A.7, Qhd.saw-3A.8, Qhd.saw-3A.9, Qhd.saw-4A.4, and Qhd.saw-4D.3) were consistently detected across varying P conditions. The additive effects of these major QTLs showed that favorable alleles significantly delayed HD. There was a clear trend of increasing HD delay as the number of favorable alleles increased. Among them, Qhd.saw-3A.8, Qhd.saw-3A.9, and Qhd.saw-4D.3 were identified as novel QTLs with no prior reports of HD QTLs/genes in their respective intervals. Candidate gene analysis highlighted seven highly expressed genes related to Ca2+ transport, hormone signaling, glycosylation, and zinc finger proteins, likely involved in HD regulation. This research elucidates the genetic basis of wheat HD under P stress, providing critical insights for breeding high-yield, stable wheat varieties suited to low-P environments.

Investigating the genomic background of calving-related traits in Canadian Jersey cattle.

Traits related to calving have a significant impact on animal welfare and farm profitability in dairy production systems. Identifying genomic regions associated with calving traits could contribute to refining dairy cattle breeding programs and management practices in the dairy industry. Therefore, the primary objectives of this study were to estimate genetic parameters and perform genome-wide association studies (GWAS) and functional enrichment analyses for stillbirth, gestation length, calf size, and calving ease traits in North American Jersey cattle. A total of 40,503 animals with phenotypic records and 5,398 animals genotyped for 45,101 single nucleotide polymorphisms (SNPs) were included in the analyses. Genetic parameters were estimated based on animal models and Bayesian methods. The effects of SNPs were estimated using the Single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) method. The heritability (standard error) estimates ranged from 0.01 (0.01) for stillbirths (SB) in heifers to 0.11 (0.01) for gestation length (GL) in cows. The genetic correlations ranged from -0.58 (0.11) between calving ease (CE) and SB in heifers to 0.44 (0.14) between calving ease and calf size (CZ) in cows. CE showed the highest genetic correlation between heifers and cows, 0.8 (0.22) respectively. The candidate genes identified, including MTHFR, SERPINA5, IGFBP3, and ZRANB1, are involved in key biological processes and metabolic pathways related to the studied traits. Reducing environmental variation and identifying novel indicators of reproduction traits in the Jersey breed are needed given the low heritability estimates for most traits evaluated in this study. In conclusion, this study provides a characterization of the genetic background of calving-related traits in Jersey cattle. The estimates obtained can be used to improve or build selection indexes in Jersey cattle breeding programs in North America.

Genome-wide association scan reveals the reinforcing effect of nano-potassium in improving the yield and quality of salt-stressed barley via enhancing the antioxidant defense system.

Salinity is one of the major environmental factor that can greatly impact the growth, development, and productivity of barley. Our study aims to detect the natural phenotypic variation of morphological and physiological traits under both salinity and potassium nanoparticles (n-K) treatment. In addition to understanding the genetic basis of salt tolerance in barley is a critical aspect of plant breeding for stress resilience. Therefore, a foliar application of n-K was applied at the vegetative stage for 138 barley accessions to enhance salt stress resilience. Interestingly, barley accessions showed high significant increment under n-K treatment compared to saline soil. Based on genome-wide association studies (GWAS) analysis, causative alleles /reliable genomic regions were discovered underlying improved salt resilience through the application of potassium nanoparticles. On chromosome 2H, a highly significant QTN marker (A:C) was located at position 36,665,559 bp which is associated with APX, AsA, GSH, GS, WGS, and TKW under n-K treatment. Inside this region, our candidate gene is HORVU.MOREX.r3.2HG0111480 that annotated as NAC domain protein. Allelic variation detected that the accessions carrying C allele showed higher antioxidants (APX, AsA, and GSH) and barley yield traits (GS, WGS, and TKW) than the accessions carrying A allele, suggesting a positive selection of the accessions carrying C allele that could be used to develop barley varieties with improved salt stress resilience.

Genetic Diversity, Selection Signatures, and Genome-Wide Association Study Identify Candidate Genes Related to Litter Size in Hu Sheep.

Hu sheep is a renowned prolific local sheep breed in China, widely distributed across the country due to its excellent reproductive performance. Deciphering the molecular mechanisms underlying the high fecundity of Hu sheep is crucial for improving the litter size of ewes. In this study, we genotyped 830 female Hu sheep using the Illumina OvineSNP50 BeadChip and performed genetic diversity analysis, selection signature detection, and a genome-wide association study (GWAS) for litter size. Our results revealed that the Hu sheep population exhibits relatively high genetic diversity. A total of 4927 runs of homozygosity (ROH) segments were detected, with the majority (74.73%) being short in length. Different genomic inbreeding coefficients (FROH, FHOM, FGRM, and FUNI) ranged from -0.0060 to 0.0126, showing low levels of inbreeding in this population. Additionally, we identified 91 candidate genomic regions through three complementary selection signature methods, including ROH, composite likelihood ratio (CLR), and integrated haplotype score (iHS), and annotated 189 protein-coding genes. Moreover, we observed two significant SNPs related to the litter size of Hu sheep using GWAS analysis based on a repeatability model. Integrating the selection signatures and the GWAS results, we identified 15 candidate genes associated with litter size, among which BMPR1B and UNC5C were particularly noteworthy. These findings provide valuable insights for improving the reproductive performance and breeding of high-fecundity lines of Hu sheep.

Identification of quantitative trait loci for lodging and related agronomic traits in soybean (Glycine max [L.] Merr.).

BACKGROUND: Lodging, a crucial agronomic trait linked to soybean yield, poses a significant challenge in soybean production. Nevertheless, there has been less research on soybean lodging compared to other important agronomic traits, hindering progress in breeding high-yield soybeans. Our goals were to investigate lodging, pinpoint quantitative trait loci (QTL) linked to lodging, and forecast potential candidate genes linked to this trait. To achieve this, we employed a recombinant inbred line (RIL) population derived from a cross between Guizao 1 and B13 (GB) across various environments. RESULTS: The lodging score of the RIL population was found to be significantly positively correlated with flowering time, maturity time, plant height, number of main stem nodes, stem diameter, and internode length, with correlation coefficients ranging from 0.457 to 0.783. A total of 84 QTLs associated with soybean lodging and related traits were identified using the GB population. The contribution of phenotypic variance ranged from 1.26 to 66.87%, with LOD scores ranging from 2.52 to 69.22. Additionally, within these QTLs, a stable major QTL associated with lodging was newly discovered in the GB population. Out of the ten major QTLs associated with other related traits, nine of them were situated within the qLD-4-1 interval of the major lodging score locus, displaying phenotypic variations ranging from 12.10 to 66.87%. Specific alterations in gene expression were revealed through the analysis of resequencing data from the two parental lines, potentially indicating their significant roles in lodging. Subsequently, it was determined through qRT-PCR that four genes are likely to be the major genes controlling soybean lodging. CONCLUSIONS: This study's findings offer valuable insights into the genetic underpinnings of soybean lodging resistance traits. By comprehending the potential genetic factors associated with lodging, this research lays the groundwork for breeding high-yield soybeans with improved lodging resistance.

QTL analysis of sorghum grain traits based on high-density genetic map.

Sorghum grain traits are important agronomic traits directly related to yield and are key factors affecting the brewing process of distill liquor. Exploring the genes controlling those traits is of great significance for understanding the genetic mechanism of sorghum grain development. In this study, we conducted genotyping using Super-GBS technology on a recombinant inbred lines (RILs) population derived from the cross between "BTx623" and "Hongyingzi," consisting of 205 lines. The grain-related traits of the RIL population were investigated in Guiyang, Anshun in Guizhou, and Ledong in Hainan in China. By inclusive composite interval mapping (ICIM) method, a total of 47 quantitative trait locus (QTL) related to four grain traits (thousand grain weight, grain length, grain width, and length-width ratio) were identified across 10 chromosomes. Among them, 20 important QTL were repeatedly detected in multiple traits or environments and distributed on chromosomes 1 (1), 2 (2), 3 (5), 4 (5), 5 (1), 6 (2), 7 (2), 8 (1), and 9 (1). Six candidate genes were identified within the confidence interval of these QTL, and they are homologous to genes controlling rice grain development (OsMADS1, RGG2, OsNST1, SMG1, OsGRF8, and OsAP2-39). The results provide a basis for further cloning and functional verification of these candidate genes.

Identification of Salt-Stress-Responding Genes by Weighted Gene Correlation Network Analysis and Association Analysis in Wheat Leaves.

The leaf is not only the main site of photosynthesis, but also an important organ reflecting plant salt tolerance. Discovery of salt-stress-responding genes in the leaf is of great significance for the molecular improvement of salt tolerance in wheat varieties. In this study, transcriptome sequencing was conducted on the leaves of salt-tolerant wheat germplasm CH7034 seedlings at 0, 1, 6, 24, and 48 h after NaCl treatment. Based on weighted gene correlation network analysis of differentially expressed genes (DEGs) under salt stress, 12 co-expression modules were obtained, of which, 9 modules containing 4029 DEGs were related to the salt stress time-course. These DEGs were submitted to the Wheat Union database, and a total of 904,588 SNPs were retrieved from 114 wheat germplasms, distributed on 21 wheat chromosomes. Using the R language package and GAPIT program, association analysis was performed between 904,588 SNPs and leaf salt injury index of 114 wheat germplasms. The results showed that 30 single nucleotide polymorphisms (SNPs) from 15 DEGs were associated with salt tolerance. Then, nine candidate genes, including four genes (TaBAM, TaPGDH, TaGluTR, and TaAAP) encoding enzymes as well as five genes (TaB12D, TaS40, TaPPR, TaJAZ, and TaWRKY) encoding functional proteins, were identified by converting salt tolerance-related SNPs into Kompetitive Allele-Specifc PCR (KASP) markers for validation. Finally, interaction network prediction was performed on TaBAM and TaAAP, both belonging to the Turquoise module. Our results will contribute to a further understanding of the salt stress response mechanism in plant leaves and provide candidate genes and molecular markers for improving salt-tolerant wheat varieties.

Comparative Analysis of Transcriptomes Reveals Pathways and Verifies Candidate Genes for Clubroot Resistance in Brassica oleracea.

Clubroot, a soil-borne disease caused by Plasmodiophora brassicae, is one of the most destructive diseases of Brassica oleracea all over the world. However, the mechanism of clubroot resistance remains unclear. In this research, transcriptome sequencing was conducted on root samples from both resistant (R) and susceptible (S) B. oleracea plants infected by P. brassicae. Then the comparative analysis was carried out between the R and S samples at different time points during the infection stages to reveal clubroot resistance related pathways and candidate genes. Compared with 0 days after inoculation, a total of 4991 differential expressed genes were detected from the S pool, while only 2133 were found from the R pool. Gene function enrichment analysis found that the effector-triggered immunity played a major role in the R pool, while the pathogen-associated molecular pattern triggered immune response was stronger in the S pool. Simultaneously, candidate genes were identified through weighted gene co-expression network analysis, with Bol010786 (CNGC13) and Bol017921 (SD2-5) showing potential for conferring resistance to clubroot. The findings of this research provide valuable insights into the molecular mechanisms underlying clubroot resistance and present new avenues for further research aimed at enhancing the clubroot resistance of B. oleracea through breeding.

QTL Mapping and Candidate Gene Analysis for Starch-Related Traits in Tartary Buckwheat (Fagopyrum tataricum (L.) Gaertn).

Starch is the main component that determines the yield and quality of Tartary buckwheat. As a quantitative trait, using quantitative trait locus (QTL) mapping to excavate genes associated with starch-related traits is crucial for understanding the genetic mechanisms involved in starch synthesis and molecular breeding of Tartary buckwheat varieties with high-quality starch. Employing a recombinant inbred line population as research material, this study used QTL mapping to investigate the amylose, amylopectin, and total starch contents across four distinct environments. The results identified a total of 20 QTLs spanning six chromosomes, which explained 4.07% to 14.41% of the phenotypic variation. One major QTL cluster containing three stable QTLs governing both amylose and amylopectin content, qClu-4-1, was identified and located in the physical interval of 39.85-43.34 Mbp on chromosome Ft4. Within this cluster, we predicted 239 candidate genes and analyzed their SNP/InDel mutations, expression patterns, and enriched KEGG pathways. Ultimately, five key candidate genes, namely FtPinG0004897100.01, FtPinG0002636200.01, FtPinG0009329200.01, FtPinG0007371600.01, and FtPinG0005109900.01, were highlighted, which are potentially involved in starch synthesis and regulation, paving the way for further investigative studies. This study, for the first time, utilized QTL mapping to detect major QTLs controlling amylose, amylopectin, and total starch contents in Tartary buckwheat. The QTLs and candidate genes would provide valuable insights into the genetic mechanisms underlying starch synthesis and improving starch-related traits of Tartary buckwheat.

Exploring the Genetic Landscape of Vitiligo in the Pura Raza Español Horse: A Genomic Perspective.

Vitiligo is a depigmentation autoimmune disorder characterized by the progressive loss of melanocytes leading to the appearance of patchy depigmentation of the skin. The presence of vitiligo in horses is greater in those with grey coats. The aim of this study was therefore to perform a genome-wide association study (GWAS) to identify genomic regions and putative candidate loci associated with vitiligo depigmentation and susceptibility in the Pura Raza Español population. For this purpose, we performed a wssGBLUP (weighted single step genomic best linear unbiased prediction) using data from a total of 2359 animals genotyped with Affymetrix Axiom™ Equine 670 K and 1346 with Equine GeneSeek Genomic Profiler™ (GGP) Array V5. A total of 60,136 SNPs (single nucleotide polymorphisms) present on the 32 chromosomes from the consensus dataset after quality control were employed for the analysis. Vitiligo-like depigmentation was phenotyped by visual inspection of the different affected areas (eyes, mouth, nostrils) and was classified into nine categories with three degrees of severity (absent, slight, and severe). We identified one significant genomic region for vitiligo around the eyes, eight significant genomic regions for vitiligo around the mouth, and seven significant genomic regions for vitiligo around the nostrils, which explained the highest percentage of variance. These significant genomic regions contained candidate genes related to melanocytes, skin, immune system, tumour suppression, metastasis, and cutaneous carcinoma. These findings enable us to implement selective breeding strategies to decrease the incidence of vitiligo and to elucidate the genetic architecture underlying vitiligo in horses as well as the molecular mechanisms involved in the disease's development. However, further studies are needed to better understand this skin disorder in horses.

Imputation strategies for low-coverage whole-genome sequencing data and their effects on genomic prediction and genome-wide association studies in pigs.

The uncertainty resulting from missing genotypes in low-coverage whole-genome sequencing (LCWGS) data complicates genotype imputation. The aim of this study is to find out an optimal strategy for accurately imputing LCWGS data and assess its effectiveness for genomic prediction (GP) and genome-wide association study (GWAS) on economically important traits of Large White pigs. The LCWGS data of 1 423 Large White pigs were imputed using three different strategies: (1) using the high-coverage whole-genome sequencing (HCWGS) of 30 key progenitors as the reference panel (Ref_LG); (2) mixing HCWGS of key progenitors with LCWGS (Mix_HLG) and (3) self-imputation in LCWGS (Within_LG). Additionally, to compare the imputation effects of LCWGS, we also imputed SNP chip data of 1 423 Large White pigs to the whole-genome sequencing level using the reference panel consisting of key progenitors (Ref_SNP). To evaluate effects of the imputed sequencing data, we compared the accuracies of GP and statistical power of GWAS for four reproductive traits based on the chip data, sequencing data imputed from chip data and LCWGS data using an optimal strategy. The average imputation accuracies of the Within_LG, Ref_LG and Mix_HLG were 0.9893, 0.9899 and 0.9875, respectively, which were higher than that of the Ref_SNP (0.8522). Using the imputed sequencing data from LCWGS with the Ref_LG imputation strategy, the accuracies of GP for four traits improved by approximately 0.31-1.04% compared to the chip data, and by 0.7-1.05% compared to the imputed sequencing data from chip data. Furthermore, by using the sequence data imputed from LCWGS with the Ref_LG, 18 candidate genes were identified to be associated with the four reproductive traits of interest in Large White pigs: total number of piglets born - EPC2, MBD5, ORC4 and ACVR2A; number of piglets born healthy - IKBKE; total litter weight of piglets born alive - HSPA13 and CPA1; gestation length - GTF2H5, ITGAV, NFE2L2, CALCRL, ITGA4, STAT1, HOXD10, MSTN, COL5A2 and STAT4. With the exception of EPC2, ORC4, ACVR2A and MSTN, others represent novel candidates. Our findings can provide a reference for the application of LCWGS data in livestock and poultry.

Genome-wide association study for high-temperature tolerance in the Japanese flounder.

This study addresses the critical issue of high-temperature stress in Japanese flounder (Paralichthys olivaceus), a factor threatening both their survival and the growth of the aquaculture industry. The research aims to identify genetic markers associated with high-temperature tolerance, unravel the genetic regulatory mechanisms, and lay the foundation for breeding Japanese flounder with increased resistance to high temperatures. In this study, using a genome-wide association study was performed to identify single nucleotide polymorphisms (SNPs) and genes associated with high-temperature tolerance for Japanese flounder using 280 individuals with 342 311 high-quality SNPs. The traits of high-temperature tolerance were defined as the survival time and survival status of Japanese flounder at high water temperature (31℃) for 15 days cultivate. A genome-wide association study identified six loci on six chromosomes significantly correlated with survival time under high-temperature stress. Six candidate genes were successfully annotated. Additionally, 34 loci associated with survival status were identified and mapped to 15 chromosomes, with 22 candidate genes annotated. Functional analysis highlighted the potential importance of genes like traf4 and ppm1l in regulating apoptosis, impacting high-temperature tolerance in Japanese flounder. These findings provide a valuable theoretical framework for integrating molecular markers into Japanese flounder breeding programmes, serving as a molecular tool to enhance genetic traits linked to high-temperature tolerance in cultured Japanese flounder.

The genetic characterization of germplasm and identification of the litter size trait associated candidate genes in Dexin mutton and fine-wool sheep.

Xinjiang is a major province of sheep breeding in China, which plays an important role in meeting people's needs for meat products, increasing farmers' income and sustainable development of animal husbandry. However, the genetic differentiation relationship between breeds was not clear, and most sheep had low fecundity, which seriously restricted the efficient development of sheep industry. Therefore, this study used the whole genome resequencing to detect the genetic variation of Dexin mutton and fine-wool sheep, explored the selected regions and important genes of the litter size traits, analyzed the genetic mechanism of reproductive traits, and provided new insights for the high fecundity breeding of sheep. A total of 5,236.338 G genome data and 35,884,037 SNPs were obtained. Furthermore, we identified 39 selection signals spanning candidate genes, 99 genes were significantly associated related to growth, reproduction and immunity, among which, BRIP1, BMPR1B, BMP4, NGF, etc. genes, and MAKP signaling pathway, Fanconi anemia pathway and Thyroid hormone signaling pathway and other signaling pathways were significantly correlated with litter size trait. Among them, we identified NGF, TrKA and BRIP1 genes was the important genes for sheep litter size traits and the mutation frequencies of 9 SNPs in BRIP1 gene were significantly different in domestic sheep in the world. The research provided new insights for the breeding of self-cultivated meat fine-wool sheep.

Fine Mapping and Candidate Gene Analysis of Two Major Quantitative Trait Loci, qFW2.1 and qFW3.1, Controlling Fruit Weight in Pepper (Capsicum annuum).

Fruit weight is an important agronomic trait in pepper production and is closely related to yield. At present, many quantitative trait loci (QTL) related to fruit weight have been found in pepper; however, the genes affecting fruit weight remain unknown. We analyzed the fruit weight-related quantitative traits in an intraspecific Capsicum annuum cross between the cultivated species blocky-type pepper, cv. Qiemen, and the bird pepper accession, "129-1" (Capsicum annuum var. glatriusculum), which was the wild progenitor of C. annuum. Using the QTL-seq combined with the linkage-based QTL mapping approach, QTL detection was performed; and two major effects of QTL related to fruit weight, qFW2.1 and qFW3.1, were identified on chromosomes 2 and 3. The qFW2.1 maximum explained 12.28% of the phenotypic variance observed in two F2 generations, with the maximum LOD value of 11.02, respectively; meanwhile, the qFW3.1 maximum explained 15.50% of the observed phenotypic variance in the two F2 generations, with the maximum LOD value of 11.36, respectively. qFW2.1 was narrowed down to the 1.22 Mb region using homozygous recombinant screening from BC2S2 and BC2S3 populations, while qFW3.1 was narrowed down to the 4.61Mb region. According to the transcriptome results, a total of 47 and 86 differentially expressed genes (DEGs) in the candidate regions of qFW2.1 and qFW3.1 were identified. Further, 19 genes were selected for a qRT-PCR analysis based on sequence difference combined with the gene annotation. Finally, Capana02g002938 and Capana02g003021 are the most likely candidate genes for qFW2.1, and Capana03g000903 may be a candidate gene for qFW3.1. Taken together, our results identified and fine-mapped two major QTL for fruit weight in pepper that will facilitate marker-assistant breeding for the manipulation of yield in pepper.

Exploring the Molecular Mechanism of Skeletal Muscle Development in Ningxiang Pig by Weighted Gene Co-Expression Network Analysis.

With the continuous improvement in living standards, people's demand for high-quality meat is increasing. Ningxiang pig has delicious meat of high nutritional value, and is loved by consumers. However, its slow growth and low meat yield seriously restrict its efficient utilization. Gene expression is the internal driving force of life activities, so in order to fundamentally improve its growth rate, it is key to explore the molecular mechanism of skeletal muscle development in Ningxiang pigs. In this paper, Ningxiang boars were selected in four growth stages (30 days: weaning period, 90 days: nursing period, 150 days: early fattening period, and 210 days: late fattening period), and the longissimus dorsi (LD) muscle was taken from three boars in each stage. The fatty acid content, amino acid content, muscle fiber diameter density and type of LD were detected by gas chromatography, acidolysis, hematoxylin eosin (HE) staining and immunofluorescence (IF) staining. After transcription sequencing, weighted gene co-expression network analysis (WGCNA) combined with the phenotype of the LD was used to explore the key genes and signaling pathways affecting muscle development. The results showed that 10 modules were identified by WGCNA, including 5 modules related to muscle development stage, module characteristics of muscle fiber density, 5 modules characteristic of muscle fiber diameter, and a module characteristic of palmitoleic acid (C16:1) and linoleic acid (C18:2n6C). Gene ontology (GO) enrichment analysis found that 52 transcripts relating to muscle development were enriched in these modules, including 44 known genes and 8 novel genes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these genes were enriched in the auxin, estrogen and cyclic guanosine monophosphate-protein kinase G (cGMP-PKG) pathways. Twelve of these genes were transcription factors, there were interactions among 20 genes, and the interactions among 11 proteins in human, pig and mouse were stable. To sum up, through the integrated analysis of phenotype and transcriptome, this paper analyzed the key genes and possible regulatory networks of skeletal muscle development in Ningxiang pigs at various stages, to provide a reference for the in-depth study of skeletal muscle development.