BACKGROUND: The influence of linkage disequilibrium (LD), epistasis, and inbreeding on genotypic variance continues to be an important area of investigation in genetics and evolution. Although the current knowledge about biological pathways and gene networks indicates that epistasis is important in determining quantitative traits, the empirical evidence for a range of species and traits is that the genotypic variance is most additive. This has been confirmed by some recent theoretical studies. However, because these investigations assumed linkage equilibrium, considered only additive effects, or used simplified assumptions for two- and higher-order epistatic effects, the objective of this investigation was to provide additional information about the impact of LD and epistasis on genetic variances in noninbred and inbred populations, using a simulated dataset. RESULTS: In general, the most important component of the genotypic variance was additive variance. Because of positive LD values, after 10 generations of random crosses there was generally a decrease in all genetic variances and covariances, especially the nonepistatic variances. Thus, the epistatic variance/genotypic variance ratio is inversely proportional to the LD level. Increasing inbreeding increased the magnitude of the additive, additive x additive, additive x dominance, and dominance x additive variances, and decreased the dominance and dominance x dominance variances. Except for duplicate epistasis with 100% interacting genes, the epistatic variance/genotypic variance ratio was proportional to the inbreeding level. In general, the additive x additive variance was the most important component of the epistatic variance. Concerning the genetic covariances, in general, they showed lower magnitudes relative to the genetic variances and positive and negative signs. The epistatic variance/genotypic variance ratio was maximized under duplicate and dominant epistasis and minimized assuming recessive and complementary epistasis. Increasing the percentage of epistatic genes from 30 to 100% increased the epistatic variance/genotypic variance ratio by a rate of 1.3 to 12.6, especially in inbred populations. The epistatic variance/genotypic variance ratio was maximized in the noninbred and inbred populations with intermediate LD and an average allelic frequency of the dominant genes of 0.3 and in the noninbred and inbred populations with low LD and an average allelic frequency of 0.5. CONCLUSIONS: Additive variance is in general the most important component of genotypic variance. LD and inbreeding have a significant effect on the magnitude of the genetic variances and covariances. In general, the additive x additive variance is the most important component of epistatic variance. The maximization of the epistatic variance/genotypic variance ratio depends on the LD level, degree of inbreeding, epistasis type, percentage of interacting genes, and average allelic frequency.

BACKGROUND: Hi-C and its high nucleosome resolution variant Micro-C provide a window into the spatial packing of a genome in 3D within the cell. Even though both techniques do not directly depend on the binding of specific antibodies, previous work has revealed enriched interactions and domain structures around multiple chromatin marks; epigenetic modifications and transcription factor binding sites. However, the joint impact of chromatin marks in Hi-C and Micro-C interactions have not been globally characterized, which limits our understanding of 3D genome characteristics. An emerging question is whether it is possible to deduce 3D genome characteristics and interactions by integrative analysis of multiple chromatin marks and associate interactions to functionality of the interacting loci. RESULT: We come up with a probabilistic method PROBC to decompose Hi-C and Micro-C interactions by known chromatin marks. PROBC is based on convex likelihood optimization, which can directly take into account both interaction existence and nonexistence. Through PROBC, we discover histone modifications (H3K27ac, H3K9me3, H3K4me3, H3K4me1) and CTCF as particularly predictive of Hi-C and Micro-C contacts across cell types and species. Moreover, histone modifications are more effective than transcription factor binding sites in explaining the genome's 3D shape through these interactions. PROBC can successfully predict Hi-C and Micro-C interactions in given species, while it is trained on different cell types or species. For instance, it can predict missing nucleosome resolution Micro-C interactions in human ES cells trained on mouse ES cells only from these 5 chromatin marks with above 0.75 AUC. Additionally, PROBC outperforms the existing methods in predicting interactions across almost all chromosomes. CONCLUSION: Via our proposed method, we optimally decompose Hi-C interactions in terms of these chromatin marks at genome and chromosome levels. We find a subset of histone modifications and transcription factor binding sites to be predictive of both Hi-C and Micro-C interactions and TADs across human, mouse, and different cell types. Through learned models, we can predict interactions on species just from chromatin marks for which Hi-C data may be limited.

BACKGROUND: Spot form net blotch (SFNB) caused by the necrotrophic fungal pathogen Pyrenophora teres f. maculata (Ptm) is an economically important disease of barley that also infects wheat. Using genetic analysis to characterize loci in Ptm genomes associated with virulence or avirulence is an important step to identify pathogen effectors that determine compatible (virulent) or incompatible (avirulent) interactions with cereal hosts. Association mapping (AM) is a powerful tool for detecting virulence loci utilizing phenotyping and genotyping data generated for natural populations of plant pathogenic fungi. RESULTS: Restriction-site associated DNA genotyping-by-sequencing (RAD-GBS) was used to generate 4,836 single nucleotide polymorphism (SNP) markers for a natural population of 103 Ptm isolates collected from Idaho, Montana and North Dakota. Association mapping analyses were performed utilizing the genotyping and infection type data generated for each isolate when challenged on barley seedlings of thirty SFNB differential barley lines. A total of 39 marker trait associations (MTAs) were detected across the 20 barley lines corresponding to 30 quantitative trait loci (QTL); 26 novel QTL and four that were previously mapped in Ptm biparental populations. These results using diverse US isolates and barley lines showed numerous barley-Ptm genetic interactions with seven of the 30 Ptm virulence/avirulence loci falling on chromosome 3, suggesting that it is a reservoir of diverse virulence effectors. One of the loci exhibited reciprocal virulence/avirulence with one haplotype predominantly present in isolates collected from Idaho increasing virulence on barley line MXB468 and the alternative haplotype predominantly present in isolates collected from North Dakota and Montana increasing virulence on barley line CI9819. CONCLUSIONS: Association mapping provided novel insight into the host pathogen genetic interactions occurring in the barley-Ptm pathosystem. The analysis suggests that chromosome 3 of Ptm serves as an effector reservoir in concordance with previous reports for Pyrenophora teres f. teres, the causal agent of the closely related disease net form net blotch. Additionally, these analyses identified the first reported case of a reciprocal pathogen virulence locus. However, further investigation of the pathosystem is required to determine if multiple genes or alleles of the same gene are responsible for this genetic phenomenon.

BACKGROUND: Disclosure of patients' genetic information in the process of applying machine learning techniques for tumor classification hinders the privacy of personal information. Homomorphic Encryption (HE), which supports operations between encrypted data, can be used as one of the tools to perform such computation without information leakage, but it brings great challenges for directly applying general machine learning algorithms due to the limitations of operations supported by HE. In particular, non-polynomial activation functions, including softmax functions, are difficult to implement with HE and require a suitable approximation method to minimize the loss of accuracy. In the secure genome analysis competition called iDASH 2020, it is presented as a competition task that a multi-label tumor classification method that predicts the class of samples based on genetic information using HE. METHODS: We develop a secure multi-label tumor classification method using HE to ensure privacy during all the computations of the model inference process. Our solution is based on a 1-layer neural network with the softmax activation function model and uses the approximate HE scheme. We present an approximation method that enables softmax activation in the model using HE and a technique for efficiently encoding data to reduce computational costs. In addition, we propose a HE-friendly data filtering method to reduce the size of large-scale genetic data. RESULTS: We aim to analyze the dataset from The Cancer Genome Atlas (TCGA) dataset, which consists of 3,622 samples from 11 types of cancers, genetic features from 25,128 genes. Our preprocessing method reduces the number of genes to 4,096 or less and achieves a microAUC value of 0.9882 (85% accuracy) with a 1-layer shallow neural network. Using our model, we successfully compute the tumor classification inference steps on the encrypted test data in 3.75 minutes. As a result of exceptionally high microAUC values, our solution was awarded co-first place in iDASH 2020 Track 1: "Secure multi-label Tumor classification using Homomorphic Encryption". CONCLUSIONS: Our solution is the first result of implementing a neural network model with softmax activation using HE. Also, HE optimization methods presented in this work enable machine learning implementation using HE or other challenging HE applications.

BACKGROUND: Egg production is one of the most important economic traits in the poultry industry. The hypothalamic-pituitary-gonadal (HPG) axis plays an essential role in regulating reproductive activities. However, the key genes and regulatory pathways within the HPG axis dominating egg production performance remain largely unknown in ducks. RESULTS: In this study, we compared the transcriptomic profiles of the HPG-related tissues between ducks with high egg production (HEP) and low egg production (LEP) to reveal candidate genes and regulatory pathways dominating egg production. We identified 543, 759, 670, and 181 differentially expressed genes (DEGs) in the hypothalamus, pituitary, ovary stroma, and F5 follicle membrane, respectively. Gene Ontology (GO) analysis revealed that DEGs from four HPG axis-related tissues were enriched in the "cellular component" category. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that the neuroactive ligand-receptor interaction pathway was significantly enriched based on DEGs commonly identified in all four HPG axis-related tissues. Gene expression profiles and Protein-Protein Interaction (PPI) network were performed to show the regulatory relationships of the DEGs identified. Five DEGs encoding secreted proteins in the hypothalamus and pituitary have interaction with DEGs encoding targeted proteins in the ovary stroma and F5 follicle membrane, implying that they were these DEGs might play similar roles in the regulation of egg production. CONCLUSIONS: Our results revealed that neuroactive ligand-receptor interaction pathway and five key genes(VEGFC, SPARC, BMP2, THBS1, and ADAMTS15) were identified as the key signaling pathways and candidate genes within the HPG axis responsible for different egg production performance between HEP and LEP. This is the first study comparing the transcriptomic profiles of all HPG axis-related tissues in HEP and LEP using RNA-seq in ducks to the best of our knowledge. These data are helpful to enrich our understanding of the classical HPG axis regulating the egg production performance and identify candidate genes that can be used for genetic selection in ducks.

BACKGROUND: The mitochondrial (mt) genomes of 15 species of sucking lice from seven families have been studied to date. These louse species have highly dynamic, fragmented mt genomes that differ in the number of minichromosomes, the gene content, and gene order in a minichromosome between families and even between species of the same genus. RESULTS: In the present study, we analyzed the publicly available data to understand mt genome fragmentation in seal lice (family Echinophthiriidae) and gorilla louse, Pthirus gorillae (family Pthiridae), in particular the role of minichromosome split and minichromosome merger in the evolution of fragmented mt genomes. We show that 1) at least three ancestral mt minichromosomes of sucking lice have split in the lineage leading to seal lice, 2) one minichromosome ancestral to primate lice has split in the lineage to the gorilla louse, and 3) two ancestral minichromosomes of seal lice have merged in the lineage to the northern fur seal louse. Minichromosome split occurred 15-16 times in total in the lineages leading to species in six families of sucking lice investigated. In contrast, minichromosome merger occurred only four times in the lineages leading to species in three families of sucking lice. Further, three ancestral mt minichromosomes of sucking lice have split multiple times independently in different lineages of sucking lice. Our analyses of mt karyotypes and gene sequences also indicate the possibility of a host switch of crabeater seal louse to Weddell seals. CONCLUSIONS: We conclude that: 1) minichromosome split contributes more than minichromosome merger in mt genome fragmentation of sucking lice, and 2) mt karyotype comparison helps understand the phylogenetic relationships between sucking louse species.

BACKGROUND: Bees are the most important group of pollinators worldwide and their populations are declining. In natural conditions, Apis mellifera depends exclusively on food from the field to meet its physiological demands. In the period of scarcity, available resources are insufficient and artificial supplementation becomes essential for maintaining the levels of vitamins, proteins, carbohydrates, and minerals of colonies. Among these minerals, zinc is essential in all living systems, particularly for the regulation of cell division and protein synthesis, and is a component of more than 200 metalloenzymes. RESULTS: The total RNA extracted from the brain tissue of nurse bees exposed to different sources and concentrations of zinc was sequenced. A total of 1,172 genes in the treatment that received an inorganic source of zinc and 502 genes that received an organic source of zinc were found to be differentially expressed among the control group. Gene ontology enrichment showed that zinc can modulate important biological processes such as nutrient metabolism and the molting process. CONCLUSIONS: Our results indicate that zinc supplementation modulates the expression of many differentially expressed genes and plays an important role in the development of Apis mellifera bees. All the information obtained in this study can contribute to future research in the field of bee nutrigenomics.

BACKGROUND: Transcription factors (TFs) play central roles in regulating gene expression. With the rapid growth in the use of high-throughput sequencing methods, there is a need to develop a comprehensive data processing and analyzing framework for inferring influential TFs based on ChIP-seq/ATAC-seq datasets. RESULTS: Here, we introduce FindIT2 (Find Influential TFs and Targets), an R/Bioconductor package for annotating and processing high-throughput multi-omics data. FindIT2 supports a complete framework for annotating ChIP-seq/ATAC-seq peaks, identifying TF targets by the combination of ChIP-seq and RNA-seq datasets, and inferring influential TFs based on different types of data input. Moreover, benefited from the annotation framework based on Bioconductor, FindIT2 can be applied to any species with genomic annotations, which is particularly useful for the non-model species that are less well-studied. CONCLUSION: FindIT2 provides a user-friendly and flexible framework to generate results at different levels according to the richness of the annotation information of user's species. FindIT2 is compatible with all the operating systems and is released under Artistic-2.0 License. The source code and documents are freely available through Bioconductor ( https://bioconductor.org/packages/devel/bioc/html/FindIT2.html ).

BACKGROUND: Rice sheath blight, caused by Rhizoctonia solani Kühn (teleomorph: Thanatephorus cucumeris), is one of the most severe diseases in rice (Oryza sativa L.) worldwide. Studies on resistance genes and resistance mechanisms of rice sheath blight have mainly focused on indica rice. Rice sheath blight is a growing threat to rice production with the increasing planting area of japonica rice in Northeast China, and it is therefore essential to explore the mechanism of sheath blight resistance in this rice subspecies. RESULTS: In this study, RNA-seq technology was used to analyse the gene expression changes of leaf sheath at 12, 24, 36, 48, and 72 h after inoculation of the resistant cultivar 'Shennong 9819' and susceptible cultivar 'Koshihikari' with R. solani. In the early stage of R. solani infection of rice leaf sheaths, the number of differentially expressed genes (DEGs) in the inoculated leaf sheaths of resistant and susceptible cultivars showed different regularity. After inoculation, the number of DEGs in the resistant cultivar fluctuated, while the number of DEGs in the susceptible cultivar increased first and then decreased. In addition, the number of DEGs in the susceptible cultivar was always higher than that in the resistant cultivar. After inoculation with R. solani, the overall transcriptome changes corresponding to multiple biological processes, molecular functions, and cell components were observed in both resistant and susceptible cultivars. These included metabolic process, stimulus response, biological regulation, catalytic activity, binding and membrane, and they were differentially regulated. The phenylalanine metabolic pathway; tropane, piperidine, and pyridine alkaloid biosynthesis pathways; and plant hormone signal transduction were significantly enriched in the early stage of inoculation of the resistant cultivar Shennong 9819, but not in the susceptible cultivar Koshihikari. This indicates that the response of the resistant cultivar Shennong 9819 to pathogen stress was faster than that of the susceptible cultivar. The expression of plant defense response marker PR1b gene, transcription factor OsWRKY30 and OsPAL1 and OsPAL6 genes that induce plant resistance were upregulated in the resistant cultivar. These data suggest that in the early stage of rice infection by R. solani, there is a pathogen-induced defence system in resistant rice cultivars, involving the expression of PR genes, key transcription factors, PAL genes, and the enrichment of defence-related pathways. CONCLUSION: The transcriptome data revealed the molecular and biochemical differences between resistant and susceptible cultivars of rice after inoculation with R. solani, indicating that resistant cultivars have an immune response mechanism in the early stage of pathogen infection. Disease resistance is related to the overexpression of PR genes, key transcriptome factors, and PAL genes, which are potential targets for crop improvement.

