BACKGROUND: To understand the dynamics of infectious diseases, genomic epidemiology is increasingly advocated, with a need for rapid generation of genetic sequences during outbreaks for public health decision making. Here, we explore the use of metagenomic sequencing compared to specific amplicon- and capture-based sequencing, both on the Nanopore and the Illumina platform for generation of whole genomes of Usutu virus, Zika virus, West Nile virus, and Yellow Fever virus. RESULTS: We show that amplicon-based Nanopore sequencing can be used to rapidly obtain whole genome sequences in samples with a viral load up to Ct 33 and capture-based Illumina is the most sensitive method for initial virus determination. CONCLUSIONS: The choice of sequencing approach and platform is important for laboratories wishing to start whole genome sequencing. Depending on the purpose of genome sequencing the best choice can differ. The insights presented in this work and the shown differences in data characteristics can guide labs to make a well informed choice.

BACKGROUND: The genus Verbascum L. (Scrophulariaceae) is distributed in Africa, Europe, and parts of Asia, with the Mediterranean having the most species variety. Several researchers have already worked on the phylogenetic and taxonomic analysis of Verbascum by using ITS data and chloroplast genome fragments and have produced different conclusions. The taxonomy and phylogenetic relationships of this genus are unclear. RESULTS: The complete plastomes (cp) lengths for V. chaixii, V. songaricum, V. phoeniceum, V. blattaria, V. sinaiticum, V. thapsus, and V. brevipedicellatum ranged from 153,014 to 153,481 bp. The cp coded 114 unique genes comprising of 80 protein-coding genes, four ribosomal RNA (rRNA), and 30 tRNA genes. We detected variations in the repeat structures, gene expansion on the inverted repeat, and single copy (IR/SC) boundary regions. The substitution rate analysis indicated that some genes were under purifying selection pressure. Phylogenetic analysis supported the sister relationship of (Lentibulariaceae + Acanthaceae + Bignoniaceae + Verbenaceae + Pedaliaceae) and (Lamiaceae + Phyrymaceae + Orobanchaceae + Paulowniaceae + Mazaceae) in Lamiales. Within Scrophulariaceae, Verbascum was sister to Scrophularia, while Buddleja formed a monophyletic clade from (Scrophularia + Verbascum) with high bootstrap support values. The relationship of the nine species within Verbascum was highly supported. CONCLUSION: Based on the phylogenetic results, we proposed to reinstate the species status of V. brevipedicellatum (Engl.) Hub.-Mor. Additionally, three genera (Mazus, Lancea, and Dodartia) placed in the Phyrymaceae family formed a separate clade within Lamiaceae. The classification of the three genera was supported by previous studies. Thus, the current study also suggests the circumscription of these genera as documented previously to be reinstated. The divergence time of Lamiales was approximated to be 86.28 million years ago (Ma) (95% highest posterior density (HPD), 85.12-89.91 Ma). The complete plastomes sequence data of the Verbascum species will be important for understanding the Verbascum phylogenetic relationships and evolution in order Lamiales.

BACKGROUND: Acer pseudosieboldianum is a kind of excellent color-leafed plants, and well known for its red leaves in autumn. At the same time, A. pseudosieboldianum is one of the native tree species in the northeast of China, and it plays an important role in improving the lack of color-leafed plants in the north. In previous study, we found a mutant of the A. pseudosieboldianum that leaves intersect red and green in spring and summer. However, it is unclear which genes cause the color change of mutant leaves. RESULTS: In order to study the molecular mechanism of leaf color formation, we analyzed the leaves of the mutant group and the control group from A. pseudosieboldianum by RNA deep sequencing in this study. Using an Illumina sequencing platform, we obtained approximately 276,071,634 clean reads. After the sequences were filtered and assembled, the transcriptome data generated a total of 70,014 transcripts and 54,776 unigenes, of which 34,486 (62.96%) were successfully annotated in seven public databases. There were 8,609 significant DEGs identified between the control and mutant groups, including 4,897 upregulated and 3,712 downregulated genes. We identified 13 genes of DEGs for leaf color synthesis that was involved in the flavonoid pathway, 26 genes that encoded transcription factors, and eight genes associated with flavonoid transport. CONCLUSION: Our results provided comprehensive gene expression information about A. pseudosieboldianum transcriptome, and directed the further study of accumulation of anthocyanin in A. pseudosieboldianum, aiming to provide insights into leaf coloring of it through transcriptome sequencing and analysis.

BACKGROUND: Brassinosteroids (BRs) are a group of essential steroid hormones involved in diverse developmental and physiological processes in plants. The Brassinazole-resistant 1 (BZR1) transcription factors are key components of BR signaling and integrate a wide range of internal and environmental signals to coordinate plant development, growth, and resistance to abiotic and biotic stresses. Although the BZR1 family has been fully studied in Arabidopsis, celery BZR1 family genes remain largely unknown. RESULTS: Nine BZR1 genes were identified in the celery genome, and categorized into four classes based on phylogenetic and gene structure analyses. All the BZR1 proteins shared a typical bHLH (basic helix-loop-helix) domain that is highly conserved across the whole family in Arabidopsis, grape, lettuce, ginseng, and three Apiaceae species. Both duplications and losses of the BZR1 gene family were detected during the shaping of the celery genome. Whole-genome duplication (WGD) or segmental duplication contributed 55.56% of the BZR1 genes expansion, and the γ as well as celery-ω polyploidization events made a considerable contribution to the production of the BZR1 paralogs in celery. Four AgBZR1 members (AgBZR1.1, AgBZR1.3, AgBZR1.5, and AgBZR1.9), which were localized both in the nucleus and cytoplasm, exhibit transcription activation activity in yeast. AgBZR1.5 overexpression transgenic plants in Arabidopsis showed curled leaves with bent, long petioles and constitutive BR-responsive phenotypes. Furthermore, the AgBZR1 genes possessed divergent expression patterns with some overlaps in roots, petioles, and leaves, suggesting an extensive involvement of AgBZR1s in the developmental processes in celery with both functional redundancy and divergence. CONCLUSIONS: Our results not only demonstrated that AgBZR1 played a conserved role in BR signaling but also suggested that AgBZR1 might be extensively involved in plant developmental processes in celery. The findings lay the foundation for further study on the molecular mechanism of the AgBZR1s in regulating the agronomic traits and environmental adaptation of celery, and provide insights for future BR-related genetic breeding of celery and other Apiaceae crops.

