BACKGROUND: Previous studies have identified that chromosome structure plays a very important role in gene control. The transcription factor Yin Yang 1 (YY1), a multifunctional DNA binding protein, could form a dimer to mediate chromatin loops and active enhancer-promoter interactions. The deletion of YY1 or point mutations at the YY1 binding sites significantly inhibit the enhancer-promoter interactions and affect gene expression. To date, only a few computational methods are available for identifying YY1-mediated chromatin loops. RESULTS: We proposed a novel model named CapsNetYY1, which was based on capsule network architecture to identify whether a pair of YY1 motifs can form a chromatin loop. Firstly, we encode the DNA sequence using one-hot encoding method. Secondly, multi-scale convolution layer is used to extract local features of the sequence, and bidirectional gated recurrent unit is used to learn the features across time steps. Finally, capsule networks (convolution capsule layer and digital capsule layer) used to extract higher level features and recognize YY1-mediated chromatin loops. Compared with DeepYY1, the only prediction for YY1-mediated chromatin loops, our model CapsNetYY1 achieved the better performance on the independent datasets (AUC [Formula: see text]). CONCLUSION: The results indicate that CapsNetYY1 is an excellent method for identifying YY1-mediated chromatin loops. We believe that the CapsNetYY1 method will be used for predictive classification of other DNA sequences.

BACKGROUND: Disease molecular complexity requires high throughput workflows to map disease pathways through analysis of vast tissue repositories. Great progress has been made in tissue multiomics analytical technologies. To match the high throughput of these advanced analytical platforms, we have previously developed a multipurpose 96-well microplate sonicator, PIXUL, that can be used in multiple workflows to extract analytes from cultured cells and tissue fragments for various downstream molecular assays. And yet, the sample preparation devices, such as PIXUL, along with the downstream multiomics analytical capabilities have not been fully exploited to interrogate tissues because storing and sampling of such biospecimens remain, in comparison, inefficient. RESULTS: To mitigate this tissue interrogation bottleneck, we have developed a low-cost user-friendly system, CryoGrid, to catalog, cryostore and sample tissue fragments. TRIzol is widely used to isolate RNA but it is labor-intensive, hazardous, requires fume-hoods, and is an expensive reagent. Columns are also commonly used to extract RNA but they involve many steps, are prone to human errors, and are also expensive. Both TRIzol and column protocols use test tubes. We developed a microplate PIXUL-based TRIzol-free and column-free RNA isolation protocol that uses a buffer containing proteinase K (PK buffer). We have integrated the CryoGrid system with PIXUL-based PK buffer, TRIzol, and PureLink column methods to isolate RNA for gene-specific qPCR and genome-wide transcript analyses. CryoGrid-PIXUL, when integrated with either PK buffer, TRIzol or PureLink column RNA isolation protocols, yielded similar transcript profiles in frozen organs (brain, heart, kidney and liver) from a mouse model of sepsis. CONCLUSIONS: RNA isolation using the CryoGrid-PIXUL system combined with the 96-well microplate PK buffer method offers an inexpensive user-friendly high throughput workflow to study transcriptional responses in tissues in health and disease as well as in therapeutic interventions.

BACKGROUND: Single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) are the most abundant genetic variations and widely distribute across the genomes in plant. Development of SNP and InDel markers is a valuable tool for genetics and genomic research in radish (Raphanus sativus L.). RESULTS: In this study, a total of 366,679 single nucleotide polymorphisms (SNPs) and 97,973 insertion-deletion (InDel) markers were identified based on genome resequencing between 'YZH' and 'XHT'. In all, 53,343 SNPs and 4,257 InDels were detected in two cultivars by transcriptome sequencing. Among the InDel variations, 85 genomic and 15 transcriptomic InDels were newly developed and validated PCR. The 100 polymorphic InDels markers generated 207 alleles among 200 Chinese radish germplasm, with an average 2.07 of the number of alleles (Na) and with an average 0.33 of the polymorphism information content (PIC). Population structure and phylogenetic relationship revealed that the radish cultivars from northern China were clustered together and the southwest China cultivars were clustered together. RNA-Seq analysis revealed that 11,003 differentially expressed genes (DEGs) were identified between the two cultivars, of which 5,020 were upregulated and 5,983 were downregulated. In total, 145 flowering time-related DGEs were detected, most of which were involved in flowering time integrator, circadian clock/photoperiod autonomous, and vernalization pathways. In flowering time-related DGEs region, 150 transcriptomic SNPs and 9 InDels were obtained. CONCLUSIONS: The large amount of SNPs and InDels identified in this study will provide a valuable marker resource for radish genetic and genomic studies. The SNPs and InDels within flowering time-related DGEs provide fundamental insight into for dissecting molecular mechanism of bolting and flowering in radish.

BACKGROUND: Lonicera japonica Thunb. is widely used in traditional Chinese medicine. Medicinal L. japonica mainly consists of dried flower buds and partially opened flowers, thus flowers are an important quality indicator. MADS-box genes encode transcription factors that regulate flower development. However, little is known about these genes in L. japonica. RESULTS: In this study, 48 MADS-box genes were identified in L. japonica, including 20 Type-I genes (8 Mα, 2 Mβ, and 10 Mγ) and 28 Type-II genes (26 MIKCc and 2 MIKC*). The Type-I and Type-II genes differed significantly in gene structure, conserved domains, protein structure, chromosomal distribution, phylogenesis, and expression pattern. Type-I genes had a simpler gene structure, lacked the K domain, had low protein structure conservation, were tandemly distributed on the chromosomes, had more frequent lineage-specific duplications, and were expressed at low levels. In contrast, Type-II genes had a more complex gene structure; contained conserved M, I, K, and C domains; had highly conserved protein structure; and were expressed at high levels throughout the flowering period. Eleven floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in L. japonica. By integrating expression pattern and protein interaction data for these genes, we developed a possible model for floral organ identity determination. CONCLUSION: This study genome-widely identified and characterized the MADS-box gene family in L. japonica. Eleven floral homeotic MADS-box genes were identified and a possible model for floral organ identity determination was also developed. This study contributes to our understanding of the MADS-box gene family and its possible involvement in floral organ development in L. japonica.

BACKGROUND: Long non-coding RNAs (lncRNAs) are defined as transcribed molecules longer than 200 nucleotides with little to no protein-coding potential. LncRNAs can regulate gene expression of nearby genes (cis-acting) or genes located on other chromosomes (trans-acting). Several methodologies have been developed to capture lncRNAs associated with chromatin at a genome-wide level. Analysis of RNA-DNA contacts can be combined with epigenetic and RNA-seq data to define potential lncRNAs involved in the regulation of gene expression. RESULTS: We performed Chromatin Associated RNA sequencing (ChAR-seq) in Anolis carolinensis to obtain the genome-wide map of the associations that RNA molecules have with chromatin. We analyzed the frequency of DNA contacts for different classes of RNAs and were able to define cis- and trans-acting lncRNAs. We integrated the ChAR-seq map of RNA-DNA contacts with epigenetic data for the acetylation of lysine 16 on histone H4 (H4K16ac), a mark connected to actively transcribed chromatin in lizards. We successfully identified three trans-acting lncRNAs significantly associated with the H4K16ac signal, which are likely involved in the regulation of gene expression in A. carolinensis. CONCLUSIONS: We show that the ChAR-seq method is a powerful tool to explore the RNA-DNA map of interactions. Moreover, in combination with epigenetic data, ChAR-seq can be applied in non-model species to establish potential roles for predicted lncRNAs that lack functional annotations.

BACKGROUND: Morphological and traditional genetic studies of the young Pliocene genus Hyles have led to the understanding that despite its importance for taxonomy, phenotypic similarity of wing patterns does not correlate with phylogenetic relationship. To gain insights into various aspects of speciation in the Spurge Hawkmoth (Hyles euphorbiae), we assembled a chromosome-level genome and investigated some of its characteristics. RESULTS: The genome of a male H. euphorbiae was sequenced using PacBio and Hi-C data, yielding a 504 Mb assembly (scaffold N50 of 18.2 Mb) with 99.9% of data represented by the 29 largest scaffolds forming the haploid chromosome set. Consistent with this, FISH analysis of the karyotype revealed n = 29 chromosomes and a WZ/ZZ (female/male) sex chromosome system. Estimates of chromosome length based on the karyotype image provided an additional quality metric of assembled chromosome size. Rescaffolding the published male H. vespertilio genome resulted in a high-quality assembly (651 Mb, scaffold N50 of 22 Mb) with 98% of sequence data in the 29 chromosomes. The larger genome size of H. vespertilio (average 1C DNA value of 562 Mb) was accompanied by a proportional increase in repeats from 45% in H. euphorbiae (measured as 472 Mb) to almost 55% in H. vespertilio. Several wing pattern genes were found on the same chromosomes in the two species, with varying amounts and positions of repetitive elements and inversions possibly corrupting their function. CONCLUSIONS: Our two-fold comparative genomics approach revealed high gene synteny of the Hyles genomes to other Sphingidae and high correspondence to intact Merian elements, the ancestral linkage groups of Lepidoptera, with the exception of three simple fusion events. We propose a standardized approach for genome taxonomy using nucleotide homology via scaffold chaining as the primary tool combined with Oxford plots based on Merian elements to infer and visualize directionality of chromosomal rearrangements. The identification of wing pattern genes promises future understanding of the evolution of forewing patterns in the genus Hyles, although further sequencing data from more individuals are needed. The genomic data obtained provide additional reliable references for further comparative studies in hawkmoths (Sphingidae).

BACKGROUND: The sugarcane aphid (SCA; Melanaphis sacchari) has emerged as a key pest on sorghum in the United States that feeds from the phloem tissue, drains nutrients, and inflicts physical damage to plants. Previously, it has been shown that SCA reproduction was low and high on sorghum SC265 and SC1345 plants, respectively, compared to RTx430, an elite sorghum male parental line (reference line). In this study, we focused on identifying the defense-related genes that confer resistance to SCA at early and late time points in sorghum plants with varied levels of SCA resistance. RESULTS: We used RNA-sequencing approach to identify the global transcriptomic responses to aphid infestation on RTx430, SC265, and SC1345 plants at early time points 6, 24, and 48 h post infestation (hpi) and after extended period of SCA feeding for 7 days. Aphid feeding on the SCA-resistant line upregulated the expression of 3827 and 2076 genes at early and late time points, respectively, which was relatively higher compared to RTx430 and SC1345 plants. Co-expression network analysis revealed that aphid infestation modulates sorghum defenses by regulating genes corresponding to phenylpropanoid metabolic pathways, secondary metabolic process, oxidoreductase activity, phytohormones, sugar metabolism and cell wall-related genes. There were 187 genes that were highly expressed during the early time of aphid infestation in the SCA-resistant line, including genes encoding leucine-rich repeat (LRR) proteins, ethylene response factors, cell wall-related, pathogenesis-related proteins, and disease resistance-responsive dirigent-like proteins. At 7 days post infestation (dpi), 173 genes had elevated expression levels in the SCA-resistant line and were involved in sucrose metabolism, callose formation, phospholipid metabolism, and proteinase inhibitors. CONCLUSIONS: In summary, our results indicate that the SCA-resistant line is better adapted to activate early defense signaling mechanisms in response to SCA infestation because of the rapid activation of the defense mechanisms by regulating genes involved in monolignol biosynthesis pathway, oxidoreductase activity, biosynthesis of phytohormones, and cell wall composition. This study offers further insights to better understand sorghum defenses against aphid herbivory.

BACKGROUND: Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource. RESULTS: We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue. CONCLUSIONS: By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

BACKGROUND: Biocontrol is a key technology for the control of pest species. Microctonus parasitoid wasps (Hymenoptera: Braconidae) have been released in Aotearoa New Zealand as biocontrol agents, targeting three different pest weevil species. Despite their value as biocontrol agents, no genome assemblies are currently available for these Microctonus wasps, limiting investigations into key biological differences between the different species and strains. METHODS AND FINDINGS: Here we present high-quality genomes for Microctonus hyperodae and Microctonus aethiopoides, assembled with short read sequencing and Hi-C scaffolding. These assemblies have total lengths of 106.7 Mb for M. hyperodae and 129.2 Mb for M. aethiopoides, with scaffold N50 values of 9 Mb and 23 Mb respectively. With these assemblies we investigated differences in reproductive mechanisms, and association with viruses between Microctonus wasps. Meiosis-specific genes are conserved in asexual Microctonus, with in-situ hybridisation validating expression of one of these genes in the ovaries of asexual Microctonus aethiopoides. This implies asexual reproduction in these Microctonus wasps involves meiosis, with the potential for sexual reproduction maintained. Investigation of viral gene content revealed candidate genes that may be involved in virus-like particle production in M. aethiopoides, as well as a novel virus infecting M. hyperodae, for which a complete genome was assembled. CONCLUSION AND SIGNIFICANCE: These are the first published genomes for Microctonus wasps which have been deployed as biocontrol agents, in Aotearoa New Zealand. These assemblies will be valuable resources for continued investigation and monitoring of these biocontrol systems. Understanding the biology underpinning Microctonus biocontrol is crucial if we are to maintain its efficacy, or in the case of M. hyperodae to understand what may have influenced the significant decline of biocontrol efficacy. The potential for sexual reproduction in asexual Microctonus is significant given that empirical modelling suggests this asexual reproduction is likely to have contributed to biocontrol decline. Furthermore the identification of a novel virus in M. hyperodae highlights a previously unknown aspect of this biocontrol system, which may contribute to premature mortality of the host pest. These findings have potential to be exploited in future in attempt to increase the effectiveness of M. hyperodae biocontrol.

BACKGROUND: DNA methylation is one of the best characterized epigenetic modifications in the mammalian nuclear genome and is known to play a significant role in various biological processes. Nonetheless, the presence of 5-methylcytosine (5mC) in mitochondrial DNA remains controversial, as data ranging from the lack of 5mC to very extensive 5mC have been reported. RESULTS: By conducting comprehensive bioinformatic analyses of both published and our own data, we reveal that previous observations of extensive and strand-biased mtDNA-5mC are likely artifacts due to a combination of factors including inefficient bisulfite conversion, extremely low sequencing reads in the L strand, and interference from nuclear mitochondrial DNA sequences (NUMTs). To reduce false positive mtDNA-5mC signals, we establish an optimized procedure for library preparation and data analysis of bisulfite sequencing. Leveraging our modified workflow, we demonstrate an even distribution of 5mC signals across the mtDNA and an average methylation level ranging from 0.19% to 0.67% in both cell lines and primary cells, which is indistinguishable from the background noise. CONCLUSIONS: We have developed a framework for analyzing mtDNA-5mC through bisulfite sequencing, which enables us to present multiple lines of evidence for the lack of extensive 5mC in mammalian mtDNA. We assert that the data available to date do not support the reported presence of mtDNA-5mC.

