Search

The impacts of invertebrate RNA virus population dynamics on virulence and infection outcomes are poorly understood. Deformed wing virus (DWV), the main viral pathogen of honey bees, negatively impacts bee health, which can lead to colony death. Despite previous reports on the reduction of DWV diversity following the arrival of the parasitic mite Varroa destructor, the key DWV vector, we found high genetic diversity of DWV in infested United States honey bee colonies. Phylogenetic analysis showed that divergent US DWV genotypes are of monophyletic origin and were likely generated as a result of diversification after a genetic bottleneck. To investigate the population dynamics of this divergent DWV, we designed a series of novel infectious cDNA clones corresponding to coexisting DWV genotypes, thereby devising a reverse-genetics system for an invertebrate RNA virus quasispecies. Equal replication rates were observed for all clone-derived DWV variants in single infections. Surprisingly, individual clones replicated to the same high levels as their mixtures and even the parental highly diverse natural DWV population, suggesting that complementation between genotypes was not required to replicate to high levels. Mixed clone–derived infections showed a lack of strong competitive exclusion, suggesting that the DWV genotypes were adapted to coexist. Mutational and recombination events were observed across clone progeny, providing new insights into the forces that drive and constrain virus diversification. Accordingly, our results suggest that Varroa influences DWV dynamics by causing an initial selective sweep, which is followed by virus diversification fueled by negative frequency-dependent selection for new genotypes. We suggest that this selection might reflect the ability of rare lineages to evade host defenses, specifically antiviral RNA interference (RNAi). In support of this hypothesis, we show that RNAi induced against one DWV strain is less effective against an alternate strain from the same population. [ABSTRACT FROM AUTHOR]

Genome-wide association studies (GWAS) have commonly been used to identify candidate genes that control economically important traits in livestock. Our objective was to detect potential candidate genes associated mainly with muscle development traits related to dimension of hindquarter in cattle. A next generation sequencing (NGS) dataset to imputed to 12 million single nucleotide polymorphisms (SNPs) (from 1252 Simmental beef cattle) were used to search for genes affecting hindquarter traits using a linear, mixed model approach. We also used haplotype and linkage disequilibrium blocks to further support our identifications. We identified 202 significant SNPs in the bovine BTA4 chromosome region associated with width of hind leg, based on a stringent statistical threshold (p = 0.05/ effective number of SNPs identified). After exploring the region around these SNPs, we found candidate genes that were potentially related to the associated markers. More importantly, we identified a region of approximately 280 Kb on the BTA4 chromosome that harbored several muscle specific candidate genes, genes to be in a potential region for muscle development. However, we also found candidate gene SLC13A1 on BTA4, which seems to be associated with bone disorders (such as chondrodysplasia) in Simmental beef cattle. [ABSTRACT FROM AUTHOR]

Limited translational genomic research data have been reported on the application of exome sequencing and parallel gene testing for preconception carrier screening (PCS). Here, we present individual-level data from a large PCS program in which exome sequencing was routinely performed on either gamete donors (5,845) or infertile patients (8,280) undergoing in vitro fertilization (IVF) treatment without any known family history of inheritable genetic conditions. Individual-level data on pathogenic variants were used to define conditions for PCS based on criteria for severity, penetrance, inheritance pattern, and age of onset. Fetal risk was defined based on actual carrier frequency data accounting for the specific inheritance pattern (fetal disease risk, FDR). In addition, large-scale application of exome sequencing for PCS allowed a deep investigation of the incidence of medically actionable secondary findings in this population. Exome sequencing achieved remarkable clinical sensitivity for reproductive risk of highly penetrant childhood-onset disorders (1/337 conceptions) through analysis of 114 selected gene-condition pairs. A significant contribution to fetal disease risk was observed for rare (carrier rate < 1:100) and X-linked conditions (16.7% and 41.2% of total FDR, respectively). Subgroup analysis of 776 IVF couples identified 37 at increased reproductive risk (4.8%; 95% CI = 3.4–6.5). Further, two additional couples had increased risk for very rare conditions when both members of a parental pair were treated as a unit and the search was extended to the entire exome. About 2.3% of participants showed at least one pathogenic variant for genes included in the updated American College of Medical Genetics and Genomics v2.0 list of secondary findings. Gamete donors and IVF couples showed similar carrier burden for both carrier screening and secondary findings, indicating no causal relationship to fertility. These translational research data will facilitate development of more effective PCS strategies that maximize clinical sensitivity with minimal counterproductive effects. [ABSTRACT FROM AUTHOR]

MicroRNAs are conserved, endogenous small RNAs with critical post-transcriptional regulatory functions throughout eukaryota, including prominent roles in development and disease. Despite much effort, microRNA annotations still contain errors and are incomplete due especially to challenges related to identifying valid miRs that have small numbers of reads, to properly locating hairpin precursors and to balancing precision and recall. Here, we present miRWoods, which solves these challenges using a duplex-focused precursor detection method and stacked random forests with specialized layers to detect mature and precursor microRNAs, and has been tuned to optimize the harmonic mean of precision and recall. We trained and tuned our discovery pipeline on data sets from the well-annotated human genome, and evaluated its performance on data from mouse. Compared to existing approaches, miRWoods better identifies precursor spans, and can balance sensitivity and specificity for an overall greater prediction accuracy, recalling an average of 10% more annotated microRNAs, and correctly predicts substantially more microRNAs with only one read. We apply this method to the under-annotated genomes of Felis catus (domestic cat) and Bos taurus (cow). We identified hundreds of novel microRNAs in small RNA sequencing data sets from muscle and skin from cat, from 10 tissues from cow and also from human and mouse cells. Our novel predictions include a microRNA in an intron of tyrosine kinase 2 (TYK2) that is present in both cat and cow, as well as a family of mirtrons with two instances in the human genome. Our predictions support a more expanded miR-2284 family in the bovine genome, a larger mir-548 family in the human genome, and a larger let-7 family in the feline genome. [ABSTRACT FROM AUTHOR]

Two novel and related pathogenic variants of the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene were structurally and functionally characterized. These alterations have not been previously described in literature. Two patients with diagnosis of Cystic Fibrosis (CF) based on the presence of one mutated allele, p.Phe508del, pathological sweat test and clinical symptoms were studied. To complete the genotypes of both patients, an extensive genetic and functional analysis of the CFTR gene was performed. Extensive genetic characterization confirmed the presence of p.Phe508del pathogenic variant and revealed, in both patients, the presence of an insertion of part of intron 10 in intron 9 of the CFTR gene, within the (TG)m repeat, with a variable poly-T stretch. The molecular lesions resulted to be very similar in both patients, with only a difference in the number of T in the poly-T stretch. The functional characterization at RNA level revealed a complete anomalous splicing, without exon 10, from the allele with the insertion of both patients. Consequently, the alleles with the insertions are expected not to contribute to the formation of a functional CFTR protein. Molecular and functional features of these alterations are compatible with the definition of novel CF-causing variants of the CFTR gene. This also allowed the completion of the genetic characterization of both patients. [ABSTRACT FROM AUTHOR]

Lutzomyia longipalpis is the main vector of Leishmania infantum, the causative agent of visceral leishmaniasis in the Neotropical realm. Its taxonomic status has been widely discussed once it encompasses a complex of species. The knowledge about the genetic structure of insect vector populations helps the elucidation of components and interactions of the disease ecoepidemiology. Thus, the objective of this study was to genotypically analyze populations of the Lu. longipalpis complex from a macrogeographic perspective using Next Generation Sequencing. Polymorphism analysis of three molecular markers was used to access the levels of population genetic structure among nine different populations of sand flies. Illumina Amplicon Sequencing Protocol® was used to identify possible polymorphic sites. The library was sequenced on paired-end Illumina MiSeq platform. Significant macrogeographical population differentiation was observed among Lu. longipalpis populations via PCA and DAPC analyses. Our results revealed that populations of Lu. longipalpis from the nine municipalities were grouped into three clusters. In addition, it was observed that the levels of Lu. longipalpis population structure could be associated with distance isolation. This new sequencing method allowed us to study different molecular markers after a single sequencing run, and to evaluate population and inter-species differences on a macrogeographic scale. [ABSTRACT FROM AUTHOR]