BACKGROUND: Shrimp have the ability to accommodate viruses in long term, persistent infections without signs of disease. Endogenous viral elements (EVE) play a role in this process probably via production of negative-sense Piwi-interacting RNA (piRNA)-like fragments. These bind with Piwi proteins to dampen viral replication via the RNA interference (RNAi) pathway. We searched a genome sequence (GenBank record JABERT000000000) of the giant tiger shrimp (Penaeus monodon for the presence of EVE related to a shrimp parvovirus originally named infectious hypodermal and hematopoietic necrosis virus (IHHNV). RESULTS: The shrimp genome sequence contained three piRNA-like gene clusters containing scrambled IHHNV EVE. Two clusters were located distant from one another in pseudochromosome 35 (PC35). Both PC35 clusters contained multiple sequences with high homology (99%) to GenBank records DQ228358 and EU675312 that were both called "non-infectious IHHNV Type A" (IHHNV-A) when originally discovered. However, our results and those from a recent Australian P. monodon genome assembly indicate that the relevant GenBank records for IHHNV-A are sequence-assembly artifacts derived from scrambled and fragmental IHHNV-EVE. Although the EVE in the two PC35 clusters showed high homology only to IHHNV-A, the clusters were separate and distinct with respect to the arrangement (i.e., order and reading direction) and proportional content of the IHHNV-A GenBank records. We conjecture that these 2 clusters may constitute independent allele-like clusters on a pair of homologous chromosomes. The third EVE cluster was found in pseudochromosome 7 (PC7). It contained EVE with high homology (99%) only to GenBank record AF218266 with the potential to protect shrimp against current types of infectious IHHNV. One disadvantage was that some EVE in PC7 can give false positive PCR test results for infectious IHHNV. CONCLUSIONS: Our results suggested the possibility of viral-type specificity in EVE clusters. Specificity is important because whole EVE clusters for one viral type would be transmitted to offspring as collective hereditary units. This would be advantageous if one or more of the EVE within the cluster were protective against the disease caused by the cognate virus. It would also facilitate gene editing for removal of non-protective EVE clusters or for transfer of protective EVE clusters to genetically improve existing shrimp breeding stocks that might lack them.

BACKGROUND: Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes. RESULTS: Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation. CONCLUSIONS: The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes.

Erratum in BMC Genomics. 2022 Aug 25;23(1):615.

BACKGROUND: Type 2C protein phosphatase (PP2C) is a negative regulator of ABA signaling pathway, which plays important roles in stress signal transduction in plants. However, little research on the PP2C genes family of cucumber (Cucumis sativus L.), as an important economic vegetable, has been conducted. RESULTS: This study conducted a genome-wide investigation of the CsPP2C gene family. Through bioinformatics analysis, 56 CsPP2C genes were identified in cucumber. Based on phylogenetic analysis, the PP2C genes of cucumber and Arabidopsis were divided into 13 groups. Gene structure and conserved motif analysis showed that CsPP2C genes in the same group had similar gene structure and conserved domains. Collinearity analysis showed that segmental duplication events played a key role in the expansion of the cucumber PP2C genes family. In addition, the expression of CsPP2Cs under different abiotic treatments was analyzed by qRT-PCR. The results reveal that CsPP2C family genes showed different expression patterns under ABA, drought, salt, and cold treatment, and that CsPP2C3, 11-17, 23, 45, 54 and 55 responded significantly to the four stresses. By predicting the cis-elements in the promoter, we found that all CsPP2C members contained ABA response elements and drought response elements. Additionally, the expression patterns of CsPP2C genes were specific in different tissues. CONCLUSIONS: The results of this study provide a reference for the genome-wide identification of the PP2C gene family in other species and provide a basis for future studies on the function of PP2C genes in cucumber.

BACKGROUND: Mycobacterium tuberculosis presents several lineages each with distinct characteristics of evolutionary status, transmissibility, drug resistance, host interaction, latency, and vaccine efficacy. Whole genome sequencing (WGS) has emerged as a new diagnostic tool to reliably inform the occurrence of phylogenetic lineages of Mycobacterium tuberculosis and examine their relationship with patient demographic characteristics and multidrug-resistance development. METHODS: 191 Mycobacterium tuberculosis isolates obtained from a 2017/2018 Tanzanian drug resistance survey were sequenced on the Illumina Miseq platform at Supranational Tuberculosis Reference Laboratory in Uganda. Obtained fast-q files were imported into tools for resistance profiling and lineage inference (Kvarq v0.12.2, Mykrobe v0.8.1 and TBprofiler v3.0.5). Additionally for phylogenetic tree construction, RaxML-NG v1.0.3(25) was used to generate a maximum likelihood phylogeny with 800 bootstrap replicates. The resulting trees were plotted, annotated and visualized using ggtree v2.0.4 RESULTS: Most [172(90.0%)] of the isolates were from newly treated Pulmonary TB patients. Coinfection with HIV was observed in 33(17.3%) TB patients. Of the 191 isolates, 22(11.5%) were resistant to one or more commonly used first line anti-TB drugs (FLD), 9(4.7%) isolates were MDR-TB while 3(1.6%) were resistant to all the drugs. Of the 24 isolates with any resistance conferring mutations, 13(54.2%) and 10(41.6%) had mutations in genes associated with resistance to INH and RIF respectively. The findings also show four major lineages i.e. Lineage 3[81 (42.4%)], followed by Lineage 4 [74 (38.7%)], the Lineage 1 [23 (12.0%)] and Lineages 2 [13 (6.8%)] circulaing in Tanzania. CONCLUSION: The findings in this study show that Lineage 3 is the most prevalent lineage in Tanzania whereas drug resistant mutations were more frequent among isolates that belonged to Lineage 4.

BACKGROUND: Cyclic nucleotide-gated ion channels (CNGCs) are calcium-permeable channels that participate in a variety of biological functions, such as signaling pathways, plant development, and environmental stress and stimulus responses. Nevertheless, there have been few studies on CNGC gene family in cotton. RESULTS: In this study, a total of 114 CNGC genes were identified from the genomes of 4 cotton species. These genes clustered into 5 main groups: I, II, III, IVa, and IVb. Gene structure and protein motif analysis showed that CNGCs on the same branch were highly conserved. In addition, collinearity analysis showed that the CNGC gene family had expanded mainly by whole-genome duplication (WGD). Promoter analysis of the GhCNGCs showed that there were a large number of cis-acting elements related to abscisic acid (ABA). Combination of transcriptome data and the results of quantitative RT-PCR (qRT-PCR) analysis revealed that some GhCNGC genes were induced in response to salt and drought stress and to exogenous ABA. Virus-induced gene silencing (VIGS) experiments showed that the silencing of the GhCNGC32 and GhCNGC35 genes decreased the salt tolerance of cotton plants (TRV:00). Specifically, physiological indexes showed that the malondialdehyde (MDA) content in gene-silenced plants (TRV:GhCNGC32 and TRV:GhCNGC35) increased significantly under salt stress but that the peroxidase (POD) activity decreased. After salt stress, the expression level of ABA-related genes increased significantly, indicating that salt stress can trigger the ABA signal regulatory mechanism. CONCLUSIONS: we comprehensively analyzed CNGC genes in four cotton species, and found that GhCNGC32 and GhCNGC35 genes play an important role in cotton salt tolerance. These results laid a foundation for the subsequent study of the involvement of cotton CNGC genes in salt tolerance.

BACKGROUND: Plants synthesize metabolites to adapt to a continuously changing environment. Metabolite biosynthesis often occurs in response to the tissue-specific combinatorial developmental cues that are transcriptionally regulated. Polyphyllins are the major bioactive components in Paris species that demonstrate hemostatic, anti-inflammatory and antitumor effects and have considerable market demands. However, the mechanisms underlying polyphyllin biosynthesis and regulation during plant development have not been fully elucidated. RESULTS: Tissue samples of P. polyphylla var. yunnanensis during the four dominant developmental stages were collected and investigated using high-performance liquid chromatography and RNA sequencing. Polyphyllin concentrations in the different tissues were found to be highly dynamic across developmental stages. Specifically, decreasing trends in polyphyllin concentration were observed in the aerial vegetative tissues, whereas an increasing trend was observed in the rhizomes. Consistent with the aforementioned polyphyllin concentration trends, different patterns of spatiotemporal gene expression in the vegetative tissues were found to be closely related with polyphyllin biosynthesis. Additionally, molecular dissection of the pathway components revealed 137 candidate genes involved in the upstream pathway of polyphyllin backbone biosynthesis. Furthermore, gene co-expression network analysis revealed 74 transcription factor genes and one transporter gene associated with polyphyllin biosynthesis and allocation. CONCLUSIONS: Our findings outline the framework for understanding the biosynthesis and accumulation of polyphyllins during plant development and contribute to future research in elucidating the molecular mechanism underlying polyphyllin regulation and accumulation in P. polyphylla.

BACKGROUND: Protein-protein interaction (PPI) is very important for many biochemical processes. Therefore, accurate prediction of PPI can help us better understand the role of proteins in biochemical processes. Although there are many methods to predict PPI in biology, they are time-consuming and lack accuracy, so it is necessary to build an efficiently and accurately computational model in the field of PPI prediction. RESULTS: We present a novel sequence-based computational approach called DCSE (Double-Channel-Siamese-Ensemble) to predict potential PPI. In the encoding layer, we treat each amino acid as a word, and map it into an N-dimensional vector. In the feature extraction layer, we extract features from local and global perspectives by Multilayer Convolutional Neural Network (MCN) and Multilayer Bidirectional Gated Recurrent Unit with Convolutional Neural Networks (MBC). Finally, the output of the feature extraction layer is then fed into the prediction layer to output whether the input protein pair will interact each other. The MCN and MBC are siamese and ensemble based network, which can effectively improve the performance of the model. In order to demonstrate our model's performance, we compare it with four machine learning based and three deep learning based models. The results show that our method outperforms other models in all evaluation criteria. The Accuracy, Precision, [Formula: see text], Recall and MCC of our model are 0.9303, 0.9091, 0.9268, 0.9452, 0.8609. For the other seven models, the highest Accuracy, Precision, [Formula: see text], Recall and MCC are 0.9288, 0.9243, 0.9246, 0.9250, 0.8572. We also test our model in the imbalanced dataset and transfer our model to another species. The results show our model is excellent. CONCLUSION: Our model achieves the best performance by comparing it with seven other models. NLP-based coding method has a good effect on PPI prediction task. MCN and MBC extract protein sequence features from local and global perspectives and these two feature extraction layers are based on siamese and ensemble network structures. Siamese-based network structure can keep the features consistent and ensemble based network structure can effectively improve the accuracy of the model.

BACKGROUND: Advancements in genomic sequencing continually improve personalized medicine, and recent breakthroughs generate multimodal data on a cellular level. We introduce MOSCATO, a technique for selecting features across multimodal single-cell datasets that relate to clinical outcomes. We summarize the single-cell data using tensors and perform regularized tensor regression to return clinically-associated variable sets for each 'omic' type. RESULTS: Robustness was assessed over simulations based on available single-cell simulation methods, and applicability was assessed through an example using CITE-seq data to detect genes associated with leukemia. We find that MOSCATO performs favorably in selecting network features while also shown to be applicable to real multimodal single-cell data. CONCLUSIONS: MOSCATO is a useful analytical technique for supervised feature selection in multimodal single-cell data. The flexibility of our approach enables future extensions on distributional assumptions and covariate adjustments.

BACKGROUND: There is a need to investigate mechanisms of phenotypic plasticity in marine invertebrates as negative effects of climate change, like ocean acidification, are experienced by coastal ecosystems. Environmentally-induced changes to the methylome may regulate gene expression, but methylome responses can be species- and tissue-specific. Tissue-specificity has implications for gonad tissue, as gonad-specific methylation patterns may be inherited by offspring. We used the Pacific oyster (Crassostrea gigas) - a model for understanding pH impacts on bivalve molecular physiology due to its genomic resources and importance in global aquaculture- to assess how low pH could impact the gonad methylome. Oysters were exposed to either low pH (7.31 ± 0.02) or ambient pH (7.82 ± 0.02) conditions for 7 weeks. Whole genome bisulfite sequencing was used to identify methylated regions in female oyster gonad samples. C- > T single nucleotide polymorphisms were identified and removed to ensure accurate methylation characterization. RESULTS: Analysis of gonad methylomes revealed a total of 1284 differentially methylated loci (DML) found primarily in genes, with several genes containing multiple DML. Gene ontologies for genes containing DML were involved in development and stress response, suggesting methylation may promote gonad growth homeostasis in low pH conditions. Additionally, several of these genes were associated with cytoskeletal structure regulation, metabolism, and protein ubiquitination - commonly-observed responses to ocean acidification. Comparison of these DML with other Crassostrea spp. exposed to ocean acidification demonstrates that similar pathways, but not identical genes, are impacted by methylation. CONCLUSIONS: Our work suggests DNA methylation may have a regulatory role in gonad and larval development, which would shape adult and offspring responses to low pH stress. Combined with existing molluscan methylome research, our work further supports the need for tissue- and species-specific studies to understand the potential regulatory role of DNA methylation.

BACKGROUND: Macrobrachium nipponense is an economically and ecologically important freshwater prawn that is widely farmed in China. In contrast to other species of marine shrimp, M. nipponense has a short sexual maturity period, resulting in not only high stocking densities, but also a reduced survival rate and increased risk of hypoxia. Therefore, there is an urgent need to study the molecular mechanisms underlying fast ovary maturation in this species. RESULTS: Comparative transcriptome analysis was performed using hepatopancreatic tissue from female M. nipponense across five ovarian maturation stages to explore differentially expressed genes and pathways involved in ovarian maturation. In total, 118.01 Gb of data were generated from 15 transcriptomes. Approximately 90.46% of clean reads were mapped from the M. nipponense reference genome. A comprehensive comparative analysis between successive ovarian maturation stages generated 230-5814 differentially expressed genes. Gene Ontology (GO) enrichment was highly concentrated in the "biological process" category in all four comparison groups, and mainly focused on energy synthesis and accumulation, energy decomposition and transport. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results showed that, among 20 significantly enriched KEGG pathways, nine were involved in the synthesis, degradation, and metabolism of carbohydrates, lipids, and other nutrient intermediates, suggesting that the hepatopancreas has an important role in energy supply during ovarian maturation. Furthermore, the "Insect hormone biosynthesis" pathway was found to have a dominant role in the development of the ovary from immaturity to maturity, supporting the hypothesis that ecdysteroid- and juvenile hormone-signaling pathways have an important role in hepatopancreas regulation of ovarian maturation. CONCLUSION: Taken together, this study sheds light on the role of the hepatopancreas in the molecular regulation of ovary maturation in M. nipponense. The present study provided new insights for understanding the mechanisms of reproductive regulation in crustaceans.

BACKGROUND: Balanites aegyptiaca (L.) Delile, commonly known as desert date, is a thorny evergreen tree belonging to the family Zygophyllaceae and subfamily Tribuloideae that is widespread in arid and semiarid regions. This plant is an important source of food and medicines and plays an important role in conservation strategies for restoring degraded desert ecosystems. RESULTS: In the present study, we sequenced the complete plastome of B. aegyptiaca. The chloroplast genome was 155,800 bp, with a typical four-region structure: a large single copy (LSC) region of 86,562 bp, a small single copy (SSC) region of 18,102 bp, and inverted repeat regions (IRa and IRb) of 25,568 bp each. The GC content was 35.5%. The chloroplast genome of B. aegyptiaca contains 107 genes, 75 of which coding proteins, 28 coding tRNA, and 4 coding rRNA. We did not observe a large loss in plastid genes or a reduction in the genome size in B. aegyptiaca, as found previously in some species belonging to the family Zygophyllaceae. However, we noticed a divergence in the location of certain genes at the IR-LSC and IR-SSC boundaries and loss of ndh genes relative to other species. Furthermore, the phylogenetic tree constructed from the complete chloroplast genome data broadly supported the taxonomic classification of B. aegyptiaca as belonging to the Zygophyllaceae family. The plastome of B. aegyptiaca was found to be rich in single sequence repeats (SSRs), with a total of 240 SSRs. CONCLUSIONS: The genomic data available from this study could be useful for developing molecular markers to evaluate population structure, investigate genetic variation, and improve production programs for B. aegyptiaca. Furthermore, the current data will support future investigation of the evolution of the family Zygophyllaceae.

The pathophysiological effects of a number of metabolic and age-related disorders can be prevented to some extent by exercise and increased physical activity. However, the molecular mechanisms that contribute to the beneficial effects of muscle activity remain poorly explored. Availability of a fast, inexpensive, and genetically tractable model system for muscle activity and exercise will allow the rapid identification and characterization of molecular mechanisms that mediate the beneficial effects of exercise. Here, we report the development and characterization of an optogenetically-inducible muscle contraction (OMC) model in Drosophila larvae that we used to study acute exercise-like physiological responses. To characterize muscle-specific transcriptional responses to acute exercise, we performed bulk mRNA-sequencing, revealing striking similarities between acute exercise-induced genes in flies and those previously identified in humans. Our larval muscle contraction model opens a path for rapid identification and characterization of exercise-induced factors.

BACKGROUND: Selection of optimal computational strategies for analyzing metagenomics data is a decisive step in determining the microbial composition of a sample, and this procedure is complex because of the numerous tools currently available. The aim of this research was to summarize the results of crowdsourced sbv IMPROVER Microbiomics Challenge designed to evaluate the performance of off-the-shelf metagenomics software as well as to investigate the robustness of these results by the extended post-challenge analysis. In total 21 off-the-shelf taxonomic metagenome profiling pipelines were benchmarked for their capacity to identify the microbiome composition at various taxon levels across 104 shotgun metagenomics datasets of bacterial genomes (representative of various microbiome samples) from public databases. Performance was determined by comparing predicted taxonomy profiles with the gold standard. RESULTS: Most taxonomic profilers performed homogeneously well at the phylum level but generated intermediate and heterogeneous scores at the genus and species levels, respectively. kmer-based pipelines using Kraken with and without Bracken or using CLARK-S performed best overall, but they exhibited lower precision than the two marker-gene-based methods MetaPhlAn and mOTU. Filtering out the 1% least abundance species-which were not reliably predicted-helped increase the performance of most profilers by increasing precision but at the cost of recall. However, the use of adaptive filtering thresholds determined from the sample's Shannon index increased the performance of most kmer-based profilers while mitigating the tradeoff between precision and recall. CONCLUSIONS: kmer-based metagenomic pipelines using Kraken/Bracken or CLARK-S performed most robustly across a large variety of microbiome datasets. Removing non-reliably predicted low-abundance species by using diversity-dependent adaptive filtering thresholds further enhanced the performance of these tools. This work demonstrates the applicability of computational pipelines for accurately determining taxonomic profiles in clinical and environmental contexts and exemplifies the power of crowdsourcing for unbiased evaluation.

BACKGROUND: Equid gammaherpesvirus 2 (EHV2) is a gammaherpesvirus with a widespread distribution in horse populations globally. Although its pathogenic significance can be unclear in most cases of infection, EHV2 infection can cause upper respiratory tract disease in foals. Co-infection of different strains of EHV2 in an individual horse is common. Small regions of the EHV2 genome have shown considerable genetic heterogeneity. This could suggest genomic recombination between different strains of EHV2, similar to the extensive recombination networks that have been demonstrated for some alphaherpesviruses. This study examined natural recombination and genome diversity of EHV2 field isolates. RESULTS: Whole genome sequencing analysis of 18 EHV2 isolates, along with analysis of two publicly available EHV2 genomes, revealed variation in genomes sizes (from 173.7 to 184.8 kbp), guanine plus cytosine content (from 56.7 to 57.8%) and the size of the terminal repeat regions (from 17,196 to 17,551 bp). The nucleotide sequence identity between the genomes ranged from 86.2 to 99.7%. The estimated average inter-strain nucleotide diversity between the 20 EHV2 genomes was 2.9%. Individual gene sequences showed varying levels of nucleotide diversity and ranged between 0 and 38.1%. The ratio of nonsynonymous substitutions, Ka, to synonymous substitutions, Ks, (Ka/Ks) suggests that over 50% of EHV2 genes are undergoing diversifying selection. Recombination analyses of the 20 EHV2 genome sequences using the recombination detection program (RDP4) and SplitsTree revealed evidence of viral recombination. CONCLUSIONS: Analysis of the 18 new EHV2 genomes alongside the 2 previously sequenced genomes revealed a high degree of genetic diversity and extensive recombination networks. Herpesvirus genome diversification and virus evolution can be driven by recombination, and our findings are consistent with recombination being a key mechanism by which EHV2 genomes may vary and evolve.

BACKGROUND: RNA preparations contaminated with genomic DNA (gDNA) are frequently disregarded by RNA-seq studies. Such contamination may generate false results; however, their effect on the outcomes of RNA-seq analyses is unknown. To address this gap in our knowledge, here we added different concentrations of gDNA to total RNA preparations and subjected them to RNA-seq analysis. RESULTS: We found that the contaminating gDNA altered the quantification of transcripts at relatively high concentrations. Differentially expressed genes (DEGs) resulting from gDNA contamination may therefore contribute to higher rates of false enrichment of pathways compared with analogous samples lacking numerous DEGs. A strategy was developed to correct gene expression levels in gDNA-contaminated RNA samples, which assessed the magnitude of contamination to improve the reliability of the results. CONCLUSIONS: Our study indicates that caution must be exercised when interpreting results associated with low-abundance transcripts. The data provided here will likely serve as a valuable resource to evaluate the influence of gDNA contamination on RNA-seq analysis, particularly related to the detection of putative novel gene elements.

BACKGROUND: In cold regions, low temperature is the main limiting factor affecting grape production. As an important breeding resource, V. amurensis Rupr. has played a crucial role in the discovery of genes which confer cold resistance in grapes. Thus far, many cold-resistance genes have been reported based on the study of V. amurensis. In order to identify more candidate genes related to cold resistance in V. amurensis, QTL mapping and RNA-seq was conducted based on the hybrid population and different cold-resistance cultivars in this study. RESULTS: In this study, highly cold-resistant grape cultivar 'Shuangyou' (SY) which belongs to V. amurensis, and cold-sensitive cultivar 'Red Globe' (RG) which belongs to Vitis vinifera L. were used to identify cold resistance genes. Cold-resistance quantitative trait locus (QTL) mapping was performed based on genetic population construction through interspecific crossing of 'Shuangyou' and 'Red Globe'. Additionally, transcriptome analysis was conducted for the dormant buds of these two cultivars at different periods. Based on transcriptome analysis and QTL mapping, many new structural genes and transcription factors which relate to V. amurensis cold resistance were discovered, including CORs (VaCOR413IM), GSTs (VaGST-APIC, VaGST-PARB, VaGSTF9 and VaGSTF13), ARFs (VaIAA27 and VaSAUR71), ERFs (VaAIL1), MYBs (VaMYBR2, VaMYBLL and VaMYB3R-1) and bHLHs (VaICE1 and VabHLH30). CONCLUSIONS: This discovery of candidate cold-resistance genes will provide an important theoretical reference for grape cold-resistance mechanisms, research, and cold-resistant grape cultivar breeding in the future.

Tetrodotoxin (TTX) is a deadly neurotoxin and usually accumulates in large amounts in the ovaries but is non-toxic or low toxic in the testis of pufferfish. The molecular mechanism underlying sexual dimorphism accumulation of TTX in ovary and testis, and the relationship between TTX accumulation with sex related genes expression remain largely unknown. The present study investigated the effects of exogenous TTX treatment on Takifugu flavidus. The results demonstrated that exogenous TTX administration significantly incresed level of TTX concentration in kidney, cholecyst, skin, liver, heart, muscle, ovary and testis of the treatment group (TG) than that of the control group (CG). Transcriptome sequencing and analysis were performed to study differential expression profiles of mRNA and piRNA after TTX administration of the ovary and testis. The results showed that compared with female control group (FCG) and male control group (MCG), TTX administration resulted in 80 and 23 piRNAs, 126 and 223 genes up and down regulated expression in female TTX-treated group (FTG), meanwhile, 286 and 223 piRNAs, 2 and 443 genes up and down regulated expression in male TTX-treated group (MTG). The female dominant genes cyp19a1, gdf9 and foxl2 were found to be up-regulated in MTG. The cyp19a1, whose corresponding target piRNA uniq_554482 was identified as down-regulated in the MTG, indicating the gene expression feminization in testis after exogenous TTX administration. The KEGG enrichment analysis revealed that differentially expressed genes (DEGs) and piRNAs (DEpiRNAs) in MTG vs MCG group were more enriched in metabolism pathways, indicating that the testis produced more metabolic pathways in response to exogenous TTX, which might be a reason for the sexual dimorphism of TTX distribution in gonads. In addition, TdT-mediated dUTP-biotin nick end labeling staining showed that significant apoptosis was detected in the MTG testis, and the role of the cell apoptotic pathways was further confirmed. Overall, our research revealed that the response of the ovary and testis to TTX administration was largely different, the ovary is more tolerant whereas the testis is more sensitive to TTX. These data will deepen our understanding on the accumulation of TTX sexual dimorphism in Takifugu.

BACKGROUND: The expression and biological functions of circular RNAs (circRNAs) in reproductive organs have been extensively reported. However, it is still unclear whether circRNAs are involved in sex change. To this end, RNA sequencing (RNA-seq) was performed in gonads at 5 sexual stages (ovary, early intersexual stage gonad, middle intersexual stage gonad, late intersexual stage gonad, and testis) of ricefield eel, and the expression profiles and potential functions of circRNAs were studied. RESULTS: Seven hundred twenty-one circRNAs were identified, and the expression levels of 10 circRNAs were verified by quantitative real-time PCR (qRT-PCR) and found to be in accordance with the RNA-seq data, suggesting that the RNA-seq data were reliable. Then, the sequence length, category, sequence composition and the relationship between the parent genes of the circRNAs were explored. A total of 147 circRNAs were differentially expressed in the sex change process, and GO and KEGG analyses revealed that some differentially expressed (such as novel_circ_0000659, novel_circ_0004005 and novel_circ_0005865) circRNAs were closely involved in sex change. Furthermore, expression pattern analysis demonstrated that both circSnd1 and foxl2 were downregulated in the process of sex change, which was contrary to mal-miR-135b. Finally, dual-luciferase reporter assay and RNA immunoprecipitation showed that circSnd1 and foxl2 can combine with mal-miR-135b and mal-miR-135c. These data revealed that circSnd1 regulates foxl2 expression in the sex change of ricefield eel by acting as a sponge of mal-miR-135b/c. CONCLUSION: Our results are the first to demonstrate that circRNAs have potential effects on sex change in ricefield eel; and circSnd1 could regulate foxl2 expression in the sex change of ricefield eel by acting as a sponge of mal-miR-135b/c. These data will be useful for enhancing our understanding of sequential hermaphroditism and sex change in ricefield eel or other teleosts.

BACKGROUND: Lysine 2-hydroxyisobutyrylation (Khib) is a novel and conserved post-translational modification (PTM). Frankliniella occidentalis are economically important agricultural pests globally and also notorious for vectoring destructive plant viruses. To better study the disease transmission mechanism of F. occidentalis, it is necessary to conduct in-depth analysis of it. So far, no Khib modification of insects has been reported. RESULTS: In this study, a proteome-wide analysis of Khib modifications in F. occidentalis was analyzed for the first time through the combination of high performance liquid chromatography fractionation technology and 2-hydroxyisobutyrylated peptide enrichment and other advanced technologies, 4093 Khib sites were identified on 1125 modified proteins. Bioinformatics and functional enrichment analyses showed that Khib-modified proteins were significantly enriched in many cell compartments and pathways, especially related to various cellular components and biological processes, and were more concentrated in ribosomes and proteasome subunits, involved in energy metabolism, protein synthesis and degradation, compared to the other nine species including Japonica rice, Homo sapiens, P. patens, Botrytis, Ustilaginoidea virens, Saccharomyces cerevisiae, T. gondii, C. albicans, and F. oxysporum. And Khib sites on virus-interacting insect proteins were discovered for the first time, such as cyclophilin and endoCP-GN. CONCLUSIONS: After three repeated experiments, we found a total of 4093 Khib sites on 1125 proteins. These modified proteins are mainly concentrated in ribosomes and proteasome subunits, and are widely involved in a variety of critical biological activities and metabolic processes of F. occidentalis. In addition, for the first time, Khib modification sites are found on the proteome of F. occidentalis, and these sites could be acted as for the virus interaction, including cyclophilin and endoCP-GN. The global map of 2-hydroxyisobutyrylation in thrips is an invaluable resource to better understand the biological processes of thrips and provide new means for disease control and mitigation of pest damage to crops.

BACKGROUND: The objective of this study was to reveal the flavonoid biosynthesis pathway in white (Z6), red (Z27) and black (HC4) seeds of the sweet sorghum (Sorghum bicolor) using metabolomics and transcriptomics, to identify different flavonoid metabolites, and to analyze the differentially expressed genes involved in flavonoid biosynthesis. RESULTS: We analyzed the metabolomics and transcriptomics data of sweet sorghum seeds. Six hundred and fifty-one metabolites including 171 flavonoids were identified in three samples. Integrated analysis of transcriptomics and metabolomics showed that 8 chalcone synthase genes (gene19114, gene19115, gene19116, gene19117, gene19118, gene19120, gene19122 and gene19123) involved in flavonoid biosynthesis, were identified and play central role in change of color. Six flavanone including homoeriodictyol, naringin, prunin, naringenin, hesperetin and pinocembrin were main reason for the color difference. CONCLUSIONS: Our results provide valuable information on the flavonoid metabolites and the candidate genes involved in the flavonoid biosynthesis pathway in sweet sorghum seeds.

BACKGROUND: Tetrastigma hemsleyanum is a valuable traditional Chinese medicinal plant widely distributed in the subtropical areas of China. It belongs to the Cayratieae tribe, family Vitaceae, and exhibited significant anti-tumor and anti-inflammatory activities. However, obvious differences were observed on the quality of T. hemsleyanum root from different regions, requiring the discrimination strategy for the geographical origins. RESULT: This study characterized five complete chloroplast (cp) genomes of T. hemsleynum samples from different regions, and conducted a comparative analysis with other representing species from family Vitaceae to reveal the structural variations, informative markers and phylogenetic relationships. The sequenced cp genomes of T. hemsleyanum exhibited a conserved quadripartite structure with full length ranging from 160,124 bp of Jiangxi Province to 160,618 bp of Zhejiang Province. We identified 112 unique genes (80 protein-coding, 28 tRNA and 4 rRNA genes) in the cp genomes of T. hemsleyanum with highly similar gene order, content and structure. The IR contraction/expansion events occurred on the junctions of ycf1, rps19 and rpl2 genes with different degrees, causing the differences of genome sizes in T. hemsleyanum and Vitaceae plants. The number of SSR markers discovered in T. hemsleyanum was 56-57, exhibiting multiple differences among the five geographic groups. Phylogenetic analysis based on conserved cp genome proteins strongly grouped the five T. hemsleyanum species into one clade, showing a sister relationship with T. planicaule. Comparative analysis of the cp genomes from T. hemsleyanum and Vitaceae revealed five highly variable spacers, including 4 intergenic regions and one protein-coding gene (ycf1). Furthermore, five mutational hotspots were observed among T. hemsleyanum cp genomes from different regions, providing data for designing DNA barcodes trnL and trnN. The combination of molecular markers of trnL and trnN clustered the T. hemsleyanum samples from different regions into four groups, thus successfully separating specimens of Sichuan and Zhejiang from other areas. CONCLUSION: Our study obtained the chloroplast genomes of T. hemsleyanum from different regions, and provided a potential molecular tracing tool for determining the geographical origins of T. hemsleyanum, as well as important insights into the molecular identification approach and and phylogeny in Tetrastigma genus and Vitaceae family.

BACKGROUND: The reduction of the chromosome number from 48 in the Great Apes to 46 in modern humans is thought to result from the end-to-end fusion of two ancestral non-human primate chromosomes forming the human chromosome 2 (HSA2). Genomic signatures of this event are the presence of inverted telomeric repeats at the HSA2 fusion site and a block of degenerate satellite sequences that mark the remnants of the ancestral centromere. It has been estimated that this fusion arose up to 4.5 million years ago (Mya). RESULTS: We have developed an enhanced algorithm for the detection and efficient counting of the locally over-represented weak-to-strong (AT to GC) substitutions. By analyzing the enrichment of these substitutions around the fusion site of HSA2 we estimated its formation time at 0.9 Mya with a 95% confidence interval of 0.4-1.5 Mya. Additionally, based on the statistics derived from our algorithm, we have reconstructed the evolutionary distances among the Great Apes (Hominoidea). CONCLUSIONS: Our results shed light on the HSA2 fusion formation and provide a novel computational alternative for the estimation of the speciation chronology.

BACKGROUND: The histidine metabolism and transport (his) genes are controlled by a variety of RNA-dependent regulatory systems among diverse taxonomic groups of bacteria including T-box riboswitches in Firmicutes and Actinobacteria and RNA attenuators in Proteobacteria. Using a comparative genomic approach, we previously identified a novel DNA-binding transcription factor (named HisR) that controls the histidine metabolism genes in diverse Gram-positive bacteria from the Firmicutes phylum. RESULTS: Here we report the identification of HisR-binding sites within the regulatory regions of the histidine metabolism and transport genes in 395 genomes representing the Bacilli, Clostridia, Negativicutes, and Tissierellia classes of Firmicutes, as well as in 97 other HisR-encoding genomes from the Actinobacteria, Proteobacteria, and Synergistetes phyla. HisR belongs to the TrpR family of transcription factors, and their predicted DNA binding motifs have a similar 20-bp palindromic structure but distinct lineage-specific consensus sequences. The predicted HisR-binding motif was validated in vitro using DNA binding assays with purified protein from the human gut bacterium Ruminococcus gnavus. To fill a knowledge gap in the regulation of histidine metabolism genes in Firmicutes genomes that lack a hisR repressor gene, we systematically searched their upstream regions for potential RNA regulatory elements. As result, we identified 158 T-box riboswitches preceding the histidine biosynthesis and/or transport genes in 129 Firmicutes genomes. Finally, novel candidate RNA attenuators were identified upstream of the histidine biosynthesis operons in six species from the Bacillus cereus group, as well as in five Eubacteriales and six Erysipelotrichales species. CONCLUSIONS: The obtained distribution of the HisR transcription factor and two RNA-mediated regulatory mechanisms for histidine metabolism genes across over 600 species of Firmicutes is discussed from functional and evolutionary points of view.

BACKGROUND: Vagococcus fluvialis is a species of lactic acid bacteria found both free-living in river and seawater and associated to hosts, such as marine sponges. This species has been greatly understudied, with no complete genome assembly available to date, which is essential for the characterisation of the mobilome. RESULTS: We sequenced and assembled de novo the complete genome sequences of five V. fluvialis isolates recovered from marine sponges. Pangenome analysis of the V. fluvialis species (total of 17 genomes) showed a high intraspecific diversity, with 45.5% of orthologous genes found to be strain specific. Despite this diversity, analyses of gene functions clustered all V. fluvialis species together and separated them from other sequenced Vagococcus species. V. fluvialis strains from different habitats were highly similar in terms of functional diversity but the sponge-isolated strains were enriched in several functions related to the marine environment. Furthermore, sponge-isolated strains carried a significantly higher number of mobile genetic elements (MGEs) compared to previously sequenced V. fluvialis strains from other environments. Sponge-isolated strains carried up to 4 circular plasmids each, including a 48-kb conjugative plasmid. Three of the five strains carried an additional circular extrachromosomal sequence, assumed to be an excised prophage as it contained mainly viral genes and lacked plasmid replication genes. Insertion sequences (ISs) were up to five times more abundant in the genomes of sponge-isolated strains compared to the others, including several IS families found exclusively in these genomes. CONCLUSIONS: Our findings highlight the dynamics and plasticity of the V. fluvialis genome. The abundance of mobile genetic elements in the genomes of sponge-isolated V. fluvialis strains suggests that the mobilome might be key to understanding the genomic signatures of symbiosis in bacteria.

BACKGROUND: The use of archived formalin-fixed paraffin-embedded (FFPE) tumor tissues has become a common practice in clinical and epidemiologic genetic research. Simultaneous extraction of DNA and RNA from FFPE tissues is appealing but can be practically challenging. Here we report our results and lessons learned from processing FFPE breast tumor tissues for a large epidemiologic study. METHODS: Qiagen AllPrep DNA/RNA FFPE kit was adapted for dual extraction using tissue punches or sections from breast tumor tissues. The yield was quantified using Qubit and fragmentation analysis by Agilent Bioanalyzer. A subset of the DNA samples were used for genome-wide DNA methylation assays and RNA samples for sequencing. The QC metrices and performance of the assays were analyzed with pre-analytical variables. RESULTS: A total of 1859 FFPE breast tumor tissues were processed. We found it critical to adjust proteinase K digestion time based on tissue volume to achieve balanced yields of DNA and RNA. Tissue punches taken from tumor-enriched regions provided the most reliable output. A median of 1475 ng DNA and 1786 ng RNA per sample was generated. The median DNA integrity number (DIN) was 3.8 and median DV200 for RNA was 33.2. Of 1294 DNA samples used in DNA methylation assays, 97% passed quality check by qPCR and 92% generated data deemed high quality. Of the 130 RNA samples with DV200 ≥ 20% used in RNA-sequencing, all but 5 generated usable transcriptomic data with a mapping rate ≥ 60%. CONCLUSIONS: Dual DNA/RNA purification using Qiagen AllPrep FFPE extraction protocol is feasible for clinical and epidemiologic studies. We recommend tissue punches as a reliable source material and fine tuning of proteinase K digestion time based on tissue volume. IMPACT: Our protocol and recommendations may be adapted by future studies for successful extraction of archived tumor tissues.

Individuals with schizophrenia (SCZ) have, on average, a 10- to 20-year shorter expected life span than the rest of the population, primarily due to cardiovascular disease comorbidity. Genome-wide association studies (GWAS) have previously been used to separately identify common variants in SCZ and cardiometabolic traits. However, genetic variants jointly influencing both traits remain to be fully characterised. To assess overlaps (if any) between the genetic architecture of SCZ and cardiometabolic traits, we used conditional false discovery rate (FDR) and local genetic correlation statistical framework analyses. A conjunctional FDR was used to identify shared genetic traits between SCZ and cardiometabolic risk factors. We identified 144 genetic variants which were shared between SCZ and body mass index (BMI), and 15 variants shared between SCZ and triglycerides (TG). Furthermore, we discovered four novel single nucleotide polymorphisms (SNPs) (rs3865350, rs9860913, rs13307 and rs9614186) and four proximate genes (DERL2, SNX4, LY75 and EFCAB6) which were shared by SCZ and BMI. We observed that the novel genetic variant rs13307 and the most proximate gene LY75 exerted potential effects on SCZ and BMI comorbidity. Also, we observed a mixture of concordant and opposite direction associations with shared genetic variants. We demonstrated a moderate to high genetic overlap between SCZ and cardiometabolic traits associated with a pattern of bidirectional associations. Our data suggested a complex interplay between metabolism-related gene pathways in SCZ pathophysiology.

DNA and RNA sequencing are widely used techniques to investigate genomic modifications and gene expression. The costs for sequencing dropped dramatically in the last decade. However, due to material and labor intense steps, the sample preparation costs could not keep up with that pace. About 80% of the total costs occur prior to sequencing during DNA/RNA extraction, enrichment steps and subsequent library preparation. In this study, we investigate the potential of pooling different organisms samples prior to DNA/RNA extraction to significantly reduce costs in preparative steps. Similar to the common procedure of ligated DNA tags to pool (c)DNA samples, sequence diversity of different organisms intrinsically provide unique sequences that allow separation of reads after sequencing. With this approach, sample pooling can occur before DNA/RNA isolation and library preparation. We show that pooled sequencing of three related bacterial organisms is possible without loss of data quality at a cost reduction of approx. 50% in DNA- and RNA-seq approaches. Furthermore, we show that this approach is highly efficient down to the level of a shared genus and is, therefore, widely applicable in sequencing facilities and companies with diverse sample pools.

BACKGROUND: Emerging transcriptome-wide high-throughput screenings reveal the landscape and functions of RNAs, such as circular RNAs (circRNAs), in human cancer. In addition, the post-transcriptional RNA internal modifications, especially N6-methyladenosine (m6A), greatly enrich the variety of RNAs metabolism. However, the m6A modification on circRNAs has yet to be addressed. RESULTS: Here, we report an epitranscriptome-wide mapping of m6A-modified circRNAs (m6A-circRNA) in oral squamous cell carcinoma (OSCC). Utilizing the data of m6A methylated RNA immunoprecipitation sequencing (MeRIP-seq) and m6A-circRNAs microarray, we found that m6A-circRNAs exhibited particular modification styles in OSCC, which was independent of m6A-mRNA. Besides, m6A modification on circRNAs frequently occurred on the long exons in the front part of the coding sequence (CDS), which was distinct from m6A-mRNA that in 3'-UTR or stop codon. CONCLUSION: In conclusion, our work preliminarily demonstrates the traits of m6A-circRNAs, which may bring enlighten for the roles of m6A-circRNAs in OSCC.

BACKGROUND: Nematodes are parasitic animals that cause over 100 billion US dollars loss in agricultural business. The whole-genomes of two Streptomyces strains, Streptomyces spectabilis KCTC9218T and Streptomyces sp. AN091965, were sequenced. Both strains produce spectinabilin, an antinematode drug. Its secondary metabolism was examined to aid the development of an efficient nematicidal drug-producing host strain. RESULTS: The whole-genome sequences of S. spectabilis KCTC9218T and Streptomyces sp. AN091965 were analyzed using PacBio and Illumina sequencing platforms, and assembled using hybrid methodology. The total contig lengths for KCTC9218T and AN091965 were 9.97 Mb and 9.84 Mb, respectively. A total of 8,374 and 8,054 protein-coding genes, as well as 39 and 45 secondary metabolite biosynthetic gene clusters were identified in KCTC9218T and AN091965, respectively. 18.4 ± 6.45 mg/L and 213.89 ± 21.30 mg/L of spectinabilin were produced by S. spectabilis KCTC9218T and Streptomyces sp. AN091965, respectively. Pine wilt disease caused by nematode was successfully prevented by lower concentration of spectinabilin injection than that of abamectin recommended by its manufacturer. Production of multiple antinematode drugs, including spectinabilin, streptorubin B, and undecylprodigiosin was observed in both strains using high-resolution liquid chromatography mass spectrometry (LC-MS) analysis. CONCLUSIONS: Whole-genome sequencing of spectinabilin-producing strains, coupled with bioinformatics and mass spectrometry analyses, revealed the production of multiple nematicidal drugs in the KCTC9218T and AN091965 strains. Especially, Streptomyces sp. AN091965 showed high production level of spectinabilin, and this study provides crucial information for the development of potential nematicidal drug producers.

BACKGROUND: Dendrobium catenatum/D. officinale (here after D. catenatum), a well-known economically important traditional medicinal herb, produces a variety of bioactive metabolites including polysaccharides, alkaloids, and flavonoids with excellent pharmacological and clinical values. Although many genes associated with the biosynthesis of medicinal components have been cloned and characterized, the biosynthetic pathway, especially the downstream and regulatory pathway of major medicinal components in the herb, is far from clear. β-glucosidases (BGLUs) comprise a diverse group of enzymes that widely exist in plants and play essential functions in cell wall modification, defense response, phytohormone signaling, secondary metabolism, herbivore resistance, and scent release by hydrolyzing β-D-glycosidic bond from a carbohydrate moiety. The recent release of the chromosome-level reference genome of D. catenatum enables the characterization of gene families. Although the genome-wide analysis of the BGLU gene family has been successfully conducted in various plants, no systematic analysis is available for the D. catenatum. We previously isolated DcBGLU2 in the BGLU family as a key regulator for polysaccharide biosynthesis in D. catenatum. Yet, the exact number of DcBGLUs in the D. catenatum genome and their possible roles in bioactive compound production deserve more attention. RESULTS: To investigate the role of BGLUs in active metabolites production, 22 BGLUs (DcBGLU1-22) of the glycoside hydrolase family 1 (GH1) were identified from D. catenatum genome. Protein prediction showed that most of the DcBGLUs were acidic and phylogenetic analysis classified the family into four distinct clusters. The sequence alignments revealed several conserved motifs among the DcBGLU proteins and analyses of the putative signal peptides and N-glycosylation site revealed that the majority of DcBGLU members dually targeted to the vacuole and/or chloroplast. Organ-specific expression profiles and specific responses to MeJA and MF23 were also determined. Furthermore, four DcBGLUs were selected to test their involvement in metabolism regulation. Overexpression of DcBGLU2, 6, 8, and 13 significantly increased contents of flavonoid, reducing-polysaccharide, alkaloid and soluble-polysaccharide, respectively. CONCLUSION: The genome-wide systematic analysis identified candidate DcBGLU genes with possible roles in medicinal metabolites production and laid a theoretical foundation for further functional characterization and molecular breeding of D. catenatum.

BACKGROUND: For over a century, scientists have studied host-pathogen interactions between the crayfish plague disease agent Aphanomyces astaci and freshwater crayfish. It has been hypothesised that North American crayfish hosts are disease-resistant due to the long-lasting coevolution with the pathogen. Similarly, the increasing number of latent infections reported in the historically sensitive European crayfish hosts seems to indicate that similar coevolutionary processes are occurring between European crayfish and A. astaci. Our current understanding of these host-pathogen interactions is largely focused on the innate immunity processes in the crayfish haemolymph and cuticle, but the molecular basis of the observed disease-resistance and susceptibility remain unclear. To understand how coevolution is shaping the host's molecular response to the pathogen, susceptible native European noble crayfish and invasive disease-resistant marbled crayfish were challenged with two A. astaci strains of different origin: a haplogroup A strain (introduced to Europe at least 50 years ago, low virulence) and a haplogroup B strain (signal crayfish in lake Tahoe, USA, high virulence). Here, we compare the gene expression profiles of the hepatopancreas, an integrated organ of crayfish immunity and metabolism. RESULTS: We characterised several novel innate immune-related gene groups in both crayfish species. Across all challenge groups, we detected 412 differentially expressed genes (DEGs) in the noble crayfish, and 257 DEGs in the marbled crayfish. In the noble crayfish, a clear immune response was detected to the haplogroup B strain, but not to the haplogroup A strain. In contrast, in the marbled crayfish we detected an immune response to the haplogroup A strain, but not to the haplogroup B strain. CONCLUSIONS: We highlight the hepatopancreas as an important hub for the synthesis of immune molecules in the response to A. astaci. A clear distinction between the innate immune response in the marbled crayfish and the noble crayfish is the capability of the marbled crayfish to mobilise a higher variety of innate immune response effectors. With this study we outline that the type and strength of the host immune response to the pathogen is strongly influenced by the coevolutionary history of the crayfish with specific A. astaci strains.

BACKGROUND: Bacillus Calmette-Guérin (BCG) refers to a group of vaccine strains with unique genetic characteristics. BCG is the only available vaccine for preventing tuberculosis (TB). Genetic and biochemical variations among the BCG vaccine strains have been considered as one of the significant parameters affecting the variable protective efficacy of the vaccine against pulmonary tuberculosis. To track genetic variations, here two vaccine strains (Danish 1331 and Pasteur 1173P2) popularly used according to the BCG World Atlas were subjected to a comparative analysis against the Mycobacterium tuberculosis H37Rv, Mycobacterium bovis AF2122/97, and Mycobacterium tuberculosis variant bovis BCG str. Pasteur 1173P2 reference genomes. Besides, the presence or absence of the experimentally verified human T cell epitopes was examined. RESULTS: Only two variants were identified in BCG Danish 1331 that have not been reported previously in any BCG strains with the complete submitted genome yet. Furthermore, we identified a DU1-like 14,577 bp region in BCG Danish 1331; The duplication which was previously seemed to be exclusive to the BCG Pasteur. We also found that 35% of the T cell epitopes are absent from both strains, and epitope sequences are more conserved than the rest of the genome. CONCLUSIONS: We provided a comprehensive catalog of single nucleotide polymorphisms (SNPs) and short insertions and deletions (indels) in BCG Danish 1331 and BCG Pasteur 1173P2. These findings may help determine the effect of genetic variations on the variable protective efficacy of BCG vaccine strains.

BACKGROUND: Functional genomics uses unbiased systematic genome-wide gene disruption or analyzes natural variations such as gene expression profiles of different tissues from multicellular organisms to link gene functions to particular phenotypes. Functional genomics approaches are of particular importance to identify large sets of genes that are specifically important for a particular biological process beyond known candidate genes, or when the process has not been studied with genetic methods before. RESULTS: Here, we present a large set of genes whose disruption interferes with the function of the odoriferous defensive stink glands of the red flour beetle Tribolium castaneum. This gene set is the result of a large-scale systematic phenotypic screen using RNA interference applied in a genome-wide forward genetics manner. In this first-pass screen, 130 genes were identified, of which 69 genes could be confirmed to cause phenotypic changes in the glands upon knock-down, which vary from necrotic tissue and irregular reservoir size to irregular color or separation of the secreted gland compounds. Gene ontology analysis revealed that many of those genes are encoding enzymes (peptidases and cytochromes P450) as well as proteins involved in membrane trafficking with an enrichment in lysosome and mineral absorption pathways. The knock-down of 13 genes caused specifically a strong reduction of para-benzoquinones in the gland reservoirs, suggesting a specific function in the synthesis of these toxic compounds. Only 14 of the 69 confirmed gland genes are differentially overexpressed in stink gland tissue and thus could have been detected in a transcriptome-based analysis. However, only one out of eight genes identified by a transcriptomics approach known to cause phenotypic changes of the glands upon knock-down was recognized by this phenotypic screen, indicating the limitation of such a non-redundant first-pass screen. CONCLUSION: Our results indicate the importance of combining diverse and independent methodologies to identify genes necessary for the function of a certain biological tissue, as the different approaches do not deliver redundant results but rather complement each other. The presented phenotypic screen together with a transcriptomics approach are now providing a set of close to hundred genes important for odoriferous defensive stink gland physiology in beetles.

BACKGROUND: Transcription factors (TFs) play important roles in plants. Among the major TFs, GATA plays a crucial role in plant development, growth, and stress responses. However, there have been few studies on the GATA gene family in foxtail millet (Setaria italica). The release of the foxtail millet reference genome presents an opportunity for the genome-wide characterization of these GATA genes. RESULTS: In this study, we identified 28 GATA genes in foxtail millet distributed on seven chromosomes. According to the classification method of GATA members in Arabidopsis, SiGATA was divided into four subfamilies, namely subfamilies I, II, III, and IV. Structural analysis of the SiGATA genes showed that subfamily III had more introns than other subfamilies, and a large number of cis-acting elements were abundant in the promoter region of the SiGATA genes. Three tandem duplications and five segmental duplications were found among SiGATA genes. Tissue-specific results showed that the SiGATA genes were mainly expressed in foxtail millet leaves, followed by peels and seeds. Many genes were significantly induced under the eight abiotic stresses, such as SiGATA10, SiGATA16, SiGATA18, and SiGATA25, which deserve further attention. CONCLUSIONS: Collectively, these findings will be helpful for further in-depth studies of the biological function of SiGATA, and will provide a reference for the future molecular breeding of foxtail millet.

BACKGROUND: The Qinba region is the transition region between Indica and Japonica varieties in China. It has a long history of Indica rice planting of more than 7000 years and is also a planting area for fine-quality Indica rice. The aims of this study are to explore different genetic markers applied to the analysis population structure, genetic diversity, selection and optimization of molecular markers of Indica rice, thus providing more information for the protection and utilization on germplasm resources of Indica rice. METHODS: Fifteen phenotypic traits, a core set of 48 SSR markers which originated protocol for identification of rice varieties-SSR marker method in agricultural industry standard of the People's Republic of China (Ministry of Agriculture of the PRC, NY/T1433-2014, Protocol for identification of rice varieties-SSR marker method, 2014), and SNPs data obtained by genotyping-by-sequencing (GBS, NlaIII and MseI digestion, referred to as SNPs-NlaIII and SNPs-MseI, respectively) for this panel of 93 samples using the Illumina HiSeq2000 sequencing platform, were employed to explore the genetic diversity and population structure of 93 samples. RESULTS: The average of coefficient of variation (CV) and diversity index (He) were 29.72% and 1.83 ranging from 3.07% to 137.43%, and from 1.45 to 2.03, respectively. The correlation coefficient between 15 phenotypic traits ranged from 0.984 to -0.604. The first four PCs accounted for 70.693% phenotypic variation based on phenotypic analysis. A total of 379 alleles were obtained using SSR markers, encompassing an average of 8.0 alleles per primer. Polymorphic bands (PPB) and polymorphism information content (PIC) was 88.65% and 0.77, respectively. The Mantel test showed that the correlation between the genetic distance matrix based on SNPs-NlaIII and SNPs-MseI was the largest (R2=0.88), and that based on 15 phenotypic traits and SSR was the smallest (R2=0.09). The 93 samples could be clustered into two subgroups by 3 types of genetic markers. Molecular variance analysis revealed that the genetic variation was 2% among populations and 98% within populations (the Nm was 0.16), Tajima's D value was 1.66, the FST between the two populations was 0.61 based on 72,824 SNPs. CONCLUSIONS: The population genetic variation explained by SNPs was larger than that explained by SSRs. The gene flow of 93 samples used in this study was larger than that of naturally self-pollinated crops, which may be caused by long-term breeding selection of Indica rice in the Qinba region. The genetic structure of the 93 samples was simple and lacked rare alleles.

BACKGROUND: Mycoplasma (M.) bovis is a major etiological agent of bovine respiratory disease, which is the most economically costly disease of cattle worldwide. Cattle disease surveillance on M. bovis is increasingly using gene-based techniques, such as multilocus sequence typing (MLST), or genome-based techniques such as core genome MLST that both require only partial genomic data. However, accurate up-to-date surveillance also demands complete, circular genomes that can be used as reference to track the evolution of the different lineages. Yet, in France, two of the main subtypes currently circulating still have no representing genome in public databases. Here, to address this gap, we provide and compare three new complete M. bovis genomes obtained from recent clinical isolates that represent major subtypes circulating in France and Europe. RESULTS: Genomes were obtained using a hybrid assembly strategy (Illumina and Nanopore) with fine-tuning of settings and inputs used in the Unicycler assembly pipeline, such as size selection of reads and quality trimming of the FASTQ files. The main characteristics and synteny of the genomes were compared. The three genomes mainly differed by their content in terms of mobile genetic elements, i.e. integrative conjugative elements (ICE) and insertion sequences (IS), a feature that impacts their structure. For instance, strain L15527, representing subtype3 (st3), harbours an exceptionally high number of ICEs, which results in a bigger-sized genome than all those previously described and could be associated with the propensity of st3 to gain and fix mutations through chromosomal transfer mechanisms. In contrast, strain F9160, of st1, is very close to the PG45 type strain isolated in 1961 in the USA, and harbours a huge number of IS. These features may be associated with an evolution towards a host-restricted state or in a "closed" host or environment reservoir until a recent re-emergence. CONCLUSIONS: Whole-genome comparison of the three French M. bovis subtypes provides valuable resources for future studies combining epidemiology, phylogenetic data, and phylodynamic methods.

BACKGROUND: Rice sheath blight, which is caused by Rhizoctonia solani, is the most destructive disease affecting rice production, but the resistance mechanism to this pathogen has not been fully elucidated. RESULTS: In this study, we selected two rice cultivars based on their resistance to the pathogen and analyzed and compared the transcriptomic profiles of two cultivars, the moderately resistant variety Gangyuan8 and the highly susceptible variety Yanfeng47, at different time points after inoculation. The comparative transcriptome profiling showed that the expression of related genes gradually increased after pathogen inoculation. The number of differentially expressed genes (DEGs) in Yanfeng47 was higher than that in Gangyuan8, and this result revealed that Yanfeng47 was more susceptible to fungal attack. At the early stage (24 and 48 h), the accumulation of resistance genes and a resistance metabolism occurred earlier in Ganguan8 than in Yanfeng47, and the resistance enrichment entries were more abundant in Ganguan8 than in Yanfeng47. CONCLUSIONS: Based on the GO and KEGG enrichment analyses at five infection stages, we concluded that phenylalanine metabolism and the jasmonic acid pathway play a crucial role in the resistance of rice to sheath blight. Through a comparative transcriptome analysis, we preliminarily analyzed the molecular mechanism responsible for resistance to sheath blight in rice, and the results lay the foundation for the development of gene mining and functional research on rice resistance to sheath blight.

BACKGROUND: Morinda officinalis How (MO) is a perennial herb distributed in tropical and subtropical regions, which known as one of the "Four Southern Herbal Medicines". The extent of genetic variability and the population structure of MO are presently little understood. Here, nine morphological traits, six chemical components and Single nucleotide polymorphism (SNP) markers were used in integrative research of MO germplasm variation among 88 individuals collected from ten populations across four geographical provinces of China. RESULTS: Both phenotype and chemical composition have significant genetic variation, and there is a certain correlation between them such as root diameter and the nystose content, as well as geographical distribution. The principal component analysis (PCA) showed the leaf length, leaf width, nystose, 1F-furanosaccharide nystose, and the section color were the major contributors to diversity. The cluster analysis based on phenotypic and oligosaccharide data distinguished three significant groups, which was consistent with the result of a corresponding analysis with 228,615 SNP markers, and importantly, they all showed a significant correlation with geographical origin. However, there was little similarity between two cluster results. The Shannon's information index (I) varied from 0.17 to 0.53 with a mean of 0.37, suggesting a high level of genetic diversity in MO populations, which mainly existed among individuals within populations, accounting for 99.66% of the total according to the analysis of molecular variance (AMOVA) results. Each population also maintains the connection because of certain gene communication, so that the genetic differentiation between populations was not very significant. The STRUCTURE software was used to analyse the population structure and the result showed that 88 accessions were clustered into three groups, and 67% of them were pure type, which was also confirmed through PCA. CONCLUSIONS: The comprehensive study of phenotypic, chemical and molecular markers will provide valuable information for future breeding plans and understanding the phylogenetic relationship of MO population.

It is currently assumed that around 100 million years ago, the common ancestor to the Fabales, Fagales, Rosales and Cucurbitales in Gondwana, developed a root nodule symbiosis with a nitrogen-fixing bacterium. The symbiotic trait evolved first in Frankia cluster-2; thus, strains belonging to this cluster are the best extant representatives of this original symbiont. Most cluster-2 strains could not be cultured to date, except for Frankia coriariae, and therefore many aspects of the symbiosis are still elusive. Based on phylogenetics of cluster-2 metagenome-assembled genomes (MAGs), it has been shown that the genomes of strains originating in Eurasia are highly conserved. These MAGs are more closely related to Frankia cluster-2 in North America than to the single genome available thus far from the southern hemisphere, i.e., from Papua New Guinea.To unravel more biodiversity within Frankia cluster-2 and predict routes of dispersal from Gondwana, we sequenced and analysed the MAGs of Frankia cluster-2 from Coriaria japonica and Coriaria intermedia growing in Japan, Taiwan and the Philippines. Phylogenetic analyses indicate there is a clear split within Frankia cluster-2, separating a continental from an island lineage. Presumably, these lineages already diverged in Gondwana.Based on fossil data on the host plants, we propose that these two lineages dispersed via at least two routes. While the continental lineage reached Eurasia together with their host plants via the Indian subcontinent, the island lineage spread towards Japan with an unknown host plant.

BACKGROUND: The development of the ovaries is an important factor that affects egg production performance in geese. Ovarian development is regulated by genes that are expressed dynamically and stage-specifically. The transcriptome profile analysis on ovarian tissues of goose at different egg laying stages could provide an important basis for screening and identifying key genes regulating ovarian development. RESULTS: In this study, 4 ovary tissues at each breeding period of pre-laying (PP), laying (LP), and ceased-laying period (CP), respectively, with significant morphology difference, were used for RNA extraction and mRNAs, lncRNAs, and miRNAs comparison in Yili geese. CeRNA regulatory network was constructed for key genes screening. A total of 337, 1136, and 525 differentially expressed DE mRNAs, 466, 925, and 742 DE lncRNAs and 258, 1131 and 909 DE miRNAs were identified between PP and LP, between CP and LP, and between CP and PP groups, respectively. Functional enrichment analysis showed that the differentially expressed mRNAs and non-coding RNA target genes were mainly involved in the cell process, cytokine-cytokine receptor interaction, phagosome, calcium signaling pathway, steroid biosynthesis and ECM-receptor interaction. Differential genes and non-coding RNAs, PDGFRB, ERBB4, LHCGR, MSTRG.129094.34, MSTRG.3524.1 and gga-miR-145-5p, related to reproduction and ovarian development were highly enriched. Furthermore, lncRNA-miRNA-mRNA regulatory networks related to ovary development were constructed. CONCLUSIONS: Our study found dramatic transcriptomic differences in ovaries of Yili geese at different egg-laying stages, and a differential lncRNA-miRNA-mRNA regulatory network related to cell proliferation, differentiation and apoptosis and involved in stromal follicle development were established and preliminarily validated, which could be regarded as a key regulatory pathway of ovarian development in Yili geese.

BACKGROUND: Flavonols are the largest subgroup of flavonoids, possessing multiple functions in plants including protection against ultraviolet radiation, antimicrobial activities, and flower pigmentation together with anthocyanins. They are of agronomical and economical importance because the major off-taste component in rapeseed protein isolates is a flavonol derivative, which limits rapeseed protein use for human consumption. Flavonol production in Arabidopsis thaliana is mainly regulated by the subgroup 7 (SG7) R2R3-MYB transcription factors MYB11, MYB12, and MYB111. Recently, the SG19 MYBs MYB21, MYB24, and MYB57 were shown to regulate flavonol accumulation in pollen and stamens. The members of each subgroup are closely related, showing gene redundancy and tissue-specific expression in A. thaliana. However, the evolution of these flavonol regulators inside the Brassicaceae, especially inside the Brassiceae, which include the rapeseed crop species, is not fully understood. RESULTS: We studied the SG7 and SG19 MYBs in 44 species, including 31 species of the Brassicaceae, by phylogenetic analyses followed by synteny and gene expression analyses. Thereby we identified a deep MYB12 and MYB111 duplication inside the Brassicaceae, which likely occurred before the divergence of Brassiceae and Thelypodieae. These duplications of SG7 members were followed by the loss of MYB11 after the divergence of Eruca vesicaria from the remaining Brassiceae species. Similarly, MYB21 experienced duplication before the emergence of the Brassiceae tribe, where the gene loss of MYB24 is also proposed to have happened. The members of each subgroup revealed frequent overlapping spatio-temporal expression patterns in the Brassiceae member B. napus, which are assumed to compensate for the loss of MYB11 and MYB24 in the analysed tissues. CONCLUSIONS: We identified a duplication of MYB12, MYB111, and MYB21 inside the Brassicaceae and MYB11 and MYB24 gene loss inside the tribe Brassiceae. We propose that polyploidization events have shaped the evolution of the flavonol regulators in the Brassicaceae, especially in the Brassiceae.

BACKGROUND: Due to the prolonged usage of Bt-based biopesticides and Bt-transgenic crops worldwide, insects are continually developing resistance against Cry toxins. This resistance may occur if any mechanistic step in the insecticidal process is disrupted possibly because of the alteration in Cry-receptor binding affinity due to mutation in receptor genes. Compared to other lepidopteran insects, Cry receptor-related research has made asymmetric progress in the model insect Galleria mellonella. RESULTS: Present study describes the molecular characterization and functional analysis of five Cry toxin receptor-related genes (prohibitin, GLTP, α-amylase, ADAM and UDP-GT) and a gut repair gene (arylphorin) from the gut tissues of G. mellonella. Protein-protein docking analysis revealed that Cry1AcF putatively binds with all the five candidate proteins, suggesting their receptor-like function. These receptor-like genes were significantly overexpressed in the gut tissues of fourth-instar G. mellonella larvae upon early exposure to a sub-lethal dose of Cry1AcF toxin. However, targeted knockdown (by using bacterially-expressed dsRNAs) of these genes led to variable effect on insect susceptibility to Cry1AcF toxin. Insects pre-treated with prohibitin and α-amylase dsRNA exhibited significant reduction in Cry1AcF-induced mortality, suggesting their probable role as Cry receptor. By contrast, insects pre-treated with GLTP, ADAM and UDP-GT dsRNA exhibited no significant decline in mortality. This maybe explained by the possibility of RNAi feedback regulation (as few of the receptors belong to multigene family) or redundant role of GLTP, ADAM and UDP-GT in Cry intoxication process. CONCLUSION: Since the laboratory culture of G. mellonella develop Bt resistance quite rapidly, findings of the current investigation may provide some useful information for future Cry receptor-related research in the model insect.

BACKGROUND: Somatic copy number alterations (SCNAs) are an important class of genomic alteration in cancer. They are frequently observed in cancer samples, with studies showing that, on average, SCNAs affect 34% of a cancer cell's genome. Furthermore, SCNAs have been shown to be major drivers of tumour development and have been associated with response to therapy and prognosis. Large-scale cancer genome studies suggest that tumours are driven by somatic copy number alterations (SCNAs) or single-nucleotide variants (SNVs). Despite the frequency of SCNAs and their clinical relevance, the use of genomics assays in the clinic is biased towards targeted gene panels, which identify SNVs but provide limited scope to detect SCNAs throughout the genome. There is a need for a comparably low-cost and simple method for high-resolution SCNA profiling. RESULTS: We present conliga, a fully probabilistic method that infers SCNA profiles from a low-cost, simple, and clinically-relevant assay (FAST-SeqS). When applied to 11 high-purity oesophageal adenocarcinoma samples, we obtain good agreement (Spearman's rank correlation coefficient, rs=0.94) between conliga's inferred SCNA profiles using FAST-SeqS data (approximately £14 per sample) and those inferred by ASCAT using high-coverage WGS (gold-standard). We find that conliga outperforms CNVkit (rs=0.89), also applied to FAST-SeqS data, and is comparable to QDNAseq (rs=0.96) applied to low-coverage WGS, which is approximately four-fold more expensive, more laborious and less clinically-relevant. By performing an in silico dilution series experiment, we find that conliga is particularly suited to detecting SCNAs in low tumour purity samples. At two million reads per sample, conliga is able to detect SCNAs in all nine samples at 3% tumour purity and as low as 0.5% purity in one sample. Crucially, we show that conliga's hidden state information can be used to decide when a sample is abnormal or normal, whereas CNVkit and QDNAseq cannot provide this critical information. CONCLUSIONS: We show that conliga provides high-resolution SCNA profiles using a convenient, low-cost assay. We believe conliga makes FAST-SeqS a more clinically valuable assay as well as a useful research tool, enabling inexpensive and fast copy number profiling of pre-malignant and cancer samples.