BACKGROUND: Oxford Nanopore Technologies (ONT) offers an accessible platform for long-read sequencing, which improves the reconstruction of genomes and helps to resolve complex genomic contexts, especially in the case of metagenome analysis. To take the best advantage of long-read sequencing, DNA extraction methods must be able to isolate pure high molecular weight (HMW) DNA from complex metagenomics samples, without introducing any bias. New methods released on the market, and protocols developed at the research level, were specifically designed for this application and need to be assessed. RESULTS: In this study, with different bacterial cocktail mixes, analyzed as pure or spiked in a synthetic fecal matrix, we evaluated the performances of 6 DNA extraction methods using various cells lysis and purification techniques, from quick and easy, to more time-consuming and gentle protocols, including a portable method for on-site application. In addition to the comparison of the quality, quantity and purity of the extracted DNA, the performance obtained when doing Nanopore sequencing on a MinION flow cell was also tested. From the obtained results, the Quick-DNA HMW MagBead Kit (Zymo Research) was selected as producing the best yield of pure HMW DNA. Furthermore, this kit allowed an accurate detection, by Nanopore sequencing, of almost all the bacterial species present in a complex mock community. CONCLUSION: Amongst the 6 tested methods, the Quick-DNA HMW MagBead Kit (Zymo Research) was considered as the most suitable for Nanopore sequencing and would be recommended for bacterial metagenomics studies using this technology.

BACKGROUND: Cold-adapted archaea have diverse ecological roles in a wide range of low-temperature environments. Improving our knowledge of the genomic features that enable psychrophiles to grow in cold environments helps us to understand their adaptive responses. However, samples from typical cold regions such as the remote Arctic and Antarctic are rare, and the limited number of high-quality genomes available leaves us with little data on genomic traits that are statistically associated with cold environmental conditions. RESULTS: In this study, we examined the haloarchaeal genus Halorubrum and defined a new clade that represents six isolates from polar and deep earth environments ('PD group' hereafter). The genomic G + C content and amino acid composition of this group distinguishes it from other Halorubrum and the trends are consistent with the established genomic optimization of psychrophiles. The cold adaptation of the PD group was further supported by observations of increased flexibility of proteins encoded across the genome and the findings of a growth test. CONCLUSIONS: The PD group Halorubrum exhibited denser genome packing, which confers higher metabolic potential with constant genome size, relative to the reference group, resulting in significant differences in carbon, nitrogen and sulfur metabolic patterns. The most marked feature was the enrichment of genes involved in sulfur cycling, especially the production of sulfite from organic sulfur-containing compounds. Our study provides an updated view of the genomic traits and metabolic potential of Halorubrum and expands the range of sources of cold-adapted haloarchaea.

BACKGROUND: Along with global warming, resulting in crop production, exacerbating the global food crisis. Therefore, it is urgent to study the mechanism of plant heat resistance. However, crop resistance genes were lost due to long-term artificial domestication. By analyzing the potential heat tolerance genes and molecular mechanisms in other wild materials, more genetic resources can be provided for improving the heat tolerance of crops. Elephant grass (Pennisetum purpureum Schum.) has strong adaptability to heat stress and contains abundant heat-resistant gene resources. RESULTS: Through sequence structure analysis, a total of 36 RWP-RK members were identified in elephant grass. Functional analysis revealed their close association with heat stress. Four randomly selected RKDs (RKD1.1, RKD4.3, RKD6.6, and RKD8.1) were analyzed for expression, and the results showed upregulation under high temperature conditions, suggesting their active role in response to heat stress. The members of RWP-RK gene family (36 genes) in elephant grass were 2.4 times higher than that of related tropical crops, rice (15 genes) and sorghum (15 genes). The 36 RWPs of elephant grass contain 15 NLPs and 21 RKDs, and 73% of RWPs are related to WGD. Among them, combined with the DAP-seq results, it was found that RWP-RK gene family expansion could improve the heat adaptability of elephant grass by enhancing nitrogen use efficiency and peroxidase gene expression. CONCLUSIONS: RWP-RK gene family expansion in elephant grass is closely related to thermal adaptation evolution and speciation. The RKD subgroup showed a higher responsiveness than the NLP subgroup when exposed to high temperature stress. The promoter region of the RKD subgroup contains a significant number of MeJA and ABA responsive elements, which may contribute to their positive response to heat stress. These results provided a scientific basis for analyzing the heat adaptation mechanism of elephant grass and improving the heat tolerance of other crops.

BACKGROUND: Structural descriptions of complete genomes have elucidated evolutionary processes in angiosperms. In Cactaceae (Caryophyllales), a high structural diversity of the chloroplast genome has been identified within and among genera. In this study, we assembled the first mitochondrial genome (mtDNA) for the short-globose cactus Mammillaria huitzilopochtli. For comparative purposes, we used the published genomes of 19 different angiosperms and the gymnosperm Cycas taitungensis as an external group for phylogenetic issues. RESULTS: The mtDNA of M. huitzilopochtli was assembled into one linear chromosome of 2,052,004 bp, in which 65 genes were annotated. These genes account for 57,606 bp including 34 protein-coding genes (PCGs), 27 tRNAs, and three rRNAs. In the non-coding sequences, repeats were abundant, with a total of 4,550 (179,215 bp). In addition, five complete genes (psaC and four tRNAs) of chloroplast origin were documented. Negative selection was estimated for most (23) of the PCGs. The phylogenetic tree showed a topology consistent with previous analyses based on the chloroplast genome. CONCLUSIONS: The number and type of genes contained in the mtDNA of M. huitzilopochtli were similar to those reported in 19 other angiosperm species, regardless of their phylogenetic relationships. Although other Caryophyllids exhibit strong differences in structural arrangement and total size of mtDNA, these differences do not result in an increase in the typical number and types of genes found in M. huitzilopochtli. We concluded that the total size of mtDNA in angiosperms increases by the lengthening of the non-coding sequences rather than a significant gain of coding genes.

BACKGROUND: Helcococcus ovis (H. ovis) is an emerging bacterial pathogen that commonly causes opportunistic respiratory, mammary, and uterine infections across mammalian hosts. This study applied long- and short-read whole genome sequencing technologies to identify virulence factors in five H. ovis isolates with low, medium, and high virulence phenotypes. RESULTS: The resulting assemblies contained one circular chromosome ranging from 1,744,566 to 1,850,083 bp in length and had a mean GC content of 27.6%. Phylogenetic and nucleotide identity analyses found low virulence strain KG38 to be part of a clade that forms an outgroup apart from the rest of the H. ovis taxon. Assembling the first complete genomes of the species revealed major genomic rearrangements in KG38. One to six prophage regions were identified in each genome. A novel pathogenicity island was found exclusively in the two high virulence strains (KG37 and KG104), along with two hypothetical transmembrane proteins designated as putative VFs. Finally, three zinc ABC transporters and three Type-II/IV secretion systems were identified as possible virulence determinants in this species. The low virulence strain KG38 has fewer intact paralogs of these operons in its genome compared to the higher virulence isolates, which strongly suggests a role in virulence. This strain is also missing four putative virulence factors (VFs) found in other isolates associated with adherence (collagen adhesin precursor), immune evasion (choline-binding protein A and a PspA-like hypothetical protein) and cell wall synthesis (glycerol-3-phosphate cytidylyltransferase). CONCLUSIONS: In this study, we assembled reference-quality complete genomes for five H. ovis strains to identify putative virulence factors. Phylogenetic analyses of H. ovis isolates revealed the presence of a clade representing a potentially novel species within the genus Helcococcus. A novel pathogenicity island and two hypothetical transmembrane proteins were found exclusively in high-virulence strains. The identification of Zinc ABC transporters and Type-II/IV secretion systems as possible virulence determinants, along with the differences in operon content between the low and high virulence isolates, strongly suggests they also play a role in the bacterium's pathogenicity. Taken together, these findings are a valuable first step toward deciphering the pathogenesis of H. ovis infections.

BACKGROUND: Blueberries (Vaccinium corymbosum) are regarded as "superfoods" attributed to large amounts of anthocyanins, a group of flavonoid metabolites, which provide pigmentation in plant and beneficial effects for human health. MYB transcription factor is one of vital components in the regulation of plant secondary metabolism, which occupies a dominant position in the regulatory network of anthocyanin biosynthesis. However, the role of MYB family in blueberry responding to anthocyanin biosynthesis remains elusive. RESULTS: In this study, we conducted a comprehensive analysis of VcMYBs in blueberry based on the genome data, including phylogenetic relationship, conserved motifs, identification of differentially expressed MYB genes during fruit development and their expression profiling, etc. A total of 437 unique MYB sequences with two SANT domains were identified in blueberry, which were divided into 3 phylogenetic trees. Noticeably, there are many trigenic and tetragenic VcMYBs pairs with more than 95% identity to each other. Meanwhile, the transcript accumulations of VcMYBs were surveyed underlying blueberry fruit development, and they showed diverse expression patterns, suggesting various functional roles in fruit ripening. More importantly, distinct transcript profiles between skin and pulp of ripe fruit were observed for several VcMYBs, such as VcMYB437, implying the potential roles in anthocyanin biosynthesis. CONCLUSIONS: Totally, 437 VcMYBs were identified and characterized. Subsequently, their transcriptional patterns were explored during fruit development and fruit tissues (skin and pulp) closely related to anthocyanin biosynthesis. These genome-wide data and findings will contribute to demonstrating the functional roles of VcMYBs and their regulatory mechanisms for anthocyanins production and accumulation in blueberry in the future study.

BACKGROUND: As an important reproductive hormone, melatonin plays an important role in regulating the reproductive activities of sheep and other mammals. Hu sheep is a breed favoring for meat, with prolific traits. In order to explore the relationship between melatonin and reproductive function of Hu sheep, 7,694,759 SNPs were screened out through the whole genome sequencing analysis from high and low melatonin production Hu sheep. RESULTS: A total of 68,673 SNPs, involving in 1126 genes, were identified by ED association analysis. Correlation analysis of SNPs of AANAT/ASMT gene and MTNR1A/MTNR1B gene were carried out. The melatonin level of CG genotype 7,981,372 of AANAT, GA genotype 7,981,866 of ASMT and GG genotype 17,355,171 of MTNR1A were higher than the average melatonin level of 1.64 ng/mL. High melatonin Hu sheep appear to have better multiple reproductive performance. CONCLUSIONS: By using different methods, three SNPs which are associated with high melatonin production trait have been identified in Hu sheep. These 3 SNPs are located in melatonin synthetase AANAT/ASMT and receptor MTNR1A, respectively. Considering the positive association between melatonin production and reproductive performance in ruminants, these three SNPs can be served as the potential molecular markers for breading Hu sheep with the desirable reproductive traits.

BACKGROUND: The emergence and wide spread of carbapenemase-producing Enterobacteriaceae (CPE) poses a growing threat to global public health. However, clinically derived carbapenemase-producing Citrobacter causing multiple infections has rarely been investigated. Here we first report the isolation and comparative genomics of two blaNDM-5 carrying Citrobacter freundii (C. freundii) isolates from a patient with bloodstream and urinary tract infections. RESULTS: Antimicrobial susceptibility testing showed that both blaNDM-5 carrying C. freundii isolates were multidrug-resistant. Positive modified carbapenem inactivation method (mCIM) and EDTA-carbapenem inactivation method (eCIM) results suggested metallo-carbapenemase production. PCR and sequencing confirmed that both metallo-carbapenemase producers were blaNDM-5 positive. Genotyping and comparative genomics analyses revealed that both isolates exhibited a high level of genetic similarity. Plasmid analysis confirmed that the blaNDM-5 resistance gene is located on IncX3 plasmid with a length of 46,161 bp, and could successfully be transferred to the recipient Escherichia coli EC600 strain. A conserved structure sequence (ISAba125-IS5-blaNDM-5-trpF-IS26-umuD-ISKox3) was found in the upstream and downstream of the blaNDM-5 gene. CONCLUSIONS: The data presented in this study showed that the conjugative blaNDM-5 plasmid possesses a certain ability to horizontal transfer. The dissemination of NDM-5-producing C. freundii isolates should be of close concern in future clinical surveillance. To our knowledge, this is the first study to characterize C. freundii strains carrying the blaNDM-5 gene from one single patient with multiple infections.

BACKGROUND: Laccases are green biocatalysts with wide industrial applications. The study of efficient and specific laccase producers remains a priority. Cerrena species have been shown to be promising basidiomycete candidates for laccase production. Although two sets of Cerrena genome data have been publicly published, no comprehensive bioinformatics study of laccase gene family in C. unicolor has been reported, particularly concerning the analysis of their three-dimensional (3D) structures and molecular docking to substrates, like ABTS and aflatoxin B1 (AFB1). RESULTS: In this study, we conducted a comprehensive genome-wide analysis of laccase gene family in C. unicolor 87613. We identified eighteen laccase genes (CuLacs) and classified them into three clades using phylogenetic analysis. We characterized these laccases, including their location in contig 5,6,9,12,15,19,26,27, gene structures of different exon-intron arrangements, molecular weight ranging from 47.89 to 141.41 kDa, acidic pI value, 5-15 conserved protein motifs, signaling peptide of extracellular secretion (harbored by 13 CuLacs) and others. In addition, the analysis of cis-acting element in laccase promoters indicated that the transcription response of CuLac gene family was regulatable and complex under different environmental cues. Furthermore, analysis of transcription pattern revealed that CuLac8, 12 and CuLac2, 13 were the predominant laccases in response to copper ions or oxidative stress, respectively. Finally, we focused on the 3D structure analysis of CuLac proteins. Seven laccases with extra transmembrane domains or special sequences were particularly interesting. Predicted structures of each CuLac protein with or without these extra sequences showed altered interacting amino acid residues and binding sites, leading to varied affinities to both ABTS and AFB1. As far as we know, it is the first time to discuss the influence of the extra sequence on laccase's affinity to substrates. CONCLUSIONS: Our findings provide robust genetic data for a better understanding of the laccase gene family in C. unicolor 87613, and create a foundation for the molecular redesign of CuLac proteins to enhance their industrial applications.

BACKGROUND: The Mongolian gazelle (Procapra gutturosa) population has shown a considerable range of contractions and local extinctions over the last century, owing to habitat fragmentation and poaching. A thorough understanding of the genetic diversity and structure of Mongolian gazelle populations in fragmented habitats is critical for planning effective conservation strategies. RESULT: In this study, we used eight microsatellite loci and mitochondrial cytochrome b (Cytb) to compare the levels of genetic diversity and genetic structure of Mongolian gazelle populations in the Hulun Lake National Nature Reserve (HLH) with those in the China-Mongolia border area (BJ). The results showed that the nucleotide diversity and observed heterozygosity of the HLH population were lower than those of the BJ population. Moreover, the HLH and BJ populations showed genetic differentiation. We concluded that the HLH population had lower genetic diversity and a distinct genetic structure compared with the BJ population. CONCLUSION: The genetic diversity of fragmented Mongolian gazelle populations, can be improved by protecting these populations while reinforcing their gene exchange with other populations. For example, attempts can be made to introduce new individuals with higher genetic diversity from other populations to reduce inbreeding.

