Deciphering mycobiota and its functional dynamics in root hairs of Rhododendron campanulatum D. Don through Next-gen sequencing

The Himalayas provide unique opportunities for the extension of shrubs beyond the upper limit of the tree. However, little is known about the limitation of the biotic factors belowground of shrub growth at these cruising altitudes. To fill this gap, the present study deals with the documentation of root-associated microbiota with their predicted functional profiles and interactions in the host Rhododendron campanulatum, a krummholz species. While processing 12 root samples of R. campanulatum from the sites using Omics we could identify 134 root-associated fungal species belonging to 104 genera, 74 families, 39 orders, 17 classes, and 5 phyla. The root-associated microbiota members of Ascomycota were unambiguously dominant followed by Basidiomycota. Using FUNGuild, we reported that symbiotroph and pathotroph as abundant trophic modes. Furthermore, FUNGuild revealed the dominant prevalence of the saptroptroph guild followed by plant pathogens and wood saprotrophs. Alpha diversity was significantly different at the sites. The heatmap dendrogram showed the correlation between various soil nutrients and some fungal species. The study paves the way for a more in-depth exploration of unidentified root fungal symbionts, their interactions and their probable functional roles, which may serve as an important factor for the growth and conservation of these high-altitude ericaceous plants.

Using team-based precision medicine to advance understanding of rare genetic brain disorders

We describe a multidisciplinary teamwork approach known as "Operation IDD Gene Team" developed by the Rose F. Kennedy Intellectual and Developmental Disabilities Research Center (RFK IDDRC) at the Albert Einstein College of Medicine. This initiative brings families affected by rare genetic diseases that cause intellectual and developmental disability together with physicians, basic scientists, and their trainees. At team meetings, family members share their child's medical and personal history, physicians describe the broader clinical consequences of the condition, and scientists provide accessible tutorials focused on the fundamental biology of relevant genes. When appropriate, possible treatment approaches are also discussed. The outcomes of team meetings have been overwhelmingly positive, with families not only expressing deep gratitude, but also becoming empowered to establish foundations dedicated to their child's specific condition. Physicians, and in particular the scientists and their trainees, have gained a deeper understanding of challenges faced by affected families, broadening their perspective on how their research can extend beyond the laboratory. Remarkably, research by the scientists following the Gene Team meetings have often included focus on the actual gene variants exhibited by the participating children. As these investigations progress and newly created foundations expand their efforts, national as well as international collaborations are forged. These developments emphasize the importance of rare diseases as windows into previously unexplored molecular and cellular processes, which can offer fresh insights into both normal function as well as more common diseases. Elucidating the mechanisms of and treatments for rare and ultra-rare diseases thus has benefits for all involved-families, physicians, and scientists and their trainees, as well as the broader medical community. While the RFK IDDRC's Operation IDD Gene Team program has focused on intellectual disabilities affecting children, we believe it has the potential to be applied to rare genetic diseases impacting individuals of any age and encompassing a wide variety of developmental disorders affecting multiple organ systems.

Tissue-specific expansion of Zika virus isogenic variants drive disease pathogenesis

BACKGROUND

The Asian lineage Zika virus (ZIKV) emerged as a public health emergency in 2016 causing severe neurological pathologies with no apparent historical correlate to the mild, disease-causing innocuous member of the mosquito-borne flavivirus genus that was discovered in Africa in 1947. Replication error rate of RNA viruses combined with viral protein/RNA structural plasticity can lead to evolution of virus-induced pathogenicity that is critical to identify and validate.

METHODS

Infection studies in cells and A129 interferon alpha/beta receptor deficient mice with ZIKV French Polynesian H/PF/2013 clinical isolate, plaque-purified isogenic clone derivatives as well as infectious cDNA clone derived wild-type and site-specific mutant viruses, were employed together with Next-Generation Sequencing (NGS) to pin-point the contributions of specific viral variants in neurovirulence recapitulated in our ZIKV mouse model.

FINDINGS