Dehydration responsive element binding protein (DREB) is a significant transcription factor class known to be implicated in abiotic stresses. In this study, we systematically conducted a genome-wide identification and expression analysis of the DREB gene family, including gene structures, evolutionary relationships, chromosome distribution, conserved domains, and expression patterns. A total of 65 DREB family gene members were identified in Chinese cabbage (Brassica rapa L.) and were classified into five subgroups based on phylogenetic analysis. Through analysis of the conserved domains of BrDREB family genes, only one exon existed in the gene structure. Through the analysis of cis-acting elements, these genes were mainly involved in hormone regulation and adversity stress. In order to identify the function of BrDREB2B, overexpressed transgenic Arabidopsis was constructed. After different stress treatments, the germination rate, root growth, survival rate, and various plant physiological indicators were measured. The results showed that transgenic Arabidopsis thaliana plants overexpressing BrDREB2B exhibited enhanced tolerance to salt, heat and drought stresses. Taken together, our results are the first to report the BrDREB2B gene response to drought and heat stresses in Chinese cabbage and provide a basis for further studies to determine the function of BrDREBs in response to abiotic stresses.

BACKGROUND: Longissimus muscle samples were collected from lambs exposed in utero to mycotoxins [E-, endophyte-free tall fescue seed without ergot alkaloids (negative control) or E + , endophyte-infected tall fescue seed containing ergot alkaloids] during mid-gestation (MID; E + /E-) or late-gestation (LATE; E-/E +) harvested at two developmental stages (FETAL, gestational d133) or (MAT, near maturity, 250 d of age; n = 3/treatment/developmental stage). Muscle samples were examined to determine the impact of in utero mycotoxin exposure on skeletal muscle fiber hypertrophy and the miRNA profile at FETAL and MAT. RESULTS: Longissimus weight was greater (P < 0.05) in E + /E- lambs compared to E-/E + lambs at MAT; however, FETAL longissimus weight did not differ (P > 0.10) between fescue treatments. Type I fiber cross sectional area was larger (P < 0.10) for E + /E- than E-/E + at MAT but did not differ (P > 0.10) between fescue treatments at FETAL. Type II fiber area was larger (P < 0.05) at MAT in E + /E- compared to E-/E + but did not differ (P < 0.05) between fescue treatments at FETAL. Cross-sectional Type I and Type II longissimus muscle fiber area increased (P < 0.05) from FETAL to MAT by 6.86-fold and 10.83-fold, respectively. The ratio of Type II:Type I muscle fibers was lower (P = 0.04) at MAT compared to FETAL. There were 120 miRNA differentially expressed (q < 0.05) between FETAL and MAT. Maternal fescue treatment did not alter (q > 0.05) expression of miRNAs in the longissimus muscle. miR-133, -29a, -22-3p, and -410-3p were identified as highly significant with a log2 fold change > 4. In vitro satellite cell cultures showed that selected miRNAs (miR-22-3p, 29a, 27a, and 133a) are differentially regulated during proliferation and differentiation indicating a role of miRNA in muscle hypertrophy. CONCLUSIONS: Exposure to mycotoxins did not alter fiber type but had long-term impacts on postnatal muscle hypertrophy and cross-sectional area. The miRNA profile of the longissimus was not altered by Maternal mycotoxin exposure at FETAL or MAT. Developmental age altered the miRNA transcriptome and mRNA expression of known genes related to muscle growth. These results indicate that Maternal exposure to E + fescue seed during LATE gestation can alter postnatal muscle hypertrophy in sheep; however, these changes are not regulated by the miRNA transcriptome of the longissimus muscle.