BACKGROUND: Cadmium (Cd) flows into the ocean with industrial and agricultural pollution and significantly affects the growth and development of economic cephalopods such as Sepia esculenta, Amphioctopus fangsiao, and Loligo japonica. As of now, the reasons why Cd affects the growth and development of S. esculenta are not yet clear. RESULTS: In this study, transcriptome and four oxidation and toxicity indicators are used to analyze the toxicological mechanism of Cd-exposed S. esculenta larvae. Indicator results indicate that Cd induces oxidative stress and metal toxicity. Functional enrichment analysis results suggest that larval ion transport, cell adhesion, and some digestion and absorption processes are inhibited, and the cell function is damaged. Comprehensive analysis of protein-protein interaction network and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was used to explore S. esculenta larval toxicological mechanisms, and we find that among the 20 identified key genes, 14 genes are associated with neurotoxicity. Most of them are down-regulated and enriched to the neuroactive ligand-receptor interaction signaling pathway, suggesting that larval nervous system might be destroyed, and the growth, development, and movement process are significantly affected after Cd exposure. CONCLUSIONS: S. esculenta larvae suffered severe oxidative damage after Cd exposure, which may inhibit digestion and absorption functions, and disrupt the stability of the nervous system. Our results lay a function for understanding larval toxicological mechanisms exposed to heavy metals, promoting the development of invertebrate environmental toxicology, and providing theoretical support for S. esculenta artificial culture.

The avian eggshell is formed in the uterus. Changes in uterine function may have a significant effect on eggshell quality. To identify the vital genes impacting uterine functional maintenance in the chicken, uteri in three different periods (22W, 31W, 51W) were selected for RNA sequencing and bioinformatics analysis. In our study, 520, 706 and 736 differentially expressed genes (DEGs) were respectively detected in the W31 vs W22 group, W51 vs W31 group and W51 vs W22 group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated DEGs were enriched in the extracellular matrix, extracellular region part, extracellular region, extracellular matrix structural constituent, ECM receptor interaction, collagen-containing extracellular matrix and collagen trimer in the uterus (P < 0.05). Protein-protein interaction analysis revealed that FN1, LOX, THBS2, COL1A1, COL1A2, COL5A1, COL5A2, POSTN, MMP13, VANGL2, RAD54B, SPP1, SDC1, BTC, ANGPTL3 might be key candidate genes for uterine functional maintenance in chicken. This study discovered dominant genes and pathways which enhanced our knowledge of chicken uterine functional maintenance.

BACKGROUND: COVID-19 waves caused by specific SARS-CoV-2 variants have occurred globally at different times. We focused on Omicron variants to understand the genomic diversity and phylogenetic relatedness of SARS-CoV-2 strains in various regions of Pakistan. METHODS: We studied 276,525 COVID-19 cases and 1,031 genomes sequenced from December 2021 to August 2022. Sequences were analyzed and visualized using phylogenetic trees. RESULTS: The highest case numbers and deaths were recorded in Sindh and Punjab, the most populous provinces in Pakistan. Omicron variants comprised 93% of all genomes, with BA.2 (32.6%) and BA.5 (38.4%) predominating. The first Omicron wave was associated with the sequential identification of BA.1 in Sindh, then Islamabad Capital Territory, Punjab, Khyber Pakhtunkhwa (KP), Azad Jammu Kashmir (AJK), Gilgit-Baltistan (GB) and Balochistan. Phylogenetic analysis revealed Sindh to be the source of BA.1 and BA.2 introductions into Punjab and Balochistan during early 2022. BA.4 was first introduced in AJK and BA.5 in Punjab. Most recent common ancestor (MRCA) analysis revealed relatedness between the earliest BA.1 genome from Sindh with Balochistan, AJK, Punjab and ICT, and that of first BA.1 from Punjab with strains from KPK and GB. CONCLUSIONS: Phylogenetic analysis provides insights into the introduction and transmission dynamics of the Omicron variant in Pakistan, identifying Sindh as a hotspot for viral dissemination. Such data linked with public health efforts can help limit surges of new infections.

BACKGROUND: Pectate lyase (PL, EC 4.2.2.2), as an endo-acting depolymerizing enzyme, cleaves α-1,4-glycosidic linkages in esterified pectin and involves a broad range of cell wall modifications. However, the knowledge concerning the genome-wide analysis of the PL gene family in Fragaria vesca has not been thoroughly elucidated. RESULTS: In this study, sixteen PLs members in F. vesca were identified based on a genome-wide investigation. Substantial divergences existed among FvePLs in gene duplication, cis-acting elements, and tissue expression patterns. Four clusters were classified according to phylogenetic analysis. FvePL6, 8 and 13 in cluster II significantly contributed to the significant expansions during evolution by comparing orthologous PL genes from Malus domestica, Solanum lycopersicum, Arabidopsis thaliana, and Fragaria×ananassa. The cis-acting elements implicated in the abscisic acid signaling pathway were abundant in the regions of FvePLs promoters. The RNA-seq data and in situ hybridization revealed that FvePL1, 4, and 7 exhibited maximum expression in fruits at twenty days after pollination, whereas FvePL8 and FvePL13 were preferentially and prominently expressed in mature anthers and pollens. Additionally, the co-expression networks displayed that FvePLs had tight correlations with transcription factors and genes implicated in plant development, abiotic/biotic stresses, ions/Ca2+, and hormones, suggesting the potential roles of FvePLs during strawberry development. Besides, histological observations suggested that FvePL1, 4 and 7 enhanced cell division and expansion of the cortex, thus negatively influencing fruit firmness. Finally, FvePL1-RNAi reduced leaf size, altered petal architectures, disrupted normal pollen development, and rendered partial male sterility. CONCLUSION: These results provide valuable information for characterizing the evolution, expansion, expression patterns and functional analysis, which help to understand the molecular mechanisms of the FvePLs in the development of strawberries.

BACKGROUND: CD-1 is an outbred mouse stock that is frequently used in toxicology, pharmacology, and fundamental biomedical research. Although inbred strains are typically better suited for such studies due to minimal genetic variability, outbred stocks confer practical advantages over inbred strains, such as improved breeding performance and low cost. Knowledge of the full genetic variability of CD-1 would make it more useful in toxicology, pharmacology, and fundamental biomedical research. RESULTS: We performed deep genomic DNA sequencing of CD-1 mice and used the data to identify genome-wide SNPs, indels, and germline transposable elements relative to the mm10 reference genome. We used multiple genome-wide sequencing data types and previously published CD-1 SNPs to validate our called variants. We used the called variants to construct a strain-specific CD-1 reference genome, which we show can improve mappability and reduce experimental biases from genome-wide sequencing data derived from CD-1 mice. Based on previously published ChIP-seq and ATAC-seq data, we find evidence that genetic variation between CD-1 mice can lead to alterations in transcription factor binding. We also identified a number of variants in the coding region of genes which could have effects on translation of genes. CONCLUSIONS: We have identified millions of previously unidentified CD-1 variants with the potential to confound studies involving CD-1. We used the identified variants to construct a CD-1-specific reference genome, which can improve accuracy and reduce bias when aligning genomics data derived from CD-1 mice.

BACKGROUND: Arthropods are the largest group in the animal kingdom and are morphologically characterized by heterorhythmic segments. Brachyuran decapod crustaceans undergo brachyurization metamorphosis in the early developmental process, characterized by a reduced abdomen that is folded beneath the cephalothorax and inserted between the pereiopods or in a special cavity. As the main cause of major alterations in the evolution of animal body plans, Hox genes encode transcription factors and are involved in bilaterian anterior-posterior axis patterning. RESULTS: We found eight Hox genes (labial, proboscipedia, Deformed, zerknüllt, Sex combs reduced, Antennapedia, Ultrabithorax, fushi tarazu, abdominal-A and Abdominal-B) in Eriocheir sinensis. The phylogenetic topology of 13 arthropod Hox genes was closely related to traditional taxonomic groupings. Genome collinearity analysis was performed using genomic data and chromosomal location data of E. sinensis and Portunus trituratus. We found that their chromosomes were highly collinear, and there was a corresponding collinear relationship between the three Hox genes (lab, ftz and Abd-B). The mRNA expression levels of Scr and Antp fluctuated significantly in different developmental stages of E. sinensis, especially in the brachyurization stages. Evolutionary analysis indicated the presence of positively selected sites in Ubx. CONCLUSIONS: In this study, we used genome-wide analysis to identify and analyze all members of the Hox genes in E. sinensis. Our data will contribute to a better understanding of Hox genes in E. sinensis and provide useful molecular evolutionary information for further investigation on their roles in the brachyurization of crabs.

BACKGROUND: Fatty acids composition in poultry muscle is directly related to its tenderness, flavour, and juiciness, whereas its genetic mechanisms have not been elucidated. In this study, the genetic structure and key regulatory genes of the breast muscle fatty acid composition of local Chinese chicken, Gushi-Anka F2 resource population by integrating genome-wide association study (GWAS) and weighted gene co-expression network analysis (WGCNA) strategies. GWAS was performed based on 323,306 single nucleotide polymorphisms (SNPs) obtained by genotyping by sequencing (GBS) method and 721 chickens from the Gushi-Anka F2 resource population with highly variable fatty acid composition traits in the breast muscle. And then, according to the transcriptome data of the candidate genes that were obtained and phenotypic data of fatty acid composition traits in breast muscle of Gushi chickens at 14, 22, and 30 weeks of age, we conducted a WGCNA. RESULTS: A total of 128 suggestive significantly associated SNPs for 11 fatty acid composition traits were identified and mapped on chromosomes (Chr) 2, 3, 4, 5, 13, 17, 21, and 27. Of these, the two most significant SNPs were Chr13:5,100,140 (P = 4.56423e-10) and Chr13:5,100,173 (P = 4.56423e-10), which explained 5.6% of the phenotypic variation in polyunsaturated fatty acids (PUFA). In addition, six fatty acid composition traits, including C20:1, C22:6, saturated fatty acid (SFA), unsaturated fatty acids (UFA), PUFA, and average chain length (ACL), were located in the same QTL intervals on Chr13. We obtained 505 genes by scanning the linkage disequilibrium (LD) regions of all significant SNPs and performed a WGCNA based on the transcriptome data of the above 505 genes. Combining two strategies, 9 hub genes (ENO1, ADH1, ASAH1, ADH1C, PIK3CD, WISP1, AKT1, PANK3, and C1QTNF2) were finally identified, which could be the potential candidate genes regulating fatty acid composition traits in chicken breast muscle. CONCLUSION: The results of this study deepen our understanding of the genetic mechanisms underlying the regulation of fatty acid composition traits, which is helpful in the design of breeding strategies for the subsequent improvement of fatty acid composition in poultry muscle.

BACKGROUND: Plant organs grow in a certain direction and organ twisted growth, a rare and distinctive trait, is associated with internal structure changes and special genes. The twisted branch mutant of Chinese jujube jujube, an important fruit tree native to China and introduced to nearly 50 countries, provides new typical materials for exploration of plant twisted growth. RESULTS: In this study, the cytological characteristics and related genes of twisted branches in Chinese jujube were revealed by microscopy observation and transcriptome analysis. The unique coexistence of primary and secondary structures appeared in the twisted parts of branches, and special structures such as collateral bundle, cortical bundles, and internal phloem were formed. Ninety differentially expressed genes of 'Dongzao' and its twisted mutant were observed, in which ZjTBL43, ZjFLA11, ZjFLA12 and ZjIQD1 were selected as candidate genes. ZjTBL43 was homologous to AtTBL43 in Arabidopsis, which was involved in the synthesis and deposition of cellular secondary wall cellulose. The attbl43 mutant showed significant inflorescence stem bending growth. The transgenic lines of attbl43 with overexpression of ZjTBL43 were phenotypically normal.The branch twisted growth may be caused by mutations in ZjTBL43 in Chinese jujube. AtIQD10, AtFLA11 and AtFLA12 were homologous to ZjIQD1, ZjFLA11 and ZjFLA12, respectively. However, the phenotype of their function defect mutants was normal. CONCLUSION: In summary, these findings will provide new insights into the plant organ twisted growth and a reference for investigation of controlling mechanisms of plant growth direction.

BACKGROUND: Growing evidence has shown that gut microbiome composition is associated with breast cancer (BC), but the causality remains unknown. We aimed to investigate the link between BC prognosis and the gut microbiome at various oestrogen receptor (ER) statuses. METHODS: We performed a genome-wide association study (GWAS) to analyse the gut microbiome of BC patients, the dataset for which was collected by the Breast Cancer Association Consortium (BCAC). The analysis was executed mainly via inverse variance weighting (IVW); the Mendelian randomization (MR) results were verified by heterogeneity tests, sensitivity analysis, and pleiotropy analysis. RESULTS: Our findings identified nine causal relationships between the gut microbiome and total BC cases, with ten and nine causal relationships between the gut microbiome and ER-negative (ER-) and ER-positive (ER+) BC, respectively. The family Ruminococcaceae and genus Parabacteroides were most apparent among the three categories. Moreover, the genus Desulfovibrio was expressed in ER- BC and total BC, whereas the genera Sellimonas, Adlercreutzia and Rikenellaceae appeared in the relationship between ER + BC and total BC. CONCLUSION: Our MR inquiry confirmed that the gut microbiota is causally related to BC. This further explains the link between specific bacteria for prognosis of BC at different ER statuses. Considering that potential weak instrument bias impacts the findings and that the results are limited to European females due to data constraints, further validation is crucial.

Despite recent efforts to increase diversity in genome-wide association studies (GWASs), most loci currently associated with kidney function are still limited to European ancestry due to the underlying sample selection bias in available GWASs. We set out to identify susceptibility loci associated with estimated glomerular filtration rate (eGFRcrea) in 80027 individuals of African-ancestry from the UK Biobank (UKBB), Million Veteran Program (MVP), and Chronic Kidney Disease genetics (CKDGen) consortia.We identified 8 lead SNPs, 7 of which were previously associated with eGFR in other populations. We identified one novel variant, rs77408001 which is an intronic variant mapped to the ELN gene. We validated three previously reported loci at GATM-SPATA5L1, SLC15A5 and AGPAT3. Fine-mapping analysis identified variants rs77121243 and rs201602445 as having a 99.9% posterior probability of being causal. Our results warrant designing bigger studies within individuals of African ancestry to gain new insights into the pathogenesis of Chronic Kidney Disease (CKD), and identify genomic variants unique to this ancestry that may influence renal function and disease.