Alternative splicing diversifies mRNA transcripts in human cells. While the spliceosome pairs exons with a high degree of accuracy, the rates of rare aberrant and non-canonical pre-mRNA splicing have not been evaluated at the nucleotide level to determine the quantity and identity of these events across splice junctions. Using ultra-deep sequencing the frequency of aberrant and non-canonical splicing events for three splice junctions flanking exon 7 of SMN1 were determined at single nucleotide resolution. After correction for background noise introduced by PCR amplification and sequencing steps, pre-mRNA splicing was shown to maintain a low overall rate of aberrant and non-canonically spliced events. Several previously unannotated splicing events across 3 exon|intron junctions in SMN1 were identified. Mutations within SMN exon 7 were shown to affect splicing fidelity by modulating RNA secondary structures, by altering the binding site of regulatory proteins and by changing the 5’ splice site strength. Mutations also create a truncated SMN1 exon 7 through the introduction of a de novo non-canonical 5’ splice site. The results from the ultra-deep sequencing approach highlight the impressive fidelity of pre-mRNA splicing and demonstrate that the immediate sequence context around splice sites is the main driving force behind non-canonical splice site pairing. [ABSTRACT FROM AUTHOR]

Characterization of Human Endogenous Retrovirus (HERV) expression within the transcriptomic landscape using RNA-seq is complicated by uncertainty in fragment assignment because of sequence similarity. We present Telescope, a computational software tool that provides accurate estimation of transposable element expression (retrotranscriptome) resolved to specific genomic locations. Telescope directly addresses uncertainty in fragment assignment by reassigning ambiguously mapped fragments to the most probable source transcript as determined within a Bayesian statistical model. We demonstrate the utility of our approach through single locus analysis of HERV expression in 13 ENCODE cell types. When examined at this resolution, we find that the magnitude and breadth of the retrotranscriptome can be vastly different among cell types. Furthermore, our approach is robust to differences in sequencing technology, and demonstrates that the retrotranscriptome has potential to be used for cell type identification. We compared our tool with other approaches for quantifying TE expression, and found that Telescope has the greatest resolution, as it estimates expression at specific TE insertions rather than at the TE subfamily level. Telescope performs highly accurate quantification of the retrotranscriptomic landscape in RNA-seq experiments, revealing a differential complexity in the transposable element biology of complex systems not previously observed. Telescope is available at . [ABSTRACT FROM AUTHOR]

Telomeres are nucleo-protein structures that protect the ends of eukaryotic chromosomes. They are not completely synthesized during DNA replication and are elongated by specific mechanisms. The structure of the telomeres and the elongation mechanism have not been determined in Dictyostelium discoideum. This organism presents extrachromosomal palindromic elements containing two copies of the rDNA, also present at the end of the chromosomes. In this article the structure of the terminal region of the rDNA is shown to consist of repetitions of the A(G)n sequence where the number of Gs is variable. These repeats extend as a 3’ single stranded region. The G-rich region is preceded by four tandem repetitions of two different DNA motifs. D. discoideum telomere reverse transcriptase homologous protein (TERTHP) presented RNase-sensitive enzymatic activity and was required to maintain telomere structure since terthp-mutant strains presented reorganizations of the DNA terminal regions. These modifications were different in several terthp-mutants and changed with their prolonged culture and subcloning. However, the terthp gene is not essential for D. discoideum proliferation. Telomeres could be maintained in terthp-mutant strains by homologous recombination mechanisms such as ALT (Alternative Lengthening of Telomeres) or HAATI (heterochromatin amplification-mediated and telomerase-independent). In agreement with this hypothesis, the expression of mRNAs coding for several proteins involved in homologous recombination was induced in terthp-mutant strains. Extrachromosomal rDNA could serve as substrate in these DNA homologous recombination reactions. [ABSTRACT FROM AUTHOR]

HIV-1 replicates via a low-fidelity polymerase with a high mutation rate; strong conservation of individual nucleotides is highly indicative of the presence of critical structural or functional properties. Identifying such conservation can reveal novel insights into viral behaviour. We analysed 3651 publicly available sequences for the presence of nucleic acid conservation beyond that required by amino acid constraints, using a novel scale-free method that identifies regions of outlying score together with a codon scoring algorithm. Sequences with outlying score were further analysed using an algorithm for producing local RNA folds whilst accounting for alignment properties. 11 different conserved regions were identified, some corresponding to well-known cis-acting functions of the HIV-1 genome but also others whose conservation has not previously been noted. We identify rational causes for many of these, including cis functions, possible additional reading frame usage, a plausible mechanism by which the central polypurine tract primes second-strand DNA synthesis and a conformational stabilising function of a region at the 5′ end of env. [ABSTRACT FROM AUTHOR]

The wealth of information deliverable from transcriptome sequencing (RNA-seq) is significant, however current applications for variant detection still remain a challenge due to the complexity of the transcriptome. Given the ability of RNA-seq to reveal active regions of the genome, detection of RNA-seq SNPs can prove valuable in understanding the phenotypic diversity between populations. Thus, we present a novel computational workflow named VAP (Variant Analysis Pipeline) that takes advantage of multiple RNA-seq splice aware aligners to call SNPs in non-human models using RNA-seq data only. We applied VAP to RNA-seq from a highly inbred chicken line and achieved high accuracy when compared with the matching whole genome sequencing (WGS) data. Over 65% of WGS coding variants were identified from RNA-seq. Further, our results discovered SNPs resulting from post transcriptional modifications, such as RNA editing, which may reveal potentially functional variation that would have otherwise been missed in genomic data. Even with the limitation in detecting variants in expressed regions only, our method proves to be a reliable alternative for SNP identification using RNA-seq data. The source code and user manuals are available at . [ABSTRACT FROM AUTHOR]

The diagnosis of Marfan spectrum includes a large number of clinical criteria. Although the identification of pathogenic variants contributes to the diagnostic process, its value to the prediction of clinical outcomes is still limited. An important novelty of the present study is represented by the statistical approach adopted to investigate genotype-phenotype correlation. The analysis has been improved considering the extended genetic information obtained by Next Generation Sequencing (NGS) and combining the effects of both rare and common genetic variants in an inclusive model. To this aim a cohort of 181 patients were analyzed with a NGS panel including 11 genes associated with Marfan spectrum. The genotype-phenotype correlation was also investigated considering the possibility to predict presence of a pathological mutation in Marfan syndrome (MFS) main genes based only on the analysis of phenotypic traits. Results obtained indicate that information about clinical traits can be summarized in a new variable that resulted significantly associated with the probability to find a pathological mutation in MFS main genes. This is important since the choice of the genetic test is often influenced by the phenotypic characterization of patients. Moreover, both rare and common variants were found to significantly contribute to clinical spectrum and their combination allowed to increase the percentage of phenotype variability that could be explained based on genetic factors. Results highlight the opportunity to take advantage of the overall genetic information obtained by NGS data to have a better clinical classification of patients. [ABSTRACT FROM AUTHOR]

Manure from animals that have been treated with antibiotics is often used to fertilize agricultural soils and its application has previously been shown to enrich for genes associated with antibiotic resistance in agroecosystems. To investigate the magnitude of this effect, we designed a column experiment simulating manure-treated agricultural soil that utilizes artificial subsurface drainage to determine the duration and extent which this type of manure fertilization impacts the set of genes associated with antibiotic resistance in drainage water. We classified ARGs in manure-treated drainage effluent water by its source of origin. Overall, we found that 61% and 7% of the total abundance of ARGs found in drainage water samples could be attributed to manure enrichment and manure addition, respectively. Among these ARGs, we identified 75 genes unique to manure that persisted in both soil and drainage water throughout a drainage season typical of the Upper Midwestern United States. While most of these genes gradually decreased in abundance over time, the IS6100-associated tet(33) gene accrued. These results demonstrate the influence of manure applications on the composition of the resistome observed in agricultural drainage water and highlight the importance of anthropogenic ARGs in the environment. [ABSTRACT FROM AUTHOR]

The human gut microbiota is well established as an important factor in health and disease. Fecal sample microbiota are often analyzed as a proxy for gut microbiota, and characterized with respect to their composition profiles. Modern approaches employ whole genome shotgun next-generation sequencing as the basis for these analyses. Sequencing depth as well as choice of next-generation sequencing data analysis method constitute two main interacting methodological factors for such an approach. In this study, we used 200 million sequence read pairs from one fecal sample for comparing different taxonomy classification methods, using default and custom-made reference databases, at different sequencing depths. A mock community data set with known composition was used for validating the classification methods. Results suggest that sequencing beyond 60 million read pairs does not seem to improve classification. The phylogeny prediction pattern, when using the default databases and the consensus database, appeared to be similar for all three methods. Moreover, these methods predicted rather different species. We conclude that the choice of sequencing depth and classification method has important implications for taxonomy composition prediction. A multi-method-consensus approach for robust gut microbiota NGS analysis is recommended. [ABSTRACT FROM AUTHOR]