BACKGROUND: Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively. RESULTS: These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development. CONCLUSIONS: This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

BACKGROUND: Monochasma savatieri Franch. ex Maxim is a medicinally valuable herb. However, the collection and protection of the wild germplasm resources of M. savatieri are still insufficient, and their genetic diversity and population structure have been poorly studied. RESULTS: We collected and examined 46 M. savatieri individuals from Fujian, Hunan, Jiangxi, and Zhejiang provinces for genetic diversity and population structure, using 33 newly developed expressed sequence tag-simple sequence repeat (EST-SSR) markers. Applying these markers, we detected a total of 208 alleles, with an average of 6.303 alleles per locus. The polymorphic information content varied from 0.138 to 0.884 (average: 0.668), indicating a high level of polymorphism. At the population level, there was a low degree of genetic diversity among populations (I = 0.535, He = 0.342), with Zhejiang individuals showing the highest genetic diversity among the four populations (Fst = 0.497), which indicated little gene flow within the M. savatieri populations (Nm = 0.253). Mantel test analysis revealed a significant positive correlation between geographical and genetic distance among populations (R2 = 0.3304, p < 0.05), and structure and principal coordinate analyses supported classification of populations into three clusters, which was consistent with the findings of cluster analysis. CONCLUSIONS: As a rare medicinal plants, the protection of M. savatieri does not look optimistic, and accordingly, protective efforts should be beefed up on the natural wild populations. This study provided novel tools and insights for designing effective collection and conservation strategies for M. savatieri.