This paper aims to explore the role of circRNA expression profiles and circRNA-associated ceRNA networks in the regulation of myogenesis in the longissimus dorsi of cattle breeds surviving under subtropical conditions in southern China by RNA sequencing and bioinformatics analysis. It also aims to provide comprehensive understanding of the differences in muscle fibers in subtropical cattle breeds and to expand the knowledge of the molecular networks that regulate myogenesis. With regard to meat quality indicators, results showed that the longissimus dorsi of LQC had lower pH (P < 0.0001), lower redness (P < 0.01), lower shear force (P < 0.05), and higher brightness (P < 0.05) than the longissimus dorsi of LFC. With regard to muscle fiber characteristics, the longissimus dorsi of LQC had a smaller diameter (P < 0.0001) and higher density of muscle fibers (P < 0.05). The analysis results show that the function of many circRNA-targeted mRNAs was related to myogenesis and metabolic regulation. Furthermore, in the analysis of the function of circRNA source genes, we hypothesized that btacirc_00497 and btacirc_034497 may regulate the function and type of myofibrils by affecting the expression of MYH6, MYH7, and NEB through competitive linear splicing.

BACKGROUND: Alfalfa (Medicago sativa) is the most widely planted legume forage and one of the most economically valuable crops in the world. The periodic changes in its growth and development and abiotic stress determine its yield and economic benefits. Auxin controls many aspects of alfalfa growth by regulating gene expression, including organ differentiation and stress response. Auxin response factors (ARF) are transcription factors that play an essential role in auxin signal transduction and regulate the expression of auxin-responsive genes. However, the function of ARF transcription factors is unclear in autotetraploid-cultivated alfalfa. RESULT: A total of 81 ARF were identified in the alfalfa genome in this study. Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed, identifying that ARF genes are mainly involved in transcriptional regulation and plant hormone signal transduction pathways. Phylogenetic analysis showed that MsARF was divided into four clades: I, II, III, and IV, each containing 52, 13, 7, and 9 genes, respectively. The promoter region of the MsARF gene contained stress-related elements, such as ABRE, TC-rich repeats, MBS, LTR. Proteins encoded by 50 ARF genes were localized in the nucleus without guide peptides, signal peptides, or transmembrane structures, indicating that most MsARF genes are not secreted or transported but only function in the nucleus. Protein structure analysis revealed that the secondary and tertiary structures of the 81 MsARF genes varied. Chromosomal localization analysis showed 81 MsARF genes were unevenly distributed on 25 chromosomes, with the highest distribution on chromosome 5. Furthermore, 14 segmental duplications and two sets of tandem repeats were identified. Expression analysis indicated that the MsARF was differentially expressed in different tissues and under various abiotic stressors. The quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that the expression profiles of 23 MsARF genes were specific to abiotic stresses such as drought, salt, high temperature, and low temperature, as well as tissue-specific and closely related to the duration of stress. CONCLUSION: This study identified MsARF in the cultivated alfalfa genome based on the autotetraploid level, which GO, KEGG analysis, phylogenetic analysis, sequence characteristics, and expression pattern analysis further confirmed. Together, these findings provide clues for further investigation of MsARF functional verification and molecular breeding of alfalfa. This study provides a novel approach to systematically identify and characterize ARF transcription factors in autotetraploid cultivated alfalfa, revealing 23 MsARF genes significantly involved in response to various stresses.

BACKGROUND: Chromatin-associated phase separation proteins establish various biomolecular condensates via liquid-liquid phase separation (LLPS), which regulates vital biological processes spatially and temporally. However, the widely used methods to characterize phase separation proteins are still based on low-throughput experiments, which consume time and could not be used to explore protein LLPS properties in bulk. RESULTS: By combining gradient 1,6-hexanediol (1,6-HD) elution and quantitative proteomics, we developed chromatin enriching hexanediol separation coupled with liquid chromatography-mass spectrometry (CHS-MS) to explore the LLPS properties of different chromatin-associated proteins (CAPs). First, we found that CAPs were enriched more effectively in the 1,6-HD treatment group than in the isotonic solution treatment group. Further analysis showed that the 1,6-HD treatment group could effectively enrich CAPs prone to LLPS. Finally, we compared the representative proteins eluted by different gradients of 1,6-HD and found that the representative proteins of the 2% 1,6-HD treatment group had the highest percentage of IDRs and LCDs, whereas the 10% 1,6-HD treatment group had the opposite trend. CONCLUSION: This study provides a convenient high-throughput experimental method called CHS-MS. This method can efficiently enrich proteins prone to LLPS and can be extended to explore LLPS properties of CAPs in different biological systems.

BACKGROUND: Immune traits are considered to serve as potential biomarkers for pig's health. Medium to high heritabilities have been observed for some of the immune traits suggesting genetic variability of these phenotypes. Consideration of previously established genetic correlations between immune traits can be used to identify pleiotropic genetic markers. Therefore, genome-wide association study (GWAS) approaches are required to explore the joint genetic foundation for health biomarkers. Usually, GWAS explores phenotypes in a univariate (uv), trait-by-trait manner. Besides two uv GWAS methods, four multivariate (mv) GWAS approaches were applied on combinations out of 22 immune traits for Landrace (LR) and Large White (LW) pig lines. RESULTS: In total 433 (LR: 351, LW: 82) associations were identified with the uv approach implemented in PLINK and a Bayesian linear regression uv approach (BIMBAM) software. Single Nucleotide Polymorphisms (SNPs) that were identified with both uv approaches (n = 32) were mostly associated with immune traits such as haptoglobin, red blood cell characteristics and cytokines, and were located in protein-coding genes. Mv GWAS approaches detected 647 associations for different mv immune trait combinations which were summarized to 133 Quantitative Trait Loci (QTL). SNPs for different trait combinations (n = 66) were detected with more than one mv method. Most of these SNPs are associated with red blood cell related immune trait combinations. Functional annotation of these QTL revealed 453 immune-relevant protein-coding genes. With uv methods shared markers were not observed between the breeds, whereas mv approaches were able to detect two conjoint SNPs for LR and LW. Due to unmapped positions for these markers, their functional annotation was not clarified. CONCLUSIONS: This study evaluated the joint genetic background of immune traits in LR and LW piglets through the application of various uv and mv GWAS approaches. In comparison to uv methods, mv methodologies identified more significant associations, which might reflect the pleiotropic background of the immune system more accurately. In genetic research of complex traits, the SNP effects are generally small. Furthermore, one genetic variant can affect several correlated immune traits at the same time, termed pleiotropy. As mv GWAS methods consider strong dependencies among traits, the power to detect SNPs can be boosted. Both methods revealed immune-relevant potential candidate genes. Our results indicate that one single test is not able to detect all the different types of genetic effects in the most powerful manner and therefore, the methods should be applied complementary.

BACKGROUND: DNA methylation plays an important role in life processes by affecting gene expression, but it is still unclear how DNA methylation is controlled and how it regulates gene transcription under high temperature stress conditions in Strongylocentrotus intermedius. The potential link between DNA methylation variation and gene expression changes in response to heat stress in S. intermedius was investigated by MethylRAD-seq and RNA-seq analysis. We screened DNA methylation driver genes in order to comprehensively elucidate the regulatory mechanism of its high temperature adaptation at the DNA/RNA level. RESULTS: The results revealed that high temperature stress significantly affected not only the DNA methylation and transcriptome levels of S. intermedius (P < 0.05), but also growth. MethylRAD-seq analysis revealed 12,129 CG differential methylation sites and 966 CWG differential methylation sites, and identified a total of 189 differentially CG methylated genes and 148 differentially CWG methylated genes. Based on KEGG enrichment analysis, differentially expressed genes (DEGs) are mostly enriched in energy and cell division, immune, and neurological damage pathways. Further RNA-seq analysis identified a total of 1968 DEGs, of which 813 genes were upregulated and 1155 genes were downregulated. Based on the joint MethylRAD-seq and RNA-seq analysis, metabolic processes such as glycosaminoglycan degradation, oxidative phosphorylation, apoptosis, glutathione metabolism, thermogenesis, and lysosomes are regulated by DNA methylation. CONCLUSIONS: High temperature affected the DNA methylation and expression levels of genes such as MOAP-1, GGT1 and RDH8, which in turn affects the metabolism of HPSE, Cox, glutathione, and retinol, thereby suppressing the immune, energy metabolism, and antioxidant functions of the organism and finally manifesting as stunted growth. In summary, the observations in the present study improve our understanding of the molecular mechanism of the response to high temperature stress in sea urchin.

BACKGROUND: Soybean is one of the most important oil crops in the world. The domestication of wild soybean has resulted in significant changes in the seed oil content and seed size of cultivated soybeans. To better understand the molecular mechanisms of seed formation and oil content accumulation, WDD01514 (E1), ZYD00463 (E2), and two extreme progenies (E23 and E171) derived from RILs were used for weighted gene coexpression network analysis (WGCNA) combined with transcriptome analysis. RESULTS: In this study, both seed weight and oil content in E1 and E171 were significantly higher than those in E2 and E23, and 20 DAF and 30 DAF may be key stages of soybean seed oil content accumulation and weight increase. Pathways such as "Photosynthesis", "Carbon metabolism", and "Fatty acid metabolism", were involved in oil content accumulation and grain formation between wild and cultivated soybeans at 20 and 30 DAF according to RNA-seq analysis. A total of 121 oil content accumulation and 189 seed formation candidate genes were screened from differentially expressed genes. WGCNA identified six modules related to seed oil content and seed weight, and 76 candidate genes were screened from modules and network. Among them, 16 genes were used for qRT-PCR and tissue specific expression pattern analysis, and their expression-levels in 33-wild and 23-cultivated soybean varieties were subjected to correlation analysis; some key genes were verified as likely to be involved in oil content accumulation and grain formation. CONCLUSIONS: Overall, these results contribute to an understanding of seed lipid metabolism and seed size during seed development, and identify potential functional genes for improving soybean yield and seed oil quantity.

Peanut (Arachis hypogaea L.) is an important oilseed crop worldwide. Improving its yield is crucial for sustainable peanut production to meet increasing food and industrial requirements. Deciphering the genetic control underlying peanut kernel weight and size, which are essential components of peanut yield, would facilitate high-yield breeding. A high-density single nucleotide polymorphism (SNP)-based linkage map was constructed using a recombinant inbred lines (RIL) population derived from a cross between the variety Yuanza9102 and a germplasm accession wt09-0023. Kernel weight and size quantitative trait loci (QTLs) were co-localized to a 0.16 Mb interval on Arahy07 using inclusive composite interval mapping (ICIM). Analysis of SNP, and Insertion or Deletion (INDEL) markers in the QTL interval revealed a gene encoding a pentatricopeptide repeat (PPR) superfamily protein as a candidate closely linked with kernel weight and size in cultivated peanut. Examination of the PPR gene family indicated a high degree of collinearity of PPR genes between A. hypogaea and its diploid progenitors, Arachis duranensis and Arachis ipaensis. The candidate PPR gene, Arahy.JX1V6X, displayed a constitutive expression pattern in developing seeds. These findings lay a foundation for further fine mapping of QTLs related to kernel weight and size, as well as validation of candidate genes in cultivated peanut.

BACKGROUND: Liver metastasis is the major challenge in the treatment for malignant tumors. Genomic profiling is increasingly used in the diagnosis, treatment and prediction of prognosis in malignancies. In this study, we constructed a gene mutation-based risk model to predict the survival of liver metastases. METHOD: We identified the gene mutations associated with survival and constructed the risk model in the training cohort including 800 patients with liver metastases from Memorial Sloan-Kettering Cancer Center (MSKCC) dataset. Other 794 patients with liver metastases were collected from 4 cohorts for validation. Furthermore, the analyses of tumor microenvironment (TME) and somatic mutations were performed on 51 patients with breast cancer liver metastases (BCLM) who had both somatic mutation data and RNA-sequencing data. RESULTS: A gene mutation-based risk model involved 10 genes was constructed to divide patients with liver metastases into the high- and low-risk groups. Patients in the low-risk group had a longer survival time compared to those in the high-risk group, which was observed in both training and validation cohorts. The analyses of TME in BCLM showed that the low-risk group exhibited more immune infiltration than the high-risk group. Furthermore, the mutation signatures of the high-risk group were completely different from those of the low-risk group in patients with BCLM. CONCLUSIONS: The gene mutation-based risk model constructed in our study exhibited the reliable ability of predicting the prognosis in liver metastases. The difference of TME and somatic mutations among BCLM patients with different risk score can guide the further research and treatment decisions for liver metastases.

BACKGROUND: As the characteristic functional component in ginger, gingerols possess several health-promoting properties. Long non-coding RNAs (lncRNAs) act as crucial regulators of diverse biological processes. However, lncRNAs in ginger are not yet identified so far, and their potential roles in gingerol biosynthesis are still unknown. In this study, metabolomic and transcriptomic analyses were performed in three main ginger cultivars (leshanhuangjiang, tonglingbaijiang, and yujiang 1 hao) in China to understand the potential roles of the specific lncRNAs in gingerol accumulation. RESULTS: A total of 744 metabolites were monitored by metabolomics analysis, which were divided into eleven categories. Among them, the largest group phenolic acid category contained 143 metabolites, including 21 gingerol derivatives. Of which, three gingerol analogs, [8]-shogaol, [10]-gingerol, and [12]-shogaol, accumulated significantly. Moreover, 16,346 lncRNAs, including 2,513, 1,225, and 2,884 differentially expressed (DE) lncRNA genes (DELs), were identified in all three comparisons by transcriptomic analysis. Gene ontology enrichment (GO) analysis showed that the DELs mainly enriched in the secondary metabolite biosynthetic process, response to plant hormones, and phenol-containing compound metabolic process. Correlation analysis revealed that the expression levels of 11 DE gingerol biosynthesis enzyme genes (GBEGs) and 190 transcription factor genes (TF genes), such as MYB1, ERF100, WRKY40, etc. were strongly correlation coefficient with the contents of the three gingerol analogs. Furthermore, 7 and 111 upstream cis-acting lncRNAs, 1,200 and 2,225 upstream trans-acting lncRNAs corresponding to the GBEGs and TF genes were identified, respectively. Interestingly, 1,184 DELs might function as common upstream regulators to these GBEGs and TFs genes, such as LNC_008452, LNC_006109, LNC_004340, etc. Furthermore, protein-protein interaction networks (PPI) analysis indicated that three TF proteins, MYB4, MYB43, and WRKY70 might interact with four GBEG proteins (PAL1, PAL2, PAL3, and 4CL-4). CONCLUSION: Based on these findings, we for the first time worldwide proposed a putative regulatory cascade of lncRNAs, TFs genes, and GBEGs involved in controlling of gingerol biosynthesis. These results not only provide novel insights into the lncRNAs involved in gingerol metabolism, but also lay a foundation for future in-depth studies of the related molecular mechanism.