There are a few well-known indigenous breeds of Chinese rabbits in Sichuan and Fujian provinces, for which the genetic diversity and population structure have been poorly investigated. In the present study, we successfully employed the restriction-site-associated DNA sequencing (RAD-seq) approach to comprehensively discover genome-wide SNPs of 104 rabbits from four Chinese indigenous breeds: 30 Sichuan White, 34 Tianfu Black, 32 Fujian Yellow and eight Fujian Black. A total of 7,055,440 SNPs were initially obtained, from which 113,973 high-confidence SNPs (read depth ≥ 3, calling rate = 100% and biallelic SNPs) were selected to study the genetic diversity and population structure. The mean polymorphism information content (PIC) and nucleotide diversity (π) of each breed slightly varied with ranging from 0.2000 to 0.2281 and from 0.2678 to 0.2902, respectively. On the whole, Fujian Yellow rabbits showed the highest genetic diversity, which was followed by Tianfu Black and Sichuan White rabbits. The principal component analysis (PCA) revealed that the four breeds were clearly distinguishable. Our results first reveal the genetic differences among these four rabbit breeds in the Sichuan and Fujian provinces and also provide a high-confidence set of genome-wide SNPs for Chinese indigenous rabbits that could be employed for gene linkage and association analyses in the future. [ABSTRACT FROM AUTHOR]

In this study, we describe the development of a PCR-free whole exome sequencing method. Using this method, 2 μg DNA was sufficient for library preparation for whole exome sequencing. Furthermore, the method is simple and makes use of a commercial kit, with additional step of concentrating the captured library by ethanol precipitation. The accuracy of the PCR-free method was found to be equivalent to that of unique molecular identifier-corrected analysis method, which is the commonly used method to detect rare mutations. Thus, the PCR-free whole exome sequencing method is cost-effective as well as efficient in detecting rare mutations. [ABSTRACT FROM AUTHOR]

Plasmodium vivax invasion of human erythrocytes depends on the Duffy Binding Protein (PvDBP) which interacts with the Duffy antigen. PvDBP copy number has been recently shown to vary between P. vivax isolates in Sub-Saharan Africa. However, the extent of PvDBP copy number variation, the type of PvDBP multiplications, as well as its significance across broad samples are still unclear. We determined the prevalence and type of PvDBP duplications, as well as PvDBP copy number variation among 178 Ethiopian P. vivax isolates using a PCR-based diagnostic method, a novel quantitative real-time PCR assay and whole genome sequencing. For the 145 symptomatic samples, PvDBP duplications were detected in 95 isolates, of which 81 had the Cambodian and 14 Malagasy-type PvDBP duplications. PvDBP varied from 1 to >4 copies. Isolates with multiple PvDBP copies were found to be higher in symptomatic than asymptomatic infections. For the 33 asymptomatic samples, PvDBP was detected with two copies in two of the isolates, and both were the Cambodian-type PvDBP duplication. PvDBP copy number in Duffy-negative heterozygotes was not significantly different from that in Duffy-positives, providing no support for the hypothesis that increased copy number is a specific association with Duffy-negativity, although the number of Duffy-negatives was small and further sampling is required to test this association thoroughly. [ABSTRACT FROM AUTHOR]

Background: Tyrosine kinase inhibitors targeted to anaplastic lymphoma kinase (ALK) have been demonstrated to be effective for lung cancer patients with an ALK fusion gene. Application of liquid biopsy, i.e., detection and quantitation of the fusion product in plasma cell-free DNA (cfDNA), could improve clinical practice. To detect ALK fusions, because fusion breakpoints occur somewhere in intron 19 of the ALK gene, sequencing of the entire intron is required to locate breakpoints. Results: We constructed a target sequencing system using an adapter and a set of primers that cover the entire ALK intron 19. This system can amplify fragments, including breakpoints, regardless of fusion partners. The data analysis pipeline firstly detected fusions by alignment to selected target sequences, and then quantitated the fusion alleles aligning to the identified breakpoint sequences. Performance was validated using 20 cfDNA samples from ALK-positive non-small cell lung cancer patients and samples from 10 healthy volunteers. Sensitivity and specificity were 50 and 100%, respectively. Conclusions: We demonstrated that PCR-based target sequencing using a tiling primer set and two-step mapping/alignment quantitatively detected ALK fusions in cfDNA from lung cancer patients. The system offers an alternative to existing approaches based on hybridization capture. [ABSTRACT FROM AUTHOR]

Background: More than 11,000 laboratories and companies developed their own next-generation sequencing (NGS) for screening and diagnosis of various diseases including cancer. Although inconsistencies of mutation calls as high as 43% in databases such as GDSC (Genomics of Drug Sensitivity in Cancer) and CCLE (Cancer Cell Line Encyclopedia) have been reported, not many studies on the reasons for the inconsistencies have been published. Methods: Targeted-NGS analysis of 151 genes in 35 cell lines common to GDSC and CCLE was performed, and the results were compared with those from GDSC and CCLE wherein whole-exome- or highly-multiplex NGS were employed. Results: In the comparison, GDSC and CCLE had a high rate (40–45%) of false-negative (FN) errors which would lead to high rate of inconsistent mutation calls, suggesting that highly-multiplex NGS may have high rate of FN errors. We also posited the possibility that targeted NGS, especially for the detection of low-level cancer cells in cancer tissues might suffer significant FN errors. Conclusion: FN errors may be the most important errors in NGS testing for cancer; their evaluation in laboratory-developed NGS tests is needed. [ABSTRACT FROM AUTHOR]

We investigated the feasibility of using next-generation sequencing (NGS) technique using molecular barcoding technology to detect MYD88 L265P mutation in unselected peripheral blood mononuclear cells (PBMCs) in 52 patients with Waldenström’s macroglobulinemia [1] and 21 patients with IgM-monoclonal gammopathy of undetermined significance (MGUS). The NGS technique successfully detected the MYD88 L265P in unselected PBMCs at a sensitivity of 0.02%, which was ×5 higher than that of AS-PCR. All the results between paired BM and PB samples from 2 IgM MGUS and 4 untreated WM patients matched completely. MYD88 L265P mutation was detected in 14/21 (66.7%), 14/19 (73.7%), and 10/33 (30.3%) with the median mutant allele burden of 0.36% (range, 0.06–2.85%), 0.48% (range, 0.02–32.3%), and 0.16% (range, 0.02–33.8%), in IgM-MGUS, untreated WM, and previously treated WM, respectively. Multiple linear regression analysis identified an absolute peripheral lymphocyte count as the positive predictor of PB mutant allele burden (R2 = 0,72, P<0.0001). Our non-invasive, simple NGS method has the potential to detect MYD88 L265P mutations in PBMCs of IgM MGUS and WM patients, which may especially utilized for monitoring minimal residual tumor burden after treatment. [ABSTRACT FROM AUTHOR]

Background: Genomic data have become major resources to understand complex mechanisms at fine-scale temporal and spatial resolution in functional and evolutionary genetic studies, including human diseases, such as cancers. Recently, a large number of whole genomes of evolving populations of yeast (Saccharomyces cerevisiae W303 strain) were sequenced in a time-dependent manner to identify temporal evolutionary patterns. For this type of study, a chromosome-level sequence assembly of the strain or population at time zero is required to compare with the genomes derived later. However, there is no fully automated computational approach in experimental evolution studies to establish the chromosome-level genome assembly using unique features of sequencing data. Methods and results: In this study, we developed a new software pipeline, the integrative meta-assembly pipeline (IMAP), to build chromosome-level genome sequence assemblies by generating and combining multiple initial assemblies using three de novo assemblers from short-read sequencing data. We significantly improved the continuity and accuracy of the genome assembly using a large collection of sequencing data and hybrid assembly approaches. We validated our pipeline by generating chromosome-level assemblies of yeast strains W303 and SK1, and compared our results with assemblies built using long-read sequencing and various assembly evaluation metrics. We also constructed chromosome-level sequence assemblies of S. cerevisiae strain Sigma1278b, and three commonly used fungal strains: Aspergillus nidulans A713, Neurospora crassa 73, and Thielavia terrestris CBS 492.74, for which long-read sequencing data are not yet available. Finally, we examined the effect of IMAP parameters, such as reference and resolution, on the quality of the final assembly of the yeast strains W303 and SK1. Conclusions: We developed a cost-effective pipeline to generate chromosome-level sequence assemblies using only short-read sequencing data. Our pipeline combines the strengths of reference-guided and meta-assembly approaches. Our pipeline is available online at including a Docker image, as well as a Perl script, to help users install the IMAP package, including several prerequisite programs. Users can use IMAP to easily build the chromosome-level assembly for the genome of their interest. [ABSTRACT FROM AUTHOR]