BACKGROUND: The grass carp has great economic value and occupies an important evolutionary position. Genomic information regarding this species could help better understand its rapid growth rate as well as its unique body plan and environmental adaptation. RESULTS: We assembled the chromosome-level grass carp genome using the PacBio sequencing and chromosome structure capture technique. The final genome assembly has a total length of 893.2 Mb with a contig N50 of 19.3 Mb and a scaffold N50 of 35.7 Mb. About 99.85% of the assembled contigs were anchored into 24 chromosomes. Based on the prediction, this genome contained 30,342 protein-coding genes and 43.26% repetitive sequences. Furthermore, we determined that the large genome size can be attributed to the DNA-mediated transposable elements which accounted for 58.9% of the repetitive sequences in grass carp. We identified that the grass carp has only 24 pairs of chromosomes due to the fusion of two ancestral chromosomes. Enrichment analyses of significantly expanded and positively selected genes reflected evolutionary adaptation of grass carp to the feeding habits. We also detected the loss of conserved non-coding regulatory elements associated with the development of the immune system, nervous system, and digestive system, which may be critical for grass carp herbivorous traits. CONCLUSIONS: The high-quality reference genome reported here provides a valuable resource for the genetic improvement and molecular-guided breeding of the grass carp.

BACKGROUND: Salivary gland (SMG) degeneration and dysfunction are common symptoms that occur after sex hormone deprivation, but the underlying mechanisms remain largely unknown. Additionally, immunocastration, which causes drop of sex hormones, has been developed as an alternative to surgical castration, however whether it exerts similar effects as surgical castration on the salivary glands is unknown. Through histological and RNA-seq analysis, we assessed changes in morphology and transcriptome of SMG in response to immunocastration (IM) versus surgical castration (bilateral orchiectomy, ORC). RESULTS: Compared to entire males (EM), ORC caused severe degeneration of SMG in rats, as evidenced by both decreased (P < 0.01) SMG weight and organ index, and by decreased (P < 0.01) quantity of SMG acini and ducts. IM had minimal effects (P > 0.05) on SMG weight and organ index, but it still caused degeneration (P < 0.05) of the acini and ducts. Even though, the quantity of both SMG acini and ducts was much higher (P < 0.001) in IM than in ORC. Functional enrichment analysis of the common regulated genes by ORC/IM revealed disrupted epithelial cell development, angiogenesis, anatomical structure morphogenesis and enhanced cell death are associated with SMG degeneration in deprivation of androgens. Integrated data analysis shown that there existed a selective hyperfunction of SMG ribosome and mitochondrion in ORC but not in IM, which might be associated with more severe degeneration of SMG in ORC than in IM. CONCLUSIONS: Our findings suggested that both surgical castration and immunocastration caused SMG degeneration by disrupting epithelial cell development, angiogenesis, anatomical structure morphogenesis and enhancing cell death. But, surgical castration selectively induced hyperfunction of SMG ribosome and mitochondrion, thus causing more severe degeneration of SMG than immunocastration.

BACKGROUND: Shiga toxin-producing Escherichia coli (STEC) O157:H7 is a pathogen known to reside in cattle feedlots. This retrospective study examined 181 STEC O157:H7 strains collected over 23 years from a closed-system feedlot. All strains were subjected to short-read sequencing, with a subset of 36 also subjected to long-read sequencing. RESULTS: Over 96% of the strains fell into four phylogenetically distinct clades. Clade membership was associated with multiple factors including stx composition and the alleles of a well-characterized polymorphism (tir 255 T > A). Small plasmids (2.7 to 40 kb) were found to be primarily clade specific. Within each clade, chromosomal rearrangements were observed along with a core phageome and clade specific phages. Across both core and mobile elements of the genome, multiple SNP alleles were in complete linkage disequilibrium across all strains within specific clades. Clade evolutionary rates varied between 0.9 and 2.8 SNP/genome/year with two tir A allele clades having the lowest evolutionary rates. Investigation into possible causes of the differing rates was not conclusive but revealed a synonymous based mutation in the DNA polymerase III of the fastest evolving clade. Phylogenetic trees generated through our bioinformatic pipeline versus the NCBI's pathogen detection project were similar, with the two tir A allele clades matching individual NCBI SNP clusters, and the two tir T allele clades assigned to multiple closely-related SNP clusters. CONCLUSIONS: In one ecological niche, a diverse STEC O157:H7 population exhibited different rates of evolution that associated with SNP alleles in linkage disequilibrium in the core genome and mobile elements, including tir 255 T > A.

BACKGROUND: Immortalized cell lines are widely used model systems whose genomes are often highly rearranged and polyploid. However, their genome structure is seldom deciphered and is thus not accounted for during analyses. We therefore used linked short- and long-read sequencing to perform haplotype-level reconstruction of the genome of a Drosophila melanogaster cell line (S2-DRSC) with a complex genome structure. RESULTS: Using a custom implementation (that is designed to use ultra-long reads in complex genomes with nested rearrangements) to call structural variants (SVs), we found that the most common SV was repetitive sequence insertion or deletion (> 80% of SVs), with Gypsy retrotransposon insertions dominating. The second most common SV was local sequence duplication. SNPs and other SVs were rarer, but several large chromosomal translocations and mitochondrial genome insertions were observed. Haplotypes were highly similar at the nucleotide level but structurally very different. Insertion SVs existed at various haplotype frequencies and were unlinked on chromosomes, demonstrating that haplotypes have different structures and suggesting the existence of a mechanism that allows SVs to propagate across haplotypes. Finally, using public short-read data, we found that transposable element insertions and local duplications are common in other D. melanogaster cell lines. CONCLUSIONS: The S2-DRSC cell line evolved through retrotransposon activity and vast local sequence duplications, that we hypothesize were the products of DNA re-replication events. Additionally, mutations can propagate across haplotypes (possibly explained by mitotic recombination), which enables fine-tuning of mutational impact and prevents accumulation of deleterious events, an inherent problem of clonal reproduction. We conclude that traditional linear homozygous genome representation conceals the complexity when dealing with rearranged and heterozygous clonal cells.

BACKGROUND: Global climate change together with growing desertification is leading to increased dust emissions to the atmosphere, drawing attention to possible impacts on marine ecosystems receiving dust deposition. Since microorganisms play important roles in maintaining marine homeostasis through nutrient cycling and carbon flow, detrimental changes in the composition of marine microbiota in response to increased dust input could negatively impact marine health, particularly so in seas located within the Global Dust Belt. Due to its strategic location between two deserts and unique characteristics, the Red Sea provides an attractive semi-enclosed "megacosm" to examine the impacts of large dust deposition on the vastly diverse microbiota in its exceptionally warm oligotrophic waters. RESULTS: We used culture-independent metagenomic approaches to assess temporal changes in the Red Sea microbiota in response to two severe sandstorms, one originated in the Nubian Desert in the summer 2016 and a second one originated in the Libyan Desert in the spring 2017. Despite differences in sandstorm origin and meteorological conditions, both sandstorms shifted bacterial and Archaeal groups in a similar mode. In particular, the relative abundance of autotrophic bacteria declined while those of heterotrophic bacteria, particularly Bacteroidetes, and Archaea increased. The changes peaked within six days from the start of sandstorms, and the community recovered the original assemblage within one month. CONCLUSION: Our results suggest that increased dust emission with expanding desertification could lead to undesirable impacts in ocean function, enhancing heterotrophic processes while reducing autotrophic ones, thereby affecting the marine food web in seas receiving dust deposition.

BACKGROUND: Helictotrichon virescens is a perennial grass that is primarily distributed in high altitude areas of 2000 ~ 4500 m. It is widely cultivated in the Qinghai-Tibet Plateau of China, strongly resistant to cold, and an essential part of the wild herbs in this region. However, the molecular mechanism of the response of H. virescens to low temperature stress and the key regulatory genes for specific biological processes are poorly understood. RESULTS: Physiological and transcriptome analyses were used to study the cold stress response mechanism in H virescens. During the low temperature stress period, the content of chlorophyll a and b decreased more and more with the delay of the treatment time. Among them, the difference between the controls was not significant, and the difference between the control and the treatment was significant. At the same time, the expression of related differential genes was up-regulated during low temperature treatment. In addition, the plant circadian pathway is crucial for their response to cold stress. The expression of differentially expressed genes that encode LHY and HY5 were strongly up-regulated during cold stress. CONCLUSIONS: This study should help to fully understand how H. virescens responds to low temperatures. It answers pertinent questions in the response of perennial herbs to cold stress, i.e., how light and low temperature signals integrate to regulate plant circadian rhythms and Decrease of content of chlorophylls (which can be also accompanied with decrease of total quantity of reaction centers) leads to an increase in photosynthetic damage.

BACKGROUND: Cutaneous squamous cell carcinoma (cSCC) is the second most common type of skin cancer, the prognosis for patients with metastatic cSCC remains relatively poor. Thus, there is an urgent need to identify new diagnostic, prognostic, and therapeutic targets and pathways in cSCC. RESULTS: It detected a total of 37,507 lncRNA probes and 32,825 mRNA probes and found 3593 differentially expressed lncRNAs and 3236 differentially expressed mRNAs. It has been found that mRNAs ACY3, NR1D1, MZB1 has co-expression relationship with six lncRNAs, GXYLT1P3, LINC00348, LOC101928131, A-33-p3340852, A-21-p0003442 and LOC644838. CONCLUSIONS: The aim of this study is to identify cSCC-specific lncRNAs and indicated that six unstudied lncRNAs may serve an important role in endoplasmic reticulum stress apoptosis, autophagy and the progression of cSCC by modulating ACY3, NR1D1 and MZB1.

BACKGROUND: Assisted Reproductive Technologies (ART) use can increase the risk of congenital overgrowth syndromes, such as large offspring syndrome (LOS) in ruminants. Epigenetic variations are known to influence gene expression and differentially methylated regions (DMRs) were previously determined to be associated with LOS in cattle. We observed DMRs overlapping tRNA clusters which could affect tRNA abundance and be associated with tissue specificity or overgrowth. Variations in tRNA expression have been identified in several disease pathways suggesting an important role in the regulation of biological processes. Understanding the role of tRNA expression in cattle offers an opportunity to reveal mechanisms of regulation at the translational level. We analyzed tRNA expression in the skeletal muscle and liver tissues of day 105 artificial insemination-conceived, ART-conceived with a normal body weight, and ART-conceived bovine fetuses with a body weight above the 97th percentile compared to Control-AI. RESULTS: Despite the centrality of tRNAs to translation, in silico predictions have revealed dramatic differences in the number of tRNA genes between humans and cattle (597 vs 1,659). Consistent with reports in human, only a fraction of predicted tRNA genes are expressed. We detected the expression of 474 and 487 bovine tRNA genes in the muscle and liver with the remainder being unexpressed. 193 and 198 unique tRNA sequences were expressed in all treatment groups within muscle and liver respectively. In addition, an average of 193 tRNA sequences were expressed within the same treatment group in different tissues. Some tRNA isodecoders were differentially expressed between treatment groups. In the skeletal muscle and liver, we categorized 11 tRNA isoacceptors with undetected expression as well as an isodecoder that was unexpressed in the liver (SerGGA). Our results identified variation in the proportion of tRNA gene copies expressed between tissues and differences in the highest contributing tRNA anticodon within an amino acid family due to treatment and tissue type. Out of all amino acid families, roughly half of the most highly expressed tRNA isoacceptors correlated to their most frequent codon in the bovine genome. CONCLUSION: Although the number of bovine tRNA genes is nearly triple of that of the tRNA genes in human, there is a shared occurrence of transcriptionally inactive tRNA genes in both species. We detected differential expression of tRNA genes as well as tissue- and treatment- specific tRNA transcripts with unique sequence variations that could modulate translation during protein homeostasis or cellular stress, and give rise to regulatory products targeting genes related to overgrowth in the skeletal muscle and/or tumor development in the liver of LOS individuals. While the absence of certain isodecoders may be relieved by wobble base pairing, missing tRNA species could increase the likelihood of mistranslation or mRNA degradation.

Microtubule associated proteins (MAPs), defined as proteins that bind microtubules but are not molecular motors or severing enzymes, play a key role in regulating microtubule stability in neurons. Existing studies of the evolutionary relationships between these proteins are limited to genomic data from a small number of species. We therefore used a large collection of publicly available reference-quality eukaryotic proteomes to carry out a phylogenetic analysis of microtubule associated proteins in both vertebrates and invertebrates. Complete or near-complete reference quality proteomes were obtained from Uniprot. Microtubule associated proteins were identified using InterProtScan, aligned using MUSCLE and then phylogenetic trees constructed using the WAG algorithm. We identified 889 proteins with tubulin binding domains, of which 663 were in eukaryotes, including 168 vertebrates and 64 invertebrates. The vertebrate proteins separated into three families, resembling human MAP 2, MAP4 and MAPT, respectively, while invertebrate MAPs clustered separately. We found significant variation in number of microtubule associated proteins and number of microtubule binding domains between taxa, with fish and mollusks having an unexpectedly high number of MAPs and binding domains, respectively. Our findings represent a novel analysis of the evolution of microtubule associated proteins based on publicly available proteomics data sets. We were able to confirm the phylogeny of MAPs identified based on more limited genomic analyses, and in addition, derived several novel insights on the structure and function of MAPs.

BACKGROUND: In biological systems, metabolomics can not only contribute to the discovery of metabolic signatures for disease diagnosis, but is very helpful to illustrate the underlying molecular disease-causing mechanism. Therefore, identification of disease-related metabolites is of great significance for comprehensively understanding the pathogenesis of diseases and improving clinical medicine. RESULTS: In the paper, we propose a disease and literature driven metabolism prediction model (DLMPM) to identify the potential associations between metabolites and diseases based on latent factor model. We build the disease glossary with disease terms from different databases and an association matrix based on the mapping between diseases and metabolites. The similarity of diseases and metabolites is used to complete the association matrix. Finally, we predict potential associations between metabolites and diseases based on the matrix decomposition method. In total, 1,406 direct associations between diseases and metabolites are found. There are 119,206 unknown associations between diseases and metabolites predicted with a coverage rate of 80.88%. Subsequently, we extract training sets and testing sets based on data increment from the database of disease-related metabolites and assess the performance of DLMPM on 19 diseases. As a result, DLMPM is proven to be successful in predicting potential metabolic signatures for human diseases with an average AUC value of 82.33%. CONCLUSION: In this paper, a computational model is proposed for exploring metabolite-disease pairs and has good performance in predicting potential metabolites related to diseases through adequate validation. The results show that DLMPM has a better performance in prioritizing candidate diseases-related metabolites compared with the previous methods and would be helpful for researchers to reveal more information about human diseases.

BACKGROUND: Of the > 2600 Salmonella serovars, Salmonella enterica serovar I 4,[5],12:i:- (serovar I 4,[5],12:i:-) has emerged as one of the most common causes of human salmonellosis and the most frequent multidrug-resistant (MDR; resistance to ≥3 antimicrobial classes) nontyphoidal Salmonella serovar in the U.S. Serovar I 4,[5],12:i:- isolates have been described globally with resistance to ampicillin, streptomycin, sulfisoxazole, and tetracycline (R-type ASSuT) and an integrative and conjugative element with multi-metal tolerance named Salmonella Genomic Island 4 (SGI-4). RESULTS: We analyzed 13,612 serovar I 4,[5],12:i:- strain sequences available in the NCBI Pathogen Detection database to determine global distribution, animal sources, presence of SGI-4, occurrence of R-type ASSuT, frequency of antimicrobial resistance (AMR), and potential transmission clusters. Genome sequences for serovar I 4,[5],12:i:- strains represented 30 countries from 5 continents (North America, Europe, Asia, Oceania, and South America), but sequences from the United States (59%) and the United Kingdom (28%) were dominant. The metal tolerance island SGI-4 and the R-type ASSuT were present in 71 and 55% of serovar I 4,[5],12:i:- strain sequences, respectively. Sixty-five percent of strain sequences were MDR which correlates to serovar I 4,[5],12:i:- being the most frequent MDR serovar. The distribution of serovar I 4,[5],12:i:- strain sequences in the NCBI Pathogen Detection database suggests that swine-associated strain sequences were the most frequent food-animal source and were significantly more likely to contain the metal tolerance island SGI-4 and genes for MDR compared to all other animal-associated isolate sequences. CONCLUSIONS: Our study illustrates how analysis of genomic sequences from the NCBI Pathogen Detection database can be utilized to identify the prevalence of genetic features such as antimicrobial resistance, metal tolerance, and virulence genes that may be responsible for the successful emergence of bacterial foodborne pathogens.