WD40 transcription factors (TFs) constitute a large gene family in eukaryotes, playing diverse roles in cellular processes. However, their functions in the major ornamental plant, Rhododendron simsii, remain poorly understood. In this study, we identified 258 WD40 proteins in the R. simsii genome, which exhibited an uneven distribution across chromosomes. Based on domain compositions and phylogenetic analysis, we classified these 258 RsWD40 proteins into 42 subfamilies and 47 clusters. Comparative genomic analysis suggested that the expansion of the WD40 gene family predates the divergence of green algae and higher plants, indicating an ancient origin. Furthermore, by analyzing the duplication patterns of RsWD40 genes, we found that transposed duplication played a major role in their expansion. Notably, the majority of RsWD40 gene duplication pairs underwent purifying selection during evolution. Synteny analysis identified significant orthologous gene pairs between R. simsii and Arabidopsis thaliana, Oryza sativa, Vitis vinifera, and Malus domestica. We also investigated potential candidate genes involved in anthocyanin biosynthesis during different flower development stages in R. simsii using RNA-seq data. Specifically, we identified 10 candidate genes during the bud stage and 7 candidate genes during the full bloom stage. GO enrichment analysis of these candidate genes revealed the potential involvement of the ubiquitination process in anthocyanin biosynthesis. Overall, our findings provide a valuable foundation for further investigation and functional analysis of WD40 genes, as well as research on the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species.

BACKGROUND: Telomerase reverse transcriptase (TERT) plays a crucial role in synthesizing telomeric repeats that safeguard chromosomes from damage and fusion, thereby maintaining genome stability. Mutations in the TERT gene can lead to a deviation in gene expression, impaired enzyme activity, and, as a result, abnormal telomere shortening. Genetic markers of productivity traits in livestock can be developed based on the TERT gene polymorphism for use in marker-associated selection (MAS). In this study, a bioinformatic-based approach is proposed to evaluate the effect of missense single-nucleotide polymorphisms (SNPs) in the pig TERT gene on enzyme function and structure, with the prospect of developing genetic markers. RESULTS: A comparative analysis of the coding and amino acid sequences of the pig TERT was performed with corresponding sequences of other species. The distribution of polymorphisms in the pig TERT gene, with respect to the enzyme's structural-functional domains, was established. A three-dimensional model of the pig TERT structure was obtained through homological modeling. The potential impact of each of the 23 missense SNPs in the pig TERT gene on telomerase function and stability was assessed using predictive bioinformatic tools utilizing data on the amino acid sequence and structure of pig TERT. CONCLUSIONS: According to bioinformatic analysis of 23 missense SNPs of the pig TERT gene, a predictive effect of rs789641834 (TEN domain), rs706045634 (TEN domain), rs325294961 (TRBD domain) and rs705602819 (RTD domain) on the structural and functional parameters of the enzyme was established. These SNPs hold the potential to serve as genetic markers of productivity traits. Therefore, the possibility of their application in MAS should be further evaluated in associative analysis studies.

BACKGROUND: Nutrition has not only an impact on the general wellbeing of an animal but can also affect reproductive processes. In cattle, feeding regimes can influence the age of puberty onset and alter gonadal development. We analyzed effects of different milk replacer (MR) feeding regimes during rearing on ovarian physiology with specific emphasis on the numbers as well as gene expression characteristics of granulosa cells (GCs) at the age of puberty onset. Two groups of calves received either 10% or 20% of bodyweight MR per day during their first 8 weeks. After weaning, both groups were fed the same mixed ration ad libitum until slaughter at 8 months. RESULTS: Animals of the 20% feeding group had a significantly higher body weight, but the proportion of animals having a corpus luteum at the time of slaughter was not different between groups, suggesting a similar onset of puberty. Calves of the 10% group showed a constant GC count regardless of the number of follicles (r = 0.23) whereas in the 20% group increasing numbers of GCs were detected with a higher follicle count (r = 0.71). As a first effort to find a possible molecular explanation for this unexpected limitation of GC numbers in the 10% group, we comparatively analyzed GC transcriptomes in both diet groups. The mRNA microarray analysis revealed a total of 557 differentially expressed genes comparing both groups (fold change > |1.5| and p < 0.05). OAS1X, MX2 and OAS1Z were among the top downregulated genes in the 20% vs. the 10% group, whereas top upregulated genes comprised BOLA and XCL1. All of these genes are known to be regulated by interferon. Subsequent signaling pathway analysis revealed the involvement of several immune response mechanisms in accordance with a number of interferons as upstream regulators. CONCLUSIONS: The results indicate that the plane of MR feeding in early life has an impact on the number and physiology of GCs later in life. This might influence the overall reproductive life initiated by the onset of puberty in cattle. In addition, the observed alterations in GCs of calves fed less MR might be a consequence of interferon regulated immunological pathways.

BACKGROUND: Babesia caballi is an intraerythrocytic parasite from the phylum Apicomplexa, capable of infecting equids and causing equine piroplasmosis. However, since there is limited genome information available on B. caballi, molecular mechanisms involved in host specificity and pathogenicity of this species have not been fully elucidated yet. RESULTS: Genomic DNA from a B. caballi subclone was purified and sequenced using both Illumina and Nanopore technologies. The resulting assembled sequence consisted of nine contigs with a size of 12.9 Mbp, rendering a total of 5,910 protein-coding genes. The phylogenetic tree of Apicomplexan species was reconstructed using 263 orthologous genes. We identified 481 ves1-like genes and named "ves1c". In contrast, expansion of the major facilitator superfamily (mfs) observed in closely related B. bigemina and B. ovata species was not found in B. caballi. A set of repetitive units containing an open reading frame with a size of 297 bp was also identified. CONCLUSIONS: We present a chromosome-level genome assembly of B. caballi. Our genomic data may contribute to estimating gene expansion events involving multigene families and exploring the evolution of species from this genus.

BACKGROUND: The natural products, metabolites, of gut microbes are crucial effect factors on diseases. Comprehensive identification and annotation of relationships among disease, metabolites, and microbes can provide efficient and targeted solutions towards understanding the mechanism of complex disease and development of new markers and drugs. RESULTS: We developed Gut Microbial Metabolite Association with Disease (GMMAD), a manually curated database of associations among human diseases, gut microbes, and metabolites of gut microbes. Here, this initial release (i) contains 3,836 disease-microbe associations and 879,263 microbe-metabolite associations, which were extracted from literatures and available resources and then experienced our manual curation; (ii) defines an association strength score and a confidence score. With these two scores, GMMAD predicted 220,690 disease-metabolite associations, where the metabolites all belong to the gut microbes. We think that the positive effective (with both scores higher than suggested thresholds) associations will help identify disease marker and understand the pathogenic mechanism from the sense of gut microbes. The negative effective associations would be taken as biomarkers and have the potential as drug candidates. Literature proofs supported our proposal with experimental consistence; (iii) provides a user-friendly web interface that allows users to browse, search, and download information on associations among diseases, metabolites, and microbes. The resource is freely available at http://guolab.whu.edu.cn/GMMAD . CONCLUSIONS: As the online-available unique resource for gut microbial metabolite-disease associations, GMMAD is helpful for researchers to explore mechanisms of disease- metabolite-microbe and screen the drug and marker candidates for different diseases.

BACKGROUND: Wulong geese (Anser cygnoides orientalis) are known for their excellent egg-laying performance. However, they show considerable population differences in egg-laying behavior. This study combined genome-wide selection signal analysis with transcriptome analysis (RNA-seq) to identify the genes related to high egg production in Wulong geese. RESULTS: A total of 132 selected genomic regions were screened using genome-wide selection signal analysis, and 130 genes related to high egg production were annotated in these regions. These selected genes were enriched in pathways related to egg production, including oocyte meiosis, the estrogen signaling pathway, the oxytocin signaling pathway, and progesterone-mediated oocyte maturation. Furthermore, a total of 890 differentially expressed genes (DEGs), including 340 up-regulated and 550 down-regulated genes, were identified by RNA-seq. Two genes - GCG and FAP - were common to the list of selected genes and DEGs. A non-synonymous single nucleotide polymorphism was identified in an exon of FAP. CONCLUSIONS: Based on genome-wide selection signal analysis and transcriptome data, GCG and FAP were identified as candidate genes associated with high egg production in Wulong geese. These findings could promote the breeding of Wulong geese with high egg production abilities and provide a theoretical basis for exploring the mechanisms of reproductive regulation in poultry.

BACKGROUND: Sepsis is a high mortality disease which seriously threatens human life and health, for which the pathogenetic mechanism still unclear. There is increasing evidence showed that immune and inflammation responses are key players in the development of sepsis pathology. LncRNAs, which act as ceRNAs, have critical roles in various diseases. However, the regulatory roles of ceRNA in the immunopathogenesis of sepsis have not yet been elucidated. RESULTS: In this study, we aimed to identify immune biomarkers associated with sepsis. We first generated a global immune-associated ceRNA (IMCE) network based on data describing interactions pairs of gene-miRNA and miRNA-lncRNA. Afterward, we excavated a dysregulated sepsis immune-associated ceRNA (SPIMC) network from the global IMCE network by means of a multi-step computational approach. Functional enrichment indicated that lncRNAs in SPIMC network have pivotal roles in the immune mechanism underlying sepsis. Subsequently, we identified module and hub genes (CD4 and STAT4) via construction of a sepsis immune-related PPI network. Then, we identified hub genes based on the modular structure of PPI network and generated a ceRNA subnetwork to analyze key lncRNAs associated with sepsis. Finally, 6 lncRNAs (LINC00265, LINC00893, NDUFA6-AS1, NOP14-AS1, PRKCQ-AS1 and ZNF674-AS1) that identified as immune biomarkers of sepsis. Moreover, the CIBERSORT algorithm and the infiltration of circulating immune cells types were performed to identify the inflammatory state of sepsis. Correlation analyses between immune cells and sepsis immune biomarkers showed that the LINC00265 was strongly positive correlated with the macrophages M2 (r = 0.77). CONCLUSION: Collectively, these results may suggest that these lncRNAs (LINC00265, LINC00893, NDUFA6-AS1, NOP14-AS1, PRKCQ-AS1 and ZNF674-AS1) played important roles in the immune pathogenesis of sepsis and provide potential therapeutic targets for further researches on immune therapy treatment in patients with sepsis.

BACKGROUND: Because of its social nature, the honeybee is regularly exposed to environmental toxicants such as heavy metals and xenobiotics. These toxicants are known to exert strong selective pressure on the gut microbiome's structure and diversity. For example, resistant microbial members are more likely to dominate in maintaining a stable microbiome, which is critical for bee health. Therefore, the aim of this study was to examine the Enterococcus faecium strains isolated from bee guts for their in vitro growth and tolerability to diverse heavy metals and xenobiotics. An additional aim was to analyze the genomes of E. faecium isolates to assess the molecular bases of resistance and compare them with E. faecium species isolated from other environmental sources. RESULTS: The E. faecium bee isolates were able to tolerate high levels (up to 200 mg/L) of toxicants, including cadmium, zinc, benzoate, phenol and hexane. Moreover, the isolates could tolerate toluene and copper at up to 100 mg/L. The genome of E. faecium Am5, isolated from the larval stage of Apis mellifera gut, was about 2.7 Mb in size, had a GC content of 37.9% and 2,827 predicted coding sequences. Overall, the Am5 genome features were comparable with previously sequenced bee-gut isolates, E. faecium Am1, Bee9, SM21, and H7. The genomes of the bee isolates provided insight into the observed heavy metal tolerance. For example, heavy metal tolerance and/or regulation genes were present, including czcD (cobalt/zinc/cadmium resistance), cadA (exporting ATPase), cutC (cytoplasmic copper homeostasis) and zur (zinc uptake regulation). Additionally, genes associated with nine KEGG xenobiotic biodegradation pathways were detected, including γ-hexachlorocyclohexane, benzoate, biphenyl, bisphenol A, tetrachloroethene, 1,4-dichlorobenzene, ethylbenzene, trinitrotoluene and caprolactam. Interestingly, a comparative genomics study demonstrated the conservation of toxicant resistance genes across a variety of E. faecium counterparts isolated from other environmental sources such as non-human mammals, humans, avians, and marine animals. CONCLUSIONS: Honeybee gut-derived E. faecium strains can tolerate a variety of heavy metals. Moreover, their genomes encode many xenobiotic biodegradation pathways. Further research is required to examine E. faecium strains potential to boost host resistance to environmental toxins.

BACKGROUND: Heat shock protein 20 (HSP20) is a member of the heat stress-related protein family, which plays critical roles in plant growth, development, and response to abiotic stresses. Although many HSP20 genes have been associated with heat stress in numerous types of plants, little is known about the details of the HSP20 gene family in Coix. To investigate the mechanisms of the ClHSP20 response to heat and drought stresses, the ClHSP20 gene family in Coix was identified and characterized based on genome-wide analysis. RESULTS: A total of 32 putative ClHSP20 genes were identified and characterized in Coix. Phylogenetic analysis indicated that ClHSP20s were grouped into 11 subfamilies. The duplicated event analysis demonstrated that tandem duplication and segment duplication events played crucial roles in promoting the expansion of the ClHSP20 gene family. Synteny analysis showed that Coix shared the highest homology in 36 HSP20 gene pairs with wheat, followed by 22, 19, 15, and 15 homologous gene pairs with maize, sorghum, barley, and rice, respectively. The expression profile analysis showed that almost all ClHSP20 genes had different expression levels in at least one tissue. Furthermore, 22 of the 32 ClHSP20 genes responded to heat stress, with 11 ClHSP20 genes being significantly upregulated and 11 ClHSP20 genes being significantly downregulated. Furthermore, 13 of the 32 ClHSP20 genes responded to drought stress, with 6 ClHSP20 genes being significantly upregulated and 5 ClHSP20 genes being significantly downregulated. CONCLUSIONS: Thirty-two ClHSP20 genes were identified and characterized in the genome of Coix. Tandem and segmental duplication were identified as having caused the expansion of the ClHSP20 gene family. The expression patterns of the ClHSP20 genes suggested that they play a critical role in growth, development, and response to heat and drought stress. The current study provides a theoretical basis for further research on ClHSP20s and will facilitate the functional characterization of ClHSP20 genes.

BACKGROUND: Reproduction is a key feature of the sustainability of a species and thus represents an important component in livestock genetic improvement programs. Most reproductive traits are lowly heritable. In order to gain a better understanding of the underlying genetic basis of these traits, a genome-wide association was conducted for age at first calving (AFC), first inter-calving period (ICP) and scrotal circumference (SC) within the South African Bonsmara breed. Phenotypes and genotypes (120,692 single nucleotide polymorphisms (SNPs) post editing) were available on 7,128 South African Bonsmara cattle; the association analyses were undertaken using linear mixed models. RESULTS: Genomic restricted maximum likelihood analysis of the 7,128 SA Bonsmara cattle yielded genomic heritability's of 0.183 (SE = 0.021) for AFC, 0.207 (SE = 0.022) for ICP and 0.209 (SE = 0.019) for SC. A total of 16, 23 and 51 suggestive (P ≤ 4 × 10-6) SNPs were associated with AFC, ICP and SC, while 11, 11 and 44 significant (P ≤ 4 × 10-7) SNPs were associated with AFC, ICP and SC respectively. A total of 11 quantitative trait loci (QTL) and 11 candidate genes were co-located with these associated SNPs for AFC, with 10 QTL harbouring 11 candidate genes for ICP and 41 QTL containing 40 candidate genes for SC. The QTL identified were close to genes previously associated with carcass, fertility, growth and milk-related traits. The biological pathways influenced by these genes include carbohydrate catabolic processes, cellular development, iron homeostasis, lipid metabolism and storage, immune response, ovarian follicle development and the regulation of DNA transcription and RNA translation. CONCLUSIONS: This was the first attempt to study the underlying polymorphisms associated with reproduction in South African beef cattle. Genes previously reported in cattle breeds for numerous traits bar AFC, ICP or SC were detected in this study. Over 20 different genes have not been previously reported in beef cattle populations and may have been associated due to the unique genetic composite background of the SA Bonsmara breed.