The anticariogenic activity of the extract of Glycyrrhiza uralensis (licorice) has been well documented. We recently developed an herbal lollipop containing licorice extracts with Glycyrrhizol A, the compound displaying strong antimicrobial activity against Streptococcus mutans. Preliminary testing showed that the herbal lollipop reduced salivary S. mutans counts in vivo. In this study, we aimed to further test the efficacy of this herbal lollipop for reducing salivary S. mutans levels, and investigate its impact on salivary microbiome. Using a well-established in vitro oral microbiome model, we showed that licorice extract displays targeted killing against S. mutans without affecting the biodiversity of the community. In vivo study corroborated in vitro findings, showing for high caries-risk children aged 3–6 with salivary S. mutans levels >5x105 cells/ml, daily use of 2 licorice-containing lollipops for 3 weeks significantly reduced salivary S. mutans levels compared to the control group. Salivary microbiome analysis showed either no change or even increase in phylogenetic diversity of the oral community following herbal lollipop usage. Although further study with longer term observation is needed, these results suggest that use of licorice extract-containing lollipops can be as a simple and effective way to reduce the risk of dental caries in children. [ABSTRACT FROM AUTHOR]

The analyses of large volumes of metagenomic data extracted from aggregate populations of microscopic organisms residing on and in the human body are advancing contemporary understandings of the integrated participation of microbes in human health and disease. Next generation sequencing technology facilitates said analyses in terms of diversity, community composition, and differential abundance by filtering and binning microbial 16S rRNA genes extracted from human tissues into operational taxonomic units. However, current statistical tools restrict study designs to investigations of limited numbers of host characteristics mediated by limited numbers of samples potentially yielding a loss of relevant information. This paper presents a Bayesian hierarchical negative binomial model as an efficient technique capable of compensating for multivariable sets including tens or hundreds of host characteristics as covariates further expanding analyses of human microbiome count data. Simulation studies reveal that the Bayesian hierarchical negative binomial model provides a desirable strategy by often outperforming three competing negative binomial model in terms of type I error while simultaneously maintaining consistent power. An application of the Bayesian hierarchical negative binomial model using subsets of the open data published by the American Gut Project demonstrates an ability to identify operational taxonomic units significantly differentiable among persons diagnosed by a medical professional with either inflammatory bowel disease or irritable bowel syndrome that are consistent with contemporary gastrointestinal literature. [ABSTRACT FROM AUTHOR]

Angelica dahurica is a widely grown plant species with multiple uses, especially in the medical field. However, the frequent introduction of A. dahurica to new areas has made it difficult to distinguish between varieties. Simple sequence repeats (SSRs) detected based on transcriptome analyses are very useful for constructing genetic maps and analyzing genetic diversity. They are also relevant for the molecular marker-assisted breeding of A. dahurica. We identified 33,724 genic SSR loci based on transcriptome sequencing data. A total of 114 primer pairs were designed for the SSR loci and were tested for their specificity and diversity. Ten SSR loci in untranslated regions were ultimately selected. Subsequently, 56 A. dahurica ecotypes collected from different regions were analyzed. The SSR loci comprised 2–8 alleles, with a mean of 5.2 alleles per locus. The polymorphic information content value and Shannon’s information index were 0.6274–0.2702 (average of 0.4091) and 1.3040–0.5618 (average of 0.8475), respectively. Thus, the 10 novel SSRs identified in this study were almost in accordance with Harvey-Weinberg equilibrium and will be useful for analyzing A. dahurica genetic relationships. The results of this study confirm the potential value of transcriptome databases for the development of new SSR markers. [ABSTRACT FROM AUTHOR]

Babesiosis is considered an emerging disease because its incidence has significantly increased in the last 30 years, providing evidence of the expanding range of this rare but potentially life-threatening zoonotic disease. Babesia divergens is a causative agent of babesiosis in humans and cattle in Europe. The recently sequenced genome of B. divergens revealed over 3,741 protein coding-genes and the 10.7-Mb high-quality draft become the first reference tool to study the genome structure of B. divergens. Now, by exploiting this sequence data and using new computational tools and assembly strategies, we have significantly improved the quality of the B. divergens genome. The new assembly shows better continuity and has a higher correspondence to B. bovis chromosomes. Moreover, we present a differential expression analysis using RNA sequencing of the two different stages of the asexual lifecycle of B. divergens: the free merozoite capable of invading erythrocytes and the intraerythrocytic parasite stage that remains within the erythrocyte until egress. Comparison of mRNA levels of both stages identified 1,441 differentially expressed genes. From these, around half were upregulated and the other half downregulated in the intraerythrocytic stage. Orthogonal validation by real-time quantitative reverse transcription PCR confirmed the differential expression. A moderately increased expression level of genes, putatively involved in the invasion and egress processes, were revealed in the intraerythrocytic stage compared with the free merozoite. On the basis of these results and in the absence of molecular models of invasion and egress for B. divergens, we have proposed the identified genes as putative molecular players in the invasion and egress processes. Our results contribute to an understanding of key parasitic strategies and pathogenesis and could be a valuable genomic resource to exploit for the design of diagnostic methods, drugs and vaccines to improve the control of babesiosis. [ABSTRACT FROM AUTHOR]

Long-read sequencing and novel long-range assays have revolutionized de novo genome assembly by automating the reconstruction of reference-quality genomes. In particular, Hi-C sequencing is becoming an economical method for generating chromosome-scale scaffolds. Despite its increasing popularity, there are limited open-source tools available. Errors, particularly inversions and fusions across chromosomes, remain higher than alternate scaffolding technologies. We present a novel open-source Hi-C scaffolder that does not require an a priori estimate of chromosome number and minimizes errors by scaffolding with the assistance of an assembly graph. We demonstrate higher accuracy than the state-of-the-art methods across a variety of Hi-C library preparations and input assembly sizes. The Python and C++ code for our method is openly available at . [ABSTRACT FROM AUTHOR]

The Long interspersed nuclear element 1 (LINE-1) is a primary source of genetic variation in humans and other mammals. Despite its importance, LINE-1 activity remains difficult to study because of its highly repetitive nature. Here, we developed and validated a method called TeXP to gauge LINE-1 activity accurately. TeXP builds mappability signatures from LINE-1 subfamilies to deconvolve the effect of pervasive transcription from autonomous LINE-1 activity. In particular, it apportions the multiple reads aligned to the many LINE-1 instances in the genome into these two categories. Using our method, we evaluated well-established cell lines, cell-line compartments and healthy tissues and found that the vast majority (91.7%) of transcriptome reads overlapping LINE-1 derive from pervasive transcription. We validated TeXP by independently estimating the levels of LINE-1 autonomous transcription using ddPCR, finding high concordance. Next, we applied our method to comprehensively measure LINE-1 activity across healthy somatic cells, while backing out the effect of pervasive transcription. Unexpectedly, we found that LINE-1 activity is present in many normal somatic cells. This finding contrasts with earlier studies showing that LINE-1 has limited activity in healthy somatic tissues, except for neuroprogenitor cells. Interestingly, we found that the amount of LINE-1 activity was associated with the with the amount of cell turnover, with tissues with low cell turnover rates (e.g. the adult central nervous system) showing lower LINE-1 activity. Altogether, our results show how accounting for pervasive transcription is critical to accurately quantify the activity of highly repetitive regions of the human genome. [ABSTRACT FROM AUTHOR]

There is a need to investigate and better understand the inherited risk of cancer to ensure that clinical applications provide more accurate assessments and management strategies. Developing research-based next-generation sequencing gene panels that not only target (present-day) clinically actionable susceptibility genes but also genes that currently lack sufficient evidence for risk as well as candidate genes, such as those in DNA repair pathways, can help aid this effort. Therefore, gene panel B.O.P. (reast, varian, and rostate) was developed to evaluate the genetic risk of breast, ovarian and/or prostate cancer, and this manuscript serves as an introduction to B.O.P. and highlights its initial analytical validity assessment. B.O.P targets 87 genes that have been suggested, predicted, or clinically proven to be associated with breast, ovarian, and/or prostate cancer risk using Agilent Technologies Haloplex probes. The probes were designed for 100 base pair reads on an Illumina platform and target both coding and non-coding exons as well as 10 intronic base pairs flanking the intron-exon boundaries. The initial B.O.P screening involved 43 individuals from the Alabama Hereditary Cancer Cohort, and an average sequencing depth of 809X was obtained. Upon variant filtering and validation, true positives had an average sequencing depth of 659X and allele balance of 0.51. The average false positive sequencing depth was 34X and allele balance was 0.33. Although low sequencing depth was not always indicative of a false positive, high sequencing depths (>100X) signified a true positive. Furthermore, sensitivity and specificity of BRCA1/2 were calculated to be 100% and 92.3%, respectively. Overall, this screening enabled the establishment of criteria to alleviate future validation efforts and strongly supports the use of B.O.P. to further elucidate hereditary cancer susceptibility. Ultimately, continued B.O.P. screening will provide insights toward the genetic risk of and overlap between breast, ovarian, and/or prostate cancer. [ABSTRACT FROM AUTHOR]

