Accurate variant detection in the coding regions of the human genome is a key requirement for molecular diagnostics of Mendelian disorders. Efficiency of variant discovery from next-generation ...sequencing (NGS) data depends on multiple factors, including reproducible coverage biases of NGS methods and the performance of read alignment and variant calling software. Although variant caller benchmarks are published constantly, no previous publications have leveraged the full extent of available gold standard whole-genome (WGS) and whole-exome (WES) sequencing datasets.
In this work, we systematically evaluated the performance of 4 popular short read aligners (Bowtie2, BWA, Isaac, and Novoalign) and 9 novel and well-established variant calling and filtering methods (Clair3, DeepVariant, Octopus, GATK, FreeBayes, and Strelka2) using a set of 14 "gold standard" WES and WGS datasets available from Genome In A Bottle (GIAB) consortium. Additionally, we have indirectly evaluated each pipeline's performance using a set of 6 non-GIAB samples of African and Russian ethnicity. In our benchmark, Bowtie2 performed significantly worse than other aligners, suggesting it should not be used for medical variant calling. When other aligners were considered, the accuracy of variant discovery mostly depended on the variant caller and not the read aligner. Among the tested variant callers, DeepVariant consistently showed the best performance and the highest robustness. Other actively developed tools, such as Clair3, Octopus, and Strelka2, also performed well, although their efficiency had greater dependence on the quality and type of the input data. We have also compared the consistency of variant calls in GIAB and non-GIAB samples. With few important caveats, best-performing tools have shown little evidence of overfitting.
The results show surprisingly large differences in the performance of cutting-edge tools even in high confidence regions of the coding genome. This highlights the importance of regular benchmarking of quickly evolving tools and pipelines. We also discuss the need for a more diverse set of gold standard genomes that would include samples of African, Hispanic, or mixed ancestry. Additionally, there is also a need for better variant caller assessment in the repetitive regions of the coding genome.
Many studies on long chain non-coding RNAs (lncRNAs) are published in recent years. But the roles of lncRNAs in aortic dissection (AD) are still unclear and should be further examined. The present ...work focused on determining the molecular mechanisms underlying lncRNAs regulation in aortic dissection on the basis of the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network. This study collected the lncRNAs (GSE52093), mRNAs (GSE52093) and miRNAs (GSE92427) expression data within human tissue samples with aortic dissection group and normal group based on Gene Expression Omnibus (GEO) database. This study identified three differentially expressed lncRNAs (DELs), 19 differentially expressed miRNAs (DEmiRs) and 1046 differentially expressed mRNAs (DEGs) identified regarding aortic dissection. Furthermore, we constructed a lncRNA-miRNA-mRNA network through three lncRNAs (including two with up-regulation and one with down-regulation), five miRNAs (five with up-regulation), as well as 211 mRNAs (including 103 with up-regulation and 108 with down-regulation). Simultaneously, we conducted functional enrichment and pathway analyses on genes within the as-constructed ceRNA network. According to our PPI/ceRNA network and functional enrichment analysis results, four critical genes were found (E2F2, IGF1R, BDNF and PPP2R1B). In addition, E2F2 level was possibly modulated via lncRNA FAM87A-hsa-miR-31-5p/hsa-miR-7-5p or lncRNA C9orf106-hsa-miR-7-5p. The expression of IGF1R may be regulated by lncRNA FAM87A-hsa-miR-16-5p/hsa-miR-7-5p or lncRNA C9orf106-hsa-miR-7-5p. In conclusion, the ceRNA interaction axis we identified is a potentially critical target for treating AD. Our results shed more lights on the possible pathogenic mechanism in AD using a lncRNA-associated ceRNA network.
Emergence of antibiotic resistance is a global public health concern. The relationships between antibiotic use, the gut community composition, normal physiology and metabolism, and individual and ...public health are still being defined. Shifts in composition of bacteria, antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) after antibiotic treatment are not well-understood.
This project used next-generation sequencing, custom-built metagenomics pipeline and differential abundance analysis to study the effect of antibiotic monotherapy on resistome and taxonomic composition in the gut of Balb/c mice infected with E. coli via transurethral catheterization to investigate the evolution and emergence of antibiotic resistance.
There is a longitudinal decrease of gut microbiota diversity after antibiotic treatment. Various ARGs are enriched within the gut microbiota despite an overall reduction of the diversity and total amount of bacteria after antibiotic treatment. Sometimes treatment with a specific class of antibiotics selected for ARGs that resist antibiotics of a completely different class (e.g. treatment of ciprofloxacin or fosfomycin selected for cepA that resists ampicillin). Relative abundance of some MGEs increased substantially after antibiotic treatment (e.g. transposases in the ciprofloxacin group).
Antibiotic treatment caused a remarkable reduction in diversity of gut bacterial microbiota but enrichment of certain types of ARGs and MGEs. These results demonstrate an emergence of cross-resistance as well as a profound change in the gut resistome following oral treatment of antibiotics.
The potato genome sequence derived from the Solanum tuberosum Group Phureja clone DM1-3 516 R44 provides unparalleled insight into the genome composition and organisation of this important crop. A ...key class of genes that comprises the vast majority of plant resistance (R) genes contains a nucleotide-binding and leucine-rich repeat domain, and is collectively known as NB-LRRs.
As part of an effort to accelerate the process of functional R gene isolation, we performed an amino acid motif based search of the annotated potato genome and identified 438 NB-LRR type genes among the ~39,000 potato gene models. Of the predicted genes, 77 contain an N-terminal toll/interleukin 1 receptor (TIR)-like domain, and 107 of the remaining 361 non-TIR genes contain an N-terminal coiled-coil (CC) domain. Physical map positions were established for 370 predicted NB-LRR genes across all 12 potato chromosomes. The majority of NB-LRRs are physically organised within 63 identified clusters, of which 50 are homogeneous in that they contain NB-LRRs derived from a recent common ancestor.
By establishing the phylogenetic and positional relationship of potato NB-LRRs, our analysis offers significant insight into the evolution of potato R genes. Furthermore, the data provide a blueprint for future efforts to identify and more rapidly clone functional NB-LRR genes from Solanum species.
Tartary buckwheat (Fagopyrum tataricum Gaertn.) is a widely cultivated medicinal and edible crop with excellent economic and nutritional value. The development of tartary buckwheat seeds is a very ...complex process involving many expression-dependent physiological changes and regulation of a large number of genes and phytohormones. In recent years, the gene regulatory network governing the physiological changes occurring during seed development have received little attention.
Here, we characterized the seed development of tartary buckwheat using light and electron microscopy and measured phytohormone and nutrient accumulation by using high performance liquid chromatography (HPLC) and by profiling the expression of key genes using RNA sequencing with the support of the tartary buckwheat genome. We first divided the development of tartary buckwheat seed into five stages that include complex changes in development, morphology, physiology and phytohormone levels. At the same time, the contents of phytohormones (gibberellin, indole-3-acetic acid, abscisic acid, and zeatin) and nutrients (rutin, starch, total proteins and soluble sugars) at five stages were determined, and their accumulation patterns in the development of tartary buckwheat seeds were analyzed. Second, gene expression patterns of tartary buckwheat samples were compared during three seed developmental stages (13, 19, and 25 days postanthesis, DPA), and 9 765 differentially expressed genes (DEGs) were identified. We analyzed the overlapping DEGs in different sample combinations and measured 665 DEGs in the three samples. Furthermore, expression patterns of DEGs related to phytohormones, flavonoids, starch, and storage proteins were analyzed. Third, we noted the correlation between the trait (physiological changes, nutrient changes) and metabolites during seed development, and discussed the key genes that might be involved in the synthesis and degradation of each of them.
We provided abundant genomic resources for tartary buckwheat and Polygonaceae communities and revealed novel molecular insights into the correlations between the physiological changes and seed development of tartary buckwheat.
The phytopathogenic bacterium Xylella fastidiosa was thought to be restricted to the Americas where it infects and kills numerous hosts. Its detection worldwide has been blooming since 2013 in Europe ...and Asia. Genetically diverse, this species is divided into six subspecies but genetic traits governing this classification are poorly understood.
SkIf (Specific k-mers Identification) was designed and exploited for comparative genomics on a dataset of 46 X. fastidiosa genomes, including seven newly sequenced individuals. It was helpful to quickly check the synonymy between strains from different collections. SkIf identified specific SNPs within 16S rRNA sequences that can be employed for predicting the distribution of Xylella through data mining. Applied to inter- and intra-subspecies analyses, it identified specific k-mers in genes affiliated to differential gene ontologies. Chemotaxis-related genes more prevalently possess specific k-mers in genomes from subspecies fastidiosa, morus and sandyi taken as a whole group. In the subspecies pauca increased abundance of specific k-mers was found in genes associated with the bacterial cell wall/envelope/plasma membrane. Most often, the k-mer specificity occurred in core genes with non-synonymous SNPs in their sequences in genomes of the other subspecies, suggesting putative impact in the protein functions. The presence of two integrative and conjugative elements (ICEs) was identified, one chromosomic and an entire plasmid in a single strain of X. fastidiosa subsp. pauca. Finally, a revised taxonomy of X. fastidiosa into three major clades defined by the subspecies pauca (clade I), multiplex (clade II) and the combination of fastidiosa, morus and sandyi (clade III) was strongly supported by k-mers specifically associated with these subspecies.
SkIf is a robust and rapid software, freely available, that can be dedicated to the comparison of sequence datasets and is applicable to any field of research. Applied to X. fastidiosa, an emerging pathogen in Europe, it provided an important resource to mine for identifying genetic markers of subspecies to optimize the strategies attempted to limit the pathogen dissemination in novel areas.
Pearl millet (Pennisetum glaucum) is a cereal crop that possesses the ability to withstand drought, salinity and high temperature stresses. The NAC NAM (No Apical Meristem), ATAF1 (Arabidopsis ...thaliana Activation Factor 1), and CUC2 (Cup-shaped Cotyledon) transcription factor family is one of the largest transcription factor families in plants. NAC family members are known to regulate plant growth and abiotic stress response. Currently, no reports are available on the functions of the NAC family in pearl millet.
Our genome-wide analysis found 151 NAC transcription factor genes (PgNACs) in the pearl millet genome. Thirty-eight and 76 PgNACs were found to be segmental and dispersed duplicated respectively. Phylogenetic analysis divided these NAC transcription factors into 11 groups (A-K). Three PgNACs (- 073, - 29, and - 151) were found to be membrane-associated transcription factors. Seventeen other conserved motifs were found in PgNACs. Based on the similarity of PgNACs to NAC proteins in other species, the functions of PgNACs were predicted. In total, 88 microRNA target sites were predicted in 59 PgNACs. A previously performed transcriptome analysis suggests that the expression of 30 and 42 PgNACs are affected by salinity stress and drought stress, respectively. The expression of 36 randomly selected PgNACs were examined by quantitative reverse transcription-PCR. Many of these genes showed diverse salt- and drought-responsive expression patterns in roots and leaves. These results confirm that PgNACs are potentially involved in regulating abiotic stress tolerance in pearl millet.
The pearl millet genome contains 151 NAC transcription factor genes that can be classified into 11 groups. Many of these genes are either upregulated or downregulated by either salinity or drought stress and may therefore contribute to establishing stress tolerance in pearl millet.
DNA methylation is an epigenetic regulatory form that plays an important role in regulating the gene expression and the tissues development.. However, DNA methylation regulators involved in sheep ...muscle development remain unclear. To explore the functional importance of genome-scale DNA methylation during sheep muscle growth, this study systematically investigated the genome-wide DNA methylation profiles at key stages of Hu sheep developmental (fetus and adult) using deep whole-genome bisulfite sequencing (WGBS).
Our study found that the expression levels of DNA methyltransferase (DNMT)-related genes were lower in fetal muscle than in the muscle of adults. The methylation levels in the CG context were higher than those in the CHG and CHH contexts, and methylation levels were highest in introns, followed by exons and downstream regions. Subsequently, we identified 48,491, 17, and 135 differentially methylated regions (DMRs) in the CG, CHG, and CHH sequence contexts and 11,522 differentially methylated genes (DMGs). The results of bisulfite sequencing PCR (BSP) correlated well with the WGBS-Seq data. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation analysis revealed that some DMGs were involved in regulating skeletal muscle development and fatty acid metabolism. By combining the WGBS-Seq and previous RNA-Seq data, a total of 159 overlap genes were obtained between differentially expressed genes (DEGs) and DMGs (FPKM > 10 and fold change > 4). Finally, we found that 9 DMGs were likely to be involved in muscle growth and metabolism of Hu sheep.
We systemically studied the global DNA methylation patterns of fetal and adult muscle development in Hu sheep, which provided new insights into a better understanding of the epigenetic regulation of sheep muscle development.
Mutations in minor spliceosome components such as U12 snRNA (cerebellar ataxia) and U4atac snRNA (microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1)) result in tissue-specific symptoms. ...Given that the minor spliceosome is ubiquitously expressed, we hypothesized that these restricted phenotypes might be caused by the tissue-specific regulation of the minor spliceosome targets, i.e. minor intron-containing genes (MIGs). The current model of inefficient splicing is thought to apply to the regulation of the ~ 500 MIGs identified in the U12DB. However this database was created more than 10 years ago. Therefore, we first wanted to revisit the classification of minor introns in light of the most recent reference genome. We then sought to address specificity of MIG expression, minor intron retention, and alternative splicing (AS) across mouse and human tissues.
We employed position-weight matrices to obtain a comprehensive updated list of minor introns, consisting of 722 mouse and 770 human minor introns. These can be found in the Minor Intron DataBase (MIDB). Besides identification of 99% of the minor introns found in the U12DB, we also discovered ~ 150 new MIGs. We then analyzed the RNAseq data from eleven different mouse tissues, which revealed tissue-specific MIG expression and minor intron retention. Additionally, many minor introns were efficiently spliced compared to their flanking major introns. Finally, we identified several novel AS events across minor introns in both mouse and human, which were also tissue-dependent. Bioinformatics analysis revealed that several of the AS events could result in the production of novel tissue-specific proteins. Moreover, like the major introns, we found that these AS events were more prevalent in long minor introns, while retention was favoured in shorter introns.
Here we show that minor intron splicing and AS across minor introns is a highly organised process that might be regulated in coordination with the major spliceosome in a tissue-specific manner. We have provided a framework to further study the impact of the minor spliceosome and the regulation of MIG expression. These findings may shed light on the mechanism underlying tissue-specific phenotypes in diseases associated with minor spliceosome inactivation. MIDB can be accessed at https://midb.pnb.uconn.edu .
Plant natural products have been co-opted for millennia by humans for various uses such as flavor, fragrances, and medicines. These compounds often are only produced in relatively low amounts and are ...difficult to chemically synthesize, limiting access. While elucidation of the underlying biosynthetic processes might help alleviate these issues (e.g., via metabolic engineering), investigation of this is hindered by the low levels of relevant gene expression and expansion of the corresponding enzymatic gene families. However, the often-inducible nature of such metabolic processes enables selection of those genes whose expression pattern indicates a role in production of the targeted natural product.
Here, we combine metabolomics and transcriptomics to investigate the inducible biosynthesis of the bioactive diterpenoid tanshinones from the Chinese medicinal herb, Salvia miltiorrhiza (Danshen). Untargeted metabolomics investigation of elicited hairy root cultures indicated that tanshinone production was a dominant component of the metabolic response, increasing at later time points. A transcriptomic approach was applied to not only define a comprehensive transcriptome (comprised of 20,972 non-redundant genes), but also its response to induction, revealing 6,358 genes that exhibited differential expression, with significant enrichment for up-regulation of genes involved in stress, stimulus and immune response processes. Consistent with our metabolomics analysis, there appears to be a slower but more sustained increased in transcript levels of known genes from diterpenoid and, more specifically, tanshinone biosynthesis. Among the co-regulated genes were 70 transcription factors and 8 cytochromes P450, providing targets for future investigation.
Our results indicate a biphasic response of Danshen terpenoid metabolism to elicitation, with early induction of sesqui- and tri- terpenoid biosynthesis, followed by later and more sustained production of the diterpenoid tanshinones. Our data provides a firm foundation for further elucidation of tanshinone and other inducible natural product metabolism in Danshen.