BACKGROUND: Numerous circular RNAs (circRNAs) have been recently identified in porcine tissues and cell types. Nevertheless, their significance in porcine spleen development is yet unelucidated. Herein, we reported an extensive overlook of circRNA expression profile during spleen development in Meishan pigs. RESULTS: Overall, 39,641 circRNAs were identified from 6,914 host genes. Among them, many circRNAs are up- or down-regulated at different time points of pig spleen development. Using WGCNA analysis, we revealed two essential modules for protein-coding genes and circRNAs. Subsequent correlation analysis revealed 67 circRNAs/co-expressed genes that participated in immnue-associated networks. Furthermore, a competing endogenous RNA (ceRNA) network analysis of circRNAs revealed that 12 circRNAs modulated CD226, MBD2, SAMD3, SIT1, SRP14, SYTL3 gene expressions via acting as miRNA sponges. Moreover, the circRNA_21767/miR-202-3p axis regulated SIT1 expression in a ceRNA manner, which is critical for the immune-based regulation of spleen development in Meishan pigs. CONCLUSIONS: Overall, our results demonstrated that the circRNAs were differentially expressed during different stages of porcine spleen development, meanwhile the circRNAs interacted with immune-related genes in a ceRNA-based fashion. Moreover, we presented biomedical researchers with RNAseqTools, a user-friendly and powerful software for the visualization of transcriptome profile data.

BACKGROUND: Tilapia is one of the most essential farmed fishes in the world. It is a tropical and subtropical freshwater fish well adapted to warm water but sensitive to cold weather. Extreme cold weather could cause severe stress and mass mortalities in tilapia. The present study was carried out to investigate the effects of cold stress on the up-regulation of antifreeze protein (AFP) genes in Nile tilapia (Oreochromis niloticus). Two treatment groups of fish were investigated (5 replicates of 15 fish for each group in fibreglass tanks/70 L each): 1) a control group; the fish were acclimated to lab conditions for two weeks and the water temperature was maintained at 25 °C during the whole experimental period with feeding on a commercial diet (30% crude protein). 2) Cold stress group; the same conditions as the control group except for the temperature. Initially, the temperature was decreased by one degree every 12 h. The fish started showing death symptoms when the water temperature reached 6-8 °C. In this stage the tissue (muscle) samples were taken from both groups. The immune response of fish exposed to cold stress was detected and characterized using Differential Display-PCR (DD-PCR). RESULTS: The results indicated that nine different up-regulation genes were detected in the cold-stressed fish compared to the control group. These genes are Integrin-alpha-2 (ITGA-2), Gap junction gamma-1 protein-like (GJC1), WD repeat-containing protein 59 isoform X2 (WDRP59), NUAK family SNF1-like kinase, G-protein coupled receptor-176 (GPR-176), Actin cytoskeleton-regulatory complex protein pan1-like (PAN-1), Whirlin protein (WHRN), Suppressor of tumorigenicity 7 protein isoform X2 (ST7P) and ATP-binding cassette sub-family A member 1-like isoform X2 (ABCA1). The antifreeze gene type-II amplification using a specific PCR product of 600 bp, followed by cloning and sequencing analysis revealed that the identified gene is antifreeze type-II, with similarity ranging from 70 to 95%. The in-vitro transcribed gene induced an antifreeze protein with a molecular size of 22 kDa. The antifreeze gene, ITGA-2 and the WD repeat protein belong to the lectin family (sugar-protein). CONCLUSIONS: In conclusion, under cold stress, Nile tilapia express many defence genes, an antifreeze gene consisting of one open reading frame of approximately 0.6 kbp.

The genus Sophora (Fabaceae) includes medicinal plants that have been used in East Asian countries since antiquity. Sophora flavescens is a perennial herb indigenous to China, India, Japan, Korea, and Russia. Its dried roots have antioxidant, anti-inflammatory, antibacterial, apoptosis-modulating, and antitumor efficacy. The congeneric S. koreensis is endemic to Korea and its genome is less than half the size of that of S. flavescens. Nevertheless, this discrepancy can be used to assemble and validate the S. flavescens genome. A comparative genomic study of the two genomes can disclose the recent evolutionary divergence of the polymorphic phenotypic profiles of these species. Here, we used the PacBio sequencing platform to sequence and assemble the S. koreensis and S. flavescens genomes. We inferred that it was mainly small-scale duplication that occurred in S. flavescens. A KEGG analysis revealed pathways that might regulate the pharmacologically important secondary metabolites in S. flavescens and S. koreensis. The genome assemblies of Sophora spp. could be used in comparative genomics and data mining for various plant natural products.

BACKGROUND: The glyoxalase system includes glyoxalase I (GLXI), glyoxalase II (GLXII) and glyoxalase III (GLXIII), which are responsible for methylglyoxal (MG) detoxification and involved in abiotic stress responses such as drought, salinity and heavy metal. RESULTS: In this study, a total of 620 GLX family genes were identified from 21 different plant species. The results of evolutionary analysis showed that GLX genes exist in all species from lower plants to higher plants, inferring that GLX genes might be important for plants, and GLXI and GLXII account for the majority. In addition, motif showed an expanding trend in the process of evolution. The analysis of cis-acting elements in 21 different plant species showed that the promoter region of the GLX genes were rich in phytohormones and biotic and abiotic stress-related elements, indicating that GLX genes can participate in a variety of life processes. In cotton, GLXs could be divided into two groups and most GLXIs distributed in group I, GLXIIs and GLXIIIs mainly belonged to group II, indicating that there are more similarities between GLXII and GLXIII in cotton evolution. The transcriptome data analysis and quantitative real-time PCR analysis (qRT-PCR) show that some members of GLX family would respond to high temperature treatment in G.hirsutum. The protein interaction network of GLXs in G.hirsutum implied that most members can participate in various life processes through protein interactions. CONCLUSIONS: The results elucidated the evolutionary history of GLX family genes in plants and lay the foundation for their functions analysis in cotton.

BACKGROUND: Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection. RESULTS: A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively. CONCLUSIONS: The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.

BACKGROUND: Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group. RESULTS: We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy. CONCLUSIONS: The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

BACKGROUND: The nutritional status during early life can have enduring effects on an animal's metabolism, although the mechanisms underlying these long-term effects are still unclear. Epigenetic modifications are considered a prime candidate mechanism for encoding early-life nutritional memories during this critical developmental period. However, the extent to which these epigenetic changes occur and persist over time remains uncertain, in part due to challenges associated with directly stimulating the fetus with specific nutrients in viviparous mammalian systems. RESULTS: In this study, we used medaka as an oviparous vertebrate model to establish an early-life high-fat diet (HFD) model. Larvae were fed with HFD from the hatching stages (one week after fertilization) for six weeks, followed by normal chow (NC) for eight weeks until the adult stage. We examined the changes in the transcriptomic and epigenetic state of the liver over this period. We found that HFD induces simple liver steatosis, accompanied by drastic changes in the hepatic transcriptome, chromatin accessibility, and histone modifications, especially in metabolic genes. These changes were largely reversed after the long-term NC, demonstrating the high plasticity of the epigenetic state in hepatocytes. However, we found a certain number of genomic loci showing non-reversible epigenetic changes, especially around genes related to cell signaling, liver fibrosis, and hepatocellular carcinoma, implying persistent changes in the cellular state of the liver triggered by early-life HFD feeding. CONCLUSION: In summary, our data show that early-life HFD feeding triggers both reversible and persistent epigenetic changes in medaka hepatocytes. Our data provide novel insights into the epigenetic mechanism of nutritional programming and a comprehensive atlas of the long-term epigenetic state in an early-life HFD model of non-mammalian vertebrates.

BACKGROUND: Drought stress is a prevalent abiotic stress that significantly hinders the growth and development of plants. According to studies, β-aminobutyric acid (BABA) can influence the ABA pathway through the AtIBI1 receptor gene to enhance cold resistance in Arabidopsis. However, the Aspartate tRNA-synthetase (AspRS) gene family, which acts as the receptor for BABA, has not yet been investigated in poplar. Particularly, it is uncertain how the AspRS gene family (PtrIBIs)r can resist drought stress after administering various concentrations of BABA to poplar. RESULTS: In this study, we have identified 12 AspRS family genes and noted that poplar acquired four PtrIBI pairs through whole genome duplication (WGD). We conducted cis-action element analysis and found a significant number of stress-related action elements on different PtrIBI genes promoters. The expression of most PtrIBI genes was up-regulated under beetle and mechanical damage stresses, indicating their potential role in responding to leaf damage stress. Our results suggest that a 50 mM BABA treatment can alleviate the damage caused by drought stress in plants. Additionally, via transcriptome sequencing, we observed that the partial up-regulation of BABA receptor genes, PtrIBI2/4/6/8/11, in poplars after drought treatment. We hypothesize that poplar responds to drought stress through the BABA-PtrIBIs-PtrVOZ coordinated ABA signaling pathway. Our research provides molecular evidence for understanding how plants respond to drought stress through external application of BABA. CONCLUSIONS: In summary, our study conducted genome-wide analysis of the AspRS family of P. trichocarpa and identified 12 PtrIBI genes. We utilized genomics and bioinformatics to determine various characteristics of PtrIBIs such as chromosomal localization, evolutionary tree, gene structure, gene doubling, promoter cis-elements, and expression profiles. Our study found that certain PtrIBI genes are regulated by drought, beetle, and mechanical damage implying their crucial role in enhancing poplar stress tolerance. Additionally, we observed that external application of low concentrations of BABA increased plant drought resistance under drought stress. Through the BABA-PtrIBIs-PtrVOZ signaling module, poplar plants were able to transduce ABA signaling and regulate their response to drought stress. These results suggest that the PtrIBI genes in poplar have the potential to improve drought tolerance in plants through the topical application of low concentrations of BABA.

BACKGROUND: The evolution of pregnancy-specific glycoprotein (PSG) genes within the CEA gene family of primates correlates with the evolution of hemochorial placentation about 45 Myr ago. Thus, we hypothesized that hemochorial placentation with intimate contact between fetal cells and maternal immune cells favors the evolution and expansion of PSGs. With only a few exceptions, all rodents have hemochorial placentas thus the question arises whether Psgs evolved in all rodent genera. RESULTS: In the analysis of 94 rodent species from 4 suborders, we identified Psg genes only in the suborder Myomorpha in three families (characteristic species in brackets), namely Muridae (mouse), Cricetidae (hamster) and Nesomyidae (giant pouched rat). All Psgs are located, as previously described for mouse and rat, in a region of the genome separated from the Cea gene family locus by several megabases, further referred to as the rodent Psg locus. In the suborders Castorimorpha (beaver), Hystricognatha (guinea pig) and Sciuromorpha (squirrel), neither Psg genes nor so called CEA-related cell adhesion molecule (Ceacam) genes were found in the Psg locus. There was even no evidence for the existence of Psgs in any other genomic region. In contrast to the Psg-harboring rodent species, which do not have activating CEACAMs, we were able to identify Ceacam genes encoding activating CEACAMs in all other rodents studied. In the Psg locus, there are genes encoding three structurally distinct CEACAM/PSGs: (i) CEACAMs composed of one N- and one A2-type domain (CEACAM9, CEACAM15), (ii) composed of two N domains (CEACAM11-CEACAM14) and (iii) composed of three to eight N domains and one A2 domain (PSGs). All of them were found to be secreted glycoproteins preferentially expressed by trophoblast cells, thus they should be considered as PSGs. CONCLUSION: In rodents Psg genes evolved only recently in the suborder Myomorpha shortly upon their most recent common ancestor (MRCA) has coopted the retroviral genes syncytin-A and syncytin-B which enabled the evolution of the three-layered trophoblast. The expansion of Psgs is limited to the Psg locus most likely after a translocation of a CEA-related gene - possibly encoding an ITAM harboring CEACAM. According to the expression pattern two waves of gene amplification occurred, coding for structurally different PSGs.

The H9c2 myoblast cell line, isolated from the left ventricular tissue of rat, is currently used in vitro as a mimetic for skeletal and cardiac muscle due to its biochemical, morphological, and electrical/hormonal signaling properties. During culture, H9c2 cells acquire a myotube phenotype, where a critical component is the inclusion of retinoic acid (RA). The results from some authors on H9c2 suggested that thousands of genes respond to RA stimuli, while others report hundreds of genes responding to RA over different cell types. In this article, using a more appropriate experimental design, we first confirm the H9c2 cardiac phenotype with and without RA and report transcriptomic and physiological changes regarding calcium handling, bioenergetics, and other biological concepts. Interestingly, of the 2360 genes showing a transcriptional change, 622 genes were statistically associated with the RA response. Of these genes, only 305 were RA-specific, and the rest also showed a culture-time component. Thus, the major expression changes (from 74 to 87%) were indeed due to culture conditions over time. Unexpectedly, only a few components of the retinol pathway in KEGG responded to RA. Our results show the role of RA in the H9c2 cultures impacting the interpretation using H9c2 as an in vitro model.

BACKGROUND: Isatis tinctoria Linnaeus and Isatis indigotica Fortune are very inconsistent in their morphological characteristics, but the Flora of China treats them as the same species. In this work, a new technology that differs from conventional barcodes is developed to prove that they are different species and to clarify their classification. RESULTS AND METHODS: I. indigotica was indistinguishable from I. tinctoria when using ITS2. CPGAVAS2 was used to construct the chloroplast genomes. MAFFT and DnaSP were used to calculate nucleotide polymorphism, the chloroplast genomes of the two have high diversity in the rpl32 ~ trnL-UAG short region. When using this region as a mini barcode, it was found that there are obvious differences in the base numbers of I. tinctoria and different ploidy I. indigotica were found, but diploid and tetraploid I. indigotica had the same number of bases. Moreover, the reconstruction of the maximum likelihood (ML) tree, utilizing the mini-barcode, demonstrated that I. tinctoria and both diploid and tetraploid I. indigotica are located on distinct branches. The genome size of tetraploid I. indigotica was approximately 643.773 MB, the heterozygosity rate was approximately 0.98%, and the repeat sequence content was approximately 90.43%. This species has a highly heterozygous, extremely repetitive genome. CONCLUSION: A new method was established to differentiate between I. indigotica and I. tinctoria. Furthermore, this approach provides a reference and basis for the directional breeding of Isatis.