BACKGROUND: The growth and development of muscle stem cells (MuSCs) are significant events known to affect muscle plasticity, disease, meat production, and meat quality, which involves the types and functions of mRNA and non-coding RNA. Here, MuSCs were cultured from Guangxi fetal cattle. RNA sequencing was used to analyze the RNA expression of mRNA and non-coding RNAs during the cell proliferation and differentiation phases. RESULTS: Two thousand one hundred forty-eight mRNAs and 888 non-coding RNAs were differentially expressed between cell proliferation and differentiation phases, including 113 miRNAs, 662 lncRNAs, and 113 circRNAs. RT-qPCR verified the differential expression levels of mRNAs and non-coding RNAs, and the differentially expressed circUBE2Q2 was subsequently characterized. Expression profile analysis revealed that circUBE2Q2 was abundant in muscle tissues and intramuscular fat. The expression of cricUBE2Q2 was also significantly upregulated during MuSCs myogenic differentiation and SVFs adipogenic differentiation and decreased with age in cattle muscle tissue. Finally, the molecular mechanism of circUBE2Q2 regulating MuSCs function that affects skeletal muscle development was investigated. The results showed that circUBE2Q2 could serve as a sponge for miR-133a, significantly promoting differentiation and apoptosis of cultured MuSCs, and inhibiting proliferation of MuSCs. CONCLUSIONS: CircUBE2Q2 is associated with muscle growth and development and induces MuSCs myogenic differentiation through sponging miR-133a. This study will provide new clues for the mechanisms by which mRNAs and non-coding RNAs regulate skeletal muscle growth and development, affecting muscle quality and diseases.

BACKGROUND: The plant-specific IQ67-domain (IQD) gene family plays an important role in plant development and stress responses. However, little is known about the IQD family in common wheat (Triticum aestivum L), an agriculturally important crop that provides more than 20% of the calories and protein consumed in the modern human diet. RESULTS: We identified 125 IQDs in the wheat genome and divided them into four subgroups by phylogenetic analysis. The IQDs belonging to the same subgroup had similar exon-intron structure and conserved motif composition. Polyploidization contributed significantly to the expansion of IQD genes in wheat. Characterization of the expression profile of these genes revealed that a few T. aestivum (Ta)IQDs showed high tissue-specificity. The stress-induced expression pattern also revealed a potential role of TaIQDs in environmental adaptation, as TaIQD-2A-2, TaIQD-3A-9 and TaIQD-1A-7 were significantly induced by cold, drought and heat stresses, and could be candidates for future functional characterization. In addition, IQD genes in the A, B and D subgenomes displayed an asymmetric evolutionary pattern, as evidenced by their different gain or loss of member genes, expression levels and nucleotide diversity. CONCLUSIONS: This study elucidated the potential biological functions and evolutionary relationships of the IQD gene family in wheat and revealed the divergent fates of IQD genes during polyploidization.

BACKGROUND: Wound healing is a representative phenomenon of potato tubers subjected to mechanical injuries. Our previous results found that benzo-(1,2,3)-thiadiazole-7-carbothioic acid S-methyl ester (BTH) promoted the wound healing of potato tubers. However, the molecular mechanism related to inducible wound healing remains unknown. RESULTS: Transcriptomic evaluation of healing tissues from potato tubers at three stages, namely, 0 d (nonhealing), 5 d (wounded tubers healed for 5 d) and 5 d (BTH-treated tubers healed for 5 d) using RNA-Seq and differentially expressed genes (DEGs) analysis showed that more than 515 million high-quality reads were generated and a total of 7665 DEGs were enriched, and 16 of these DEGs were selected by qRT-PCR analysis to further confirm the RNA sequencing data. Gene ontology (GO) enrichment analysis indicated that the most highly DEGs were involved in metabolic and cellular processes, and KEGG enrichment analysis indicated that a large number of DEGs were associated with plant hormones, starch and sugar metabolism, fatty acid metabolism, phenylpropanoid biosynthesis and terpenoid skeleton biosynthesis. Furthermore, a few candidate transcription factors, including MYB, NAC and WRKY, and genes related to Ca2+-mediated signal transduction were also found to be differentially expressed during wound healing. Most of these enriched DEGs were upregulated after BTH treatment. CONCLUSION: This comparative expression profile provided useful resources for studies of the molecular mechanism via these promising candidates involved in natural or elicitor-induced wound healing in potato tubers.

BACKGROUND: Development of new pan-genome analysis tools is important, as the pangenome of a microbial species has become an important method to define the diversity of a selected taxon, most commonly a species, in the last years. This enables comparison of strains from different ecological niches and can be used to define the functional potential in a bacterial population. It gives us a much better view of microbial genomics than can be gained from singular genomes which after all are just single representatives of a much more varied population. RESULTS: We present Panakeia, a tool which strives to be easy to use and providing a detailed view of the pangenome structure which can efficiently be utilised for discovery, or further in-depth analysis, of features of interest. It analyses synteny and multiple structural patterns of the pangenome, giving insights into the biological diversity and evolution of the studied taxon. Panakeia hence provides both broad and detailed information on the structure of a pangenome, for diverse and highly clonal populations of bacteria. CONCLUSIONS: Previously published pangenome tools often reduce the information to a presence/absence matrix of unconnected genes or generate massive hard to interpret output graphs. However, Panakeia includes synteny and structural information and presents it in a way that can readily be used for further analysis. Panakeia can be downloaded at https://github.com/BioSina/Panakeia together with a detailed User Guide.

BACKGROUND: The characteristics of muscle fibers determine the growth and meat quality of poultry. In this study, we performed a weighted gene co-expression network analysis (WGCNA) on the muscle fiber characteristics and transcriptome profile of the breast muscle tissue of Gushi chicken at 6, 14, 22, and 30 weeks. RESULTS: A total of 27 coexpressed biological functional modules were identified, of which the midnight blue module had the strongest correlation with muscle fiber and diameter. In addition, 7 hub genes were found from the midnight blue module, including LC8 dynein light chain 2 (DYNLL2). Combined with miRNA transcriptome data, miR-148a-3p was found to be a potential target miRNA of DYNLL2. Experiments on chicken primary myoblasts (CPMs) demonstrated that miR-148a-3p promotes the expression of myosin heavy chain (MYHC) protein by targeting DYNLL2, proving that it can promote differentiation of myoblasts. CONCLUSIONS: This study proved that the hub gene DYNLL2 and its target miR-148-3p are important regulators in chicken myogenesis. These results provide novel insights for understanding the molecular regulation mechanisms related to the development of chicken breast muscle.

BACKGROUND: Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development. RESULTS: Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination. CONCLUSIONS: We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology.

BACKGROUND: Cool temperature egg storage prior to incubation is a common practice in the broiler industry; however, prolonged egg storage causes increased embryonic mortality and decreased hatchability and growth in surviving chicks. Exposing eggs to short periods of incubation during egg storage (SPIDES) reduces the adverse consequences of prolonged storage. SPIDES increases blastodermal cell viability by reducing apoptosis, though the counteracting mechanisms are unclear. To define the impact of prolonged storage and SPIDES, transcriptome analysis compared gene expression from blastoderms isolated from eggs exposed to the following treatments: control (CR, stored at 17 °C for 4 days), prolonged storage (NSR, stored at 17 °C for 21 days), SPIDES (SR, stored at 17 °C for 21 days with SPIDES), and incubated control (C2, stored at 17 °C for 4 days followed by incubation to HH (Hamburger-Hamilton) stage 2, used as the ideal standard development) (n = 3/group). Data analysis was performed using the CLC Genomics Workbench platform. Functional annotation was performed using DAVID and QIAGEN Ingenuity Pathway Analysis. RESULTS: In total, 4726 DEGs (differentially expressed genes) were identified across all experimental group comparisons (q < 0.05, FPKM> 20, |fold change| > 1.5). DEGs common across experimental comparisons were involved in cellular homeostasis and cytoskeletal protein binding. The NSR group exhibited activation of ubiquitination, apoptotic, and cell senescence processes. The SR group showed activation of cell viability, division, and metabolic processes. Through comparison analysis, cellular respiration, tRNA charging, cell cycle control, and HMBG1 signaling pathways were significantly impacted by treatment and potential regulatory roles for ribosomal protein L23a (RPL23A) and MYC proto-oncogene, BHLH transcription factor (MYC) were identified. CONCLUSIONS: Prolonged egg storage (NSR) resulted in enriched cell stress and death pathways; while SPIDES (SR) resulted in enriched basic cell and anti-apoptotic pathways. New insights into DNA repair mechanisms, RNA processing, shifts in metabolism, and chromatin dynamics in relation to egg storage treatment were obtained through this study. Although egg storage protocols have been examined through targeted gene expression approaches, this study provided a global view of the extensive molecular networks affected by prolonged storage and SPIDES and helped to identify potential upstream regulators for future experiments to optimize egg storage parameters.

BACKGROUND: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused global disruption of human health and activity. Being able to trace the early outbreak of SARS-CoV-2 within a locality can inform public health measures and provide insights to contain or prevent viral transmission. Investigation of the transmission history requires efficient sequencing methods and analytic strategies, which can be generally useful in the study of viral outbreaks. METHODS: The County of Los Angeles (hereafter, LA County) sustained a large outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To learn about the transmission history, we carried out surveillance viral genome sequencing to determine 142 viral genomes from unique patients seeking care at the University of California, Los Angeles (UCLA) Health System. 86 of these genomes were from samples collected before April 19, 2020. RESULTS: We found that the early outbreak in LA County, as in other international air travel hubs, was seeded by multiple introductions of strains from Asia and Europe. We identified a USA-specific strain, B.1.43, which was found predominantly in California and Washington State. While samples from LA County carried the ancestral B.1.43 genome, viral genomes from neighboring counties in California and from counties in Washington State carried additional mutations, suggesting a potential origin of B.1.43 in Southern California. We quantified the transmission rate of SARS-CoV-2 over time, and found evidence that the public health measures put in place in LA County to control the virus were effective at preventing transmission, but might have been undermined by the many introductions of SARS-CoV-2 into the region. CONCLUSION: Our work demonstrates that genome sequencing can be a powerful tool for investigating outbreaks and informing the public health response. Our results reinforce the critical need for the USA to have coordinated inter-state responses to the pandemic.

BACKGROUND: Colorectal cancer (CRC) is the third most common cancer worldwide and is influenced by environmental and genetic factors. Although numerous genetic loci for CRC have been identified, the overall understanding of the genetic factors is yet to be elucidated. We sought to discover new genes involved in CRC applying genetic association analysis and functional study. RESULTS: We conducted exome array analysis on 194 CRC and 600 control subjects for discovering new candidate CRC genes. Fisher's exact test detected one exome-wide significant functional locus for CRC on SMCO1 (P < 10-6) and two suggestive functional loci on HLA-C and NUTM1 (10-6 ≤ P < 10-4). To evaluate the biological role of three candidate CRC genes, the differential expression of these genes between CRC and non-cancer colorectal cells was analyzed using qRT-PCR and publicly available gene expression data. Of three genes, HLA-C consistently revealed the significant down-regulation in CRC cells. In addition, we detected a reduction in cell viability in the HLA-C overexpression CRC cell line, implying the functional relevance of HLA-C in CRC. To understand the underlying mechanism exerted by HLA-C in CRC development, we conducted RNA sequencing analyses of HLA-C overexpression CRC cells and non-cancer colorectal cells. Pathway analysis detected that significantly down-regulated genes in HLA-C overexpression CRC cells were highly enriched in cancer-related signaling pathways such as JAK/STAT, ErbB, and Hedgehog signaling pathways. CONCLUSIONS: Exome array CRC case-control analysis followed by functional validation demonstrated that HLA-C likely exerts its influence on CRC development via cancer-related signaling pathways.

BACKGROUND: Allotetraploid oilseed rape (Brassica napus L.) is an important worldwide oil-producing crop. The origin of rapeseed is still undetermined due to the lack of wild resources. Despite certain genetic architecture and phylogenetic studies have been done focus on large group of Brassica nuclear genomes, the organelle genomes information under global pattern is largely unknown, which provide unique material for phylogenetic studies of B. napus. Here, based on de novo assemblies of 1,579 B. napus accessions collected globally, we constructed the chloroplast and mitochondrial pan-genomes of B. napus, and investigated the genetic diversity, phylogenetic relationships of B. napus, B. rapa and B. oleracea. RESULTS: Based on mitotype-specific markers and mitotype-variant ORFs, four main cytoplasmic haplotypes were identified in our groups corresponding the nap, pol, ole, and cam mitotypes, among which the structure of chloroplast genomes was more conserved without any rearrangement than mitochondrial genomes. A total of 2,092 variants were detected in chloroplast genomes, whereas only 326 in mitochondrial genomes, indicating that chloroplast genomes exhibited a higher level of single-base polymorphism than mitochondrial genomes. Based on whole-genome variants diversity analysis, eleven genetic difference regions among different cytoplasmic haplotypes were identified on chloroplast genomes. The phylogenetic tree incorporating accessions of the B. rapa, B. oleracea, natural and synthetic populations of B. napus revealed multiple origins of B. napus cytoplasm. The cam-type and pol-type were both derived from B. rapa, while the ole-type was originated from B. oleracea. Notably, the nap-type cytoplasm was identified in both the B. rapa population and the synthetic B. napus, suggesting that B. rapa might be the maternal ancestor of nap-type B. napus. CONCLUSIONS: The phylogenetic results provide novel insights into the organelle genomic evolution of Brassica species. The natural rapeseeds contained at least four cytoplastic haplotypes, of which the predominant nap-type might be originated from B. rapa. Besides, the organelle pan-genomes and the overall variation data offered useful resources for analysis of cytoplasmic inheritance related agronomical important traits of rapeseed, which can substantially facilitate the cultivation and improvement of rapeseed varieties.

BACKGROUND: The cohesin complex is essential for proper chromosome structure and gene expression. Defects in cohesin subunits and regulators cause changes in cohesin complex dynamics and thereby alter three-dimensional genome organization. However, the molecular mechanisms that drive cohesin localization and function remain poorly understood. RESULTS: In this study, we observe that loss of WIZ causes changes to cohesin localization that are distinct from loss of the known WIZ binding partner G9a. Whereas loss of WIZ uniformly increases cohesin levels on chromatin at known binding sites and leads to new, ectopic cohesin binding sites, loss of G9a does not. Ectopic cohesin binding on chromatin after the loss of WIZ occurs at regions that are enriched for activating histone modifications and transcription factors motifs. Furthermore, loss of WIZ causes changes in cohesin localization that are distinct from those observed by loss of WAPL, the canonical cohesin unloading factor. CONCLUSIONS: The evidence presented here suggests that WIZ can function independently from its previously identified role with G9a and GLP in heterochromatin formation. Furthermore, while WIZ limits the levels and localization pattern of cohesin across the genome, it appears to function independently of WAPL-mediated cohesin unloading.