Glaesserella parasuis is the cause of Glӓsser’s disease in pigs and is a significant contributor to post-weaning mortality in the swine industry. Prevention of G. parasuis disease relies primarily on bacterin vaccines, which have shown good homologous protection and variable heterologous protection. Bacterin production involves large scale growth of the bacteria and proteins produced during the proliferation phase of production become important antigens that stimulate the immune response. In order to evaluate genes activated during G. parasuis growth on different media substrates, the transcriptome of broth and agar grown G. parasuis strain 29755 were sequenced and compared. The transcription of most purported virulence genes were comparable between broth and agar grown G. parasuis; however, four virulence-associated genes, including ompA and vapD, had elevated expression under agar growth, while six virulence-associate genes had elevated expression during broth growth, including several protease genes. Additionally, there were metabolic shifts toward increased protein and lipid production and increased cellular division in broth grown G. parasuis. The results contribute to the understanding of how growth substrate alters gene transcription and protein expression, which may impact vaccine efficacy if immunogens important to the protective immune response are not produced under specific in vitro conditions. While the results of this work are unable to fully elucidate which growth medium presents a transcriptome more representative of in vivo samples or best suited for bacterin production, it forms a foundation that can be used for future comparisons and provides a better understanding of the metabolic differences in broth and agar grown bacteria. [ABSTRACT FROM AUTHOR]

While sequencing ancient DNA (aDNA) from archaeological material is now commonplace, very few attempts to sequence ancient transcriptomes have been made, even from typically stable deposition environments such as permafrost. This is presumably due to assumptions that RNA completely degrades relatively quickly, particularly when dealing with autolytic, nuclease-rich mammalian tissues. However, given the recent successes in sequencing ancient RNA (aRNA) from various sources including plants and animals, we suspect that these assumptions may be incorrect or exaggerated. To challenge the underlying dogma, we generated shotgun RNA data from sources that might normally be dismissed for such study. Here, we present aRNA data generated from two historical wolf skins, and permafrost-preserved liver tissue of a 14,300-year-old Pleistocene canid. Not only is the latter the oldest RNA ever to be sequenced, but it also shows evidence of biologically relevant tissue specificity and close similarity to equivalent data derived from modern-day control tissue. Other hallmarks of RNA sequencing (RNA-seq) data such as exon-exon junction presence and high endogenous ribosomal RNA (rRNA) content confirms our data’s authenticity. By performing independent technical library replicates using two high-throughput sequencing platforms, we show not only that aRNA can survive for extended periods in mammalian tissues but also that it has potential for tissue identification. aRNA also has possible further potential, such as identifying in vivo genome activity and adaptation, when sequenced using this technology. [ABSTRACT FROM AUTHOR]

Whole genome sequencing can provide essential public health information. However, it is now known that widely used short-read methods have the potential to miss some randomly-distributed segments of genomes. This can prevent phages, plasmids, and virulence factors from being detected or properly identified. Here, we compared assemblies of three complete Shiga toxin-producing Escherichia coli (STEC) O26:H11/H- genomes from two different sequence types (ST21 and 29), each acquired using the Nextera XT MiSeq, MinION nanopore-based sequencing, and Pacific Biosciences (PacBio) sequencing. Each closed genome consisted of a single chromosome, approximately 5.7 Mb for CFSAN027343, 5.6 Mb for CFSAN027346, and 5.4 MB for CFSAN027350. However, short-read whole genome sequencing (WGS) using Nextera XT MiSeq failed to identify some virulence genes in plasmids and on the chromosome, both of which were detected using the long-read platforms. Results from long-read MinION and PacBio allowed us to identify differences in plasmid content: a single 88 kb plasmid in CFSAN027343; a 157kb plasmid in CFSAN027350; and two plasmids in CFSAN027346 (one 95 Kb, one 72 Kb). These data enabled rapid characterization of the virulome, detection of antimicrobial genes, and composition/location of Stx phages. Taken together, positive correlations between the two long-read methods for determining plasmids, virulome, antimicrobial resistance genes, and phage composition support MinION sequencing as one accurate and economical option for closing STEC genomes and identifying specific virulence markers. [ABSTRACT FROM AUTHOR]

Hepatitis B infection is a world-wide public health burden causing serious liver complications. Previous studies suggest that hepatitis B integration into the human genome plays a crucial role in triggering oncogenic process and may also constitutively produce viral antigens. Despite the progress in HBV biology and sequencing technology, our fundamental understanding of how many hepatocytes in the liver actually carry viral integrations is still lacking. Herein we provide evidence that the HBV virus integrates with a lower-bound frequency of 0.84 per diploid genome in hepatitis B positive hepatocellular cancer patients. Moreover, we calculate that integrated viral DNA generates ~80% of the HBsAg transcripts in these patients. These results underscore the need to re-evaluate the clinical end-point and treatment strategies for chronic hepatitis B patients. [ABSTRACT FROM AUTHOR]

CD36 is a platelet membrane glycoprotein whose engagement with oxidized low-density lipoprotein (oxLDL) results in platelet activation. The CD36 gene has been associated with platelet count, platelet volume, as well as lipid levels and CVD risk by genome-wide association studies. Platelet CD36 expression levels have been shown to be associated with both the platelet oxLDL response and an elevated risk of thrombo-embolism. Several genomic variants have been identified as associated with platelet CD36 levels, however none have been conclusively demonstrated to be causative. We screened 81 expression quantitative trait loci (eQTL) single nucleotide polymorphisms (SNPs) associated with platelet CD36 expression by a Massively Parallel Reporter Assay (MPRA) and analyzed the results with a novel Bayesian statistical method. Ten eQTLs located 13kb to 55kb upstream of the CD36 transcriptional start site of transcript ENST00000309881 and 49kb to 92kb upstream of transcript ENST00000447544, demonstrated significant transcription shifts between their minor and major allele in the MPRA assay. Of these, rs2366739 and rs1194196, separated by only 20bp, were confirmed by luciferase assay to alter transcriptional regulation. In addition, electromobility shift assays demonstrated differential DNA:protein complex formation between the two alleles of this locus. Furthermore, deletion of the genomic locus by CRISPR/Cas9 in K562 and Meg-01 cells results in upregulation of CD36 transcription. These data indicate that we have identified a variant that regulates expression of CD36, which in turn affects platelet function. To assess the clinical relevance of our findings we used the PhenoScanner tool, which aggregates large scale GWAS findings; the results reinforce the clinical relevance of our variants and the utility of the MPRA assay. The study demonstrates a generalizable paradigm for functional testing of genetic variants to inform mechanistic studies, support patient management and develop precision therapies. [ABSTRACT FROM AUTHOR]

Pathogen monitoring, detection and removal are essential to public health and outbreak management. Systems are in place for monitoring the microbial load of hospitals and public health facilities with strategies to mitigate pathogen spread. However, no such strategies are in place for ambulances, which are tasked with transporting at-risk individuals in immunocompromised states. As standard culturing techniques require a laboratory setting, and are time consuming and labour intensive, our approach was designed to be portable, inexpensive and easy to use based on the MinION third-generation sequencing platform from Oxford Nanopore Technologies. We developed a transferable sampling-to-analysis pipeline to characterize the microbial community in emergency medical service vehicles. Our approach identified over sixty-eight organisms in ambulances to the genera level, with a proportion of these being connected with health-care associated infections, such as Clostridium spp. and Staphylococcus spp. We also monitored the microbiome of different locations across three ambulances over time, and examined the dynamic community of microorganisms found in emergency medical service vehicles. Observed differences identified hot spots, which may require heightened monitoring and extensive cleaning. Through metagenomics analysis it is also possible to identify how microorganisms spread between patients and colonize an ambulance over time. The sequencing results aid in the development of practices to mitigate disease spread, while also providing a useful tool for outbreak prediction through ongoing analysis of the ambulance microbiome to identify new and emerging pathogens. Overall, this pipeline allows for the tracking and monitoring of pathogenic microorganisms of epidemiological interest, including those related to health-care associated infections. [ABSTRACT FROM AUTHOR]