BACKGROUND: Plants possess mitochondrial genomes that are large and complex compared to animals. Despite their size, plant mitochondrial genomes do not contain significantly more genes than their animal counterparts. Studies into the sequence and structure of plant mitochondrial genomes heavily imply that the main mechanism driving replication of plant mtDNA, and offer valuable insights into plant evolution, energy production, and environmental adaptation. RESULTS: This study presents the first comprehensive analysis of Agrostis stolonifera's mitochondrial genome, characterized by a branched structure comprising three contiguous chromosomes, totaling 560,800 bp with a GC content of 44.07%. Annotations reveal 33 unique protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The predominant codons for alanine and glutamine are GCU and CAA, respectively, while cysteine and phenylalanine exhibit weaker codon usage biases. The mitogenome contains 73, 34, and 23 simple sequence repeats (SSRs) on chromosomes 1, 2, and 3, respectively. Chromosome 1 exhibits the most frequent A-repeat monomeric SSR, whereas chromosome 2 displays the most common U-repeat monomeric SSR. DNA transformation analysis identifies 48 homologous fragments between the mitogenome and chloroplast genome, representing 3.41% of the mitogenome's total length. The PREP suite detects 460 C-U RNA editing events across 33 mitochondrial PCGs, with the highest count in the ccmFn gene and the lowest in the rps7 gene. Phylogenetic analysis confirms A. stolonifera's placement within the Pooideae subfamily, showing a close relationship to Lolium perenne, consistent with the APG IV classification system. Numerous homologous co-linear blocks are observed in A. stolonifera's mitogenomes and those of related species, while certain regions lack homology. CONCLUSIONS: The unique features and complexities of the A. stolonifera mitochondrial genome, along with its similarities and differences to related species, provide valuable insights into plant evolution, energy production, and environmental adaptation. The findings from this study significantly contribute to the growing body of knowledge on plant mitochondrial genomes and their role in plant biology.

BACKGROUND: Phloem protein 2 (PP2) proteins play a vital role in the Phloem-based defense (PBD) and participate in many abiotic and biotic stress. However, research on PP2 proteins in cotton is still lacking. RESULTS: A total of 25, 23, 43, and 47 PP2 genes were comprehensively identified and characterized in G.arboretum, G.raimondii, G.barbadense, and G.hirsutum. The whole genome duplication (WGD) and allopolyploidization events play essential roles in the expansion of PP2 genes. The promoter regions of GhPP2 genes contain many cis-acting elements related to abiotic stress and the weighted gene co-expression network analysis (WGCNA) analysis displayed that GhPP2s could be related to salt stress. The qRT-PCR assays further confirmed that GhPP2-33 could be dramatically upregulated during the salt treatment. And the virus-induced gene silencing (VIGS) experiment proved that the silencing of GhPP2-33 could decrease salt tolerance. CONCLUSIONS: The results in this study not only offer new perspectives for understanding the evolution of PP2 genes in cotton but also further explore their function under salt stress.

BACKGROUND: Numerous studies in the past have expanded our understanding of the genetic differences of global distributed cyanobacteria that originated around billions of years ago, however, unraveling how gene gain and loss drive the genetic evolution of cyanobacterial species, and the trade-off of these evolutionary forces are still the central but poorly understood issues. RESULTS: To delineate the contribution of gene flow in mediating the hereditary differentiation and shaping the microbial evolution, a global genome-wide study of bloom-forming cyanobacterium, Microcystis aeruginosa species complex, provided robust evidence for genetic diversity, reflected by enormous variation in gene repertoire among various strains. Mathematical extrapolation showed an 'open' microbial pan-genome of M. aeruginosa species, since novel genes were predicted to be introduced after new genomes were sequenced. Identification of numerous horizontal gene transfer's signatures in genome regions of interest suggested that genome expansion via transformation and phage-mediated transduction across bacterial lineage as an evolutionary route may contribute to the differentiation of Microcystis functions (e.g., carbohydrate metabolism, amino acid metabolism, and energy metabolism). Meanwhile, the selective loss of some dispensable genes at the cost of metabolic versatility is as a mean of adaptive evolution that has the potential to increase the biological fitness. CONCLUSIONS: Now that the recruitment of novel genes was accompanied by a parallel loss of some other ones, a trade-off in gene content may drive the divergent differentiation of M. aeruginosa genomes. Our study provides a genetic framework for the evolution of M. aeruginosa species and illustrates their possible evolutionary patterns.

BACKGROUND: Folic acid is a water-soluble B vitamin (B9), which is closely related to the body's immune and other metabolic pathways. The folic acid synthesized by rumen microbes has been unable to meet the needs of high-yielding dairy cows. The incidence rate of subclinical mastitis in dairy herds worldwide ranged between 25%~65% with no obvious symptoms, but it significantly causes a decrease in lactation and milk quality. Therefore, this study aims at exploring the effects of folic acid supplementation on the expression profile of lncRNAs, exploring the molecular mechanism by which lncRNAs regulate immunity in subclinical mastitic dairy cows. RESULTS: The analysis identified a total of 4384 lncRNA transcripts. Subsequently, differentially expressed lncRNAs in the comparison of two groups (SF vs. SC, HF vs. HC) were identified to be 84 and 55 respectively. Furthermore, the weighted gene co-expression network analysis (WGCNA) and the KEGG enrichment analysis result showed that folic acid supplementation affects inflammation and immune response-related pathways. The two groups have few pathways in common. One important lncRNA MSTRG.11108.1 and its target genes (ICAM1, CCL3, CCL4, etc.) were involved in immune-related pathways. Finally, through integrated analysis of lncRNAs with GWAS data and animal QTL database, we found that differential lncRNA and its target genes could be significantly enriched in SNPs and QTLs related to somatic cell count (SCC) and mastitis, such as MSTRG.11108.1 and its target gene ICAM1, CXCL3, GRO1. CONCLUSIONS: For subclinical mastitic cows, folic acid supplementation can significantly affect the expression of immune-related pathway genes such as ICAM1 by regulating lncRNAs MSTRG.11108.1, thereby affecting related immune phenotypes. Our findings laid a ground foundation for theoretical and practical application for feeding folic acid supplementation in subclinical mastitic cows.

as one of the oldest cultivated crops in the world, quinoa has been widely valued for its rich nutritional value and green health. In this study, 22 CDK genes (CqCDK01-CqCDK22) were identified from quinoa genome using bioinformatics method. The number of amino acids was 173-811, the molecular weight was 19,554.89 Da-91,375.70 Da, and the isoelectric point was 4.57-9.77. The phylogenetic tree divided 21 CqCDK genes into six subfamilies, the gene structure showed that 12 (54.5%) CqCDK genes (CqCDK03, CqCDK04, CqCDK05, CqCDK06, CqCDK07, CqCDK11, CqCDK14, CqCDK16, CqCDK18, CqCDK19, CqCDK20 and CqCDK21) had UTR regions at 5' and 3' ends. Each CDK protein had different motifs (3-9 motifs), but the genes with the same motifs were located in the same branch. Promoter analysis revealed 41 cis-regulatory elements related to plant hormones, abiotic stresses, tissue-specific expression and photoresponse. The results of real-time fluorescence quantitative analysis showed that the expression level of some CDK genes was higher under drought and salt stress, which indicated that CDK genes could help plants to resist adverse environmental effects. Subcellular localization showed that CqCDK15 gene was localized to the nucleus and cytoplasm, and transgenic plants overexpressing CqCDK15 gene showed higher drought and salt tolerance compared to the controls. Therefore, CDK genes are closely related to quinoa stress resistance. In this study, the main functions of quinoa CDK gene family and its expression level in different tissues and organs were analyzed in detail, which provided some theoretical support for quinoa stress-resistant breeding. Meanwhile, this study has important implications for further understanding the function of the CDK gene family in quinoa and our understanding of the CDK family in vascular plant.

BACKGROUND: In breast cancer, as in all cancers, genetic and epigenetic deregulations can result in out-of-context expressions of a set of normally silent tissue-specific genes. The activation of some of these genes in various cancers empowers tumours cells with new properties and drives enhanced proliferation and metastatic activity, leading to a poor survival prognosis. RESULTS: In this work, we undertook an unprecedented systematic and unbiased analysis of out-of-context activations of a specific set of tissue-specific genes from testis, placenta and embryonic stem cells, not expressed in normal breast tissue as a source of novel prognostic biomarkers. To this end, we combined a strict machine learning framework of transcriptomic data analysis, and successfully created a new robust tool, validated in several independent datasets, which is able to identify patients with a high risk of relapse. This unbiased approach allowed us to identify a panel of five biomarkers, DNMT3B, EXO1, MCM10, CENPF and CENPE, that are robustly and significantly associated with disease-free survival prognosis in breast cancer. Based on these findings, we created a new Gene Expression Classifier (GEC) that stratifies patients. Additionally, thanks to the identified GEC, we were able to paint the specific molecular portraits of the particularly aggressive tumours, which show characteristics of male germ cells, with a particular metabolic gene signature, associated with an enrichment in pro-metastatic and pro-proliferation gene expression. CONCLUSIONS: The GEC classifier is able to reliably identify patients with a high risk of relapse at early stages of the disease. We especially recommend to use the GEC tool for patients with the luminal-A molecular subtype of breast cancer, generally considered of a favourable disease-free survival prognosis, to detect the fraction of patients undergoing a high risk of relapse.

BACKGROUND: Approximately 4-8% of the world suffers from a rare disease. Rare diseases are often difficult to diagnose, and many do not have approved therapies. Genetic sequencing has the potential to shorten the current diagnostic process, increase mechanistic understanding, and facilitate research on therapeutic approaches but is limited by the difficulty of novel variant pathogenicity interpretation and the communication of known causative variants. It is unknown how many published rare disease variants are currently accessible in the public domain. RESULTS: This study investigated the translation of knowledge of variants reported in published manuscripts to publicly accessible variant databases. Variants, symptoms, biochemical assay results, and protein function from literature on the SLC6A8 gene associated with X-linked Creatine Transporter Deficiency (CTD) were curated and reported as a highly annotated dataset of variants with clinical context and functional details. Variants were harmonized, their availability in existing variant databases was analyzed and pathogenicity assignments were compared with impact algorithm predictions. 24% of the pathogenic variants found in PubMed articles were not captured in any database used in this analysis while only 65% of the published variants received an accurate pathogenicity prediction from at least one impact prediction algorithm. CONCLUSIONS: Despite being published in the literature, pathogenicity data on patient variants may remain inaccessible for genetic diagnosis, therapeutic target identification, mechanistic understanding, or hypothesis generation. Clinical and functional details presented in the literature are important to make pathogenicity assessments. Impact predictions remain imperfect but are improving, especially for single nucleotide exonic variants, however such predictions are less accurate or unavailable for intronic and multi-nucleotide variants. Developing text mining workflows that use natural language processing for identifying diseases, genes and variants, along with impact prediction algorithms and integrating with details on clinical phenotypes and functional assessments might be a promising approach to scale literature mining of variants and assigning correct pathogenicity. The curated variants list created by this effort includes context details to improve any such efforts on variant curation for rare diseases.

BACKGROUND: Acute heat stress could induce high mortality and cause huge economic losses in the poultry industry. Although many studies have revealed heat stress-induced injuries of multiple tissues, the main target tissue and molecular mechanism of death under acute heat stress was largely unknown. This study systematically compared the transcriptome data of five main visceral tissues in chickens to reveal the response of multiple tissues to acute heat stress and determine the main target tissue of acute heat stress, further revealing the injuries of main target tissue and their potential mechanism by combing pathological section and qRT-PCR technologies. RESULTS: The transcriptome data of five visceral tissues revealed that acute heat stress broadly caused inflammatory response and damaged tissues metabolic homeostasis. Among the five tested visceral tissues, the number of differentially expressed genes in the lung was the highest, and their fold changes were the greatest, indicating that the lung was the main target tissue of acute heat stress. The results of pathological section revealed severe inflammation, emphysema and pulmonary hemorrhage in the lung under acute heat stress. Our study found that some pro-inflammatory genes, including CNTFR, FURIN, CCR6, LIFR and IL20RA, were significantly up-regulated both in the heat-stress and heat-death groups, and their fold changes in the heat-death group were significantly greater than that in the heat-stress group. We also found an anti-inflammatory gene, AvBD9, exhibiting an extremely high expression in the heat-stress group but a low expression in the heat-death group. CONCLUSIONS: Our study found that acute heat stress caused multiple tissue injuries broadly and the lung was the main target tissue of acute heat stress in chicken. Acute heat stress caused a severe inflammatory response, emphysema, and pulmonary haemorrhage, The severe inflammatory response in the heat-death group was related to the up-regulation of pro-inflammatory genes and down-regulation of anti-inflammatory genes.

In recent years, some common themes in the development of sex-specific traits in different animal lineages have started to emerge since the discovery of the Dmrt (doublesex-mab3-related transcription factor gene) genes. Bivalves are characterized by a diversity of sexual systems, including simultaneous hermaphroditism, sequential hermaphroditism, and strict gonochorism. However, to date, no research has focused on the genome-wide characterization and analysis of Dmrt genes in bivalves. In this study, the identification and analysis of Dmrt genes in 15 bivalves were performed using bioinformatics methods. A total of 55 Dmrt genes were retrieved in the studied bivalve genomes. The number of Dmrt genes in different species ranged from 3 to 5. The phylogenetic tree showed that Dmrt genes in bivalves can be subdivided into 5 classes: the Dmrt2-like class, Dmrt3-like class, Dmrt4/5-like class, Dsx-like class, and scallop-specific Dmrt class. The Ka/Ks ratios suggested that all Dmrt classes underwent purifying selection pressure. Furthermore, the spatiotemporal expression of Dmrt genes in four bivalve species suggested that different Dmrt genes may have different functions, and scallop-specific Dmrt genes may play a key role in sex determination/differentiation. In general, this study provides a molecular basis for in-depth examination of the functions of Dmrt genes and phylogenomic analyses in bivalves.

BACKGROUND: Lifespan extension has independently evolved several times during mammalian evolution, leading to the emergence of a group of long-lived animals. Though mammalian/mechanistic target of rapamycin (mTOR) signaling pathway is shown as a central regulator of lifespan and aging, the underlying influence of mTOR pathway on the evolution of lifespan in mammals is not well understood. RESULTS: Here, we performed evolution analyses of 72 genes involved in the mTOR network across 48 mammals to explore the underlying mechanism of lifespan extension. We identified a total of 20 genes with significant evolution signals unique to long-lived species, including 12 positively selected genes, four convergent evolution genes, and five longevity associated genes whose evolution rate related to the maximum lifespan (MLS). Of these genes, four positively selected genes, two convergent evolution genes and one longevity-associated gene were involved in the autophagy response and aging-related diseases, while eight genes were known as cancer genes, indicating the long-lived species might have evolved effective regulation mechanisms of autophagy and cancer to extend lifespan. CONCLUSION: Our study revealed genes with significant evolutionary signals unique to long-lived species, which provided new insight into the lifespan extension of mammals and might bring new strategies to extend human lifespan.

