RNA-Seq Reveals Pathways Responsible for Meat Quality Characteristic Differences between Two Yunnan Indigenous Chicken Breeds and Commercial Broilers.

Poultry is a source of meat that is in great demand in the world. The quality of meat is an imperative point for shoppers. To explore the genes controlling meat quality characteristics, the growth and meat quality traits and muscle transcriptome of two indigenous Yunnan chicken breeds, Wuding chickens (WDs) and Daweishan mini chickens (MCs), were compared with Cobb broilers (CBs). The growth and meat quality characteristics of these two indigenous breeds were found to differ from CB. In particular, the crude fat (CF), inosine monophosphate content, amino acid (AA), and total fatty acid (TFA) content of WDs were significantly higher than those of CBs and MCs. In addition, it was found that MC pectoralis had 420 differentially expressed genes (DEGs) relative to CBs, and WDs had 217 DEGs relative to CBs. Among them, 105 DEGs were shared. The results of 10 selected genes were also confirmed by qPCR. The differentially expressed genes were six enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) biological pathways including lysosomes, phagosomes, PPAR signaling pathways, cell adhesion molecules, cytokine-cytokine receptor interaction, and phagosome sphingolipid metabolism. Interestingly, four genes (LPL, GK, SCD, and FABP7) in the PPAR signal pathway related to fatty acid (FA) metabolism were elevated in WD muscles, which may account for higher CF, inosine monophosphate content, and AA and FA contents, key factors affecting meat quality. This work laid the foundation for improving the meat quality of Yunnan indigenous chickens, especially WD. In future molecular breeding, the genes in this study can be used as molecular screening markers and applied to the molecular breeding of chicken quality characteristics.

Association analysis of transcriptome and quasi-targeted metabolomics reveals the regulation mechanism underlying broiler muscle tissue development at different levels of dietary guanidinoacetic acid.

The development and characteristics of muscle fibers in broilers are critical determinants that influence their growth performance, as well as serve as essential prerequisites for the production of high-quality chicken meat. Guanidinoacetic acid (GAA) is a crucial endogenous substance in animal creatine synthesis, and its utilization as a feed additive has been demonstrated the capabilities to enhance animal performance, optimize muscle yield, and augment carcass quality. The objective of this study was to investigate the regulation and molecular mechanism underlying muscle development in broilers at different levels of GAA via multiple omics analysis. The 90 Cobb broilers, aged 1 day, were randomly allocated into three treatments consisting of five replicates of six chickens each. The control group was provided with a basal diet, while the Normal GAA and High GAA groups received a basal diet supplemented with 1.2 g/kg and 3.6 g/kg of GAA, respectively. After a feeding period of 42 days, the pectoralis muscles were collected for histomorphological observation, transcriptome and metabolomic analysis. The results demonstrated that the addition of 1.2 g/kg GAA in the diet led to an augmentation in muscle fiber diameter and up-regulation of IGF1, IHH, ASB2, and ANKRD2 gene expression. However, a high dose of 3.6 g/kg GAA in the diet potentially reversed the beneficial effects on chicken breast development by excessively activating the TGF-ß signaling pathway and reducing nucleotide metabolite content. These findings would provide a theoretical foundation for enhancing the performance and meat quality of broilers by incorporating GAA as a feed additive.

Annotations of four high-quality indigenous chicken genomes identify more than one thousand missing genes in subtelomeric regions and micro-chromosomes with high G/C contents.

BACKGROUND: Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches. RESULTS: We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought. CONCLUSION: The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.

Analysis of microRNA Expression Profiles in Broiler Muscle Tissues by Feeding Different Levels of Guanidinoacetic Acid.

The aim of this study was to explore the molecular mechanisms through which different levels of GAA affect chicken muscle development by influencing miRNA expression, to lay a theoretical foundation for the identification of key functional small RNAs related to poultry muscle development, and to provide new insights into the regulatory mechanisms of GAA on muscle development and meat quality in broilers. It provides a new theoretical basis for using GAA as a feed additive to improve feed performance. Small RNA sequencing technology was utilized to obtain the expression profiles of miRNA in the broiler pectoral muscle fed with different levels of GAA (0 g/kg, 1.2 g/kg and 3.6 g/kg). An analysis of differentially expressed miRNAs revealed 90 such miRNAs in the three combination comparisons, with gga-miR-130b-5p exhibiting significant differences across all three combinations. Furthermore, three of the differentially expressed miRNAs were performed by RT-qPCR verification, yielding results consistent with those obtained from small RNA sequencing. Target gene prediction, as well as the GO and KEGG enrichment analysis of differentially expressed miRNAs, indicated their involvement in muscle cell differentiation and other processes, particularly those associated with the MAPK signaling pathway. This study has, thus, provided valuable insights and resources for the further exploration of the miRNA molecular mechanism underlying the influence of guanidine acetic acid on broiler muscle development. Combined with previous studies and small RNA sequencing, adding 1.2 g/kg GAA to the diet can better promote the muscle development of broilers.