BACKGROUND: Maize kernel row number (KRN) is one of the most important yield traits and has changed greatly during maize domestication and selection. Elucidating the genetic basis of KRN will be helpful to improve grain yield in maize. RESULTS: Here, we measured KRN in four environments using a nested association mapping (NAM) population named HNAU-NAM1 with 1,617 recombinant inbred lines (RILs) that were derived from 12 maize inbred lines with a common parent, GEMS41. Then, five consensus quantitative trait loci (QTLs) distributing on four chromosomes were identified in at least three environments along with the best linear unbiased prediction (BLUP) values by the joint linkage mapping (JLM) method. These QTLs were further validated by the separate linkage mapping (SLM) and genome-wide association study (GWAS) methods. Three KRN genes cloned through the QTL assay were found in three of the five consensus QTLs, including qKRN1.1, qKRN2.1 and qKRN4.1. Two new QTLs of KRN, qKRN4.2 and qKRN9.1, were also identified. On the basis of public RNA-seq and genome annotation data, five genes highly expressed in ear tissue were considered candidate genes contributing to KRN. CONCLUSIONS: This study carried out a comprehensive analysis of the genetic architecture of KRN by using a new NAM population under multiple environments. The present results provide solid information for understanding the genetic components underlying KRN and candidate genes in qKRN4.2 and qKRN9.1. Single-nucleotide polymorphisms (SNPs) closely linked to qKRN4.2 and qKRN9.1 could be used to improve inbred yield during molecular breeding in maize.