NGS analysis of the low-passage inoculum virus as well as mouse serum, brain and testis derived virus, revealed specific enrichment in the mouse brain that were not found in the other tissues. Specifically, non-structural (NS) protein 2A variant at position 117 along with changes in NS1 and NS4B were uniquely associated with the mouse brain isolate. Mutational analysis of these variants in cDNA infectious clones identified the NS2A A117V as the lethal pathogenic determinant with potential epistatic contribution of NS1 and NS4B variants in ZIKV brain penetrance.

INTERPRETATION

Our findings confirm that viral subpopulations drive ZIKV neuropathogenicity and identify specific sequence variants that expand in the mouse brain that associates with this phenotype which can serve as predictors of severe epidemics.

FUNDING

Duke-NUS Khoo Post-doctoral Fellowship Award 2020 (KWKC) and National Medical Research Council of Singapore grants MOH-000524 (OFIRG) (SW) and MOH-OFIRG20nov-0002 (SGV).

Next-generation sequencing of cerebrospinal fluid for clinical molecular diagnostics in pediatric, adolescent and young adult brain tumor patients

BACKGROUND

Safe sampling of central nervous system tumor tissue for diagnostic purposes may be difficult if not impossible, especially in pediatric patients, and an unmet need exists to develop less invasive diagnostic tests.

METHODS

We report our clinical experience with minimally invasive molecular diagnostics using a clinically validated assay for sequencing of cerebrospinal fluid (CSF) cell-free DNA (cfDNA). All CSF samples were collected as part of clinical care, and results reported to both clinicians and patients/families.

RESULTS

We analyzed 64 CSF samples from 45 pediatric, adolescent and young adult (AYA) patients (pediatric = 25; AYA = 20) with primary and recurrent brain tumors across 12 histopathological subtypes including high-grade glioma (n = 10), medulloblastoma (n = 10), pineoblastoma (n = 5), low-grade glioma (n = 4), diffuse leptomeningeal glioneuronal tumor (DLGNT) (n = 4), retinoblastoma (n = 4), ependymoma (n = 3), and other (n = 5). Somatic alterations were detected in 30/64 samples (46.9%) and in at least one sample per unique patient in 21/45 patients (46.6%). CSF cfDNA positivity was strongly associated with the presence of disseminated disease at the time of collection (81.5% of samples from patients with disseminated disease were positive). No association was seen between CSF cfDNA positivity and the timing of CSF collection during the patient's disease course.

CONCLUSIONS

We identified three general categories where CSF cfDNA testing provided additional relevant diagnostic, prognostic, and/or therapeutic information, impacting clinical assessment and decision making: (1) diagnosis and/or identification of actionable alterations; (2) monitor response to therapy; and (3) tracking tumor evolution. Our findings support broader implementation of clinical CSF cfDNA testing in this population to improve care.

Detection and molecular characterization of astro and bocaviruses in dogs in Minnesota

Canine astrovirus (CAstV) and canine bocavirus (CBoV) are involved in cases of mild, and sometimes severe, gastroenteritis in dogs. Fecal samples from two dead dogs with gastroenteritis were received at the University of Minnesota Veterinary Diagnostic Laboratory to determine the cause of death. Small round viruses of 20-35 nm diameter were observed by negative contrast electron microscopy. The samples were subjected to Illumina MiSeq sequencing. Both samples were strongly positive for CAstV; all viral reads were related to CAstV. In addition, sample number 1 had a few reads of CBoV. Two complete sequences of CAstV were identified (6625 and 6627 nt in length) with 95% nt identity. RT-PCR and PCR were used to confirm CAstV and CBoV infections in successive samples of canine gastroenteritis. Sanger sequencing was done on nucleic acids from positive samples. Of a total of ten samples, CAstV and CBoV infections were confirmed in six and eight animals, respectively. Four animals had mixed infection with both viruses. All sequences of ORF1b gene of CAstVs showed closest clusters in phylogenetic tree with 96-100% nucleotide and amino acids identity. On the other hand, identity between VP2 gene of different CBoV strains in this study ranged from 93%- 100%. All strains were located close to each other except the divergent MT078234 strain, which was arranged in a separate branch and was closer to reference strain JN648103/USA/2010. This study highlights the importance of electron microscopy and next generation sequencing for early detection and characterization of viruses associated with dog gastroenteritis.

Enteric Viruses Associated with Mid-growth Turkey Enteritis