BACKGROUND: Gram-negative bacteria are important pathogens in cattle, causing severe infectious diseases, including mastitis. Lipopolysaccharides (LPS) are components of the outer membrane of Gram-negative bacteria and crucial mediators of chronic inflammation in cattle. LPS modulations of bovine immune responses have been studied before. However, the single-cell transcriptomic and chromatin accessibility analyses of bovine peripheral blood mononuclear cells (PBMCs) and their responses to LPS stimulation were never reported. RESULTS: We performed single-cell RNA sequencing (scRNA-seq) and single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) in bovine PBMCs before and after LPS treatment and demonstrated that seven major cell types, which included CD4 T cells, CD8 T cells, and B cells, monocytes, natural killer cells, innate lymphoid cells, and dendritic cells. Bioinformatic analyses indicated that LPS could increase PBMC cell cycle progression, cellular differentiation, and chromatin accessibility. Gene analyses further showed significant changes in differential expression, transcription factor binding site, gene ontology, and regulatory interactions during the PBMC responses to LPS. Consistent with the findings of previous studies, LPS induced activation of monocytes and dendritic cells, likely through their upregulated TLR4 receptor. NF-κB was observed to be activated by LPS and an increased transcription of an array of pro-inflammatory cytokines, in agreement that NF-κB is an LPS-responsive regulator of innate immune responses. In addition, by integrating LPS-induced differentially expressed genes (DEGs) with large-scale GWAS of 45 complex traits in Holstein, we detected trait-relevant cell types. We found that selected DEGs were significantly associated with immune-relevant health, milk production, and body conformation traits. CONCLUSION: This study provided the first scRNAseq and scATAC-seq data for cattle PBMCs and their responses to the LPS stimulation to the best of our knowledge. These results should also serve as valuable resources for the future study of the bovine immune system and open the door for discoveries about immune cell roles in complex traits like mastitis at single-cell resolution.

BACKGROUND: Drought has become the major abiotic stress that causes losses in rice yields and consequently is one of the main environmental factors threatening food security. Long non-coding RNA (lncRNA) is known to play an important role in plant response to drought stress, while the mechanisms of competing endogenous RNA (ceRNA) in drought resistance in upland rice have been rarely reported. RESULTS: In our study, a total of 191 lncRNAs, 2115 mRNAs and 32 miRNAs (microRNAs) were found by strand-specific sequencing and small RNA sequencing to be differentially expressed in drought-stressed rice. Functional analysis of results indicate that they play important roles in hormone signal transduction, chlorophyll synthesis, protein synthesis and other pathways. Construction of a ceRNA network revealed that MSTRG.28732.3 may interact with miR171 in the chlorophyll biosynthesis pathway and affect the ability of plants to withstand drought stress by regulating Os02g0662700, Os02g0663100 and Os06g0105350. The accuracy of the regulatory network was verified by qRT-PCR. CONCLUSION: Our results provide a theoretical basis for future studies on the potential function of lncRNA in plant drought resistance, and they provide new genetic resources for drought-resistant rice breeding.

BACKGROUND: Human granulocytic anaplasmosis is a tick-borne zoonotic disease caused by Anaplasma phagocytophilum. Coinfections with A. phagocytophilum and other tick-borne pathogens are reported frequently, whereas the relationship between A. phagocytophilum and flea-borne Yersnia pestis is rarely concerned. RESULTS: A. phagocytophilum and Yersnia pestis were discovered within a Marmota himalayana found dead in the environment, as determined by 16S ribosomal rRNA sequencing. Comparative genomic analyses of marmot-derived A. phagocytophilum isolate demonstrated its similarities and a geographic isolation from other global strains. The 16S rRNA gene and GroEL amino acid sequence identity rates between marmot-derived A. phagocytophilum (JAHLEX000000000) and reference strain HZ (CP000235.1) are 99.73% (1490/1494) and 99.82% (549/550), respectively. 16S rRNA and groESL gene screenings show that A. phagocytophilum is widely distributed in marmots; the bacterium was more common in marmots found dead (24.59%, 15/61) than in captured marmots (19.21%, 29/151). We found a higher Y. pestis isolation rate in dead marmots harboring A. phagocytophilum than in those without it (2 = 4.047, p < 0.05). Marmot-derived A. phagocytophilum was able to live in L929 cells and BALB/c mice but did not propagate well. CONCLUSIONS: In this study, A. phagocytophilum was identified for the first time in Marmota himalayana, a predominant Yersinia pestis host. Our results provide initial evidence for M. himalayana being a reservoir for A. phagocytophilum; moreover, we found with the presence of A. phagocytophilum, marmots may be more vulnerable to plague. Humans are at risk for co-infection with both pathogens by exposure to such marmots.

BACKGROUND: Poleroviruses, such as turnip yellows virus (TuYV), are plant viruses strictly transmitted by aphids in a persistent and circulative manner. Acquisition of either virus particles or plant material altered by virus infection is expected to induce gene expression deregulation in aphids which may ultimately alter their behavior. RESULTS: By conducting an RNA-Seq analysis on viruliferous aphids fed either on TuYV-infected plants or on an artificial medium containing purified virus particles, we identified several hundreds of genes deregulated in Myzus persicae, despite non-replication of the virus in the vector. Only a few genes linked to receptor activities and/or vesicular transport were common between the two modes of acquisition with, however, a low level of deregulation. Behavioral studies on aphids after virus acquisition showed that M. persicae locomotion behavior was affected by feeding on TuYV-infected plants, but not by feeding on the artificial medium containing the purified virus particles. Consistent with this, genes potentially involved in aphid behavior were deregulated in aphids fed on infected plants, but not on the artificial medium. CONCLUSIONS: These data show that TuYV particles acquisition alone is associated with a moderate deregulation of a few genes, while higher gene deregulation is associated with aphid ingestion of phloem from TuYV-infected plants. Our data are also in favor of a major role of infected plant components on aphid behavior.

Erratum in BMC Genomics. 2022 May 19;23(1):384.

BACKGROUND: The pro-neural transcription factor ASCL1 is a master regulator of neurogenesis and a key factor necessary for the reprogramming of permissive cell types to neurons. Endogenously, ASCL1 expression is often associated with neuroblast stem-ness. Moreover, ASCL1-mediated reprogramming of fibroblasts to differentiated neurons is commonly achieved using artificially high levels of ASCL1 protein, where ASCL1 acts as an "on-target" pioneer factor. However, the genome-wide effects of enhancing ASCL1 activity in a permissive neurogenic environment has not been thoroughly investigated. Here, we overexpressed ASCL1 in the neuronally-permissive context of neuroblastoma (NB) cells where modest endogenous ASCL1 supports the neuroblast programme. RESULTS: Increasing ASCL1 in neuroblastoma cells both enhances binding at existing ASCL1 sites and also leads to creation of numerous additional, lower affinity binding sites. These extensive genome-wide changes in ASCL1 binding result in significant reprogramming of the NB transcriptome, redirecting it from a proliferative neuroblastic state towards one favouring neuronal differentiation. Mechanistically, ASCL1-mediated cell cycle exit and differentiation can be increased further by preventing its multi-site phosphorylation, which is associated with additional changes in genome-wide binding and gene activation profiles. CONCLUSIONS: Our findings show that enhancing ASCL1 activity in a neurogenic environment both increases binding at endogenous ASCL1 sites and also results in additional binding to new low affinity sites that favours neuronal differentiation over the proliferating neuroblast programme supported by the endogenous protein. These findings have important implications for controlling processes of neurogenesis in cancer and cellular reprogramming.

BACKGROUND: Ophiorrhiza pumila (Rubiaceae) is capable of producing camptothecin (CPT), one monoterpene indole alkaloid extensively employed in the treatment of multiple cancers. Transcription factors (TFs) GATA are a group of transcription regulators involved in plant development and metabolism, and show the feature of binding to the GATA motif within the promoters of target genes. However, GATA TFs have not been characterized in O. pumila. RESULT: In this study, a total of 18 GATA genes classified into four subfamilies were identified, which randomly distributed on 11 chromosomes of O. pumila. Synteny analysis of GATA genes between O. pumila and other plant species such as Arabidopsis thaliana, Oryza sativa, Glycine max, Solanum lycopersicum, Vitis vinifera, and Catharanthus roseus genomes were analyzed. Tissue expression pattern revealed that OpGATA1 and OpGATA18 were found to be correlated with ASA, MK, CPR and GPPS, which were highly expressed in leaves. OpGATA7, showed high expression in roots as most of the CPT biosynthetic pathway genes did, suggesting that these OpGATAs may be potential candidates regulating CPT biosynthesis in O. pumila. CONCLUSIONS: In this study, we systematically analyzed the OpGATA TFs, and provided insights into the involvement of OpGATA TFs from O. pumila in CPT biosynthesis.

BACKGROUND: Oculocutaneous albinism (OCA) is a group of heterogeneous genetic diseases characterized by a reduction or complete lack of pigmentation in the hair, skin, and eyes. It is associated with reduced visual acuity, nystagmus, photophobia, and strabismus. OCA type 1 (OCA1) and type 2 (OCA2) are caused by mutations in the tyrosinase (TYR) and OCA2 genes, which are responsible for most cases of OCA. The present study aimed to identify the mutational spectra of 18 southwest Chinese probands with OCA. RESULTS: We used a skin disease-targeted panel to sequence more than 400 genes, including 23 genes (TYR, OCA2, AP3B1, BLOC1S3, BLOC1S6, C10orf11, DTNBP1, FRMD7, GPR143, HPS1, HPS3, HPS4, HPS5, HPS6, LYST, MC1R, MITF, MLPH, MYO5A, RAB27A, SLC24A5, SLC45A2, TYRP1) associated with syndromic and non-syndromic albinism. The targeted panel was applied to 18 patients from southwest China, nine (50%) patients were diagnosed with OCA1, and nine (50%) were diagnosed with OCA2. Our data indicate that OCA1 and OCA2, the most common subtypes, probably have the same prevalence in southwest China. In total, we identified 26 variants in TYR and OCA2 from 18 OCA cases using the NGS technology, including 24 variants presented in the Human Gene Mutation Database Professional (HGMD) and two novel variants, c.559_560insCATTATTATGTGTCAAATTATCCCC in TYR and c.1514 T > C in OCA2, which have not been previously reported. According to the American College of Medical Genetics and Genomics (ACMG) classification, c.559_560insCATTATTATGTGTCAAATTATCCCC (p.G190Cfs*12) is classified as a pathogenic variant, and c.1514 T > C (p.F505S) is evaluated as a likely pathogenic variant. CONCLUSIONS: Two novel variants were identified which will expand the mutational spectra of TYR and OCA2. The results of the present study may have implications for genetic counseling, carrier screening, and clinical management of the disease.

BACKGROUND: Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic. RESULTS: There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits. CONCLUSIONS: These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle.

BACKGROUND: Chicken provides humans with a large amount of animal protein every year, in which skeletal muscle plays a leading role. The embryonic skeletal muscle development determines the number of muscle fibers and will affect the muscle production of chickens. CircRNAs are involved in a variety of important biological processes, including muscle development. However, studies on circRNAs in the chicken embryo muscle development are still lacking. RESULTS: In the study, we collected chicken leg muscles at 14 and 20-day embryo ages both in the fast- and slow-growing groups for RNA-seq. We identified 245 and 440 differentially expressed (DE) circRNAs in the comparison group F14vsF20 and S14vsS20 respectively. GO enrichment analysis for the host genes of DE circRNAs showed that biological process (BP) terms in the top 20 related to growth in F14vsF20 were found such as positive regulation of transcription involved in G1/S phase of mitotic cell cycle, multicellular organismal macromolecule metabolic process, and multicellular organismal metabolic process. In group S14vsS20, we also found some BP terms associated with growth in the top 20 including actomyosin structure organization, actin cytoskeleton organization and myofibril assembly. A total of 7 significantly enriched pathways were obtained, containing Adherens junction and Tight junction. Further analysis of those pathways found three crucial host genes MYH9, YBX3, IGF1R in both fast- and slow-growing groups, three important host genes CTNNA3, AFDN and CREBBP only in the fast-growing group, and six host genes FGFR2, ACTN2, COL1A2, CDC42, DOCK1 and MYL3 only in the slow-growing group. In addition, circRNA-miRNA network also revealed some key regulation pairs such as novel_circ_0007646-miR-1625-5p, novel_circ_0007646-miR-1680-5p, novel_circ_0008913-miR-148b-5p, novel_circ_0008906-miR-148b-5p and novel_circ_0001640-miR-1759-3p. CONCLUSIONS: Comprehensive analysis of circRNAs and their targets would contribute to a better understanding of the molecular mechanisms in poultry skeletal muscle and it also plays an important guiding role in the next research.

The cosmopolitan Thalassionema species are often dominant components of the plankton diatom flora and sediment diatom assemblages in all but the Polar regions, making important ecological contribution to primary productivity. Historical studies concentrated on their indicative function for the marine environment based primarily on morphological features and essentially ignored their genomic information, hindering in-depth investigation on Thalassionema biodiversity. In this project, we constructed the complete chloroplast genomes (cpDNAs) of seven Thalassionema strains representing three different species, which were also the first cpDNAs constructed for any species in the order Thalassionematales that includes 35 reported species and varieties. The sizes of these Thalassionema cpDNAs, which showed typical quadripartite structures, varied from 124,127 bp to 140,121 bp. Comparative analysis revealed that Thalassionema cpDNAs possess conserved gene content inter-species and intra-species, along with several gene losses and transfers. Besides, their cpDNAs also have expanded inverted repeat regions (IRs) and preserve large intergenic spacers compared to other diatom cpDNAs. In addition, substantial genome rearrangements were discovered not only among different Thalassionema species but also among strains of a same species T. frauenfeldii, suggesting much higher diversity than previous reports. In addition to confirming the phylogenetic position of Thalassionema species, this study also estimated their emergence time at approximately 38 Mya. The availability of the Thalassionema species cpDNAs not only helps understand the Thalassionema species, but also facilitates phylogenetic analysis of diatoms.

BACKGROUND: Auxin responsive factor (ARF) family is one of core components in auxin signalling pathway, which governs diverse developmental processes and stress responses. Blueberry is an economically important berry-bearing crop and prefers to acidic soil. However, the understandings of ARF family has not yet been reported in blueberry. RESULTS: In the present study, 60 ARF genes (VcARF) were identified in blueberry, and they showed diverse gene structures and motif compositions among the groups and similar within each group in the phylogenetic tree. Noticeably, 9 digenic, 5 trigenic and 6 tetragenic VcARF pairs exhibited more than 95% identity to each other. Computational analysis indicated that 23 VcARFs harbored the miRNA responsive element (MRE) of miR160 or miR167 like other plant ARF genes. Interestingly, the MRE of miR156d/h-3p was observed in the 5'UTR of 3 VcARFs, suggesting a potentially novel post-transcriptional control. Furthermore, the transcript accumulations of VcARFs were investigated during fruit development, and three categories of transcript profiles were observed, implying different functional roles. Meanwhile, the expressions of VcARFs to different pH conditions (pH4.5 and pH6.5) were surveyed in pH-sensitive and tolerant blueberry species, and a number of VcARFs showed different transcript accumulations. More importantly, distinct transcriptional response to pH stress (pH6.5) were observed for several VcARFs (such as VcARF6s and VcARF19-3/19-4) between pH-sensitive and tolerant species, suggesting their potential roles in adaption to pH stress. CONCLUSIONS: Sixty VcARF genes were identified and characterized, and their transcript profiles were surveyed during fruit development and in response to pH stress. These findings will contribute to future research for eliciting the functional roles of VcARFs and regulatory mechanisms, especially fruit development and adaption to pH stress.

BACKGROUND: Sex chromosomes have evolved numerous times, as revealed by recent genomic studies. However, large gaps in our knowledge of sex chromosome diversity across the tree of life remain. Filling these gaps, through the study of novel species, is crucial for improved understanding of why and how sex chromosomes evolve. Characterization of sex chromosomes in already well-studied organisms is also important to avoid misinterpretations of population genomic patterns caused by undetected sex chromosome variation. RESULTS: Here we present findZX, an automated Snakemake-based computational pipeline for detecting and visualizing sex chromosomes through differences in genome coverage and heterozygosity between any number of males and females. A main feature of the pipeline is the option to perform a genome coordinate liftover to a reference genome of another species. This allows users to inspect sex-linked regions over larger contiguous chromosome regions, while also providing important between-species synteny information. To demonstrate its effectiveness, we applied findZX to publicly available genomic data from species belonging to widely different taxonomic groups (mammals, birds, reptiles, and fish), with sex chromosome systems of different ages, sizes, and levels of differentiation. We also demonstrate that the liftover method is robust over large phylogenetic distances (> 80 million years of evolution). CONCLUSIONS: With findZX we provide an easy-to-use and highly effective tool for identification of sex chromosomes. The pipeline is compatible with both Linux and MacOS systems, and scalable to suit different computational platforms.