BACKGROUND: Carcass backfat thickness (BFT), carcass lean percentage (CLP) and carcass fat percentage (CFP) are important to the commercial pig industry. Nevertheless, the genetic architecture of BFT, CLP and CFP is still elusive. Here, we performed a genome-wide association study (GWAS) based on specific-locus amplified fragment sequencing (SLAF-seq) to analyze seven fatness-related traits, including five BFTs, CLP, and CFP on 223 four-way crossbred pigs. RESULTS: A total of 227, 921 highly consistent single nucleotide polymorphisms (SNPs) evenly distributed throughout the genome were used to perform GWAS. Using the mixed linear model (MLM), a total of 20 SNP loci significantly related to these traits were identified on ten Sus scrofa chromosomes (SSC), of which 10 SNPs were located in previously reported quantitative trait loci (QTL) regions. On SSC7, two SNPs (SSC7:29,503,670 and rs1112937671) for average backfat thickness (ABFT) exceeded 1% and 10% Bonferroni genome-wide significance levels, respectively. These two SNP loci were located within an intron region of the COL21A1 gene, which was a protein-coding gene that played an important role in the porcine backfat deposition by affecting extracellular matrix (ECM) remodeling. In addition, based on the other three significant SNPs on SSC7, five candidate genes, ZNF184, ZNF391, HMGA1, GRM4 and NUDT3 were proposed to influence BFT. On SSC9, two SNPs for backfat thickness at 6-7 ribs (67RBFT) and one SNP for CLP were in the same locus region (19 kb interval). These three SNPs were located in the PGM2L1 gene, which encoded a protein that played an indispensable role in glycogen metabolism, glycolysis and gluconeogenesis as a key enzyme. Finally, one significant SNP on SSC14 for CLP was located within the PLBD2 gene, which participated in the lipid catabolic process. CONCLUSIONS: A total of two regions on SSC7 and SSC9 and eight potential candidate genes were found for fatness-related traits in pigs. The results of this GWAS based on SLAF-seq will greatly advance our understanding of the genetic architecture of BFT, CLP, and CFP traits. These identified SNP loci and candidate genes might serve as a biological basis for improving the important fatness-related traits of pigs.

BACKGROUND: Hypertension-induced cardiac hypertrophy is one of the most common pre-conditions that accompanies heart failure. This study aimed to identify the key pathogenic genes in the disease process. METHODS: GSE18224 was re-analyzed and differentially expressed genes (DEGs) were obtained. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out. Networks of transcription factor (TF)-mRNA, microRNA (miRNA)-mRNA and Protein-Protein interaction (PPI) were constructed, and a key module was further screened out from PPI network. GSE36074 dataset and our transverse aortic constriction (TAC) mouse model were used to validate gene expression in the module. Finally, the correlation between the genes and biomarkers of cardiac hypertrophy were evaluated. RESULTS: Totally, there were 348 DEGs in GSE18224, which were mainly enriched in biological processes including collagen fibril organization, cellular response to transforming growth factor-beta stimulus and were involved in ECM-receptor interaction and Oxytocin signaling pathway. There were 387 miRNAs targeted by 257 DEGs, while 177 TFs targeted 71 DEGs. The PPI network contained 222 nodes and 770 edges, with 18 genes screened out into the module. After validation, 8 genes, which were also significantly upregulated in the GSE36074 dataset, were selected from the 18 DEGs. 2 of the 8 DEGs, including Eln and Tgfb3 were significantly upregulated in our mouse model of myocardial hypertrophy. Finally, the expression of Eln and Tgfb3 were found to be positively correlated with the level of the disease biomarkers. CONCLUSIONS: Upregulated key genes Eln and Tgfb3 were positively correlated with the severity of cardiac hypertrophy, which may provide potential therapeutic targets for the disease.

BACKGROUND: Heart failure (HF) is a prevalent cause of mortality and morbidity. The molecular drivers of HF are still largely unknown. RESULTS: We aimed to identify circulating proteins causally associated with HF by leveraging genome-wide genetic association data for HF including 47,309 cases and 930,014 controls. We performed two-sample Mendelian randomization (MR) with multiple cis instruments as well as network and enrichment analysis using data from blood protein quantitative trait loci (pQTL) (2,965 blood proteins) measured in 3,301 individuals. Nineteen blood proteins were causally associated with HF, were not subject to reverse causality and were enriched in ligand-receptor and glycosylation molecules. Network pathway analysis of the blood proteins showed enrichment in NF-kappa B, TGF beta, lipid in atherosclerosis and fluid shear stress. Cross-phenotype analysis of HF identified genetic overlap with cardiovascular drugs, myocardial infarction, parental longevity and low-density cholesterol. Multi-trait MR identified causal associations between HF-associated blood proteins and cardiovascular outcomes. Multivariable MR showed that association of BAG3, MIF and APOA5 with HF were mediated by the blood pressure and coronary artery disease. According to the directional effect and biological action, 7 blood proteins are targets of existing drugs or are tractable for the development of novel therapeutics. Among the pathways, sialyl Lewis x and the activin type II receptor are potential druggable candidates. CONCLUSIONS: Integrative MR analyses of the blood proteins identified causally-associated proteins with HF and revealed pleiotropy of the blood proteome with cardiovascular risk factors. Some of the proteins or pathway related mechanisms could be targeted as novel treatment approach in HF.

BACKGROUND: The cardiovascular developmental process is a tightly regulated network involving multiple genes. The current understanding of the molecular mechanism behind cardiovascular development is insufficient and requires further research. RESULTS: Transcriptome sequencing of three developmental stages in zebrafish embryos was performed and revealed three key cardiovascular developmental stages. Then, the differentially expressed genes (DEGs) involved in cardiovascular development were screened out. The three developmental stages were 18 (T1), 24 (T2), and 42 h post fertilization (hpf) (T3), and the three stages were confirmed by detecting differences in expression between cardiomyocyte and endothelial marker genes (cmlc2, fli1) using in situ hybridization, which represents the characteristics of cardiovascular development. Thousands of DEGs were identified using transcriptome analysis. Of them, 2605 DEGs were in T1-vs-T2, including 2003 up-regulated and 602 down-regulated genes, 6446 DEGs were in T1-vs-T3, consisting of 4608 up-regulated and 1838 down-regulated genes, and 3275 DEGs were in T2-vs-T3, including 2420 up-regulated and 855 down-regulated genes. There were 644 common DEGs and 167 common five-fold higher differentially expressed genes (HDEGs) identified, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID). Significant differences was observed in the levels of gene expression among different developmental stages in multiple GO terms and KEGG pathways, such as cell migration to the midline involved in heart development, cardiovascular system development, circulatory system process for biological processes of GO terms; and cardiac muscle contraction, adrenergic signaling in cardiomyocytes for KEGG pathways. These results demonstrated that these three stages were important period for the development of the cardiovascular system. Lastly, we used quantitative real-time PCR (qPCR) to validate the reliability of RNA-sequencing by selecting 21 DEGs. CONCLUSIONS: These results demonstrated that these three stages represented the important periods for cardiovascular system development of zebrafish and some candidate genes was obtained and provided a solid foundation for additional functional studies of the DEGs.

BACKGROUND: In the life history of Daphnia, the reproductive mode of parthenogenesis and sexual reproduction alternate in aquatic ecosystem, which are often affected by environmental and genetic factors. Recently, the sex-biased genes are of great significance for clarifying the origin and evolution of reproductive transformation and the molecular regulation mechanism of sex determination in Daphnia. Although some genes on reproductive transition of Daphnia had been researched, molecular mechanism on the maintenance of sexually dimorphic phenotypes of Daphnia are still not well known, including differentially expressed genes in different life-history stages. RESULTS: In this study, four life-history stages of Daphnia sinensis, juvenile female (JF), parthenogenetic female (PF), sexual female (SF) and male (M), were performed for transcriptome, and male-biased genes were screened. A total of 110437 transcripts were obtained and assembled into 22996 unigenes. In the four life-history stages (JF, PF, SF and M), the number of unique unigenes is respectively 2863, 445, 437 and 586, and the number of common unigenes is 9708. The differentially expressed genes (DEGs) between male and other three female stages (M vs JF, M vs PF and M vs SF) were 4570, 4358 and 2855, respectively. GO gene enrichment analysis showed that the up-regulated genes in male were mainly enriched in hydrolase activity and peptidase activity. Thirty-six genes in male were significantly higher expression than in the three female stages, including one Doublesex (Dsx) gene, one laminin gene, five trypsin genes and one serine protease genes, and one chitin synthase gene and two chitinase genes. CONCLUSIONS: Our results showed that thirty-six candidate genes may be as the male-biased genes involving in the maintenance of sexually dimorphic phenotypes. This work will provide a reference for further exploring the functional genes related to sex differentiation in Daphnia species. Moreover, according to previous investigations, we thought that the expression level of functional genes may be related to the life-history stages of organisms, and may be also affected by different Daphnia species.

BACKGROUND: Porcine Epidemic Diarrhea Virus (PEDV) is a coronavirus that seriously affects the swine industry. MicroRNAs and long noncoding RNAs are two relevant non-coding RNAs (ncRNAs) class and play crucial roles in a variety of physiological processes. Increased evidence indicates a complex interaction between mRNA and ncRNA. However, our understanding of the function of ncRNA involved in host-PEDV interaction is limited. RESULTS: A total of 1,197 mRNA transcripts, 539 lncRNA transcripts, and 208 miRNA transcripts were differentially regulated at 24 h and 48 h post-infection. Gene ontology (GO) and KEGG pathway enrichment analysis showed that DE mRNAs and DE lncRNAs were mainly involved in biosynthesis, innate immunity, and lipid metabolism. Moreover, we constructed a miRNA-mRNA-pathway network using bioinformatics, including 12 DE mRNAs, 120 DE miRNAs, and 11 pathways. Finally, the target genes of DE miRNAs were screened by bioinformatics, and we constructed immune-related lncRNA-miRNA-mRNA ceRNA networks. Then, the selected DE genes were validated by qRT-PCR, which were consistent with the results from RNA-Seq data. CONCLUSIONS: This study provides the comprehensive analysis of the expression profiles of mRNAs, lncRNAs, and miRNAs during PEDV infection. We characterize the ceRNA networks which can provide new insights into the pathogenesis of PEDV.