Since August 2014, the University of Minnesota Veterinary Diagnostic Laboratory has received cases of turkey enteritis that are clinically different from previously described cases of poult enteritis syndrome and light turkey syndrome. The birds develop dark green and extremely foul-smelling diarrhea starting at 8-10 wk of age, which may last up to 15-16 wk of age. The affected turkey flocks show poor uniformity, and feed conversion and market weights are reduced. Multiple-age farms are affected more often than the single-age farms. Morbidity varies from flock to flock and in some cases reaches 100%. At necropsy, undigested feed with increased mucus is observed in the intestines along with prominent mucosal congestion and/or hemorrhage. Microscopically, lymphocytic infiltrates expand the villi in duodenum and jejunum to form lymphoid follicles, which are often accompanied by heterophils. Next generation sequencing (Illumina Miseq) on a pool of feces from affected birds identified genetic sequences of viruses belonging to Astroviridae, Reoviridae, Picornaviridae, Picobirnaviridae, and Adenoviridae. On testing pools of fecal samples from apparently healthy (16 pools) and affected birds (30 pools), there was a higher viral load in the feces of affected birds. Picobirnavirus was detected only in the affected birds; 20 of 30 pools (66.7%) were positive. These results indicate that a high viral load of turkey picobirnavirus alone, or in association with novel picornaviruses, may be a cause of this new type of turkey enteritis.

Genome-wide identification and characterization of InDels and SNPs in Glycine max and Glycine soja for contrasting seed permeability traits

BACKGROUND

Water permeability governed by seed coat is a major facet of seed crops, especially soybean, whose seeds lack physiological dormancy and experience rapid deterioration in seed viability under prolonged storage. Moreover, the physiological and chemical characteristics of soybean seeds are known to vary with seed coat color. Thus, to underpin the genes controlling water permeability in soybean seeds, we carried out an in-depth characterization of the associated genomic variation.

RESULTS

In the present study, we have analyzed genomic variation between cultivated soybean and its wild progenitor with implications on seed permeability, a trait related to seed storability. Whole genome resequencing of G.max and G. soja, identified SNPs and InDels which were further characterized on the basis of their genomic location and impact on gene expression. Chromosomal density distribution of the variation was assessed across the genome and genes carrying SNPs and InDels were characterized into different metabolic pathways. Seed hardiness is a complex trait that is affected by the allelic constitution of a genetic locus as well as by a tricky web of plant hormone interactions. Seven genes that hold a probable role in the determination of seed permeability were selected and their expression differences at different stages of water imbibition were analyzed. Variant interaction network derived 205 downstream interacting partners of 7 genes confirmed their role in seed related traits. Interestingly, genes encoding for Type I- Inositol polyphosphate 5 phosphatase1 and E3 Ubiquitin ligase could differentiate parental genotypes, revealed protein conformational deformations and were found to segregate among RILs in coherence with their permeability scores. The 2 identified genes, thus showed a preliminary association with the desirable permeability characteristics.

CONCLUSION

In the light of above outcomes, 2 genes were identified that revealed preliminary, but a relevant association with soybean seed permeability trait and hence could serve as a primary material for understanding the molecular pathways controlling seed permeability traits in soybean.

Comparison of two bioinformatics tools used to characterize the microbial diversity and predictive functional attributes of microbial mats from Lake Obersee, Antarctica

In this study, using NextGen sequencing of the collective 16S rRNA genes obtained from two sets of samples collected from Lake Obersee, Antarctica, we compared and contrasted two bioinformatics tools, PICRUSt and Tax4Fun. We then developed an R script to assess the taxonomic and predictive functional profiles of the microbial communities within the samples. Taxa such as Pseudoxanthomonas, Planctomycetaceae, Cyanobacteria Subsection III, Nitrosomonadaceae, Leptothrix, and Rhodobacter were exclusively identified by Tax4Fun that uses SILVA database; whereas PICRUSt that uses Greengenes database uniquely identified Pirellulaceae, Gemmatimonadetes A1-B1, Pseudanabaena, Salinibacterium and Sinobacteraceae. Predictive functional profiling of the microbial communities using Tax4Fun and PICRUSt separately revealed common metabolic capabilities, while also showing specific functional IDs not shared between the two approaches. Combining these functional predictions using a customized R script revealed a more inclusive metabolic profile, such as hydrolases, oxidoreductases, transferases; enzymes involved in carbohydrate and amino acid metabolisms; and membrane transport proteins known for nutrient uptake from the surrounding environment. Our results present the first molecular-phylogenetic characterization and predictive functional profiles of the microbial mat communities in Lake Obersee, while demonstrating the efficacy of combining both the taxonomic assignment information and functional IDs using the R script created in this study for a more streamlined evaluation of predictive functional profiles of microbial communities.