BACKGROUND: Most crop seeds are F1 hybrids. Seed providers and plant breeders must be confident that the seed supplied to growers is of known, and uniform, genetic makeup. This requires maintenance of pure genotypes of the parental lines and testing to ensure the genetic purity of the F1 seed. Traditionally, seed purity has been assessed with a grow-out test (GOT) in the field, a time consuming and costly venture. Early in the last decade, seed testing with molecular markers was introduced as a replacement for GOT, and Kompetitive allele specific PCR (KASP) markers were recognized as promising tools for genetic testing of seeds. However, the markers available at that time could be inaccurate and applicable to only a small number of accessions or varieties due to the limited genetic information and reference genomes available. RESULTS: We identified 4,925,742 SNPs in 50 accessions of the Brasscia rapa core collection. From these, we identified 2,925 SNPs as accession-specific, considering properties of flanking region harboring accession-specific SNPs and genic region conservation among accessions by the Next Generation Sequencing (NGS) analysis. In total, 100 accession-specific markers were developed as accession-specific KASP markers. Based on the results of our validation experiments, the accession-specific markers successfully distinguised individuals from the mixed population including 50 target accessions from B. rapa core collection and the outgroup. Additionally, the marker set we developed here discriminated F1 hybrids and their parental lines with distinct clusters. CONCLUSIONS: This study provides efficient methods for developing KASP markers to distinguish individuals from the mixture comprised of breeding lines and germplasms from the resequencing data of Chinese cabbage (Brassica rapa spp. pekinensis).

BACKGROUND: Ginsenoside, as the main active substance in ginseng, has the function of treating various diseases. However, the ginsenosides content of cultivated ginseng is obviously affected by the growth years, but the molecular mechanism is not clear. In addition, there are significant differences in morphology and physiology between wild ginseng and cultivated ginseng, and the effect of growth years on ginsenoside synthesis not yet understood in wild ginseng. RESULTS: Transcriptome sequencing on the roots, stems and leaves of cultivated ginseng and wild ginseng with different growth years was performed in this study, exploring the effect of growth years on gene expression in ginseng. The number of differentially expressed genes (DEGs) from comparison groups in cultivated ginseng was higher than that in wild ginseng. The result of weighted gene co-expression network analysis (WGCNA) showed that growth years significantly affected the gene expression of Mitogen-activated protein kinases (MAPK) signaling pathway and terpenoid backbone biosynthesis pathway in cultivated ginseng, but had no effects in wild ginseng. Furthermore, the growth years had significant effects on the genes related to ginsenoside synthesis in cultivated ginseng, and the effects were different in the roots, stems and leaves. However, it had little influence on the expression of genes related to ginsenoside synthesis in wild ginseng. Growth years might affect the expression of genes for ginsenoside synthesis by influencing the expression of these transcription factors (TFs), like my elob lastosis (MYB), NAM, ATAF1 and 2, and CUC2 (NAC), APETALA2/ethylene-responsive factor (AP2/ERF), basic helix-loop-helix (bHLH) and WRKY, etc., thereby affecting the content of ginsenosides. CONCLUSIONS: This study complemented the gaps in the genetic information of wild ginseng in different growth periods and helped to clarify the potential mechanisms of the effect of growth years on the physiological state in wild ginseng and cultivated ginseng, which also provided a new insight into the mechanism of ginsenoside regulation.

BACKGROUND: Structural variants (SVs) play a crucial role in gene regulation, trait association, and disease in humans. SV genotyping has been extensively applied in genomics research and clinical diagnosis. Although a growing number of SV genotyping methods for long reads have been developed, a comprehensive performance assessment of these methods has yet to be done. RESULTS: Based on one simulated and three real SV datasets, we performed an in-depth evaluation of five SV genotyping methods, including cuteSV, LRcaller, Sniffles, SVJedi, and VaPoR. The results show that for insertions and deletions, cuteSV and LRcaller have similar F1 scores (cuteSV, insertions: 0.69-0.90, deletions: 0.77-0.90 and LRcaller, insertions: 0.67-0.87, deletions: 0.74-0.91) and are superior to other methods. For duplications, inversions, and translocations, LRcaller yields the most accurate genotyping results (0.84, 0.68, and 0.47, respectively). When genotyping SVs located in tandem repeat region or with imprecise breakpoints, cuteSV (insertions and deletions) and LRcaller (duplications, inversions, and translocations) are better than other methods. In addition, we observed a decrease in F1 scores when the SV size increased. Finally, our analyses suggest that the F1 scores of these methods reach the point of diminishing returns at 20× depth of coverage. CONCLUSIONS: We present an in-depth benchmark study of long-read SV genotyping methods. Our results highlight the advantages and disadvantages of each genotyping method, which provide practical guidance for optimal application selection and prospective directions for tool improvement.

BACKGROUND: In some vertebrate species, gene-environment interactions can determine sex, driving bipotential gonads to differentiate into either ovaries or testes. In the central bearded dragon (Pogona vitticeps), the genetic influence of sex chromosomes (ZZ/ZW) can be overridden by high incubation temperatures, causing ZZ male to female sex reversal. Previous research showed ovotestes, a rare gonadal phenotype with traits of both sexes, develop during sex reversal, leading to the hypothesis that sex reversal relies on high temperature feminisation to outcompete the male genetic cue. To test this, we conducted temperature switching experiments at key developmental stages, and analysed the effect on gonadal phenotypes using histology and transcriptomics. RESULTS: We found sexual fate is more strongly influenced by the ZZ genotype than temperature. Any exposure to low temperatures (28 °C) caused testes differentiation, whereas sex reversal required longer exposure to high temperatures. We revealed ovotestes exist along a spectrum of femaleness to male-ness at the transcriptional level. We found inter-individual variation in gene expression changes following temperature switches, suggesting both genetic sensitivity to, and the timing and duration of the temperature cue influences sex reversal. CONCLUSIONS: These findings bring new insights to the mechanisms underlying sex reversal, improving our understanding of thermosensitive sex systems in vertebrates.

BACKGROUND: Previous phylogeographic studies of the lion (Panthera leo) have improved our insight into the distribution of genetic variation, as well as a revised taxonomy which now recognizes a northern (Panthera leo leo) and a southern (Panthera leo melanochaita) subspecies. However, existing whole range phylogeographic studies on lions either consist of very limited numbers of samples, or are focused on mitochondrial DNA and/or a limited set of microsatellites. The geographic extent of genetic lineages and their phylogenetic relationships remain uncertain, clouded by massive sampling gaps, sex-biased dispersal and incomplete lineage sorting. RESULTS: In this study we present results of low depth whole genome sequencing and subsequent variant calling in ten lions sampled throughout the geographic range, resulting in the discovery of >150,000 Single Nucleotide Polymorphisms (SNPs). Phylogenetic analyses revealed the same basal split between northern and southern populations, as well as four population clusters on a more local scale. Further, we designed a SNP panel, including 125 autosomal and 14 mitochondrial SNPs, which was tested on >200 lions from across their range. Results allow us to assign individuals to one of these four major clades (West & Central Africa, India, East Africa, or Southern Africa) and delineate these clades in more detail. CONCLUSIONS: The results presented here, particularly the validated SNP panel, have important applications, not only for studying populations on a local geographic scale, but also for tracing samples of unknown origin for forensic purposes, and for guiding conservation management of ex situ populations. Thus, these genomic resources not only contribute to our understanding of the evolutionary history of the lion, but may also play a crucial role in conservation efforts aimed at protecting the species in its full diversity.

BACKGROUND: The basic leucine zipper (bZIP) transcription factor (TF) is one of the largest families of transcription factors (TFs). It is widely distributed and highly conserved in animals, plants, and microorganisms. Previous studies have shown that the bZIP TF family is involved in plant growth, development, and stress responses. The bZIP family has been studied in many plants; however, there is little research on the bZIP gene family in tobacco. RESULTS: In this study, 77 bZIPs were identified in tobacco and named NtbZIP01 through to NtbZIP77. These 77 genes were then divided into eleven subfamilies according to their homology with Arabidopsis thaliana. NtbZIPs were unevenly distributed across twenty-two tobacco chromosomes, and we found sixteen pairs of segmental duplication. We further studied the collinearity between these genes and related genes of six other species. Quantitative real-time polymerase chain reaction analysis identified that expression patterns of bZIPs differed, including in different organs and under various abiotic stresses. NtbZIP49 might be important in the development of flowers and fruits; NtbZIP18 might be an important regulator in abiotic stress. CONCLUSIONS: In this study, the structures and functions of the bZIP family in tobacco were systematically explored. Many bZIPs may play vital roles in the regulation of organ development, growth, and responses to abiotic stresses. This research has great significance for the functional characterisation of the tobacco bZIP family and our understanding of the bZIP family in higher plants.

BACKGROUND: Tibial dyschondroplasia (TD) is a bone disorder in which dead chondrocytes accumulate as a result of apoptosis and non-vascularization in the tibial bone of broiler chickens. The pathogenicity of TD is under extensive research but is yet not fully understood. Several studies have linked it to apoptosis and non-vascularization in the tibial growth plate (GP). We conceived the idea to find the differentially expressed genes (DEGs) in chicken erythrocytes which vary in expression over time using a likelihood-ratio test (LRT). Thiram was used to induce TD in chickens, and then injected Ex-FABP protein at 0, 20, and 50 μg.kg-1 to evaluate its therapeutic effect on 30 screened immunity and angiogenesis-related genes using quantitative PCR (qPCR). The histopathology was also performed in TD chickens to explore the shape, circularity, arrangements of chondrocytes and blood vessels. RESULTS: Clinical lameness was observed in TD chickens, which decreased with the injection of Ex-FABP. Histopathological findings support Ex-FABP as a therapeutic agent for the morphology and vascularization of affected chondrocytes in TD chickens. qPCR results of 10 immunity (TLR2, TLR3, TLR4, TLR5, TLR7, TLR15, IL-7, MyD88, MHCII, and TRAF6) and 20 angiogenesis-related genes (ITGAV, ITGA2, ITGB2, ITGB3, ITGA5, IL1R1, TBXA2R, RPL17, F13A1, CLU, RAC2, RAP1B, GIT1, FYN, IQGAP2, PTCH1, NCOR2, VAV-like, PTPN11, MAML3) regulated when Ex-FABP is injected to TD chickens. CONCLUSION: Immunity and angiogenesis-related genes can be responsible for apoptosis of chondrocytes and vascularization in tibial GP. Injection of Ex-FABP protein to thiram induced TD chickens decrease the chondrocytes damage and improves vascularization.

BACKGROUND: Whole mitochondrial genomes are quickly becoming markers of choice for the exploration of within-species genealogical and among-species phylogenetic relationships. Most often, 'primer walking' or 'long PCR' strategies plus Sanger sequencing or low-pass whole genome sequencing using Illumina short reads are used for the assembling of mitochondrial chromosomes. In this study, we first confirmed that mitochondrial genomes can be sequenced from long reads using nanopore sequencing data exclusively. Next, we examined the accuracy of the long-reads assembled mitochondrial chromosomes when comparing them to a 'gold' standard reference mitochondrial chromosome assembled using Illumina short-reads sequencing. RESULTS: Using a specialized bioinformatics tool, we first produced a short-reads mitochondrial genome assembly for the silky shark C. falciformis with an average base coverage of 9.8x. The complete mitochondrial genome of C. falciformis was 16,705 bp in length and 934 bp shorter than a previously assembled genome (17,639 bp in length) that used bioinformatics tools not specialized for the assembly of mitochondrial chromosomes. Next, low-pass whole genome sequencing using a MinION ONT pocket-sized platform plus customized de-novo and reference-based workflows assembled and circularized a highly accurate mitochondrial genome in the silky shark Carcharhinus falciformis. Indels at the flanks of homopolymer regions explained most of the dissimilarities observed between the 'gold' standard reference mitochondrial genome (assembled using Illumina short reads) and each of the long-reads mitochondrial genome assemblies. Although not completely accurate, mitophylogenomics and barcoding analyses (using entire mitogenomes and the D-Loop/Control Region, respectively) suggest that long-reads assembled mitochondrial genomes are reliable for identifying a sequenced individual, such as C. falciformis, and separating the same individual from others belonging to closely related congeneric species. CONCLUSIONS: This study confirms that mitochondrial genomes can be sequenced from long-reads nanopore sequencing data exclusively. With further development, nanopore technology can be used to quickly test in situ mislabeling in the shark fin fishing industry and thus, improve surveillance protocols, law enforcement, and the regulation of this fishery. This study will also assist with the transferring of high-throughput sequencing technology to middle- and low-income countries so that international scientists can explore population genomics in sharks using inclusive research strategies. Lastly, we recommend assembling mitochondrial genomes using specialized assemblers instead of other assemblers developed for bacterial and/or nuclear genomes.

BACKGROUND: Over 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the SARS-CoV-2 consensus genomes obtained from the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs. RESULTS: Our results show that the MiSeq had a significantly higher number of consensus genomes classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher genome coverage and mutation counts than the GridION. CONCLUSION: Due to the low genome coverage, high number of indels, and sensitivity to SARS-CoV-2 viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for SARS-CoV-2 genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole consensus genomes.

BACKGROUND: Structural Variations (SVs) are genomic rearrangements derived from duplication, deletion, insertion, inversion, and translocation events. In the past, SVs detection was limited to cytological approaches, then to Next-Generation Sequencing (NGS) short reads and partitioned assemblies. Nowadays, technologies such as DNA long read sequencing and optical mapping have revolutionized the understanding of SVs in genomes, due to the enhancement of the power of SVs detection. This study aims to investigate performance of two techniques, 1) long-read sequencing obtained with the MinION device (Oxford Nanopore Technologies) and 2) optical mapping obtained with Saphyr device (Bionano Genomics) to detect and characterize SVs in the genomes of the two ecotypes of Arabidopsis thaliana, Columbia-0 (Col-0) and Landsberg erecta 1 (Ler-1). RESULTS: We described the SVs detected from the alignment of the best ONT assembly and DLE-1 optical maps of A. thaliana Ler-1 against the public reference genome Col-0 TAIR10.1. After filtering (SV > 1 kb), 1184 and 591 Ler-1 SVs were retained from ONT and Bionano technologies respectively. A total of 948 Ler-1 ONT SVs (80.1%) corresponded to 563 Bionano SVs (95.3%) leading to 563 common locations. The specific locations were scrutinized to assess improvement in SV detection by either technology. The ONT SVs were mostly detected near TE and gene features, and resistance genes seemed particularly impacted. CONCLUSIONS: Structural variations linked to ONT sequencing error were removed and false positives limited, with high quality Bionano SVs being conserved. When compared with the Col-0 TAIR10.1 reference genome, most of the detected SVs discovered by both technologies were found in the same locations. ONT assembly sequence leads to more specific SVs than Bionano one, the latter being more efficient to characterize large SVs. Even if both technologies are complementary approaches, ONT data appears to be more adapted to large scale populations studies, while Bionano performs better in improving assembly and describing specificity of a genome compared to a reference.