BACKGROUND: Hull colour is an important morphological marker for selection in seed production of foxtail millet. However, the molecular mechanisms underlying hull colour variation remain unknown. RESULTS: An F7 recombinant inbred line (RIL) population containing 215 lines derived from Hongjiugu × Yugu18 was used to analyze inheritance and detect the quantitative trait loci (QTL) for four hull colour traits using major gene plus polygene mixed inheritance analysis and composite interval mapping (CIM) in four environments. Genetic analysis revealed that the hull colour L* value (HCL*) was controlled by two major genes plus additive polygenes, the hull colour a* value (HCa*) was controlled by three major genes, the hull colour b* value (HCb*) was controlled by two major genes plus polygenes, and the hull colour C* value (HCC*) was controlled by four major genes. A high-density genetic linkage map covering 1227.383 cM of the foxtail millet genome, with an average interval of 0.879 cM between adjacent bin markers, was constructed using 1420 bin markers. Based on the genetic linkage map and the phenotypic data, a total of 39 QTL were detected for these four hull colour traits across four environments, each explaining 1.50%-49.20% of the phenotypic variation. Of these, six environmentally stable major QTL were co-localized to regions on chromosomes 1 and 9, playing a major role in hull colour. There were 556 annotated genes within the two QTL regions. Based on the functions of homologous genes in Arabidopsis and the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) gene annotations, five genes were predicted as candidate genes for further studies. CONCLUSIONS: This is the first study to use an inheritance model and QTL mapping to determine the genetic mechanisms of hull colour trait in foxtail millet. We identified six major environmentally stable QTL and predicted five potential candidate genes to be associated with hull colour. These results advance the current understanding of the genetic mechanisms underlying hull colour traits in foxtail millet and provide additional resources for application in genomics-assisted breeding and potential isolation and functional characterization of the candidate genes.

BACKGROUND: Trichoderma is a diverse genus of fungi that includes several species that possess biotechnological and agricultural applications, including the biocontrol of pathogenic fungi and nematodes. The mitochondrial genome of a putative strain of Trichoderma harzianum called PAR3 was analyzed after isolation from the roots of Scarlet Royal grapevine scion grafted to Freedom rootstock, located in a grapevine vineyard in Parlier, CA, USA. Here, we report the sequencing, comparative assembly, and annotation of the nuclear genome of PAR3 and confirm its identification as a strain of T. harzianum. We subsequently compared the genes found in T. harzianum PAR3 to other known T. harzianum strains. Assembly of Illumina and/or Oxford Nanopore reads by the popular long-read assemblers, Flye and Canu, and the hybrid assemblers, SPAdes and MaSuRCA, was performed and the quality of the resulting assemblies were compared to ascertain which assembler generated the highest quality draft genome assembly. RESULTS: MaSuRCA produced the most complete and high-fidelity assembly yielding a nuclear genome of 40.7 Mb comprised of 112 scaffolds. Subsequent annotation of this assembly produced 12,074 gene models and 210 tRNAs. This included 221 genes that did not have equivalent genes in other T. harzainum strains. Phylogenetic analysis of ITS, rpb2, and tef1a sequences from PAR3 and established Trichoderma spp. showed that all three sequences from PAR3 possessed more than 99% identity to those of Trichoderma harzianum, confirming that PAR3 is an isolate of Trichoderma harzianum. We also found that comparison of gene models between T. harzianum PAR3 and other T. harzianum strains resulted in the identification of significant differences in gene type and number, with 221 unique genes identified in the PAR3 strain. CONCLUSIONS: This study gives insight into the efficacy of several popular assembly platforms for assembly of fungal nuclear genomes, and found that the hybrid assembler, MaSuRCA, was the most effective program for genome assembly. The annotated draft nuclear genome and the identification of genes not found in other T. harzainum strains could be used to investigate the potential applications of T. harzianum PAR3 for biocontrol of grapevine fungal canker pathogens and as source of anti-microbial compounds.

BACKGROUND: The 1RS arm of wheat-rye 1BL.1RS translocations contains several subtelomeric tandem repeat families. To study the effect of the difference in the composition of these tandem repeats on the meiotic recombination of 1RS arms can help to enrich the genetic diversity of 1BL.1RS translocation chromosomes. RESULTS: Five wheat-rye 1BL.1RS translocation cultivars/lines were used to build two cross combinations including group 1 (20T401 × Zhou 8425B, 20T401 × Lovrin 10 and 20T401 × Chuannong 17) and group 2 (20T360-2 × Zhou 8425B, 20T360-2 × Lovrin 10 and 20T360-2 × Chuannong 17). Oligonucleotide (oligo) probes Oligo-s120.3, Oligo-TR72, and Oligo-119.2-2 produced the same signal pattern on the 1RS arms in lines 20T401 and 20T360-2, and another signal pattern in the three cultivars Zhou 8425B, Lovrin 10 and Chuannong 17. The Oligo-pSc200 signal disappeared from the 1RS arms of the line 20T401, and the signal intensity of this probe on the 1RS arms of the line 20T360-2 was weaker than that of the three cultivars. The five cultivars/lines had the same signal pattern of the probe Oligo-pSc250. The recombination rate of 1RS arms in group 1 was significantly lower than that in group 2. In the progenies from group 1, unequal meiotic recombination in the subtelomeric pSc119.2 and pSc250 tandem repeat regions, and a 1BL.1RS with inversion of 1RS segment between the pSc200 and the nucleolar organizer region were found. CONCLUSIONS: This study provides a visual tool to detect the meiotic recombination of 1RS arms. The meiotic recombination rate of 1RS arms was affected by the variation of pSc200 tandem repeat, indicating the similar composition of subtelomeric tandem repeats on these arms could increase their recombination rate. These results indicate that the 1RS subtelomeric structure will affect its recombination, and thus the localization of genes on 1RS by means of meiotic recombination might also be affected.

BACKGROUND: Rheumatoid arthritis (RA) is often accompanied by a common extra-articular manifestation known as RA-related usual interstitial pneumonia (RA-UIP), which is associated with a poor prognosis. However, the mechanism remains unclear. To identify potential mechanisms, we conducted bioinformatics analysis based on high-throughput sequencing of the Gene Expression Omnibus (GEO) database. RESULTS: Weighted gene co-expression network analysis (WGCNA) analysis identified 2 RA-positive related modules and 4 idiopathic pulmonary fibrosis (IPF)-positive related modules. A total of 553 overlapped differentially expressed genes (DEG) were obtained, of which 144 in the above modules were further analyzed. The biological process of "oxidative phosphorylation" was found to be the most relevant with both RA and IPF. Additionally, 498 up-regulated genes in lung tissues of RA-UIP were screened out and enriched by 7 clusters, of which 3 were closely related to immune regulation. The analysis of immune infiltration showed a characteristic distribution of peripheral immune cells in RA-UIP, compared with IPF-UIP in lung tissues. CONCLUSIONS: These results describe the complex molecular and functional landscape of RA-UIP, which will help illustrate the molecular pathological mechanism of RA-UIP and identify new biomarkers and therapeutic targets for RA-UIP in the future.

BACKGROUND: For Asian seabass (Lates calcarifer, Bloch 1790) cultured at sea cages various aquatic pathogens, complex environmental and stress factors are considered as leading causes of disease, causing tens of millions of dollars of annual economic losses. Over the years, we conducted farm-based challenges by exposing Asian seabass juveniles to complex natural environmental conditions. In one of these challenges, we collected a total of 1,250 fish classified as either 'sensitive' or 'robust' individuals during the 28-day observation period. RESULTS: We constructed a high-resolution linkage map with 3,089 SNPs for Asian seabass using the double digest Restriction-site Associated DNA (ddRAD) technology and a performed a search for Quantitative Trait Loci (QTL) associated with robustness. The search detected a major genome-wide significant QTL for increased robustness in pathogen-infected marine environment on linkage group 11 (ASB_LG11; 88.9 cM to 93.6 cM) with phenotypic variation explained of 81.0%. The QTL was positioned within a > 800 kb genomic region located at the tip of chromosome ASB_LG11 with two Single Nucleotide Polymorphism markers, R1-38468 and R1-61252, located near to the two ends of the QTL. When the R1-61252 marker was validated experimentally in a different mass cross population, it showed a statistically significant association with increased robustness. The majority of thirty-six potential candidate genes located within the QTL have known functions related to innate immunity, stress response or disease. By utilizing this ddRAD-based map, we detected five mis-assemblies corresponding to four chromosomes, namely ASB_LG8, ASB_LG9, ASB_LG15 and ASB_LG20, in the current Asian seabass reference genome assembly. CONCLUSION: According to our knowledge, the QTL associated with increased robustness is the first such finding from a tropical fish species. Depending on further validation in other stocks and populations, it might be potentially useful for selecting robust Asian seabass lines in selection programs.

BACKGROUND: The Pacific oyster, Crassostrea gigas, is an economically important shellfish around the world. Great efforts have been made to improve its growth rate through genetic breeding. However, the candidate marker genes, pathways, and potential lncRNAs involved in oyster growth regulation remain largely unknown. To identify genes, lncRNAs, and pathways involved in growth regulation, C. gigas spat was cultured at a low temperature (15 ℃) to yield a growth-inhibited model, which was used to conduct comparative transcriptome analysis with spat cultured at normal temperature (25 ℃). RESULTS: In total, 8627 differentially expressed genes (DEGs) and 1072 differentially expressed lncRNAs (DELs) were identified between the normal-growth oysters (cultured at 25 ℃, hereinafter referred to as NG) and slow-growth oysters (cultured at 15 ℃, hereinafter referred to as SG). Functional enrichment analysis showed that these DEGs were mostly enriched in the AMPK signaling pathway, MAPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process. LncRNAs analysis identified 265 cis-acting pairs and 618 trans-acting pairs that might participate in oyster growth regulation. The expression levels of LNC_001270, LNC_003322, LNC_011563, LNC_006260, and LNC_012905 were inducible to the culture temperature and food abundance. These lncRNAs were located at the antisense, upstream, or downstream of the SREBP1/p62, CDC42, CaM, FAS, and PIK3CA genes, respectively. Furthermore, the expression of the trans-acting lncRNAs, including XR_9000022.2, LNC_008019, LNC_015817, LNC_000838, LNC_00839, LNC_011859, LNC_007294, LNC_006429, XR_002198885.1, and XR_902224.2 was also significantly associated with the expression of genes enriched in AMPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process. CONCLUSIONS: In this study, we identified the critical growth-related genes and lncRNAs that could be utilized as candidate markers to illustrate the molecular mechanisms underlying the growth regulation of Pacific oysters.

Recent developments in plant genomics have enabled a comprehensive analysis of the medicinal potential of plants based on their gene repertoire. Genes of biosynthesis pathways can be discovered through comparative genomics and through integration of transcriptomic data. Data-driven discovery of specialized metabolites could accelerate research.

The tumor immune microenvironment (TIME) of colon cancer (CC) has been associated with extensive immune cell infiltration (IMI). Increasing evidence demonstrated that plasma cells (PC) have an extremely important role in advance of antitumor immunity. Nonetheless, there is a lack of comprehensive analyses of PC infiltration in clinical prognosis and immunotherapy in CC. This study systematically addressed the gene expression model and clinical information of CC patients. Clinical samples were obtained from the TCGA (The Cancer Genome Atlas) databases. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), GSVA, and the MAlignant Tumors using Expression data (ESTIMATE) algorithm were employed to research the potential mechanism and pathways. Immunophenoscore (IPS) was obtained to evaluate the immunotherapeutic significance of risk score. Half maximal inhibitory concentration (IC50) of chemotherapeutic medicine was predicted by employing the pRRophetic algorithm. A total of 513 CC samples (including 472 tumor samples and 41 normal samples) were collected from the TCGA-GDC database. Significant black modules and 313 candidate genes were considered PC-related genes by accessing WGCNA. Five pivotal genes were established through multiple analyses, which revealed excellent prognostic. The underlying correlation between risk score with tumor mutation burden (TMB) was further explored. In addition, the risk score was obviously correlated with various tumor immune microenvironment (TIME). Also, risk CC samples showed various signaling pathways activity and different pivotal sensitivities to administering chemotherapy. Finally, the biological roles of the CD177 gene were uncovered in CC.

BACKGROUND: The adaptive evolution of plateau indigenous animals is a current research focus. However, phenotypic adaptation is complex and may involve the interactions between multiple genes or pathways, many of which remain unclear. As a kind of livestock with important economic value, cashmere goat has a high ability of plateau adaptation, which provides us with good materials for studying the molecular regulation mechanism of animal plateau adaptation. RESULTS: In this study, 32 Jiangnan (J) and 32 Tibetan (T) cashmere goats were sequenced at an average of 10. Phylogenetic, population structure, and linkage disequilibrium analyses showed that natural selection or domestication has resulted in obvious differences in genome structure between the two breeds. Subsequently, 553 J vs. T and 608 T vs. J potential selected genes (PSGs) were screened. These PSGs showed potential relationships with various phenotypes, including myocardial development and activity (LOC106502520, ATP2A2, LOC102181869, LOC106502520, MYL2, ISL1, and LOC102181869 genes), pigmentation (MITF and KITLG genes), hair follicles/hair growth (YAP1, POGLUT1, AAK1, HES1, WNT1, PRKAA1, TNKS, WNT5A, VAX2, RSPO4, CSNK1G1, PHLPP2, CHRM2, PDGFRB, PRKAA1, MAP2K1, IRS1, LPAR1, PTEN, PRLR, IBSP, CCNE2, CHAD, ITGB7, TEK, JAK2, and FGF21 genes), and carcinogenesis (UBE2R2, PIGU, DIABLO, NOL4L, STK3, MAP4, ADGRG1, CDC25A, DSG3, LEPR, PRKAA1, IKBKB, and ABCG2 genes). Phenotypic analysis showed that Tibetan cashmere goats has finer cashmere than Jiangnan cashmere goats, which may allow cashmere goats to better adapt to the cold environment in the Tibetan plateau. Meanwhile, KRTs and KAPs expression in Jiangnan cashmere goat skin was significantly lower than in Tibetan cashmere goat. CONCLUSIONS: The mutations in these PSGs maybe closely related to the plateau adaptation ability of cashmere goats. In addition, the expression differences of KRTs and KAPs may directly determine phenotypic differences in cashmere fineness between the two breeds. In conclusion, this study provide a reference for further studying plateau adaptive mechanism in animals and goat breeding.