BACKGROUND: Quantitative real time PCR (qPCR) is a powerful tool to evaluate mRNA expression level. However, reliable qPCR results require normalization with validated reference gene(s). In this study, we investigated stable reference genes in seven tissues according to four developmental stages in minipigs. Six candidate reference genes and one target gene (ACE2) were selected and qPCR was performed. BestKeeper, geNorm, NormFinder, and delta Ct method through the RefFinder web-based tool were used to evaluate the stability of candidate reference genes. To verify the selected stable genes, relative expression of ACE2 was calculated and compared with each other. RESULTS: As a result, HPRT1 and 18S genes had lower SD value, while HMBS and GAPDH genes had higher SD value in all samples. Using statistical algorithms, HPRT1 was the most stable gene, followed by 18S, β-actin, B2M, GAPDH, and HMBS. In intestine, all candidate reference genes exhibited similar patterns of ACE2 gene expression over time, whereas in liver, lung, and kidney, gene expression pattern normalized with stable reference genes differed from those normalized with less stable genes. When normalized with the most stable genes, the expression levels of ACE2 in minipigs highly increased in intestine and kidney at PND28, which is consistent with the ACE2 expression pattern in humans. CONCLUSIONS: We suggest that HPRT1 and 18S are good choices for analyzing all these samples across the seven tissues and four developmental stages. However, this study can be a reference literature for gene expression experiments using minipig because reference gene should be validated and chosen according to experimental conditions.

BACKGROUND: Carcass traits are important in pig breeding programs for improving pork production. Understanding the genetic variants underlies complex phenotypes can help explain trait variation in pigs. In this study, we integrated a weighted single-step genome-wide association study (wssGWAS) and copy number variation (CNV) analyses to map genetic variations and genes associated with loin muscle area (LMA), loin muscle depth (LMD) and lean meat percentage (LMP) in Duroc pigs. RESULTS: Firstly, we performed a genome-wide analysis for CNV detection using GeneSeek Porcine SNP50 Bead chip data of 3770 pigs. A total of 11,100 CNVs were detected, which were aggregated by overlapping 695 CNV regions (CNVRs). Next, we investigated CNVs of pigs from the same population by whole-genome resequencing. A genome-wide analysis of 21 pigs revealed 23,856 CNVRs that were further divided into three categories (851 gain, 22,279 loss, and 726 mixed), which covered 190.8 Mb (~ 8.42%) of the pig autosomal genome. Further, the identified CNVRs were used to determine an overall validation rate of 68.5% for the CNV detection accuracy of chip data. CNVR association analyses identified one CNVR associated with LMA, one with LMD and eight with LMP after applying stringent Bonferroni correction. The wssGWAS identified eight, six and five regions explaining more than 1% of the additive genetic variance for LMA, LMD and LMP, respectively. The CNVR analyses and wssGWAS identified five common regions, of which three regions were associated with LMA and two with LMP. Four genes (DOK7, ARAP1, ELMO2 and SLC13A3) were highlighted as promising candidates according to their function. CONCLUSIONS: We determined an overall validation rate for the CNV detection accuracy of low-density chip data and constructed a genomic CNV map for Duroc pigs using resequencing, thereby proving a value genetic variation resource for pig genome research. Furthermore, our study utilized a composite genetic strategy for complex traits in pigs, which will contribute to the study for elucidating the genetic architecture that may be influenced and regulated by multiple forms of variations.

BACKGROUND: Insect metamorphosis from larvae to pupae is one of the most important stages of insect life history. Relatively comprehensive information related to gene transcription profiles during lepidopteran metamorphosis is required to understand the molecular mechanism underlying this important stage. We conducted transcriptional profiling of the brain and fat body of the cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) during its transition from last instar larva into pupa to explore the physiological processes associated with different phases of metamorphosis. RESULTS: During metamorphosis, the differences in gene expression patterns and the number of differentially expressed genes in the fat body were found to be greater than those in the brain. Each stage had a specific gene expression pattern, which contributed to different physiological changes. A decrease in juvenile hormone levels at the feeding stage is associated with increased expression levels of two genes (juvenile hormone esterase, juvenile hormone epoxide hydrolase). The expression levels of neuropeptides were highly expressed at the feeding stage and the initiation of the wandering stage and less expressed at the prepupal stage and the initiation of the pupal stage. The transcription levels of many hormone (or neuropeptide) receptors were specifically increased at the initiation of the wandering stage in comparison with other stages. The expression levels of many autophagy-related genes in the fat body were found to be gradually upregulated during metamorphosis. The activation of apoptosis was probably related to enhanced expression of many key genes (Apaf1, IAP-binding motif 1 like, cathepsins, caspases). Active proliferation might be associated with enhanced expression levels in several factors (JNK pathway: jun-D; TGF-β pathway: decapentaplegic, glass bottom boat; insulin pathway: insulin-like peptides from the fat body; Wnt pathway: wntless, TCF/Pangolin). CONCLUSIONS: This study revealed several vital physiological processes and molecular events of metamorphosis and provided valuable information for illustrating the process of insect metamorphosis from larvae to pupae.

BACKGROUND: Mitochondrial genomes are the most sequenced genomes after bacterial and fungal genomic DNA. However, little information on mitogenomes is available for multiple metazoan taxa, such as Culicoides, a globally distributed, megadiverse genus containing 1,347 species. AIM:  Generating novel mitogenomic information from single Culicoides sonorensis and C. biguttatus specimens, comparing available mitogenome mapping and de novo assembly tools, and identifying the best performing strategy and tools for Culicoides species. RESULTS: We present two novel and fully annotated mitochondrial haplotypes for two Culicoides species, C. sonorensis and C. biguttatus. We also annotated or re-annotated the only available reference mitogenome for C. sonorensis and C. arakawae. All species present a high similarity in mitogenome organization. The general gene arrangement for all Culicoides species was identical to the ancestral insect mitochondrial genome. Only short spacers were found in C. sonorensis (up to 30 bp), contrary to C. biguttatus (up to 114 bp). The mitochondrial genes ATP8, NAD2, NAD6, and LSU rRNA exhibited the highest nucleotide diversity and pairwise interspecific p genetic distance, suggesting that these genes might be suitable and complementary molecular barcodes for Culicoides identification in addition to the commonly utilized COI gene. We observed performance differences between the compared mitogenome generation strategies. The mapping strategy outperformed the de novo assembly strategy, but mapping results were partially biased in the absence of species-specific reference mitogenome. Among the utilized tools, BWA performed best for C. sonorensis while SPAdes, MEGAHIT, and MitoZ were among the best for C. biguttatus. The best-performing mitogenome annotator was MITOS2. Additionally, we were able to recover exogenous mitochondrial DNA from Bos taurus (biting midges host) from a C. biguttatus blood meal sample. CONCLUSIONS: Two novel annotated mitogenome haplotypes for C. sonorensis and C. biguttatus using High-Throughput Sequencing are presented. Current results are useful as the baseline for mitogenome reconstruction of the remaining Culicoides species from single specimens to HTS and genome annotation. Mapping to a species-specific reference mitogenome generated better results for Culicoides mitochondrial genome reconstruction than de novo assembly, while de novo assembly resulted better in the absence of a closely related reference mitogenome. These results have direct implications for molecular-based identification of these vectors of human and zoonotic diseases, setting the basis for using the whole mitochondrial genome as a marker in Culicoides identification.

BACKGROUND: High temperature induces early bolting in lettuce (Lactuca sativa L.), which affects both quality and production. However, the molecular mechanism underlying high temperature-induced bolting is still limited. RESULTS: We performed systematical analysis of morphology, transcriptome, miRNAs and methylome in lettuce under high temperature treatment. Through a comparison of RNA-Seq data between the control and the high temperature treated lettuces at different time points totally identified 2944 up-regulated genes and 2203 down-regulated genes, which cover three floral pathways including photoperiod, age and gibberellin (GA) pathways. Genome wide analysis of miRNAs and methylome during high temperature treatment indicated miRNAs and DNA methylation might play a role controlling gene expression during high temperature-induced bolting. miRNA targets included some protein kinase family proteins, which potentially play crucial roles in this process. CONCLUSIONS: Together, our results propose a possible regulation network involved in high temperature-induced bolting.

BACKGROUND: Circular RNAs (CircRNAs) play critical roles in gene expression regulation and disease development. Understanding the regulation mechanism of CircRNAs formation can help reveal the role of CircRNAs in various biological processes mentioned above. Back-splicing is important for CircRNAs formation. Back-splicing sites prediction helps uncover the mysteries of CircRNAs formation. Several methods were proposed for back-splicing sites prediction or circRNA-realted prediction tasks. Model performance was constrained by poor feature learning and using ability. RESULTS: In this study, CircCNN was proposed to predict pre-mRNA back-splicing sites. Convolution neural network and batch normalization are the main parts of CircCNN. Experimental results on three datasets show that CircCNN outperforms other baseline models. Moreover, PPM (Position Probability Matrix) features extract by CircCNN were converted as motifs. Further analysis reveals that some of motifs found by CircCNN match known motifs involved in gene expression regulation, the distribution of motif and special short sequence is important for pre-mRNA back-splicing. CONCLUSIONS: In general, the findings in this study provide a new direction for exploring CircRNA-related gene expression regulatory mechanism and identifying potential targets for complex malignant diseases. The datasets and source code of this study are freely available at: https://github.com/szhh521/CircCNN .

BACKGROUND: Sequencing viruses in many specimens is hindered by excessive background material from hosts, microbiota, and environmental organisms. Consequently, enrichment of target genomic material is necessary for practical high-throughput viral genome sequencing. Hybridization probes are widely used for enrichment in many fields, but their application to viral sequencing faces a major obstacle: it is difficult to design panels of probe oligo sequences that broadly target many viral taxa due to their rapid evolution, extensive diversity, and genetic hypervariability. To address this challenge, we created ProbeTools, a package of bioinformatic tools for generating effective viral capture panels, and for assessing coverage of target sequences by probe panel designs in silico. In this study, we validated ProbeTools by designing a panel of 3600 probes for subtyping the hypervariable haemagglutinin (HA) and neuraminidase (NA) genome segments of avian-origin influenza A viruses (AIVs). Using in silico assessment of AIV reference sequences and in vitro capture on egg-cultured viral isolates, we demonstrated effective performance by our custom AIV panel and ProbeTools' suitability for challenging viral probe design applications. RESULTS: Based on ProbeTool's in silico analysis, our panel provided broadly inclusive coverage of 14,772 HA and 11,967 NA reference sequences. For each reference sequence, we calculated the percentage of nucleotide positions covered by our panel in silico; 90% of HA and NA references sequences had at least 90.8 and 95.1% of their nucleotide positions covered respectively. We also observed effective in vitro capture on a representative collection of 23 egg-cultured AIVs that included isolates from wild birds, poultry, and humans and representatives from all HA and NA subtypes. Forty-two of forty-six HA and NA segments had over 98.3% of their nucleotide positions significantly enriched by our custom panel. These in vitro results were further used to validate ProbeTools' in silico coverage assessment algorithm; 89.2% of in silico predictions were concordant with in vitro results. CONCLUSIONS: ProbeTools generated an effective panel for subtyping AIVs that can be deployed for genomic surveillance, outbreak prevention, and pandemic preparedness. Effective probe design against hypervariable AIV targets also validated ProbeTools' design and coverage assessment algorithms, demonstrating their suitability for other challenging viral capture applications.

Exome sequencing is becoming a routine in health care, because it increases the chance of pinpointing the genetic cause of an individual patient's condition and thus making an accurate diagnosis. It is important for facilities providing genetic services to keep track of changes in the technology of exome capture in order to maximize throughput while reducing cost per sample. In this study, we focused on comparing the newly released exome probe set Agilent SureSelect Human All Exon v8 and the previous probe set v7. In preparation for higher throughput of exome sequencing using the DNBSEQ-G400, we evaluated target design, coverage statistics, and variants across these two different exome capture products. Although the target size of the v8 design has not changed much compared to the v7 design (35.24 Mb vs 35.8 Mb), the v8 probe design allows you to call more of SNVs (+ 3.06%) and indels (+ 8.49%) with the same number of raw reads per sample on the common target regions (34.84 Mb). Our results suggest that the new Agilent v8 probe set for exome sequencing yields better data quality than the current Agilent v7 set.

BACKGROUND: Patatin-like phospholipase domain containing 5 (PNPLA5) is a newly-discovered lipase. Although the PNPLA family plays critical roles in diverse biological processes, the biological functions of PNPLA5 mostly unknown. We previously found that the deletion of Pnpla5 in rats causes a variety of phenotypic abnormalities. In this study, we further explored the effects of Pnpla5 knockout (KO) on male rats. RESULTS: The body weight and testicular or epididymal tissue weight of three to six 3-month-old Pnpla5 KO or wild-type (WT) male Sprague-Dawley rats were measured. The protein expression levels were also measured via western blotting and iTRAQ (isobaric tags for relative and absolute quantitation) analyses. No significant difference between Pnpla5 KO and WT rats, regarding body weight, testicular or epididymal tissue weight, or hormone levels, were found. However, the relative testicular tissue weight of the KO (Pnpla5-/-) rats was higher (P < 0.05) than that of WT rats. Significant increases in apoptotic cells numbers (P < 0.001) and BAX and Caspase-9 expression levels were observed in the testicular tissue of Pnpla5-/- rats. Moreover, iTRAQ analysis revealed that the levels of proteins involved in steroid metabolism and wound healing were significantly decreased in Pnpla5-/- rats. CONCLUSION: This study revealed that Pnpla5 knockout induced apoptosis in rat testes. We also ascertained that Pnpla5 plays an important role in lipid metabolism, wound healing, and affects reproductive organs negatively, providing new target genes and pathways that can be analyzed to unravel the biological function of Pnpla5.

AIM: To analyze and compare the mRNA N6-methyladenosine modifications in transverse aortic constriction induced mice hearts and normal mice hearts. MATERIALS AND METHODS: Colorimetric quantification was used to probe the changes in m6A modifications in the total RNA. The expression of m6A-related enzymes was analyzed via qRT-PCR and western blotting. RNA-seq and MeRIP-seq were performed to identify genes with differences in m6A modifications or expression in the transcriptome profile. RESULTS: Compared with the control group, the TAC group exhibited higher m6A methylation levels. FTO and WTAP were downregulated after TAC, while METTL3 was significantly downregulated at the protein level. MeRIP-seq revealed that 1179 m6A peaks were upmethylated and 733 m6A peaks were downmethylated, and biological analysis of these genes exhibited a strong relationship with heart function. CONCLUSION: Our findings provide novel information regarding m6A modification and gene expression changes in cardiac hypertrophy, which may be fundamental for further research.