BACKGROUND: Inonotus obliquus is an important edible and medicinal mushroom that was shown to have many pharmacological activities in preclinical trials, including anti-inflammatory, antitumor, immunomodulatory, and antioxidant effects. However, the biosynthesis of these pharmacological components has rarely been reported. The lack of genomic information has hindered further molecular characterization of this mushroom. RESULTS: In this study, we report the genome of I. obliquus using a combined high-throughput Illumina NovaSeq with Oxford Nanopore PromethION sequencing platform. The de novo assembled 38.18 Mb I. obliquus genome was determined to harbor 12,525 predicted protein-coding genes, with 81.83% of them having detectable sequence similarities to others available in public databases. Phylogenetic analysis revealed the close evolutionary relationship of I. obliquus with Fomitiporia mediterranea and Sanghuangporus baumii in the Hymenochaetales clade. According to the distribution of reproduction-related genes, we predict that this mushroom possesses a tetrapolar heterothallic reproductive system. The I. obliquus genome was found to encode a repertoire of enzymes involved in carbohydrate metabolism, along with 135 cytochrome P450 proteins. The genome annotation revealed genes encoding key enzymes responsible for secondary metabolite biosynthesis, such as polysaccharides, polyketides, and terpenoids. Among them, we found four polyketide synthases and 20 sesquiterpenoid synthases belonging to four more types of cyclization mechanism, as well as 13 putative biosynthesis gene clusters involved in terpenoid synthesis in I. obliquus. CONCLUSIONS: To the best of our knowledge, this is the first reported genome of I. obliquus; we discussed its genome characteristics and functional annotations in detail and predicted secondary metabolic biosynthesis-related genes, which provides genomic information for future studies on its associated molecular mechanism.

BACKGROUND: Drug-resistant bacteria are important carriers of antibiotic-resistant genes (ARGs). This fact is crucial for the development of precise clinical drug treatment strategies. Long-read sequencing platforms such as the Oxford Nanopore sequencer can improve genome assembly efficiency particularly when they are combined with short-read sequencing data. RESULTS: Alcaligenes faecalis PGB1 was isolated and identified with resistance to penicillin and three other antibiotics. After being sequenced by Nanopore MinION and Illumina sequencer, its entire genome was hybrid-assembled. One chromosome and one plasmid was assembled and annotated with 4,433 genes (including 91 RNA genes). Function annotation and comparison between strains were performed. A phylogenetic analysis revealed that it was closest to A. faecalis ZD02. Resistome related sequences was explored, including ARGs, Insert sequence, phage. Two plasmid aminoglycoside genes were determined to be acquired ARGs. The main ARG category was antibiotic efflux resistance and β-lactamase (EC 3.5.2.6) of PGB1 was assigned to Class A, Subclass A1b, and Cluster LSBL3. CONCLUSIONS: The present study identified the newly isolated bacterium A. faecalis PGB1 and systematically annotated its genome sequence and ARGs.

BACKGROUND: Long noncoding RNAs (lncRNAs) are involved in physiological and pathological processes. However, no studies have been conducted on the relationship between lncRNAs and renal aging. RESULTS: First, we evaluated the histopathology of young (3-month-old) and old (24-month-old) C57BL/6J mouse kidneys. Masson trichrome staining and PAS staining showed interstitial collagen deposition and fibrosis, mesangial matrix expansion, a thicker basement membrane and renal interstitial fibrosis in old mouse kidneys. Senescence-associated β-galactosidase (SA-β-gal)-positive areas in the kidneys of old mice were significantly elevated compared to those of young mice. Then, we analyzed the differential expression of lncRNAs and mRNAs in the kidneys of young and old mouse kidneys by RNA-seq analysis. 42 known and 179 novel differentially expressed lncRNAs and 702 differential mRNAs were detected in the mouse kidney. Next, we focused on the differentially expressed mRNAs and lncRNAs by RNA-seq. GO and KEGG analyses were performed based on differentially expressed mRNAs between young and old mouse kidneys. Transregulation based on RIsearch and the correlation coefficient of mRNA-lncRNA were also calculated. The mRNA-lncRNA network was constructed by choosing a Spearman correlation coefficient > 0.9 or <-0.9. GO and KEGG pathway enrichment analyses revealed that differentially expressed mRNAs participated in aging-related pathways. A total of 10 lncRNAs and trans-regulated mRNAs were constructed. Finally, we validated the role of lncRNA Gm43360 by CCK-8, flow cytometry, western blot and SA-β-gal staining. The expression level of Adra1a was positively correlated and Csnk1a1 was negatively correlated with lncRNA Gm43360. The cell counting kit-8 (CCK-8) results showed that lncRNA Gm43360 promoted cell viability. LncRNA Gm43360 increased the percentage of S phase cells and decreased the percentage of G1 phase cells compared with the negative control. LncRNA Gm43360 decreased the expression of p53, p21 and SA-β-gal. CONCLUSIONS: LncRNA Gm43360 may play a protective role in kidney aging.

BACKGROUND: Classification of parasitic bopyrids has traditionally been based on morphological characteristics, but phylogenetic relationships have remained elusive due to limited information provided by morphological data and tendency for loss of morphological features as a result of parasitic lifestyle. Subfamily Argeiinae was separated from Bopyrinae based on morphological evidence, although the assignment of all genera has not been phylogenetically evaluated. Bopyroides hippolytes has been traditionally classified in Bopyrinae, but divergent morphological characters make this assignment questionable. To investigate the relationship of bopyrines, we sequenced the complete mitochondrial genome of B. hippolytes and four mitochondrial genes of two other Bopyrinae. RESULTS: The phylogenetic trees based on separate and combined cox1and 18S sequence data recovered Bopyridae as robustly monophyletic, but Bopyrinae as polyphyletic. Bopyroides hippolytes was a close sister to Argeia pugettensis, type species to Argeiinae. Mitochondrial phylogenomics also suggested that B. hippolytes was close to Argeiinae. We also found a novel gene order in B. hippolytes compared to other isopods. CONCLUSIONS: Bopyroides hippolytes should be excluded from the Bopyrinae and has a close affinity with Argeia pugettensis based on molecular and morphological data. The conserved syntenic blocks of mitochondrial gene order have distinctive characteristics at a subordinal level and may be helpful for understanding the higher taxonomic level relationships of Isopoda.

BACKGROUND: Deep-sea mussels living in the cold seeps with enormous biomass act as the primary consumers. They are well adapted to the extreme environment where light is absent, and hydrogen sulfide, methane, and other hydrocarbon-rich fluid seepage occur. Despite previous studies on diversity, role, evolution, and symbiosis, the changing adaptation patterns during different developmental stages of the deep-sea mussels remain largely unknown. RESULTS: The deep-sea mussels (Bathymodiolus platifrons) of two developmental stages were collected from the cold seep during the ocean voyage. The gills, mantles, and adductor muscles of these mussels were used for the Illumina sequencing. A total of 135 Gb data were obtained, and subsequently, 46,376 unigenes were generated using de-novo assembly strategy. According to the gene expression analysis, amounts of genes were most actively expressed in the gills, especially genes involved in environmental information processing. Genes encoding Toll-like receptors and sulfate transporters were up-regulated in gills, indicating that the gill acts as both intermedium and protective screen in the deep-sea mussel. Lysosomal enzymes and solute carrier responsible for nutrients absorption were up-regulated in the older mussel, while genes related to toxin resistance and autophagy were up-regulated in the younger one, suggesting that the older mussel might be in a vigorous stage while the younger mussel was still paying efforts in survival and adaptation. CONCLUSIONS: In general, our study suggested that the adaptation capacity might be formed gradually during the development of deep-sea mussels, in which the gill and the symbionts play essential roles.

BACKGROUND: Canga is the Brazilian term for the savanna-like vegetation harboring several endemic species on iron-rich rocky outcrops, usually considered for mining activities. Parkia platycephala Benth. and Stryphnodendron pulcherrimum (Willd.) Hochr. naturally occur in the cangas of Serra dos Carajás (eastern Amazonia, Brazil) and the surrounding forest, indicating high phenotypic plasticity. The morphological and physiological mechanisms of the plants' establishment in the canga environment are well studied, but the molecular adaptative responses are still unknown. To understand these adaptative responses, we aimed to identify molecular mechanisms that allow the establishment of these plants in the canga environment. RESULTS: Plants were grown in canga and forest substrates collected in the Carajás Mineral Province. RNA was extracted from pooled leaf tissue, and RNA-seq paired-end reads were assembled into representative transcriptomes for P. platycephala and S. pulcherrimum containing 31,728 and 31,311 primary transcripts, respectively. We identified both species-specific and core molecular responses in plants grown in the canga substrate using differential expression analyses. In the species-specific analysis, we identified 1,112 and 838 differentially expressed genes for P. platycephala and S. pulcherrimum, respectively. Enrichment analyses showed that unique biological processes and metabolic pathways were affected for each species. Comparative differential expression analysis was based on shared single-copy orthologs. The overall pattern of ortholog expression was species-specific. Even so, we identified almost 300 altered genes between plants in canga and forest substrates with conserved responses in the two species. The genes were functionally associated with the response to light stimulus and the circadian rhythm pathway. CONCLUSIONS: Plants possess species-specific adaptative responses to cope with the substrates. Our results also suggest that plants adapted to both canga and forest environments can adjust the circadian rhythm in a substrate-dependent manner. The circadian clock gene modulation might be a central mechanism regulating the plants' development in the canga substrate in the studied legume species. The mechanism may be shared as a common mechanism to abiotic stress compensation in other native species.

BACKGROUND: Leaf colour mutations are universally expressed at the seedling stage and are ideal materials for exploring the chlorophyll biosynthesis pathway, carotenoid metabolism and the flavonoid biosynthesis pathway in plants. RESULTS: In this research, we analysed the different degrees of albinism in apple (Malus domestica) seedlings, including white-leaf mutants (WM), piebald leaf mutants (PM), light-green leaf mutants (LM) and normal leaves (NL) using bisulfite sequencing (BS-seq) and RNA sequencing (RNA-seq). There were 61,755, 79,824, and 74,899 differentially methylated regions (DMRs) and 7566, 3660, and 3546 differentially expressed genes (DEGs) identified in the WM/NL, PM/NL and LM/NL comparisons, respectively. CONCLUSION: The analysis of the methylome and transcriptome showed that 9 DMR-associated DEGs were involved in the carotenoid metabolism and flavonoid biosynthesis pathway. The expression of different transcription factors (TFs) may also influence the chlorophyll biosynthesis pathway, carotenoid metabolism and the flavonoid biosynthesis pathway in apple leaf mutants. This study provides a new method for understanding the differences in the formation of apple seedlings with different degrees of albinism.

BACKGROUND: Dandruff is a chronic, recurring, and common scalp problem that is caused by several etiopathogeneses with complex mechanisms. Management of this condition is typically achieved via antifungal therapies. However, the precise roles played by microbiota in the development of the condition have not been elucidated. Despite their omnipresence on human scalp little is known about the co-occurrence/co-exclusion network of cutaneous microbiota. RESULTS: We characterized the scalp and hair surface bacterial and fungal communities of 95 dandruff-afflicted and healthy individuals residing in China. The degree distributions of co-occurrence/co-exclusion network in fungi-bacteria and bacteria-bacteria were higher in the healthy group (P < 0.0001), whereas the betweenness values are higher in the dandruff group (P < 0.01). Meanwhile, the co-occurrence/co-exclusion network among fungi-fungi and fungi-bacteria showed that compared to the healthy group, the dandruff group had more positive links (P < 0.0001). In addition, we observed that Malassezia slooffiae, Malassezia japonica and Malassezia furfur, were more abundant in the dandruff group than in the healthy group. These microbiota were co-exclusion by either multiple bacterial genera or Malassezia sp. in healthy group. The lactic acid bacteria on the scalp and hair surface, especially the genera Lactobacillus and Lactococcus, exhibit a negative correlation with multiple bacterial genera on the scalp and hair surface. Lactobacillus plantarum and Pediococcus lactis isolated on the healthy human scalp can inhibit the growth of Staphylococcus epidermidis in vitro. CONCLUSIONS: We showed that microbial networks on scalp and hair surface with dandruff were less integrated than their healthy counterparts, with lower node degree and more positive and stronger links which were deemed to be unstable and may be more susceptible to environmental fluctuations. Lactobacillus bacteria have extensive interactions with other bacteria or fungi in the scalp and hair surface micro-ecological network and can be used as targets for improving scalp health.

BACKGROUND: Psathyrostachys huashanica Keng has long been used as a genetic resource for improving wheat cultivar because of its genes mediating the resistance to various diseases (stripe rust, leaf rust, take-all, and powdery mildew) as well as its desirable agronomic traits. However, a high-resolution fluorescence in situ hybridization (FISH) karyotype of P. huashanica remains unavailable. RESULTS: To develop chromosome-specific FISH markers for P. huashanica, repetitive sequences, including pSc119.2, pTa535, pTa713, pAs1, (AAC)5, (CTT)12, pSc200, pTa71A-2, and Oligo-44 were used for a FISH analysis. The results indicated that the combination of pSc200, pTa71A-2 and Oligo-44 probes can clearly identify all Ns genomic chromosomes in the two P. huashanica germplasms. The homoeologous relationships between individual P. huashanica chromosomes and common wheat chromosomes were clarified by FISH painting. Marker validation analyses revealed that the combination of pSc200, pTa71A-2, and Oligo-44 for a FISH analysis can distinguish the P. huashanica Ns-genome chromosomes from wheat chromosomes, as well as all chromosomes (except 4Ns) from the chromosomes of diploid wheat relatives carrying St, E, V, I, P and R genomes. Additionally, the probes were applicable for discriminating between the P. huashanica Ns-genome chromosomes in all homologous groups and the corresponding chromosomes in Psathyrostachys juncea and most Leymus species containing the Ns genome. Furthermore, six wheat-P. huashanica chromosome addition lines (i.e., 2Ns, 3Ns, 4Ns, 7Ns chromosomes and chromosomal segments) were characterized using the newly developed FISH markers. Thus, these probes can rapidly and precisely detect P. huashanica alien chromosomes in the wheat background. CONCLUSIONS: The FISH karyotype established in this study lays a solid foundation for the efficient identification of P. huashanica chromosomes in wheat genetic improvement programs.

BACKGROUND: Advances in genome sequencing technology, particularly restriction-site associated DNA sequence (RAD-seq) and whole-genome resequencing, have greatly aided the construction of cotton interspecific genetic maps based on single nucleotide polymorphism (SNPs), Indels, and other types of markers. High-density genetic maps can improve accuracy of quantitative trait locus (QTL) mapping, narrow down location intervals, and facilitate identification of the candidate genes. RESULT: In this study, 249 individuals from an interspecific F2 population (TM-1 and Hai7124) were re-sequenced, yielding 6303 high-confidence bin markers spanning 5057.13 cM across 26 cotton chromosomes. A total of 3380 recombination hot regions RHRs were identified which unevenly distributed on the 26 chromosomes. Based on this map, 112 QTLs relating to agronomic and physiological traits from seedling to boll opening stage were identified, including 15 loci associated with 14 traits that contained genes harboring nonsynonymous SNPs. We analyzed the sequence and expression of these ten candidate genes and discovered that GhRHD3 (GH_D10G0500) may affect fiber yield while GhGPAT6 (GH_D04G1426) may affect photosynthesis efficiency. CONCLUSION: Our research illustrates the efficiency of constructing a genetic map using binmap and QTL mapping on the basis of a certain size of the early-generation population. High-density genetic map features high recombination exchanges in number and distribution. The QTLs and the candidate genes identified based on this high-density genetic map may provide important gene resources for the genetic improvement of cotton.