Metagenomics approach to the study of the gut microbiome structure and function in zebrafish Danio rerio fed with gluten formulated diet

In this study, we report the gut microbial composition and predictive functional profiles of zebrafish, Danio rerio, fed with a control formulated diet (CFD), and a gluten formulated diet (GFD) using a metagenomics approach and bioinformatics tools. The microbial communities of the GFD-fed D. rerio displayed heightened abundances of Legionellales, Rhizobiaceae, and Rhodobacter, as compared to the CFD-fed counterparts. Predicted metagenomics of microbial communities (PICRUSt) in GFD-fed D. rerio showed KEGG functional categories corresponding to bile secretion, secondary bile acid biosynthesis, and the metabolism of glycine, serine, and threonine. The CFD-fed D. rerio exhibited KEGG functional categories of bacteria-mediated cobalamin biosynthesis, which was supported by the presence of cobalamin synthesizers such as Bacteroides and Lactobacillus. Though these bacteria were absent in GFD-fed D. rerio, a comparable level of the cobalamin biosynthesis KEGG functional category was observed, which could be contributed by the compensatory enrichment of Cetobacterium. Based on these results, we conclude D. rerio to be a suitable alternative animal model for the use of a targeted metagenomics approach along with bioinformatics tools to further investigate the relationship between the gluten diet and microbiome profile in the gut ecosystem leading to gastrointestinal diseases and other undesired adverse health effects.

Adeno-associated virus gene therapy vector scAAVIGF-I for transduction of equine articular chondrocytes and RNA-seq analysis

OBJECTIVE

IGF-I is one of several anabolic factors being investigated for the treatment of osteoarthritis (OA). Due to the short biological half-life, extended administration is required for more robust cartilage healing. Here we create a self-complimentary adeno-associated virus (AAV) gene therapy vector utilizing the transgene for IGF-I.

DESIGN

Various biochemical assays were performed to investigate the cellular response to scAAVIGF-I treatment vs an scAAVGFP positive transduction control and a negative for transduction control culture. RNA-sequencing analysis was also performed to establish a differential regulation profile of scAAVIGF-I transduced chondrocytes.

RESULTS

Biochemical analyses indicated an average media IGF-I concentration of 608 ng/ml in the scAAVIGF-I transduced chondrocytes. This increase in IGF-I led to increased expression of collagen type II and aggrecan and increased protein concentrations of cellular collagen type II and media glycosaminoglycan vs both controls. RNA-seq revealed a global regulatory pattern consisting of 113 differentially regulated GO categories including those for chondrocyte and cartilage development and regulation of apoptosis.

CONCLUSIONS

This research substantiates that scAAVIGF-I gene therapy vector increased production of IGF-I to clinically relevant levels with a biological response by chondrocytes conducive to increased cartilage healing. The RNA-seq further established a set of differentially expressed genes and gene ontologies induced by the scAAVIGF-I vector while controlling for AAV infection. This dataset provides a static representation of the cellular transcriptome that, while only consisting of one time point, will allow for further gene expression analyses to compare additional cartilage healing therapeutics or a transient cellular response.

Proteogenomics for understanding oncology: recent advances and future prospects