Background: Although the tprK gene of Treponema pallidum are thought to play a critical role in the pathogenesis of syphilis, the profile of variations in tprK during the development of human syphilis infection have remained unclear. Methods/Principal findings: Through next-generation sequencing, we compared the tprK gene of 14 secondary syphilis patients with that of 14 primary syphilis patients, and the results showed an increased number of variants within the seven V regions of the tprK gene in the secondary syphilis samples. The length of the sequences within each V region also presented a 3-bp changing pattern. Interestingly, the frequencies of predominant sequences within the V regions in the secondary syphilis samples were generally decreased compared with those found in the primary syphilis samples, particularly in the V7 region, where a frequency below 60% was found in up to 57% (8/14) of all secondary samples compared with 7% (1/14) of all primary samples. Moreover, the number of minor variants distributed between frequencies of 10 and 49.9% was increased. The alignment of all amino acid sequences within each V region of the primary and secondary syphilis samples revealed that some amino acid sequences, particularly the amino acid sequences IASDGGAIKH and IASEDGSAGNLKH in V1, were highly stable. Additionally, the amino acid sequences in V6 also exhibited notable intrastrain heterogeneity and were likely to form a strain-specific pattern at the interstrain level. Conclusions: The identification of different profiles of the tprK gene in primary and secondary syphilis patients indicated that the tprK gene of T. pallidum undergoes constant variation to result in the best adaptation to the host. The highly stable peptides found in V1 are likely promising potential vaccine components. The highly heterogenetic regions (e.g., V6) could help to understand the role of tprK in immune evasion. [ABSTRACT FROM AUTHOR]

Microvascular proliferation (MVP), an aberrant vascular structure containing multilayered mitotically active endothelial- and smooth-muscle cells/pericytes, is a histopathological hallmark of glioblastoma multiforme (GBM). Although MVP tends to be associated with high-grade glioma, it has also been detected in WHO grade I pilocytic astrocytoma (PA). However, little is known about the mechanism underlying its formation. Using TP53 point mutations as a marker for tumor-derived cells, we earlier reported that MVP was partially converted from tumor cells via mesenchymal transition. In the current study we used the KIAA1549-BRAF fusion gene as a marker to assess whether MVPs in PA contained tumor-derived cells and/or phenotypically distinct tumor cells expressing vascular markers. cDNA synthesized from frozen tissue of six PA patients operated at our institute was analyzed to detect the KIAA1549-BRAF fusion gene by reverse transcription polymerase chain reaction (RT-PCR) assay. The breakpoint in the fusion gene was identified by long and accurate PCR (LA-PCR) and Sanger sequencing of genomic DNA. Distinct tumor cells and cellular components of MVP were obtained by laser microdissection. For the qualitative and quantitative detection of the KIAA1549-BRAF fusion gene we performed genomic and digital PCR assays. Fluorescence in situ hybridization (FISH) was used to assess gene fusion in cellular components of MVP. Samples from three PA patients harbored the KIAA1549 exon 15, BRAF exon 9 fusion gene. In two patient samples with abundant MVP, RT-PCR assay detected strong bands arising from the KIAA1549-BRAF fusion gene in both tumor cells and cellular components of MVP. Digital PCR showed that vis-à-vis tumor tissue, its relative expression in cellular components of MVP was 42% in one- and 76% in another sample. FISH revealed amplified signals in both tumor cells and cellular components of MVP indicative of tandem duplication. Our findings suggest that in patients with PA, some cellular components of MVP contained tumor derived cell and/or phenotypically distinct tumor cells expressing vascular markers. [ABSTRACT FROM AUTHOR]

Garden asparagus (Asparagus officinalis L.) is a popular vegetable cultivated worldwide. The secondary metabolites in its shoot are helpful for human health. We analyzed A. officinalis transcriptomes and identified differentially expressed genes (DEGs) involved in the biosynthesis of rutin and protodioscin, which are health-promoting functional compounds, and determined their association with stem color. We sequenced the complete mRNA transcriptome using the Illumina high-throughput sequencing platform in one white, three green, and one purple asparagus cultivars. A gene set was generated by de novo assembly of the transcriptome sequences and annotated using a BLASTx search. To investigate the relationship between the contents of rutin and protodioscin and their gene expression levels, rutin and protodioscin were analyzed using high-performance liquid chromatography. A secondary metabolite analysis using high-performance liquid chromatography showed that the rutin content was higher in green asparagus, while the protodioscin content was higher in white asparagus. We studied the genes associated with the biosynthesis of the rutin and protodioscin. The transcriptomes of the five cultivars generated 336 599 498 high-quality clean reads, which were assembled into 239 873 contigs with an average length of 694 bp, using the Trinity v2.4.0 program. The green and white asparagus cultivars showed 58 932 DEGs. A comparison of rutin and protodioscin biosynthesis genes revealed that 12 of the 57 genes associated with rutin and two of the 50 genes associated with protodioscin showed more than four-fold differences in expression. These DEGs might have caused a variation in the contents of these two metabolites between green and white asparagus. The present study is possibly the first to report transcriptomic gene sets in asparagus. The DEGs putatively involved in rutin and protodioscin biosynthesis might be useful for molecular engineering in asparagus. [ABSTRACT FROM AUTHOR]

Translucent flash disorder (TFD) is one of the important physiological disorders in mangosteen (Garcinia mangostana L.). TFD has symptoms such as flesh arils that become firm and appear transparent similar to watercore in apple or pear. Information on the changes of gene expression in TFD-affected tissues remain limited, and investigations into the effects of different water regimes still need to be undertaken. Through an RNA sequencing approach using the Ion Proton, 183,274 contigs with length ranging from 173–13,035 bp were constructed by de novo assembly. Functional annotation was analyzed using various public databases such as non-redundant protein NCBI, SwissProt, and Gene Ontology, and KEGG pathway. Our studies compared different water regimes to incidence and differentially expressed genes of TFD-like physiological disorders. From the differentially expressed gene (DEG) between normal air and TFD-affected aril, we identified DEG-related TFD events, which 6228 DEGs in the control condition and 3327 DEGs in under water stress treatment condition remained, and confirmed these with RT-qPCR, including sucrose synthase (SUSY), endoglucanase (GUN), xyloglucan endotransglucosylase/hydrolase (XTH), and polygalacturonase (PG) showed statistically significant. In addition, transcription factors also indicated changes in MYB, NAC and WRKY between tissues and different water regimes. [ABSTRACT FROM AUTHOR]

Objetives: The aim of this study was to identify CMV drug resistance mutations (DRM) in solid organ transplant (SOT) recipients with suspected resistance comparing next-generation sequencing (NGS) with Sanger sequencing and assessing risk factors and the clinical impact of resistance. Methods: Using Sanger sequencing as the reference method, we prospectively assessed the ability of NGS to detect CMV DRM in the UL97 and UL54 genes in a nationwide observational study from September 2013 to August 2016. Results: Among 44 patients recruited, 14 DRM were detected by Sanger in 12 patients (27%) and 20 DRM were detected by NGS, in 16 (36%). NGS confirmed all the DRM detected by Sanger. The additional six mutations detected by NGS were present in <20% of the sequenced population, being located in the UL97 gene and conferring high-level resistance to ganciclovir. The presence of DRM by NGS was associated with lung transplantation (p = 0.050), the administration of prophylaxis (p = 0.039), a higher mean time between transplantation and suspicion of resistance (p = 0.038) and longer antiviral treatment duration before suspicion (p = 0.024). However, the latter was the only factor independently associated with the presence of DRM by NGS in the multivariate analysis (OR 2.24, 95% CI 1.03 to 4.87). Conclusions: NGS showed a higher yield than Sanger sequencing for detecting CMV resistance mutations in SOT recipients. The presence of DRM detected by NGS was independently associated with longer antiviral treatment. [ABSTRACT FROM AUTHOR]