BACKGROUND: Cattle are considered to be the most desirable livestock by small scale farmers. In Africa, although comprehensive genomic studies have been carried out on cattle, the genetic variations in indigenous cattle from Nigeria have not been fully explored. In this study, genome-wide analysis based on genotyping-by-sequencing (GBS) of 193 Nigerian cattle was used to reveal new insights on the history of West African cattle and their adaptation to the tropical African environment, particularly in sub-Saharan region.  RESULTS: The GBS data were evaluated against whole-genome sequencing (WGS) data and high rate of variant concordance between the two platforms was evident with high correlated genetic distance matrices genotyped by both methods suggestive of the reliability of GBS applicability in population genetics. The genetic structure of Nigerian cattle was observed to be homogenous and unique from other African cattle populations. Selection analysis for the genomic regions harboring imprints of adaptation revealed genes associated with immune responses, growth and reproduction, efficiency of feeds utilization, and heat tolerance. Our findings depict potential convergent adaptation between African cattle, dogs and humans with adaptive genes SPRY2 and ITGB1BP1 possibly involved in common physiological activities. CONCLUSION: The study presents unique genetic patterns of Nigerian cattle which provide new insights on the history of cattle in West Africa based on their population structure and the possibility of parallel adaptation between African cattle, dogs and humans in Africa which require further investigations.

BACKGROUND: Chicken intramuscular fat (IMF) content is closely related to meat quality and performance, such as tenderness and flavor. Abdominal fat (AF) in chickens is one of the main waste products at slaughter. Excessive AF reduces feed efficiency and carcass quality. RESULTS: To analyze the differential deposition of IMF and AF in chickens, gene expression profiles in the breast muscle (BM) and AF tissues of 18 animals were analyzed by differential expression analysis and weighted co-expression network analysis. The results showed that IMF deposition in BM was associated with pyruvate and citric acid metabolism through GAPDH, LDHA, GPX1, GBE1, and other genes. In contrast, AF deposition was related to acetyl CoA and glycerol metabolism through FABP1, ELOVL6, SCD, ADIPOQ, and other genes. Carbohydrate metabolism plays an essential role in IMF deposition, and fatty acid and glycerol metabolism regulate AF deposition. CONCLUSION: This study elucidated the molecular mechanism governing IMF and AF deposition through crucial genes and signaling pathways and provided a theoretical basis for producing high-quality broilers.

BACKGROUND: Intermittent fasting (IF), the implementation of fasting periods of at least 12 consecutive hours on a daily to weekly basis, has received a lot of attention in recent years for imparting the life-prolonging and health-promoting effects of caloric restriction with no or only moderate actual restriction of caloric intake. IF is also widely practiced in the rearing of broiler breeders, the parent stock of meat-type chickens, who require strict feed restriction regimens to prevent the serious health problems associated with their intense appetites. Although intermittent fasting has been extensively used in this context to reduce feed competition and its resulting stress, the potential of IF in chickens as an alternative and complementary model to rodents has received less investigation. In both mammals and birds, the liver is a key component of the metabolic response to IF, responding to variations in energy balance. Here we use a microarray analysis to examine the liver transcriptomics of wild-type Red Jungle Fowl chickens fed either ad libitum, chronically restricted to around 70% of ad libitum daily or intermittently fasted (IF) on a 2:1 (2 days fed, 1 day fasted) schedule without actual caloric restriction. As red junglefowl are ancestral to domestic chicken breeds, these data serve as a baseline to which existing and future transcriptomic results from farmed birds such as broiler breeders can be compared. RESULTS: We find large effects of feeding regimen on liver transcriptomics, with most of the affected genes relating to energy metabolism. A cluster analysis shows that IF is associated with large and reciprocal changes in genes related to lipid and carbohydrate metabolism, but also chronic changes in genes related to amino acid metabolism (generally down-regulated) and cell cycle progression (generally up-regulated). The overall transcription pattern appears to be one of promoting high proliferative plasticity in response to fluctuations in available energy substrates. A small number of inflammation-related genes also show chronically changed expression profiles, as does one circadian rhythm gene. CONCLUSIONS: The increase in proliferative potential suggested by the gene expression changes reported here indicates that birds and mammals respond similarly to intermittent fasting practices. Our findings therefore suggest that the health benefits of periodic caloric restriction are ubiquitous and not restricted to mammals alone. Whether a common fundamental mechanism, for example involving leptin, underpins these benefits remains to be elucidated.

BACKGROUND: Measles caused by measles virus (MeV) is a highly contagious viral disease which has also been associated with complications including pneumonia, myocarditis, encephalitis, and subacute sclerosing panencephalitis. The current study isolated 33 strains belonging to 2 groups, outbreak and sporadic strains, in 13 cities of Shandong province, China from 2013 to 2019. Comparison of genetic characterization among 15 outbreak strains and 18 sporadic strains was performed by analyzing nucleotide sequences of the C-terminal region of N protein gene (N-450). RESULTS: All 33 stains belonged to genotype H1. The outbreak strains and sporadic strains distributed crossly in phylogenetic tree. Sequences alignment revealed some interesting G to A transversion which changed the amino acids on genomic sites 1317, 1422, and 1543. The nucleotide and amino acid similarities among outbreak isolates were 98-100% (0-10 nucleotide variations) and 97.7-100%, respectively; They were 97.3-100% and 96.6-100%, respectively for sporadic isolates. Evolutionary genetics analysis revealed that the mean evolution rates of outbreak and sporadic isolates were 1.26 N 10- 3 and 1.48 N 10- 3 substitutions per site per year separately, which were similar with corresponding data before 2012. Local transmission analysis suggested that there were three transmission chains in this study, two of them originated from Japan. Outbreak cases and sporadic cases emerged alternatively and were reciprocal causation on the transmission chains. CONCLUSIONS: Our study investigated the phylogeny and evolutional genetics of MeV during a 7-year surveillance, and compared epidemic and genetic characteristics of outbreak strains and sporadic strains. These results underscore the importance of evolutionary study alongside with sporadic cases in discovering and tracing possible outbreaks, especially in the stage of measles elimination.