Berberine hydrochloride is the main effective component of Coptis spp. used in Chinese herbal medicine and its underlying molecular mechanisms, responsible for inducing effects in crustacean species, are not fully understood. In this study, the molecular response of the crab Charybdis japonica to berberine hydrochloride exposure was studied using transcriptome sequencing. The survival rate, gene expression and activities of several immune enzymes were measured after berberine hydrochloride treatments, with or without injection of the pathogenic bacterium Aeromonas hydrophila. A total of 962 differentially expressed genes (464 up-regulated and 498 down-regulated) were observed during exposure to 100 mg/L of berberine hydrochloride and in the control group after 48 h. Enrichment analysis revealed that these genes are involved in metabolism, cellular processes, signal transduction and immune functions, indicating that exposure to berberine hydrochloride activated the immune complement system. This bioactive compound simultaneously activated fibrinogen beta (FGB), fibrinogen alpha (FGA), alpha-2-macroglobulin (A2M), kininogen (KNG), fibrinogen gamma chain (FGB), alpha-2-HS-glycoprotein (AHSG), caspase-8 (CASP8), cathepsin L (CTSL), adenylate cyclase 3 (Adcy3) and MMP1. Its action could significantly increase the survival rate of the crabs injected with A. hydrophila and promote the activity of LZM, Caspas8, FGA, ACP and AKP in the hepatopancreas. When A. hydrophila was added, the neutralization of 300 mg/L berberine hydrochloride maximized the activities of Caspas8, LZM, ACP and AKP. Our results provide a new understanding of the potential effects of berberine hydrochloride on the immune system mechanisms in crustaceans.

BACKGROUND: CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated proteins) systems are adaptive immune systems commonly found in prokaryotes that provide sequence-specific defense against invading mobile genetic elements (MGEs). The memory of these immunological encounters are stored in CRISPR arrays, where spacer sequences record the identity and history of past invaders. Analyzing such CRISPR arrays provide insights into the dynamics of CRISPR-Cas systems and the adaptation of their host bacteria to rapidly changing environments such as the human gut. RESULTS: In this study, we utilized 601 publicly available Bacteroides fragilis genome isolates from 12 healthy individuals, 6 of which include longitudinal observations, and 222 available B. fragilis reference genomes to update the understanding of B. fragilis CRISPR-Cas dynamics and their differential activities. Analysis of longitudinal genomic data showed that some CRISPR array structures remained relatively stable over time whereas others involved radical spacer acquisition during some periods, and diverse CRISPR arrays (associated with multiple isolates) co-existed in the same individuals with some persisted over time. Furthermore, features of CRISPR adaptation, evolution, and microdynamics were highlighted through an analysis of host-MGE network, such as modules of multiple MGEs and hosts, reflecting complex interactions between B. fragilis and its invaders mediated through the CRISPR-Cas systems. CONCLUSIONS: We made available of all annotated CRISPR-Cas systems and their target MGEs, and their interaction network as a web resource at https://omics.informatics.indiana.edu/CRISPRone/Bfragilis . We anticipate it will become an important resource for studying of B. fragilis, its CRISPR-Cas systems, and its interaction with mobile genetic elements providing insights into evolutionary dynamics that may shape the species virulence and lead to its pathogenicity.

BACKGROUND: Echinoderms are established models in experimental and developmental biology, however genomic resources are still lacking for many species. Here, we present the draft genome of Ophioderma brevispinum, an emerging model organism in the field of regenerative biology. This new genomic resource provides a reference for experimental studies of regenerative mechanisms. RESULTS: We report a de novo nuclear genome assembly for the brittle star O. brevispinum and annotation facilitated by the transcriptome assembly. The final assembly is 2.68 Gb in length and contains 146,703 predicted protein-coding gene models. We also report a mitochondrial genome for this species, which is 15,831 bp in length, and contains 13 protein-coding, 22 tRNAs, and 2 rRNAs genes, respectively. In addition, 29 genes of the Notch signaling pathway are identified to illustrate the practical utility of the assembly for studies of regeneration. CONCLUSIONS: The sequenced and annotated genome of O. brevispinum presented here provides the first such resource for an ophiuroid model species. Considering the remarkable regenerative capacity of this species, this genome will be an essential resource in future research efforts on molecular mechanisms regulating regeneration.

BACKGROUND: As a valuable medicinal plant, Rhodiola has a very long history of folk medicine used as an important adaptogen, tonic, and hemostatic. However, our knowledge of the chloroplast genome level of Rhodiola is limited. This drawback has limited studies on the identification, evolution, genetic diversity and other relevant studies on Rhodiola. RESULTS: Six Rhodiola complete chloroplast genomes were determined and compared to another Rhodiola cp genome at the genome scale. The results revealed a cp genome with a typical quadripartite and circular structure that ranged in size from 150,771 to 151,891 base pairs. High similarity of genome organization, gene number, gene order, and GC content were found among the chloroplast genomes of Rhodiola. 186 (R. wallichiana) to 200 (R. gelida) SSRs and 144 pairs of repeats were detected in the 6 Rhodiola cp genomes. Thirteen mutational hotspots for genome divergence were determined and could be used as candidate markers for phylogenetic analyses and Rhodiola species identification. The phylogenetic relationships inferred by members of Rhodiola cluster into two clades: dioecious and hermaphrodite. Our findings are helpful for understanding Rhodiola's taxonomic, phylogenetic, and evolutionary relationships. CONCLUSIONS: Comparative analysis of chloroplast genomes of Rhodiola facilitates medicinal resource conservation, phylogenetic reconstruction and biogeographical research of Rhodiola.

BACKGROUND: Genetics studies in the porcine immune system have enhanced selection practices for disease resistance phenotypes and increased the efficacy of porcine models in biomedical research; however limited functional annotation of the porcine immunome has hindered progress on both fronts. Among epigenetic mechanisms that regulate gene expression, DNA methylation is the most ubiquitous modification made to the DNA molecule and influences transcription factor binding as well as gene and phenotype expression. Human and mouse DNA methylation studies have improved mapping of regulatory elements in these species, but comparable studies in the pig have been limited in scope. RESULTS: We performed whole-genome bisulfite sequencing to assess DNA methylation patterns in nine pig immune cell populations: CD21+ and CD21- B cells, four T cell fractions (CD4+, CD8+, CD8+CD4+, and SWC6γδ+), natural killer and myeloid cells, and neutrophils. We identified 54,391 cell differentially methylated regions (cDMRs), and clustering by cDMR methylation rate grouped samples by cell lineage. 32,737 cDMRs were classified as cell lowly methylated regions (cLMRs) in at least one cell type, and cLMRs were broadly enriched in genes and regions of intermediate CpG density. We observed strong correlations between differential methylation and expression across immune cell populations, with cell-specific low methylation disproportionately impacting genes exhibiting enriched gene expression in the same cell type. Motif analysis of cLMRs revealed cell type-specific enrichment of transcription factor binding motifs, indicating that cell-specific methylation patterns may influence accessibility by trans-acting factors. Lastly, cDMRs were enriched for immune capacity GWAS SNPs, and many such overlaps occurred within genes known to influence immune cell development and function (CD8B, NDRG1). CONCLUSION: Our DNA methylation data improve functional annotation of the porcine genome through characterization of epigenomic regulatory patterns that contribute to immune cell identity and function, and increase the potential for identifying mechanistic links between genotype and phenotype.

BACKGROUND: Skeletal muscle tissue is among the largest organ systems in mammals, essential for survival and movement. Embryonic muscle development determines the quantity and quality of muscles after the birth of an individual. MicroRNAs (miRNAs) are a significant class of non-coding RNAs that bind to the 3'UTR region of mRNA to regulate gene function. Total RNA was extracted from the leg muscles of chicken embryos in different developmental stages of Chengkou Mountain Chicken and used to generate 171,407,341 clean small RNA reads. Target prediction, GO, and KEGG enrichment analyses determined the significantly enriched genes and pathways. Differential analysis determined the significantly different miRNAs between chicken embryo leg muscles at different developmental stages. Meanwhile, the weighted correlation network analysis (WGCNA) identified key modules in different developmental stages, and the hub miRNAs were screened following the KME value. RESULTS: The clean reads contained 2047 miRNAs, including 721 existing miRNAs, 1059 known miRNAs, and 267 novel miRNAs. Many genes and pathways related to muscle development were identified, including ERBB4, MEF2C, FZD4, the Wnt, Notch, and MAPK signaling pathways. The WGCNA established the greenyellow module and gga-miR-130b-5p for E12, magenta module and gga-miR-1643-5p for E16, purple module and gga-miR-12218-5p for E19, cyan module and gga-miR-132b-5p for E21. CONCLUSION: These results lay a foundation for further research on the molecular regulatory mechanism of embryonic muscle development in Chengkou mountain chicken and provide a reference for other poultry and livestock muscle development studies.

BACKGROUND: Mammals rely on the circadian clock network to regulate daily systemic metabolism and physiological activities. The liver is an important peripheral organ in mammals, and it has a unique circadian rhythm regulation process. As the only mammals that can fly, bats have attracted much research attention due to their nocturnal habits and life histories. However, few research reports exist concerning the circadian rhythms of bat liver gene expression and the relevant biological clock regulation mechanisms in the liver. RESULTS: In this study, the expression levels of liver genes of Asian particolored bats were comparatively analyzed using RNA-seq at four different time points across 24 h. A total of 996 genes were found to be rhythmic, accounting for 65% of the total number of expressed genes. The critical circadian rhythm genes Bmal1, Rev-erbα, Cry, and Ror in the liver exhibited different expression patterns throughout the day, and participated in physiological processes with rhythmic changes, including Th17 cell differentiation (ko04659), antigen processing and presentation (ko04612), the estrogen signaling pathway (ko04915), and insulin resistance (ko04931). In addition, previous studies have found that the peroxisome proliferator-activated receptor (PPAR) metabolic signaling pathway (ko03320) may play a vital role in the rhythmic regulation of the metabolic network. CONCLUSIONS: This study is the first to demonstrate diurnal changes in bat liver gene expression and related physiological processes. The results have thus further enriched our understanding of bats' biological clocks.

BACKGROUND: Saposhnikovia divaricata (Turcz.) Schischk. is a perennial herb whose dried roots are commonly used as a source of traditional medicines. To elucidate the organelle-genome-based phylogeny of Saposhnikovia species and the transfer of DNA between organelle genomes, we sequenced and characterised the mitochondrial genome (mitogenome) of S. divaricata. RESULTS: The mitogenome of S. divaricata is a circular molecule of 293,897 bp. The nucleotide composition of the mitogenome is as follows: A, 27.73%; T, 27.03%; C, 22.39%; and G, 22.85. The entire gene content is 45.24%. A total of 31 protein-coding genes, 20 tRNAs and 4 rRNAs, including one pseudogene (rpl16), were annotated in the mitogenome. Phylogenetic analysis of the organelle genomes from S. divaricata and 10 related species produced congruent phylogenetic trees. Selection pressure analysis revealed that most of the mitochondrial genes of related species are highly conserved. Moreover, 2 and 46 RNA-editing sites were found in the chloroplast genome (cpgenome) and mitogenome protein-coding regions, respectively. Finally, a comparison of the cpgenome and the mitogenome assembled from the same dataset revealed 10 mitochondrial DNA fragments with sequences similar to those in the repeat regions of the cpgenome, suggesting that the repeat regions might be transferred into the mitogenome. CONCLUSIONS: In this study, we assembled and annotated the mitogenome of S. divaricata. This study provides valuable information on the taxonomic classification and molecular evolution of members of the family Apiaceae.

BACKGROUND: OSCA (hyperosmolality-gated calcium-permeable channel) is a calcium permeable cation channel protein that plays an important role in regulating plant signal transduction. It is involved in sensing changes in extracellular osmotic potential and an increase in Ca2+ concentration. S. habrochaites is a good genetic material for crop improvement against cold, late blight, planthopper and other diseases. Till date, there is no report on OSCA in S. habrochaites. Thus, in this study, we performed a genome-wide screen to identify OSCA genes in S. habrochaites and characterized their responses to biotic and abiotic stresses. RESULTS: A total of 11 ShOSCA genes distributed on 8 chromosomes were identified. Subcellular localization analysis showed that all members of ShOSCA localized on the plasma membrane and contained multiple stress-related cis acting elements. We observed that genome-wide duplication (WGD) occurred in the genetic evolution of ShOSCA5 (Solhab04g250600) and ShOSCA11 (Solhab12g051500). In addition, repeat events play an important role in the expansion of OSCA gene family. OSCA gene family of S. habrochaites used the time lines of expression studies by qRT-PCR, do indicate OSCAs responded to biotic stress (Botrytis cinerea) and abiotic stress (drought, low temperature and abscisic acid (ABA)). Among them, the expression of ShOSCAs changed significantly under four stresses. The resistance of silencing ShOSCA3 plants to the four stresses was reduced. CONCLUSION: This study identified the OSCA gene family of S. habrochaites for the first time and analyzed ShOSCA3 has stronger resistance to low temperature, ABA and Botrytis cinerea stress. This study provides a theoretical basis for clarifying the biological function of OSCA, and lays a foundation for tomato crop improvement.

BACKGROUND: Green-fleshed radish (Raphanus sativus L.) is an economically important root vegetable of the Brassicaceae family, and chlorophyll accumulates in its root tissues. It was reported that the basic helix-loop-helix (bHLH) transcription factors play vital roles in the process of chlorophyll metabolism. Nevertheless, a comprehensive study on the bHLH gene family has not been performed in Raphanus sativus L. RESULTS: In this study, a total of 213 Raphanus sativus L. bHLH (RsbHLH) genes were screened in the radish genome, which were grouped into 22 subfamilies. 204 RsbHLH genes were unevenly distributed on nine chromosomes, and nine RsbHLH genes were located on the scaffolds. Gene structure analysis showed that 25 RsbHLH genes were intron-less. Collineation analysis revealed the syntenic orthologous bHLH gene pairs between radish and Arabidopsis thaliana/Brassica rapa/Brassica oleracea. 162 RsbHLH genes were duplicated and retained from the whole genome duplication event, indicating that the whole genome duplication contributed to the expansion of the RsbHLH gene family. RNA-seq results revealed that RsbHLH genes had a variety of expression patterns at five development stages of green-fleshed radish and white-fleshed radish. In addition, the weighted gene co-expression network analysis confirmed four RsbHLH genes closely related to chlorophyll content. CONCLUSIONS: A total of 213 RsbHLH genes were identified, and we systematically analyzed their gene structure, evolutionary and collineation relationships, conserved motifs, gene duplication, cis-regulatory elements and expression patterns. Finally, four bHLH genes closely involved in chlorophyll content were identified, which may be associated with the photosynthesis of the green-fleshed radish. The current study would provide valuable information for further functional exploration of RsbHLH genes, and facilitate clarifying the molecular mechanism underlying photosynthesis process in green-fleshed radish.