Artificial selection footprints in indigenous and commercial chicken genomes.

BACKGROUND: Although many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China. RESULTS: We found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds. CONCLUSION: We identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.

High quality assemblies of four indigenous chicken genomes and related functional data resources.

Many lines of evidence indicate that red jungle fowl (RJF) is the primary ancestor of domestic chickens. Although multiple versions of RJF (galgal2-galgal5 and GRCg6a) and commercial chickens (GRCg7b/w and Huxu) genomes have been assembled since 2004, no high-quality indigenous chicken genomes have been assembled, hampering the understanding of chicken domestication and evolution. To fill the gap, we sequenced the genomes of four indigenous chickens with distinct morphological traits in southwest China, using a combination of short, long and Hi-C reads. We assembled each genome (~1.0 Gb) into 42 chromosomes with chromosome N50 90.5-90.9 Mb, amongst the highest quality of chicken genome assemblies. To provide resources for gene annotation and functional analysis, we also sequenced transcriptomes of 10 tissues for each of the four chickens. Moreover, we corrected many mis-assemblies and assembled missing micro-chromosomes 29 and 34-39 for GRCg6a. Our assemblies, sequencing data and the correction of GRCg6a can be valuable resources for studying chicken domestication and evolution.

High-quality genome assembly of a C. crossoptilon and related functional and genetics data resources.

There are four species in the Crossoptilon genus inhibiting at from very low to very high altitudes across China, and they are in varying levels of danger of extinction. To better understand the genetic basis of adaptation to high altitudes and genetic changes due to bottleneck, we assembled the genome (~1.02 Gb) of a white eared pheasant (WT) (Crossoptilon crossoptilon) inhibiting at high altitudes (3,000~7,000 m) in northwest of Yunnan province, China, using a combination of Illumina short reads, PacBio long reads and Hi-C reads, with a contig N50 of 19.63 Mb and only six gaps. To further provide resources for gene annotation as well as functional and population genetics analyses, we sequenced transcriptomes of 20 major tissues of the WT individual and re-sequenced another 10 WT individuals and a blue eared pheasant (Crossoptilon auritum) individual inhabiting at intermediate altitudes (1,500~3,000 m). Our assembled WT genome, transcriptome data, and DNA sequencing data can be valuable resources for studying the biology, evolution and developing conservation strategies of these endangered species.

Integrated multi-omics analysis reveals gut microbiota dysbiosis and systemic disturbance in major depressive disorder.

Major depressive disorder (MDD) involves systemic changes in peripheral blood and gut microbiota, but the current understanding is incomplete. Herein, we conducted a multi-omics analysis of fecal and blood samples obtained from an observational cohort including MDD patients (n = 99) and healthy control (HC, n = 50). 16S rRNA sequencing of gut microbiota showed structural alterations in MDD, as characterized by increased Enterococcus. Metagenomics sequencing of gut microbiota showed substantial functional alterations including upregulation in the superpathway of the glyoxylate cycle and fatty acid degradation and downregulation in various metabolic pathways in MDD. Plasma metabolomics revealed decreased amino acids and bile acids, together with increased sphingolipids and cholesterol esters in MDD. Notably, metabolites involved in arginine and proline metabolism were decreased while sphingolipid metabolic pathway were increased. Mass cytometry analysis of blood immune cell subtypes showed rises in proinflammatory immune subsets and declines in anti-inflammatory immune subsets in MDD. Furthermore, our findings revealed disease severity-related factors of MDD. Interestingly, we classified MDD into two immune subtypes that were highly correlated with disease relapse. Moreover, we established discriminative signatures that differentiate MDD from HC. These findings contribute to a comprehensive understanding of the MDD pathogenesis and provide valuable resources for the discovery of biomarkers.

Molecular detection of SARS-CoV-2 being challenged by virus variation and asymptomatic infection / 药物分析学报

Coronavirus disease 2019 (COVID-19) has been a pandemic for more than a year.With the expanding second wave of the pandemic in winter,the continuous evolution of SARS-CoV-2 has brought new issues,including the significance of virus mutations in infection and the detection of asymptomatic infection.In this review,we first introduced several major SARS-CoV-2 mutations since the COVID-19 outbreak and then mentioned the widely used molecular detection techniques to diagnose COVID-19,primarily focusing on their strengths and limitations.We further discussed the effects of viral genetic variation and asymptomatic infection on the molecular detection of SARS-CoV-2 infection.The review finally sum-marized useful insights into the molecular diagnosis of COVID-19 under the special situation being challenged by virus mutation and asymptomatic infection.