BACKGROUND: Many long non-coding RNAs (lncRNAs) have been extensively identified in higher eukaryotic species. The function of lncRNAs has been reported to play important roles in diverse biological processes, including developmental regulation and behavioral plasticity. However, there are no reports of systematic characterization of long non-coding RNAs in the fire ant Solenopsis invicta. RESULTS: In this study, we performed a genome-wide analysis of lncRNAs in the brains of S. invicta from RNA-seq. In total, 1,393 novel lncRNA transcripts were identified in the fire ant. In contrast to the annotated lncRNA transcripts having at least two exons, novel lncRNAs are monoexonic transcripts with a shorter length. Besides, the transcriptome from virgin alate and dealate mated queens were analyzed and compared. The results showed 295 differentially expressed mRNA genes (DEGs) and 65 differentially expressed lncRNA genes (DELs) between virgin and mated queens, of which 17 lncRNAs were highly expressed in the virgin alates and 47 lncRNAs were highly expressed in the mated dealates. By identifying the DEL:DEG pairs with a high association in their expression (Spearman's |rho|> 0.8 and p-value < 0.01), many DELs were co-regulated with DEGs after mating. Furthermore, several remarkable lncRNAs (MSTRG.6523, MSTRG.588, and nc909) that were found to associate with particular coding genes may play important roles in the regulation of brain gene expression in reproductive transition in fire ants. CONCLUSION: This study provides the first genome-wide identification of S. invicta lncRNAs in the brains in different reproductive states. It will contribute to a fuller understanding of the transcriptional regulation underpinning reproductive changes.

BACKGROUND: Nucleosome positioning is the precise determination of the location of nucleosomes on DNA sequence. With the continuous advancement of biotechnology and computer technology, biological data is showing explosive growth. It is of practical significance to develop an efficient nucleosome positioning algorithm. Indeed, convolutional neural networks (CNN) can capture local features in DNA sequences, but ignore the order of bases. While the bidirectional recurrent neural network can make up for CNN's shortcomings in this regard and extract the long-term dependent features of DNA sequence. RESULTS: In this work, we use word vectors to represent DNA sequences and propose three new deep learning models for nucleosome positioning, and the integrative model NP_CBiR reaches a better prediction performance. The overall accuracies of NP_CBiR on H. sapiens, C. elegans, and D. melanogaster datasets are 86.18%, 89.39%, and 85.55% respectively. CONCLUSIONS: Benefited by different network structures, NP_CBiR can effectively extract local features and bases order features of DNA sequences, thus can be considered as a complementary tool for nucleosome positioning.

BACKGROUND: Crohn's disease is one of the two categories of inflammatory bowel diseases that affect the gastrointestinal tract. The heritability estimate has been reported to be 0.75. Several genes linked to Crohn's disease risk have been identified using a plethora of strategies such as linkage-based studies, candidate gene association studies, and lately through genome-wide association studies (GWAS). Nevertheless, to our knowledge, a compendium of all the genes that have been associated with CD is lacking. METHODS: We conducted functional analyses of a gene set generated from a systematic review where genes potentially related to CD found in the literature were analyzed and classified depending on the genetic evidence reported and putative biological function. For this, we retrieved and analyzed 2496 abstracts comprising 1067 human genes plus 22 publications regarding 133 genes from GWAS Catalog. Then, each gene was curated and categorized according to the type of evidence associated with Crohn's disease. RESULTS: We identified 126 genes associated with Crohn's disease risk by specific experiments. Additionally, 71 genes were recognized associated through GWAS alone, 18 to treatment response, 41 to disease complications, and 81 to related diseases. Bioinformatic analysis of the 126 genes supports their importance in Crohn's disease and highlights genes associated with specific aspects such as symptoms, drugs, and comorbidities. Importantly, most genes were not included in commercial genetic panels suggesting that Crohn's disease is genetically underdiagnosed. CONCLUSIONS: We identified a total of 126 genes from PubMed and 71 from GWAS that showed evidence of association to diagnosis, 18 to treatment response, and 41 to disease complications in Crohn's disease. This prioritized gene catalog can be explored at http://victortrevino.bioinformatics.mx/CrohnDisease .

BACKGROUND: mRNA sequencing is a powerful technique, which is used to investigate the transcriptome status of a gene of interest, such as its transcription level and splicing variants. Presently, several RNA sequencing (RNA-Seq) methods have been developed; however, the relative advantage of each method has remained unknown. Here we used three commercially available RNA-Seq library preparation kits; the traditional method (TruSeq), in addition to full-length double-stranded cDNA methods (SMARTer and TeloPrime) to investigate the advantages and disadvantages of these three approaches in transcriptome analysis. RESULTS: We observed that the number of expressed genes detected from the TeloPrime sequencing method was fewer than that obtained using the TruSeq and SMARTer. We also observed that the expression patterns between TruSeq and SMARTer correlated strongly. Alternatively, SMARTer and TeloPrime methods underestimated the expression of relatively long transcripts. Moreover, genes having low expression levels were undetected stochastically regardless of any three methods used. Furthermore, although TeloPrime detected a significantly higher proportion at the transcription start site (TSS), its coverage of the gene body was not uniform. SMARTer is proposed to be yielded for nonspecific genomic DNA amplification. In contrast, the detected splicing event number was highest in the TruSeq. The percent spliced in index (PSI) of the three methods was highly correlated. CONCLUSIONS: TruSeq detected transcripts and splicing events better than the other methods and measured expression levels of genes, in addition to splicing events accurately. However, although detected transcripts and splicing events in TeloPrime were fewer, the coverage at TSS was highest. Additionally, SMARTer was better than TeloPrime with regards to the detected number of transcripts and splicing events among the understudied full-length double-stranded cDNA methods. In conclusion, for short-read sequencing, TruSeq has relative advantages for use in transcriptome analysis.

BACKGROUND: The giant panda (Ailuropoda melanoleuca) is a threatened species endemic to China. Alopecia, characterized by thinning and broken hair, mostly occurs in breeding males. Alopecia significantly affects the health and public image of the giant panda and the cause of alopecia is unclear. RESULTS: Here, we researched gene expression profiles of four alopecia giant pandas and seven healthy giant pandas. All pandas were approximately ten years old and their blood samples collected during the breeding season. A total of 458 up-regulated DEGs and 211 down-regulated DEGs were identified. KEGG pathway enrichment identified that upregulated genes were enriched in the Notch signaling pathway and downregulated genes were enriched in ribosome, oxidative phosphorylation, and thermogenesis pathways. We obtained 28 hair growth-related DEGs, and identified three hub genes NOTCH1, SMAD3, and TGFB1 in PPI analysis. Five hair growth-related signaling pathways were identified with abnormal expression, these were Notch, Wnt, TGF-β, Mapk, and PI3K-Akt. The overexpression of NOTCH1 delays inner root sheath differentiation and results in hair shaft abnormalities. The delayed hair regression was associated with a significant decrease in the expression levels of TGFB1. CONCLUSIONS: Our data confirmed the abnormal expression of several hair-related genes and pathways and identified alopecia candidate genes in the giant panda. Results of this study provide theoretical basis for the establishment of prevention and treatment strategies for giant pandas with alopecia.

BACKGROUND: Acoustic communication is important for the survival and reproduction of anurans and masking background noise is a critical factor for their effective acoustic communication. Males of the concave-eared frog (Odorrana tormota) have evolved an ultrasonic communication capacity to avoid masking by the widespread background noise of local fast-flowing streams, whereas females exhibit no ultrasonic sensitivity. However, the molecular mechanisms underlying the high-frequency hearing differences between the sexes of O. tormota are still poorly understood. RESULTS: In this study, we sequenced the brain transcriptomes of male and female O. tormota, and compared their differential gene expression. A total of 4,605 differentially expressed genes (DEGs) between the sexes of O. tormota were identified and eleven of them were related to auditory based on the annotation and enrichment analysis. Most of these DEGs in males showed a higher expression trend than females in both quantity and expression quantity. The highly expressed genes in males were relatively concentrated in neurogenesis, signal transduction, ion transport and energy metabolism, whereas the up-expressed genes in females were mainly related to the growth and development regulation of specific auditory cells. CONCLUSIONS: The transcriptome of male and female O. tormota has been sequenced and de novo assembled, which will provide gene reference for further genomic studies. In addition, this is the first research to reveal the molecular mechanisms of sex differences in ultrasonic hearing between the sexes of O. tormota and will provide new insights into the genetic basis of the auditory adaptation in amphibians during their transition from water to land.

BACKGROUND: The scales serve as an ideal model for studying the regulatory mechanism of bone homeostasis in fish. To explore the effect of salinity acclimation on bone metabolism of juvenile rainbow trout (Oncorhynchus mykiss), three sampling time points during salinity acclimation (7D, 14D and 21D) were selected to detect variations in histological characteristics. In the histological analysis, osteoblast marker enzymes alkaline phosphatase (ALP), osteoclast marker tartrate-resistant acid phosphatase (TRAcP) and calcium salt deposit areas (Von Kossa's) were detected. Changes in calcium (Ca), phosphorus (P) and the molar mass ratio of calcium to phosphorus (Ca/P) in the scales were also detected by Inductively Coupled Plasma Mass Spectrometry (ICP-MS). In addition, the global MicroRNA (miRNA) expression profiles during salinity acclimation were examined using Illumina sequencing platform because of their important regulatory roles in teleost biological processes. RESULTS: Twelve independent miRNA libraries were constructed, a total of 664 known and 92 putative novel miRNAs were identified. A total of 290 differentially expressed (DE) miRNAs were found in clusters with significant trends in the cluster analysis, and five types of clustering patterns were obtained; 22,374 DE predicted target genes of the aforementioned 290 DE miRNAs were obtained, 5957 of which clustered in six types of clustering patterns with a significant trend. To better understand the functions of the DE miRNAs, GO and KEGG analysis was performed on the 5957 target genes, as a result, they were significantly enriched in bone metabolism related signaling pathways such as MAPK signaling pathway, Calcium signaling pathway, Wnt signaling pathway, Mineral absorption and NF-kappa B signaling pathway. Six DE miRNAs were randomly selected and their expression were verified by quantitative real-time PCR (qRT-PCR), the expression trends were consistent with the results of transcriptome sequencing. CONCLUSIONS: The DE miRNAs and DE target genes identified in this study might play an important role in regulation of bone metabolism during salinity acclimation, relative genes or pathways could serve as key candidates for further studies to elucidate molecular mechanism of teleost bone metabolism, and help performing salinity acclimation and developing marine culture of salmonid species.

BACKGROUND: Epigenome analysis relies on defined sets of genomic regions output by widely used assays such as ChIP-seq and ATAC-seq. Statistical analysis and visualization of genomic region sets is essential to answer biological questions in gene regulation. As the epigenomics community continues generating data, there will be an increasing need for software tools that can efficiently deal with more abundant and larger genomic region sets. Here, we introduce GenomicDistributions, an R package for fast and easy summarization and visualization of genomic region data. RESULTS: GenomicDistributions offers a broad selection of functions to calculate properties of genomic region sets, such as feature distances, genomic partition overlaps, and more. GenomicDistributions functions are meticulously optimized for best-in-class speed and generally outperform comparable functions in existing R packages. GenomicDistributions also offers plotting functions that produce editable ggplot objects. All GenomicDistributions functions follow a uniform naming scheme and can handle either single or multiple region set inputs. CONCLUSIONS: GenomicDistributions offers a fast and scalable tool for exploratory genomic region set analysis and visualization. GenomicDistributions excels in user-friendliness, flexibility of outputs, breadth of functions, and computational performance. GenomicDistributions is available from Bioconductor ( https://bioconductor.org/packages/release/bioc/html/GenomicDistributions.html ).

BACKGROUND: Recently genomic selection (GS) has emerged as an important tool for plant breeders to select superior genotypes. Multi-trait (MT) prediction model provides an opportunity to improve the predictive ability of expensive and labor-intensive traits. In this study, we assessed the potential use of a MT genomic prediction model by incorporating two physiological traits (canopy temperature, CT and normalized difference vegetation index, NDVI) to predict 5 complex primary traits (harvest index, HI; grain yield, GY; grain number, GN; spike partitioning index, SPI; fruiting efiiciency, FE) using two cross-validation schemes CV1 and CV2. RESULTS: In this study, we evaluated 236 wheat genotypes in two locations in 2 years. The wheat genotypes were genotyped with genotyping by sequencing approach which generated 27,466 SNPs. MT-CV2 (multi-trait cross validation 2) model improved predictive ability by 4.8 to 138.5% compared to ST-CV1(single-trait cross validation 1). However, the predictive ability of MT-CV1 was not significantly different compared to the ST-CV1 model. CONCLUSIONS: The study showed that the genomic prediction of complex traits such as HI, GN, and GY can be improved when correlated secondary traits (cheaper and easier phenotyping) are used. MT genomic selection could accelerate breeding cycles and improve genetic gain for complex traits in wheat and other crops.

BACKGROUND: Many transcription factors (TFs), such as multi zinc-finger (ZF) TFs, have multiple DNA binding domains (DBDs), and deciphering the DNA binding motifs of individual DBDs is a major challenge. One example of such a TF is CCCTC-binding factor (CTCF), a TF with eleven ZFs that plays a variety of roles in transcriptional regulation, most notably anchoring DNA loops. Previous studies found that CTCF ZFs 3-7 bind CTCF's core motif and ZFs 9-11 bind a specific upstream motif, but the motifs of ZFs 1-2 have yet to be identified. RESULTS: We developed a new approach to identifying the binding motifs of individual DBDs of a TF through analyzing chromatin immunoprecipitation sequencing (ChIP-seq) experiments in which a single DBD is mutated: we train a deep convolutional neural network to predict whether wild-type TF binding sites are preserved in the mutant TF dataset and interpret the model. We applied this approach to mouse CTCF ChIP-seq data and identified the known binding preferences of CTCF ZFs 3-11 as well as a putative GAG binding motif for ZF 1. We analyzed other CTCF datasets to provide additional evidence that ZF 1 is associated with binding at the motif we identified, and we found that the presence of the motif for ZF 1 is associated with CTCF ChIP-seq peak strength. CONCLUSIONS: Our approach can be applied to any TF for which in vivo binding data from both the wild-type and mutated versions of the TF are available, and our findings provide new potential insights binding preferences of CTCF's DBDs.

BACKGROUND: Long-term artificial selection for growth rate and lean meat rate has eventually led to meat quality deterioration. Muscle fiber type is a key factor that markedly affects meat quality. circRNAs have been reported to participate in diverse biological activities, including myofiber growth and development; thus, we herein compared porcine circRNA transcriptome between oxidative and glycolytic muscle tissues. RESULTS: Longissimus thoracis muscle tissues were obtained from Lantang and Landrace pigs at birth (LT1D and LW1D, respectively) and 90 postnatal days (LT90D and LW90D, respectively). Hematoxylin and eosin staining and quantitative real-time PCR revealed that all structural traits of the muscle showed large variations between different breeds and growth stages. In total, 329 known miRNAs and 42,081 transcript candidates were identified; 6,962 differentially expressed transcripts were found to play a key role in myogenesis by gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses. In addition, 3,352 circRNAs were identified using five predicting algorithms, and 104 circRNA candidates were differentially expressed. Integrated analysis of differentially expressed miRNAs, mRNAs, and circRNAs led to the identification of 777, 855, and 22 convincing ceRNA interactions in LT1D vs. LT90D, LW1D vs. LW90D, and LT90D vs. LW90D, respectively. Finally, we identified a circRNA candidate circKANSL1L, which showed high homology between mice and pigs, and it was found to inhibit the proliferation of C2C12 cells but promote their differentiation. CONCLUSIONS: We identified genome-wide circRNAs in 0- and 90-day-old Lantang and Landrace pigs by RNA-seq and found that circRNAs were abundant, differentially expressed, and associated with myogenesis. Our results should serve as a reference for future studies on pork quality.

BACKGROUND: Pigeon circovirus (PiCV) infections in pigeons (Columba livia) have been reported worldwide. Currently, pigeon racing is becoming increasingly popular and considered to be a national sport in China, and even, the greatest competitions of racing pigeons are taking place in China. However, there are still no epidemiologic data regarding PiCV infections among racing pigeons in China. The purpose of our study was to provide information of prevalence, genetic variation and evolution of PiCV from racing pigeons in China. RESULTS: To trace the prevalence, genetic variation and evolution of PiCV in sick and healthy racing pigeons, 622 samples were collected from 11 provinces or municipalities in China from 2016 to 2019. The results showed that the positive rate of PiCV was 19.3% (120/622) at the sample level and 59.0% (23/39) at the club level, thus suggesting that the virus was prevalent in Chinese racing pigeons. A sequence analysis revealed that the cap genes of the PiCV strains identified in our study displayed a high genetic diversity and shared nucleotide homologies of 71.9%-100% and amino acid homologies of 71.7%-100%. 28 and 36 unique amino acid substitutions were observed in the Cap and Rep proteins derived from our PiCV strains, respectively. A cladogram representation of PiCV strains phylogeny based on 90 cap gene sequences showed that the strains in this study could be further divided into seven clades (A, B, C, E, G, H, and I) and some of them were closely related to worldwide strains from different types of pigeons. A large number of recombination events (31 events) were also detected in the PiCV genomes from Chinese racing pigeons. CONCLUSIONS: These findings indicate that PiCV strains circulating in China exhibit a high genetic diversity and also contribute to information of prevalence, genetic variation and evolution of PiCV from racing pigeons in China.

BACKGROUND: The continued spread of SARS-CoV-2 and emergence of new variants with higher transmission rates and/or partial resistance to vaccines has further highlighted the need for large-scale testing and genomic surveillance. However, current diagnostic testing (e.g., PCR) and genomic surveillance methods (e.g., whole genome sequencing) are performed separately, thus limiting the detection and tracing of SARS-CoV-2 and emerging variants. RESULTS: Here, we developed DeepSARS, a high-throughput platform for simultaneous diagnostic detection and genomic surveillance of SARS-CoV-2 by the integration of molecular barcoding, targeted deep sequencing, and computational phylogenetics. DeepSARS enables highly sensitive viral detection, while also capturing genomic diversity and viral evolution. We show that DeepSARS can be rapidly adapted for identification of emerging variants, such as alpha, beta, gamma, and delta strains, and profile mutational changes at the population level. CONCLUSIONS: DeepSARS sets the foundation for quantitative diagnostics that capture viral evolution and diversity. DeepSARS uses molecular barcodes (BCs) and multiplexed targeted deep sequencing (NGS) to enable simultaneous diagnostic detection and genomic surveillance of SARS-CoV-2. Image was created using Biorender.com .

BACKGROUND: Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high. RESULTS: We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin. CONCLUSION: The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

BACKGROUND: Paspalum notatum 'Flugge' is a diploid with 20 chromosomes (2n = 20) multi-purpose subtropical herb native to South America and has a high ecological significance. It is currently widely planted in tropical and subtropical regions. Despite the gene pool of P. notatum 'Flugge' being unearthed to a large extent in the past decade, no details about the genomic information of relevant species in Paspalum have been reported. In this study, the complete genome information of P. notatum was established and annotated through sequencing and de novo assembly of its genome. RESULTS: The latest PacBio third-generation HiFi assembly and sequencing revealed that the genome size of P. notatum 'Flugge' is 541 M. The assembly result is the higher index among the genomes of the gramineous family published so far, with a contig N50 = 52Mbp, scaffold N50 = 49Mbp, and BUSCOs = 98.1%, accounting for 98.5% of the estimated genome. Genome annotation revealed 36,511 high-confidence gene models, thus providing an important resource for future molecular breeding and evolutionary research. A comparison of the genome annotation results of P. notatum 'Flugge' with other closely related species revealed that it had a close relationship with Zea mays but not close compared to Brachypodium distachyon, Setaria viridis, Oryza sativa, Puccinellia tenuiflora, Echinochloa crusgalli. An analysis of the expansion and contraction of gene families suggested that P. notatum 'Flugge' contains gene families associated with environmental resistance, increased reproductive ability, and molecular evolution, which explained its excellent agronomic traits. CONCLUSION: This study is the first to report the high-quality chromosome-scale-based genome of P. notatum 'Flugge' assembled using the latest PacBio third-generation HiFi sequencing reads. The study provides an excellent genetic resource bank for gramineous crops and invaluable perspectives regarding the evolution of gramineous plants.

BACKGROUND: VQ motif-containing (VQ) proteins are cofactors of transcriptional regulation that are widely involved in plant growth and development and respond to various stresses. The VQ gene family has been identified and characterized for many plants, but there is little research on VQ gene family proteins in wheat (Triticum aestivum L.). RESULTS: In this study, 113 TaVQ genes (40 homoeologous groups) were identified in the wheat genome. TaVQ proteins all contain the conserved motif FxxhVQxhTG, and most of the TaVQ genes do not contain introns. Phylogenetic analysis demonstrated that TaVQ proteins can be divided into 8 subgroups (I-VIII). The chromosomal location mapping analysis indicated that TaVQ genes are disproportionally distributed on 21 wheat chromosomes. Gene duplication analysis revealed that segmental duplication significantly contributes to the expansion of the TaVQ gene family. Gene expression analysis demonstrated that the expression pattern of TaVQ genes varies in different tissues. The results of quantitative real-time PCR (qRT-PCR) found that TaVQ genes displayed different expression levels under different phytohormones and abiotic stresses. The cis-elements analysis of the promoter region demonstrated that stress responses, hormone responses, growth and development, and WRKY binding elements are all widely distributed. Additionally, a potential regulatory network between TaVQ proteins and WRKY transcription factors was visualized. CONCLUSION: This study systematically analyzed the wheat TaVQ gene family, providing a reference for further functional characterization of TaVQ genes in wheat.

BACKGROUND: Circular RNA (circRNA), as a new members of noncoding RNA family, have vital functions in many biological processes by as microRNA sponges or competing endogenous RNAs (ceRNAs). However, little has been reported about the genetic mechanism of circRNAs regulation of feed efficiency in sheep. RESULTS: This study aimed to explore the expression of circRNAs in the liver of Hu sheep with High-RFI (High residual feed intake) and Low-RFI (Low residual feed intake) using transcriptome sequencing. A total of 20,729 circRNAs were identified in two groups, in which 219 circRNAs were found as significantly differentially expressed. Several circRNAs were validated by using RT-PCR, sanger sequencing and RT-qPCR methods. These results demonstrated that the RNA-seq result and expression level of circRNAs identified are reliable. Subsequently, GO and KEGG enrichment analysis of the parental genes of the differentially expressed (DE) circRNAs were mainly involved in immunity response and metabolic process. Finally, the ceRNA regulatory networks analysis showed that the target binding sites for miRNA such as novel_41, novel_115, novel_171 and oar-miR-485-3p in the identified DE cirRNAs. Importantly, two metabolic (SHISA3 and PLEKHH2) and four (RTP4, CD274, OAS1, and RFC3) immune-related target mRNAs were identified from 4 miRNAs. Association analysis showed that the polymorphism (RTP4 c.399 A > G) in the target gene RTP4 were significantly associated with RFI (P < 0.05). CONCLUSIONS: Analysis of sequencing data showed some candidate ceRNAs that may play key roles in the feed efficiency in sheep by regulating animal immune and metabolic. These results provide the basis data for further study of the biological functions of circRNAs in regulating sheep feed efficiency.

Cystic fibrosis (CF, MIM# 219,700) is an autosomal recessive disease caused by pathogenic variants within the CFTR gene. It was shown that genetic variants located in cis can affect disease severity or treatment response because of additive or epistatic effects. Studies on the prevalence of complex alleles in Russian CF patients have just begun. Aim To evaluate frequencies and genetic background of complex alleles carrying c.1521_1523delCTT (F508del) and c.1399C>T (L467F), c.2562T>G (T854=) or c.4389G>A (Q1463=) in cis; to determine clinical consequences of complex allele c.[1399C>T;1521_1523delCTT] ([L467;F508del]) in Russian CF patients. Methods Sequencing of coding regions of CFTR gene and analysis of polymorphic markers in CF patients carrying F508del variant. Comparing of clinical features in two groups patients having genotypes [L467F;F508del];[F508del] (group 1) and [F508del];[F508del] (group 2). Results Frequency of [L467F;F508del] allele linked to 2-2-21-6-17-13 haplotype was 4.42%, of [F508del;T854=;Q1463=] allele linked to haplotype 1-2-21-6-17-13 - 2.2% in F508del chromosomes. No differences in disease severity in patients carrying complex allele [L467F;F508del] and patients homozygous for F508del was found. Conclusion The frequency of complex alleles associated with F508del was at least 6.6% in Russian CF patients, which should be taken into account for the decision on optimal treatment options with CFTR modulators.

BACKGROUND: When plants are subjected to cold stress, they undergo a series of molecular and physiological changes to protect themselves from injury. Indica cultivars can usually withstand only mild cold stress in a relatively short period. Hormone-mediated defence response plays an important role in cold stress. Weighted gene co-expression network analysis (WGCNA) is a very useful tool for studying the correlation between genes, identifying modules with high phenotype correlation, and identifying Hub genes in different modules. Many studies have elucidated the molecular mechanisms of cold tolerance in different plants, but little information about the recovery process after cold stress is available. RESULTS: To understand the molecular mechanism of cold tolerance in rice, we performed comprehensive transcriptome analyses during cold treatment and recovery stage in two cultivars of near-isogenic lines (9311 and DC907). Twelve transcriptomes in two rice cultivars were determined. A total of 2509 new genes were predicted by fragment splicing and assembly, and 7506 differentially expressed genes were identified by pairwise comparison. A total of 26 modules were obtained by expression-network analysis, 12 of which were highly correlated with cold stress or recovery treatment. We further identified candidate Hub genes associated with specific modules and analysed their regulatory relationships based on coexpression data. Results showed that various plant-hormone regulatory genes acted together to protect plants from physiological damage under short-term low-temperature stress. We speculated that this may be common in rice. Under long-term cold stress, rice improved the tolerance to low-temperature stress by promoting autophagy, sugar synthesis, and metabolism. CONCLUSION: Through WGCNA analysis at the transcriptome level, we provided a potential regulatory mechanism for the cold stress and recovery of rice cultivars and identified candidate central genes. Our findings provided an important reference for the future cultivation of rice strains with good tolerance.