The concept of proteogenomics has emerged rapidly as a valuable approach to integrate mass spectrometry-derived proteomic data with genomic and transcriptomic data. It is used to harness the full potential of the former dataset in the discovery of potential biomarkers, therapeutic targets and novel proteins associated with various biological processes including diseases. Proteogenomic strategies have been successfully utilized to identify novel genes and redefine annotation of existing gene models in various genomes. In recent years, this approach has been extended to the field of cancer biology to unravel complexities in the tumor genomes and proteomes. Standard proteomics workflows employing translated cancer genomes and transcriptomes can potentially identify peptides from mutant proteins, splice variants and fusion proteins in the tumor proteome, which in addition to the currently available biomarker panels can serve as potential diagnostic and prognostic biomarkers, besides having therapeutic utility. This review focuses on the role of proteogenomics to understand cancer biology.

Hurt, tired and queasy: Specific variants in the ATPase domain of the TRAP1 mitochondrial chaperone are associated with common, chronic "functional" symptomatology including pain, fatigue and gastrointestinal dysmotility

Functional disorders are common conditions with a substantial impact on a patients' wellbeing, and can be diagnostically elusive. There are bidirectional associations between functional disorders and mitochondrial dysfunction. In this study, provided clinical information and the exon sequence of the TRAP1 mitochondrial chaperone were retrospectively reviewed with a focus on the functional categories of chronic pain, fatigue and gastrointestinal dysmotility. Very-highly conserved TRAP1 variants were identified in 73 of 930 unrelated patients. Functional symptomatology is strongly associated with specific variants in the ATPase binding pocket. In particular, the combined presence of all three functional categories is strongly associated with p.Ile253Val (OR 7.5, P = 0.0001) and with two other interacting variants (OR 18, P = 0.0005). Considering a 1-2% combined variant prevalence and high odds ratios, these variants may be an important factor in the etiology of functional symptomatology.

Intersection of toxicogenomics and high throughput screening in the Tox21 program: an NIEHS perspective

Humans are exposed to thousands of chemicals with inadequate toxicological data. Advances in computational toxicology, robotic high throughput screening (HTS), and genome-wide expression have been integrated into the Tox21 program to better predict the toxicological effects of chemicals. Tox21 is a collaboration among US government agencies initiated in 2008 that aims to shift chemical hazard assessment from traditional animal toxicology to target-specific, mechanism-based, biological observations using in vitro assays and lower organism models. HTS uses biocomputational methods for probing thousands of chemicals in in vitro assays for gene-pathway response patterns predictive of adverse human health outcomes. In 1999, NIEHS began exploring the application of toxicogenomics to toxicology and recent advances in NextGen sequencing should greatly enhance the biological content obtained from HTS platforms. We foresee an intersection of new technologies in toxicogenomics and HTS as an innovative development in Tox21. Tox21 goals, priorities, progress, and challenges will be reviewed.

Targeted genome enrichment for efficient purification of endosymbiont DNA from host DNA

Wolbachia endosymbionts are widespread in arthropods and are generally considered reproductive parasites, inducing various phenotypes including cytoplasmic incompatibility, parthenogenesis, feminization and male killing, which serve to promote their spread through populations. In contrast, Wolbachia infecting filarial nematodes that cause human diseases, including elephantiasis and river blindness, are obligate mutualists. DNA purification methods for efficient genomic sequencing of these unculturable bacteria have proven difficult using a variety of techniques. To efficiently capture endosymbiont DNA for studies that examine the biology of symbiosis, we devised a parallel strategy to an earlier array-based method by creating a set of SureSelect™ (Agilent) 120-mer target enrichment RNA oligonucleotides ("baits") for solution hybrid selection. These were designed from Wolbachia complete and partial genome sequences in GenBank and were tiled across each genomic sequence with 60 bp overlap. Baits were filtered for homology against host genomes containing Wolbachia using BLAT and sequences with significant host homology were removed from the bait pool. Filarial parasite Brugia malayi DNA was used as a test case, as the complete sequence of both Wolbachia and its host are known. DNA eluted from capture was size selected and sequencing samples were prepared using the NEBNext® Sample Preparation Kit. One-third of a 50 nt paired-end sequencing lane on the HiSeq™ 2000 (Illumina) yielded 53 million reads and the entirety of the Wolbachia genome was captured. We then used the baits to isolate more than 97.1 % of the genome of a distantly related Wolbachia strain from the crustacean Armadillidium vulgare, demonstrating that the method can be used to enrich target DNA from unculturable microbes over large evolutionary distances.