Disconnected interacting protein 2 homolog A (DIP2A) is highly expressed in nervous system and respiratory system of developing embryos. However, genes regulated by Dip2a in developing brain and lung have not been systematically studied. Transcriptome of brain and lung in embryonic 19.5 day (E19.5) were compared between wild type and Dip2a-/- mice. An average of 50 million reads per sample was mapped to the reference sequence. A total of 214 DEGs were detected in brain (82 up and 132 down) and 1900 DEGs in lung (1259 up and 641 down). GO enrichment analysis indicated that DEGs in both Brain and Lung were mainly enriched in biological processes ‘DNA-templated transcription and Transcription from RNA polymerase II promoter’, ‘multicellular organism development’, ‘cell differentiation’ and ‘apoptotic process’. In addition, COG classification showed that both were mostly involved in ‘Replication, Recombination, and Repair’, ‘Signal transduction and mechanism’, ‘Translation, Ribosomal structure and Biogenesis’ and ‘Transcription’. KEGG enrichment analysis showed that brain was mainly enriched in ‘Thyroid cancer’ pathway whereas lung in ‘Complement and Coagulation Cascades’ pathway. Transcription factor (TF) annotation analysis identified Zinc finger domain containing (ZF) proteins were mostly regulated in lung and brain. Interestingly, study identified genes Skor2, Gpr3711, Runx1, Erbb3, Frmd7, Fut10, Sox11, Hapln1, Tfap2c and Plxnb3 from brain that play important roles in neuronal cell maturation, differentiation, and survival; genes Hoxa5, Eya1, Errfi1, Sox11, Shh, Igf1, Ccbe1, Crh, Fgf9, Lama5, Pdgfra, Ptn, Rbp4 and Wnt7a from lung are important in lung development. Expression levels of the candidate genes were validated by qRT-PCR. Genome wide transcriptional analysis using wild type and Dip2a knockout mice in brain and lung at embryonic day 19.5 (E19.5) provided a genetic basis of molecular function of these genes. [ABSTRACT FROM AUTHOR]

Background: Assessment of KRAS, NRAS (RAS) and BRAF mutations is a standard in the management of patients with metastatic colorectal cancer (mCRC). Mutations could be assessed using next-generation sequencing (NGS) or real-time PCR-based assays. Times to results are 1 to 2 weeks for NGS and 1 to 3 days for real-time PCR-based assays. Using NGS can delay first-line treatment commencement and using PCR-based assays is limited by the number of possible analysed targets. The Idylla system is a real-time PCR cartridge-based assay, able to analyse hotspots mutations using one section of FFPE tumour tissue sample. To combine short delays and analysis of a large gene-panel, we propose here a laboratory workflow combining the Idylla system and NGS and compatible with FFPE samples with low tissue quantity. In this study we evaluated and validated the Idylla system for the analysis of RAS and BRAF mutations by pipetting directly DNA in the cartridge instead of FFPE section as recommended by the manufacturer. Materials and methods: DNA extracted from 29 FFPE samples from mCRC patients with NGS-characterized RAS and BRAF mutations were tested with the Idylla KRAS and the Idylla NRAS-BRAF mutation tests to assess sensitivity, specificity, reproducibility and limit of detection of each test. Results: A 100% concordance was found between NGS and Idylla results for the determination of KRAS (12/12), NRAS (12/12) and BRAF (11/11) mutations with a sensitivity and a specificity of 100%. The system showed a good reproducibility with CV inferior to 3%. LOD was reached with 2.5 ng of DNA for KRAS and NRAS mutations and 5 ng of DNA for BRAF mutations. Conclusions: The analysis of RAS and BRAF mutations using DNA pipetted directly in the cartridge of the Idylla system showed a good sensitivity, specificity, reproducibility and LOD, and can be integrated in a laboratory workflow for samples with few tissue without compromising a further complete tumour characterization using NGS. [ABSTRACT FROM AUTHOR]

Bloodstream infection (BSI) is a life-threatening condition characterized by the presence of pathogens in the blood. It is associated with increased morbidity and mortality, and has to be treated promptly as mortality increases with every hour of delayed treatment. Therefore, rapid and sensitive diagnosis of BSI is essential. The routine diagnostic method for BSI is blood culture, which can only detect culturable pathogens and takes several days to obtain results. The 16S rRNA gene is present in all bacteria and is commonly used as a target for universal bacterial detection in rapid molecular assays such as PCR. However, molecular detection of the 16S gene is hampered by the large amount of human DNA found in blood samples, making diagnostic results aspecific and less sensitive. We have optimized the selection of PCR primers targeting the 16S rRNA gene to avoid cross-reaction with human DNA background. The developed method increases specificity and sensitivity for pathogen diagnosis, and provides rapid and accurate pathogen detection for rare bacterial DNA in the presence of abundant host DNA. [ABSTRACT FROM AUTHOR]

The relationship between noncoding DNA sequence and gene expression is not well-understood. Massively parallel reporter assays (MPRAs), which quantify the regulatory activity of large libraries of DNA sequences in parallel, are a powerful approach to characterize this relationship. We present MPRA-DragoNN, a convolutional neural network (CNN)-based framework to predict and interpret the regulatory activity of DNA sequences as measured by MPRAs. While our method is generally applicable to a variety of MPRA designs, here we trained our model on the Sharpr-MPRA dataset that measures the activity of ∼500,000 constructs tiling 15,720 regulatory regions in human K562 and HepG2 cell lines. MPRA-DragoNN predictions were moderately correlated (Spearman ρ = 0.28) with measured activity and were within range of replicate concordance of the assay. State-of-the-art model interpretation methods revealed high-resolution predictive regulatory sequence features that overlapped transcription factor (TF) binding motifs. We used the model to investigate the cell type and chromatin state preferences of predictive TF motifs. We explored the ability of our model to predict the allelic effects of regulatory variants in an independent MPRA experiment and fine map putative functional SNPs in loci associated with lipid traits. Our results suggest that interpretable deep learning models trained on MPRA data have the potential to reveal meaningful patterns in regulatory DNA sequences and prioritize regulatory genetic variants, especially as larger, higher-quality datasets are produced. [ABSTRACT FROM AUTHOR]

In 2010, the rare OP354-like P[8]b rotavirus subtype was detected in children less than 2 years old in Ghana. In this follow-up study, to provide insight into the evolutionary history of the genome of Ghanaian P[8]b strains RVA/Human-wt/GHA/GHDC949/2010/G9P[8] and RVA/Human-wt/GHA/GHM0094/2010/G9P[8] detected in an infant and a 7-month old child hospitalised for acute gastroenteritis, we sequenced the complete genome using both Sanger sequencing and Illumina MiSeq technology followed by phylogenetic analysis of the near-full length sequences. Both strains possessed the Wa-like/genotype 1 constellation G9P[8]b-I1-R1-C1-M1-A1-N1-T1-E1-H1. Sequence comparison and phylogenetic inference showed that both strains were identical at the lineage level throughout the 11 genome segments. Their VP7 sequences belonged to the major sub-lineage of the G9-lineage III whereas their VP4 sequences belonged to P[8]b cluster I. The VP7 and VP4 genes of the study strains were closely related to a Senegalese G9P[8]b strain detected in 2009. In the remaining nine genome segments, both strains consistently clustered together with Wa-like RVA strains possessing either P[8]a or P[8]b mostly of African RVA origin. The introduction of a P[8]b subtype VP4 gene into the stable Wa-like strain backbone may result in strains that might propagate easily in the human population, with a potential to become an important public health concern, especially because it is not certain if the monovalent rotavirus vaccine (Rotarix) used in Ghana will be efficacious against such strains. Our analysis of the full genomes of GHM0094 and GHDC949 adds to knowledge of the genetic make-up and evolutionary dynamics of P[8]b rotavirus strains. [ABSTRACT FROM AUTHOR]

The stellate ganglia are the predominant source of sympathetic innervation to the heart. Remodeling of sympathetic nerves projecting to the heart has been observed in several cardiovascular diseases, and sympathetic dysfunction contributes to cardiac pathology. Wistar Kyoto rats are a common model for the study of cardiovascular diseases, but we lack a profile of the baseline transcriptomic and neurochemical characteristics of their cardiac sympathetic neurons. Most studies of cardiovascular disease have used male animals only, but in the future both male and female animals will be used for these types of studies; therefore, we sought to characterize the transcriptome of male and female stellate ganglia and to correlate that with catecholamine and acetylcholine content in the heart. We have generated a dataset of baseline RNA expression in male and female Wistar Kyoto rat stellate ganglia using RNA-seq, and have measured neurotransmitter levels in heart and stellate ganglia using HPLC and mass spectrometry. We identified numerous gene expression differences between male and female stellates, including genes encoding important developmental factors, receptors and neuropeptides. Female hearts had significantly higher neurotransmitter content than male hearts; however, no significant differences were detected in expression of the genes encoding neurotransmitter synthetic enzymes. Similarly, no statistically significant differences were identified between the sexes in cardiac tyrosine hydroxylase levels. [ABSTRACT FROM AUTHOR]

Background: Ischemic retinopathies (IRs) are leading causes of visual impairment. They are characterized by an initial phase of microvascular degeneration and a second phase of aberrant pre-retinal neovascularization (NV). microRNAs (miRNAs) regulate gene expression, and a number play a role in normal and pathological NV. But, post-transcriptional modulation of miRNAs in the eye during the development of IRs has not been systematically evaluated. Aims & methods: Using Next Generation Sequencing (NGS) we profiled miRNA expression in the retina and choroid during vasodegenerative and NV phases of oxygen-induced retinopathy (OIR). Results: Approximately 20% of total miRNAs exhibited altered expression (up- or down-regulation); 6% of miRNA were found highly expressed in retina and choroid of rats subjected to OIR. During OIR-induced vessel degeneration phase, miR-199a-3p, -199a-5p, -1b, -126a-3p displayed a robust decreased expression (> 85%) in the retina. While in the choroid, miR-152-3p, -142-3p, -148a-3p, -532-3p were upregulated (>200%) and miR-96-5p, -124-3p, -9a-3p, -190b-5p, -181a-1-3p, -9a-5p, -183-5p were downregulated (>70%) compared to controls. During peak pathological NV, miR-30a-5p, -30e-5p and 190b-5p were markedly reduced (>70%), and miR-30e-3p, miR-335, -30b-5p strongly augmented (by up to 300%) in the retina. Whereas in choroid, miR-let-7f-5p, miR-126a-5p and miR-101a-3p were downregulated by (>81%), and miR-125a-5p, let-7e-5p and let-7g-5p were upregulated by (>570%) during NV. Changes in miRNA observed using NGS were validated using qRT-PCR for the 24 most modulated miRNAs. In silico approach to predict miRNA target genes (using algorithms of miRSystem database) identified potential new target genes with pro-inflammatory, apoptotic and angiogenic properties. Conclusion: The present study is the first comprehensive description of retinal/choroidal miRNAs profiling in OIR (using NGS technology). Our results provide a valuable framework for the characterization and possible therapeutic potential of specific miRNAs involved in ocular IR-triggered inflammation, angiogenesis and degeneration. [ABSTRACT FROM AUTHOR]

Di(2-ethylhexyl)phthalate (DEHP) interferes with sex hormones signaling pathways (SHP). C57BL/6J mice prenatally exposed to 300 mg/kg/day DEHP develop a testicular dysgenesis syndrome (TDS) at adulthood, but similarly-exposed FVB/N mice are not affected. Here we aim to understand the reasons behind this drastic difference that should depend on the genome of the strain. In both backgrounds, pregnant female mice received per os either DEHP or corn oil vehicle and the male filiations were examined. Computer-assisted sperm analysis showed a DEHP-induced decreased sperm count and velocities in C57BL/6J. Sperm RNA sequencing experiments resulted in the identification of the 62 most differentially expressed RNAs. These RNAs, mainly regulated by hormones, produced strain-specific transcriptional responses to prenatal exposure to DEHP; a pool of RNAs was increased in FVB, another pool of RNAs was decreased in C57BL/6J. In FVB/N, analysis of non-synonymous single nucleotide polymorphisms (SNP) impacting SHP identified rs387782768 and rs29315913 respectively associated with absence of the Forkhead Box A3 (Foxa3) RNA and increased expression of estrogen receptor 1 variant 4 (NM_001302533) RNA. Analysis of the role of SNPs modifying SHP binding sites in function of strain-specific responses to DEHP revealed a DEHP-resistance allele in FVB/N containing an additional FOXA1-3 binding site at rs30973633 and four DEHP-induced beta-defensins (Defb42, Defb30, Defb47 and Defb48). A DEHP-susceptibility allele in C57BL/6J contained five SNPs (rs28279710, rs32977910, rs46648903, rs46677594 and rs48287999) affecting SHP and six genes (Svs2, Svs3b, Svs4, Svs3a, Svs6 and Svs5) epigenetically silenced by DEHP. Finally, targeted experiments confirmed increased methylation in the Svs3ab promoter with decreased SEMG2 persisting across generations, providing a molecular explanation for the transgenerational sperm velocity decrease found in C57BL/6J after DEHP exposure. We conclude that the existence of SNP-dependent mechanisms in FVB/N inbred mice may confer resistance to transgenerational endocrine disruption. [ABSTRACT FROM AUTHOR]

Avian pathogenic Escherichia coli (APEC) causes great economic loss to the poultry industry worldwide. Chicken type II pneumocytes (CP II cells) secrete surfactants and modulate lung immunity to decrease the infection of the invading pathogen. Nevertheless, the pathogenesis of CP II cells to APEC infection remains poorly understood. Therefore, we conducted global gene expression profiling of CP II cells after APEC-O78 infection to explore the host-pathogen interaction. The differentially expressed genes of CP II cells to APEC infection were characterized by RNA-seq with EB-seq algorithm. In consequence, the mRNA of 18996 genes was identified, and CP II cells responded to APEC infection with marked changes in the expression of 1390 genes. Among them, there are 803 down-regulated mRNAs and 587 up-regulated mRNAs. The KEGG prediction and Gene Ontology terms analysis revealed that the major enriched pathways were related to NF-κB signaling pathway, apoptosis pathway, tight junction, and cytokine-cytokine receptor interaction and other pathways. We adopted qRT-PCR to verify the validity of the selected gene expression. The fold induction of qPCR was similar to the RNA-seq results. These results provide a better understanding of the pathogenesis of APEC, especially apoptosis pathway involved in APEC infection. [ABSTRACT FROM AUTHOR]

The Oxford Nanopore MinION DNA sequencing device can produce large amounts of long sequences, typically several kilobases, within a few hours. This long read capacity was exploited to detect antimicrobial resistance genes (ARGs) in a large veterinary teaching hospital environment, and to assess their taxonomic origin, genetic organisation and association with mobilisation markers concurrently. Samples were collected on eight occasions between November 2016 and May 2017 (inclusive) in a longitudinal study. Nanopore sequencing was performed on total DNA extracted from the samples after a minimal enrichment step in broth. Many ARGs present in the veterinary hospital environment could potentially confer resistance to antimicrobials widely used in treating infections of companion animals, including aminoglycosides, extended-spectrum beta-lactams, sulphonamides, macrolides, and tetracyclines. High-risk ARGs, defined here as single or multiple ARGs associated with pathogenic bacterial species or with mobile genetic elements, were shared between the intensive care unit (ICU) patient cages, a dedicated laundry trolley and a floor cleaning mop-bucket. By contrast, a floor surface from an office corridor without animal contact and located outside the veterinary hospital did not contain such high-risk ARGs. Relative abundances of high-risk ARGs and co-localisation of these genes on the same sequence read were higher in the laundry trolley and mop bucket samples, compared to the ICU cages, suggesting that amplification of ARGs is likely to occur in the collection points for hospital waste. These findings have prompted the implementation of targeted intervention measures in the veterinary hospital to mitigate the risks of transferring clinically important ARGs between sites and to improve biosecurity practices in the facility. [ABSTRACT FROM AUTHOR]

Several horse breeds have been specifically selected for the ability to exhibit alternative patterns of locomotion, or gaits. A premature stop codon in the gene DMRT3 is permissive for “gaitedness” across breeds. However, this mutation is nearly fixed in both American Standardbred trotters and pacers, which perform a diagonal and lateral gait, respectively, during harness racing. This suggests that modifying alleles must influence the preferred gait at racing speeds in these populations. A genome-wide association analysis for the ability to pace was performed in 542 Standardbred horses (n = 176 pacers, n = 366 trotters) with genotype data imputed to ~74,000 single nucleotide polymorphisms (SNPs). Nineteen SNPs on nine chromosomes (ECA1, 2, 6, 9, 17, 19, 23, 25, 31) reached genome-wide significance (p < 1.44 x 10−6). Variant discovery in regions of interest was carried out via whole-genome sequencing. A set of 303 variants from 22 chromosomes with putative modifying effects on gait was genotyped in 659 Standardbreds (n = 231 pacers, n = 428 trotters) using a high-throughput assay. Random forest classification analysis resulted in an out-of-box error rate of 0.61%. A conditional inference tree algorithm containing seven SNPs predicted status as a pacer or trotter with 99.1% accuracy and subsequently performed with 99.4% accuracy in an independently sampled population of 166 Standardbreds (n = 83 pacers, n = 83 trotters). This highly accurate algorithm could be used by owners/trainers to identify Standardbred horses with the potential to race as pacers or as trotters, according to the genotype identified, prior to initiating training and would enable fine-tuning of breeding programs with designed matings. Additional work is needed to determine both the algorithm’s utility in other gaited breeds and whether any of the predictive SNPs play a physiologically functional role in the tendency to pace or tag true functional alleles. [ABSTRACT FROM AUTHOR]