Search

Abstract Background Long Adapter Single-Stranded Oligonucleotide (LASSO) probes were developed as a novel tool for massively parallel cloning of kilobase-long genomic DNA sequences. LASSO dramatically improves the capture length limit of current DNA padlock probe technology from approximately 150 bps to several kbps. High-throughput LASSO capture involves the parallel assembly of thousands of probes. However, malformed probes are indiscernible from properly formed probes using gel electrophoretic techniques. Therefore, we used next-generation sequencing (NGS) to assess the efficiency of LASSO probe assembly and how it relates to the nature of DNA capture and amplification. Additionally, we introduce a simplified single target LASSO protocol using classic molecular biology techniques for qualitative and quantitative assessment of probe specificity. Results A LASSO probe library targeting 3164 unique E. coli ORFs was assembled using two different probe assembly reaction conditions with a 40-fold difference in DNA concentration. Unique probe sequences are located within the first 50 bps of the 5′ and 3′ ends, therefore we used paired-end NGS to assess probe library quality. Properly mapped read pairs, representing correctly formed probes, accounted for 10.81 and 0.65% of total reads, corresponding to ~ 80% and ~ 20% coverage of the total probe library for the lower and higher DNA concentration conditions, respectively. Subsequently, we used single-end NGS to correlate probe assembly efficiency and capture quality. Significant enrichment of LASSO targets over non-targets was only observed for captures done using probes assembled with a lower DNA concentration. Additionally, semi-quantitative polyacrylamide gel electrophoresis revealed a ~ 10-fold signal-to-noise ratio of LASSO capture in a simplified system. Conclusions These results suggest that LASSO probe coverage for target sequences is more predictive of successful capture than probe assembly depth-enrichment. Concomitantly, these results demonstrate that DNA concentration at a critical step in the probe assembly reaction significantly impacts probe formation. Additionally, we show that a simplified LASSO capture protocol coupled to PAGE (polyacrylamide gel electrophoresis) is highly specific and more amenable to small-scale LASSO approaches, such as screening novel probes and templates.

Abstract Background Acoustic or hydrodynamic shearing, sonication and enzymatic digestion are used to fragment DNA. However, these methods have several disadvantages, such as DNA damage, difficulties in fragmentation control, irreproducibility and under-representation of some DNA segments. The DNA fragmentation tool would be a gentle enzymatic method, offering cleavage frequency high enough to eliminate DNA fragments distribution bias and allow for easy control of partial digests. Only three such frequently cleaving natural restriction endonucleases (REases) were discovered: CviJI, SetI and FaiI. Therefore, we have previously developed two artificial enzymatic specificities, cleaving DNA approximately every ~ 3-bp: TspGWI/sinefungin (SIN) and TaqII/SIN. Results In this paper we present the third developed specificity: TthHB27I/SIN(SAM) - a new genomic tool, based on Type IIS/IIC/IIG Thermus-family REases-methyltransferases (MTases). In the presence of dimethyl sulfoxide (DMSO) and S-adenosyl-L-methionine (SAM) or its analogue SIN, the 6-bp cognate TthHB27I recognition sequence 5’-CAARCA-3′ is converted into a combined 3.2–3.0-bp ‘site’ or its statistical equivalent, while a cleavage distance of 11/9 nt is retained. Protocols for various modes of limited DNA digestions were developed. Conclusions In the presence of DMSO and SAM or SIN, TthHB27I is transformed from rare 6-bp cutter to a very frequent one, approximately 3-bp. Thus, TthHB27I/SIN(SAM) comprises a new tool in the very low-represented segment of such prototype REases specificities. Moreover, this modified TthHB27I enzyme is uniquely suited for controlled DNA fragmentation, due to partial DNA cleavage, which is an inherent feature of the Thermus-family enzymes. Such tool can be used for quasi-random libraries generation as well as for other DNA manipulations, requiring high frequency cleavage and uniform distribution of cuts along DNA.

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]

We have used a transposon insertion sequencing (TIS) approach to establish the fitness landscape of the African Salmonella enterica serovar Typhimurium ST313 strain D23580, to complement our previous comparative genomic and functional transcriptomic studies. We used a genome-wide transposon library with insertions every 10 nucleotides to identify genes required for survival and growth in vitro and during infection of murine macrophages. The analysis revealed genomic regions important for fitness under two in vitro growth conditions. Overall, 724 coding genes were required for optimal growth in LB medium, and 851 coding genes were required for growth in SPI-2-inducing minimal medium. These findings were consistent with the essentiality analyses of other S. Typhimurium ST19 and S. Typhi strains. The global mutagenesis approach also identified 60 sRNAs and 413 intergenic regions required for growth in at least one in vitro growth condition. By infecting murine macrophages with the transposon library, we identified 68 genes that were required for intra-macrophage replication but did not impact fitness in vitro. None of these genes were unique to S. Typhimurium D23580, consistent with a high conservation of gene function between S. Typhimurium ST313 and ST19 and suggesting that novel virulence factors are not involved in the interaction of strain D23580 with murine macrophages. We discovered that transposon insertions rarely occurred in many pBT1 plasmid-encoded genes (36), compared with genes carried by the pSLT-BT virulence plasmid and other bacterial plasmids. The key essential protein encoded by pBT1 is a cysteinyl-tRNA synthetase, and our enzymological analysis revealed that the plasmid-encoded CysRSpBT1 had a lower ability to charge tRNA than the chromosomally-encoded CysRSchr enzyme. The presence of aminoacyl-tRNA synthetases in plasmids from a range of Gram-negative and Gram-positive bacteria suggests that plasmid-encoded essential genes are more common than had been appreciated. [ABSTRACT FROM AUTHOR]

Fin whales (Balaenoptera physalus) and blue whales (B. musculus) are the two largest species on Earth and are widely distributed across the world’s oceans. Hybrids between these species appear to be relatively widespread and have been reported in both the North Atlantic and North Pacific; they are also relatively common, and have been proposed to occur once in every thousand fin whales. However, despite known hybridization, fin and blue whales are not sibling species. Rather, the closest living relative of fin whales are humpback whales (Megaptera novaeangliae). To improve the quality of fin whale data available for analysis, we assembled and annotated a fin whale nuclear genome using in-silico mate pair libraries and previously published short-read data. Using this assembly and genomic data from a humpback, blue, and bowhead whale, we investigated whether signatures of introgression between the fin and blue whale could be found. We find no signatures of contemporary admixture in the fin and blue whale genomes, although our analyses support ancestral gene flow between the species until 2.4–1.3 Ma. We propose the following explanations for our findings; i) fin/blue whale hybridization does not occur in the populations our samples originate from, ii) contemporary hybrids are a recent phenomenon and the genetic consequences have yet to become widespread across populations, or iii) fin/blue whale hybrids are under large negative selection, preventing them from backcrossing and contributing to the parental gene pools. [ABSTRACT FROM AUTHOR]

Advances in genomics have made whole genome studies increasingly feasible across the life sciences. However, new technologies and algorithmic advances do not guarantee flawless genomic sequences or annotation. Bias, errors, and artifacts can enter at any stage of the process from library preparation to annotation. When planning an experiment that utilizes a genome sequence as the basis for the design, there are a few basic checks that, if performed, may better inform the experimental design and ideally help avoid a failed experiment or inconclusive result. [ABSTRACT FROM AUTHOR]

Unknown sequences, or gaps, are present in many published genomes across public databases. Gap filling is an important finishing step in de novo genome assembly, especially in large genomes. The gap filling problem is nontrivial and while there are many computational tools partially solving the problem, several have shortcomings as to the reliability and correctness of the output, i.e. the gap filled draft genome. SSPACE-LongRead is a scaffolding tool that utilizes long reads from multiple third-generation sequencing platforms in finding links between contigs and combining them. The long reads potentially contain sequence information to fill the gaps created in the scaffolding, but SSPACE-LongRead currently lacks this functionality. We present an automated pipeline called gapFinisher to process SSPACE-LongRead output to fill gaps after the scaffolding. gapFinisher is based on the controlled use of a previously published gap filling tool FGAP and works on all standard Linux/UNIX command lines. We compare the performance of gapFinisher against two other published gap filling tools PBJelly and GMcloser. We conclude that gapFinisher can fill gaps in draft genomes quickly and reliably. In addition, the serial design of gapFinisher makes it scale well from prokaryote genomes to larger genomes with no increase in the computational footprint. [ABSTRACT FROM AUTHOR]

Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms. [ABSTRACT FROM AUTHOR]

The confluence of deep sequencing and powerful machine learning is providing an unprecedented peek at the darkest of the dark genomic matter, the non-coding genomic regions lacking any functional annotation. While deep sequencing uncovers rare tumor variants, the heterogeneity of the disease confounds the best of machine learning (ML) algorithms. Here we set out to answer if the dark-matter of the genome encompass signals that can distinguish the fine subtypes of disease that are otherwise gnomically indistinguishable. We introduce a novel stochastic regularization, ReVeaL, that empowers ML to discriminate subtle cancer subtypes even from the same ‘cell of origin’. Analogous to heritability, implicitly defined on whole genome, we use predictability (F1 score) definable on portions of the genome. In an effort to distinguish cancer subtypes using dark-matter DNA, we applied ReVeaL to a new WGS dataset from 727 patient samples with seven forms of hematological cancers and assessed the predictivity over several genomic regions including genic, non-dark, non-coding, non-genic, and dark. ReVeaL enabled improved discrimination of cancer subtypes for all segments of the genome. The non-genic, non-coding and dark-matter had the highest F1 scores, with dark-matter having the highest level of predictability. Based on ReVeaL’s predictability of different genomic regions, dark-matter contains enough signal to significantly discriminate fine subtypes of disease. Hence, the agglomeration of rare variants, even in the hitherto unannotated and ill-understood regions of the genome, may play a substantial role in the disease etiology and deserve much more attention. [ABSTRACT FROM AUTHOR]

Genome engineering is a powerful approach to study how chromosomal architecture impacts phenotypes. However, quantifying the fitness impact of translocations independently from the confounding effect of base substitutions has so far remained challenging. We report a novel application of the CRISPR/Cas9 technology allowing to generate with high efficiency both uniquely targeted and multiple concomitant reciprocal translocations in the yeast genome. Targeted translocations are constructed by inducing two double-strand breaks on different chromosomes and forcing the trans-chromosomal repair through homologous recombination by chimerical donor DNAs. Multiple translocations are generated from the induction of several DSBs in LTR repeated sequences and promoting repair using endogenous uncut LTR copies as template. All engineered translocations are markerless and scarless. Targeted translocations are produced at base pair resolution and can be sequentially generated one after the other. Multiple translocations result in a large diversity of karyotypes and are associated in many instances with the formation of unanticipated segmental duplications. To test the phenotypic impact of translocations, we first recapitulated in a lab strain the SSU1/ECM34 translocation providing increased sulphite resistance to wine isolates. Surprisingly, the same translocation in a laboratory strain resulted in decreased sulphite resistance. However, adding the repeated sequences that are present in the SSU1 promoter of the resistant wine strain induced sulphite resistance in the lab strain, yet to a lower level than that of the wine isolate, implying that additional polymorphisms also contribute to the phenotype. These findings illustrate the advantage brought by our technique to untangle the phenotypic impacts of structural variations from confounding effects of base substitutions. Secondly, we showed that strains with multiple translocations, even those devoid of unanticipated segmental duplications, display large phenotypic diversity in a wide range of environmental conditions, showing that simply reconfiguring chromosome architecture is sufficient to provide fitness advantages in stressful growth conditions. [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]

Studying bacterial population diversity is important to understand healthcare associated infections’ epidemiology and has a significant impact on dealing with multidrug resistant bacterial outbreaks. We characterised the extended-spectrum beta-lactamase producing K. pneumoniae (ESBLp KPN) population in our hospital using mini-MLST. Then we used whole genome sequencing (WGS) to compare selected isolates belonging to the most prevalent melting types (MelTs) and the colonization/infection pair isolates collected from one patient to study the ESBLp KPN population’s genetic diversity. A total of 922 ESBLp KPN isolates collected between 7/2016 and 5/2018 were divided into 38 MelTs using mini-MLST with only 6 MelTs forming 82.8% of all isolates. For WGS, 14 isolates from the most prominent MelTs collected in the monitored period and 10 isolates belonging to the same MelTs collected in our hospital in 2014 were randomly selected. Resistome, virulome and ST were MelT specific and stable over time. A maximum of 23 SNV per core genome and 58 SNV per core and accessory genome were found. To determine the SNV relatedness cut-off values, 22 isolates representing colonization/infection pair samples obtained from 11 different patients were analysed by WGS with a maximum of 22 SNV in the core genome and 40 SNV in the core and accessory genome within pairs. The mini-MLST showed its potential for real-time epidemiology in clinical practice. However, for outbreak evaluation in a low diversity bacterial population, mini-MLST should be combined with more sensitive methods like WGS. Our findings showed there were only minimal differences within the core and accessory genome in the low diversity hospital population and gene based SNV analysis does not have enough discriminatory power to differentiate isolate relatedness. Thus, intergenic regions and mobile elements should be incorporated into the analysis scheme to increase discriminatory power. [ABSTRACT FROM AUTHOR]

Diet has been shown to influence epigenetic key players, such as DNA methylation, which can regulate the gene expression potential in both parents and offspring. Diets enriched in omega-6 and deficient in omega-3 PUFAs (low dietary omega-3/omega-6 PUFA ratio), have been associated with the promotion of pathogenesis of diseases in humans and other mammals. In this study, we investigated the impact of increased dietary intake of arachidonic acid (ARA), a physiologically important omega-6 PUFA, on 2 generations of zebrafish. Parental fish were fed either a low or a high ARA diet, while the progeny of both groups were fed the low ARA diet. We screened for DNA methylation on single base-pair resolution using reduced representation bisulfite sequencing (RRBS). The DNA methylation profiling revealed significant differences between the dietary groups in both parents and offspring. The majority of differentially methylated loci associated with high dietary ARA were found in introns and intergenic regions for both generations. Common loci between the identified differentially methylated loci in F0 and F1 livers were reported. We described overlapping gene annotations of identified methylation changes with differential expression, but based on a small number of overlaps. The present study describes the diet-associated methylation profiles across genomic regions, and it demonstrates that parental high dietary ARA modulates DNA methylation patterns in zebrafish liver. [ABSTRACT FROM AUTHOR]

Background: Schistosoma japonicum is a parasitic flatworm that causes human schistosomiasis, which is a significant cause of morbidity in China and the Philippines. A single draft genome was available for S. japonicum, yet this assembly is very fragmented and only covers 90% of the genome, which make it difficult to be applied as a reference in functional genome analysis and genes discovery. Findings: In this study, we present a high-quality assembly of the fluke S. japonicum genome by combining 20 G (~53X) long single molecule real time sequencing reads with 80 G (~ 213X) Illumina paired-end reads. This improved genome assembly is approximately 370.5 Mb, with contig and scaffold N50 length of 871.9 kb and 1.09 Mb, representing 142.4-fold and 6.2-fold improvement over the released WGS-based assembly, respectively. Additionally, our assembly captured 85.2% complete and 4.6% partial eukaryotic Benchmarking Universal Single-Copy Orthologs. Repetitive elements account for 46.80% of the genome, and 10,089 of the protein-coding genes were predicted from the improved genome, of which 96.5% have been functionally annotated. Lastly, using the improved assembly, we identified 20 significantly expanded gene families in S. japonicum, and those genes were primarily enriched in functions of proteolysis and protein glycosylation. Conclusions: Using the combination of PacBio and Illumina Sequencing technologies, we provided an improved high-quality genome of S. japonicum. This improved genome assembly, as well as the annotation, will be useful for the comparative genomics of the flukes and more importantly facilitate the molecular studies of this important parasite in the future. [ABSTRACT FROM AUTHOR]

The marker-trait association for complex traits using genotyping by sequencing (GBS) method is being widely spread in plants. The study aimed to identify significant single nucleotide polymorphism (SNP) associations for rachis length (RL), leaf area (LA) and total dry weight (TrDW) in oil palm among diverse African germplasm. The Illumina NextSeq platform has been used for SNP genotyping and retained 4031 fully informative SNPs after applying the filter criterion. These 4031 SNPs were used for genome wide association study for the above three traits. The LD decay rates of the African germplasm using GBS data of SNP is observed to be 25 Kb at 0.45 of average pair wise correlation coefficient (r2). Association mapping led to the identification of seven significant associations for three traits using MLM approach at a P value of ≤ 0.001. Three associations were identified for total dry weight, two each for leaf area index and rachis length. The qtlLA1 was found to be highly significant at a P value of 7.39E-05 (18.4% phenotypic variance) which is located on chromosome 4. Two QTLs (qtlLA2 and qtlRL1) were located on chromosome 1, which explained 11.9% and 12.4% of phenotypic variance respectively. Three QTLs for total dry weight were located on chromosome 2, 14 and 16, all-together explained 40% phenotypic variance. The results showed that the SNP-trait associations identified in the present study could be used in selection of elite oil palm germplasm for higher yields. [ABSTRACT FROM AUTHOR]

Resolving the phylogenetic relationships of closely related species using a small set of loci is challenging as sufficient information may not be captured from a limited sample of the genome. Relying on few loci can also be problematic when conflict between gene-trees arises from incomplete lineage sorting and/or ongoing hybridization, problems especially likely in recently diverged lineages. Here, we developed a method using limited genomic resources that allows identification of many low copy candidate loci from across the nuclear and chloroplast genomes, design probes for target capture and sequence the captured loci. To validate our method we present data from Eucalyptus and Melaleuca, two large and phylogenetically problematic genera within the Myrtaceae family. With one annotated genome, one transcriptome and two whole-genome shotgun sequences of one Eucalyptus and four Melaleuca species, respectively, we identified 212 loci representing 263 kbp for targeted sequence capture and sequencing. Of these, 209 were successfully tested from 47 samples across five related genera of Myrtaceae. The average percentage of reads mapped back to the reference was 57.6% with coverage of more than 20 reads per position across 83.5% of the data. The methods developed here should be applicable across a large range of taxa across all kingdoms. The core methods are very flexible, providing a platform for various genomic resource availabilities and are useful from shallow to deep phylogenies. [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]

High quality reference genomes are an important resource in genomic research projects. A consequence is that DNA fragments carrying the reference allele will be more likely to map successfully, or receive higher quality scores. This reference bias can have effects on downstream population genomic analysis when heterozygous sites are falsely considered homozygous for the reference allele. In palaeogenomic studies of human populations, mapping against the human reference genome is used to identify endogenous human sequences. Ancient DNA studies usually operate with low sequencing coverages and fragmentation of DNA molecules causes a large proportion of the sequenced fragments to be shorter than 50 bp—reducing the amount of accepted mismatches, and increasing the probability of multiple matching sites in the genome. These ancient DNA specific properties are potentially exacerbating the impact of reference bias on downstream analyses, especially since most studies of ancient human populations use pseudo-haploid data, i.e. they randomly sample only one sequencing read per site. We show that reference bias is pervasive in published ancient DNA sequence data of prehistoric humans with some differences between individual genomic regions. We illustrate that the strength of reference bias is negatively correlated with fragment length. Most genomic regions we investigated show little to no mapping bias but even a small proportion of sites with bias can impact analyses of those particular loci or slightly skew genome-wide estimates. Therefore, reference bias has the potential to cause minor but significant differences in the results of downstream analyses such as population allele sharing, heterozygosity estimates and estimates of archaic ancestry. These spurious results highlight how important it is to be aware of these technical artifacts and that we need strategies to mitigate the effect. Therefore, we suggest some post-mapping filtering strategies to resolve reference bias which help to reduce its impact substantially. [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]

Clonal expansion of HIV infected cells plays an important role in the formation and persistence of the reservoir that allows the virus to persist, in DNA form, despite effective antiretroviral therapy. We used integration site analysis to ask if there is a similar clonal expansion of SIV infected cells in macaques. We show that the distribution of HIV and SIV integration sites in vitro is similar and that both viruses preferentially integrate in many of the same genes. We obtained approximately 8000 integration sites from blood samples taken from SIV-infected macaques prior to the initiation of ART, and from blood, spleen, and lymph node samples taken at necropsy. Seven clones were identified in the pre-ART samples; one persisted for a year on ART. An additional 100 clones were found only in on-ART samples; a number of these clones were found in more than one tissue. The timing and extent of clonal expansion of SIV-infected cells in macaques and HIV-infected cells in humans is quite similar. This suggests that SIV-infected macaques represent a useful model of the clonal expansion of HIV infected cells in humans that can be used to evaluate strategies intended to control or eradicate the viral reservoir. [ABSTRACT FROM AUTHOR]

The bertha armyworm (BAW), Mamestra configurata, is a significant pest of canola (Brassica napus L. and B. rapa L.) in western North America that undergoes cyclical outbreaks every 6–8 years. During peak outbreaks millions of dollars are spent on insecticidal control and, even with control efforts, subsequent damage can result in losses worth millions of dollars. Despite the importance of this pest insect, information is lacking on the dispersal ability of BAW and the genetic variation of populations from across its geographic range which may underlie potential differences in their susceptibility to insecticides or pathogens. Here, we examined the genetic diversity of BAW populations during an outbreak across its geographic range in western North America. First, mitochondrial cytochrome oxidase 1 (CO1) barcode sequences were used to confirm species identification of insects captured in a network of pheromone traps across the range, followed by haplotype analyses. We then sequenced the BAW genome and used double-digest restriction site associated DNA sequencing, mapped to the genome, to identify 1000s of single nucleotide polymorphisms (SNP) markers. CO1 haplotype analysis identified 9 haplotypes distributed across 28 sample locations and three laboratory-reared colonies. Analysis of genotypic data from both the CO1 and SNP markers revealed little population structure across BAW’s vast range. The CO1 haplotype pattern showed a star-like phylogeny which is often associated with species whose population abundance and range has recently expanded and combined with pheromone trap data, indicates the outbreak may have originated from a single focal point in central Saskatchewan. The relatively recent introduction of canola and rapid expansion of the canola growing region across western North America, combined with the cyclical outbreaks of BAW caused by precipitous population crashes, has likely selected for a genetically homogenous BAW population adapted to this crop. [ABSTRACT FROM AUTHOR]

Multi-locus sequencing typing (MLST) is widely used to monitor the phylogeny of microbial outbreaks. However, several strains of vancomycin-resistant Enterococcus faecium (VREfm) with a missing MLST locus (pstS) have recently emerged in Australia, with a few cases also reported in England. Here, we identified similarly distinct strains circulating in two neighbouring hospitals in Scotland. Whole genome sequencing of five VREfm strains isolated from these hospitals identified four pstS-null strains in both hospitals, while the fifth was multi-locus sequence type (ST) 262, which is the first documented in the UK. All five Scottish isolates had an insertion in the tetM gene, which is associated with increased susceptibility to tetracyclines, providing no other tetracycline-resistant gene is present. Such an insertion, which encompasses a dfrG gene and two currently uncharacterised genes, was additionally identified in all tested vanA-type pstS-null VREfm strains (5 English and 68 Australian). Phylogenetic comparison with other VREfm genomes indicates that the four pstS-null Scottish isolates sequenced in this study are more closely related to pstS-null strains from Australia rather than the English pstS-null isolates. Given how rapidly such pstS-null strains have expanded in Australia, the emergence of this clone in Scotland raises concerns for a potential outbreak. [ABSTRACT FROM AUTHOR]

Precise promoter annotation is required for understanding the mechanistic basis of transcription initiation. In the context of complex genomes, such as herpesviruses where there is extensive genic overlap, identification of transcription start sites (TSSs) is particularly problematic and cannot be comprehensively accessed by standard RNA sequencing approaches. Kaposi's sarcoma-associated herpesvirus (KSHV) is an oncogenic gammaherpesvirus and the etiological agent of Kaposi’s sarcoma and the B cell lymphoma primary effusion lymphoma (PEL). Here, we leverage RNA annotation and mapping of promoters for analysis of gene expression (RAMPAGE) and define KSHV TSSs transcriptome-wide and at nucleotide resolution in two widely used models of KSHV infection, namely iSLK.219 cells and the PEL cell line TREx-BCBL1-RTA. By mapping TSSs over a 96 h time course of reactivation we confirm 48 of 50 previously identified TSSs. Moreover, we identify over 100 novel transcription start site clusters (TSCs) in each cell line. Our analyses identified cell-type specific differences in TSC positions as well as promoter strength, and defined motifs within viral core promoters. Collectively, by defining TSSs at high resolution we have greatly expanded the transcriptional landscape of the KSHV genome and identified transcriptional control mechanisms at play during KSHV lytic reactivation. [ABSTRACT FROM AUTHOR]

Background: Whole-genome sequencing (WGS) of heart transplant recipient- and donor-derived cardiac biopsies may facilitate organ matching, graft failure prediction, and immunotolerance research. The objective of this study was to determine the feasibility of WGS based on formalin-fixed paraffin-embedded endomyocardial biopsies. Methods and results: The study included serial donor- and recipient samples from patients who had undergone heart transplantation at Skane University Hospital, Lund, Sweden, between 1988 and 2009. DNA extraction and WGS were conducted. Additional WGS sequencing quality metrics and coverage were obtained with the Genome Analysis Toolkit (GATK). 455 endomyocardial samples from 37 heart transplant recipients were acquired from routine rejection monitoring and stored as formalin-fixed paraffin-embedded samples. They were analyzed after 3–26 years of storage. DNA was extracted from 114 samples and WGS was run on 85 samples. DNA extraction yielded 313 ng (IQR 96–601) for all samples. A coverage of 11.3x (IQR 9.0–15.9) was recorded for all WGS samples. Three samples stored for > 25 years yielded a coverage of > 25x. Data were generated for 1.7 billion reads per sample (IQR 1.4–2.7). A Transition/Transversion (TiTv) ratio of 2.09 ± 0.05 was calculated for all WGS samples. No associations were found among storage time, DNA yield, or sequencing quality metrics. Conclusions: The present study demonstrated the feasibility of whole-genome sequencing based on endomyocardial biopsies. This process could enable large-scale retrospective genomic studies using stored histopathological samples. [ABSTRACT FROM AUTHOR]

Malaria parasites possess the remarkable ability to maintain chronic infections that fail to elicit a protective immune response, characteristics that have stymied vaccine development and cause people living in endemic regions to remain at risk of malaria despite previous exposure to the disease. These traits stem from the tremendous antigenic diversity displayed by parasites circulating in the field. For Plasmodium falciparum, the most virulent of the human malaria parasites, this diversity is exemplified by the variant gene family called var, which encodes the major surface antigen displayed on infected red blood cells (RBCs). This gene family exhibits virtually limitless diversity when var gene repertoires from different parasite isolates are compared. Previous studies indicated that this remarkable genome plasticity results from extensive ectopic recombination between var genes during mitotic replication; however, the molecular mechanisms that direct this process to antigen-encoding loci while the rest of the genome remains relatively stable were not determined. Using targeted DNA double-strand breaks (DSBs) and long-read whole-genome sequencing, we show that a single break within an antigen-encoding region of the genome can result in a cascade of recombination events leading to the generation of multiple chimeric var genes, a process that can greatly accelerate the generation of diversity within this family. We also found that recombinations did not occur randomly, but rather high-probability, specific recombination products were observed repeatedly. These results provide a molecular basis for previously described structured rearrangements that drive diversification of this highly polymorphic gene family. [ABSTRACT FROM AUTHOR]

Virus discovery based on high-throughput sequencing relies on enrichment for virus sequences prior to library preparation to achieve a sufficient number of viral reads. In general, preparations of double-stranded RNA or total RNA preparations treated to remove rRNA are used for sequence enrichment. We used virus-specific antibodies to immunocapture virions from plant sap to conduct cDNA synthesis, followed by library preparation and HTS. For the four potato viruses PLRV, PVY, PVA and PYV, template preparation by virion immunocapture provided a simpler and less expensive method than the enrichment of total RNA by ribosomal depletion. Specific enrichment of viral sequences without an intermediate amplification step was achieved, and this high coverage of sequences across the viral genomes was important to identify rare sequence variations. Using this approach, the first complete genome sequence of a potato yellowing virus isolate (PYV, DSMZ PV-0706) was determined in this study. PYV can be confidently assigned as a distinct species in the genus Ilarvirus. [ABSTRACT FROM AUTHOR]

Pathogenic yeasts and fungi are an increasing global healthcare burden, but discovery of novel antifungal agents is slow. The mycoparasitic yeast Saccharomycopsis schoenii was recently demonstrated to be able to kill the emerging multi-drug resistant yeast pathogen Candida auris. However, the molecular mechanisms involved in the predatory activity of S. schoenii have not been explored. To this end, we de novo sequenced, assembled and annotated a draft genome of S. schoenii. Using proteomics, we confirmed that Saccharomycopsis yeasts have reassigned the CTG codon and translate CTG into serine instead of leucine. Further, we confirmed an absence of all genes from the sulfate assimilation pathway in the genome of S. schoenii, and detected the expansion of several gene families, including aspartic proteases. Using Saccharomyces cerevisiae as a model prey cell, we honed in on the timing and nutritional conditions under which S. schoenii kills prey cells. We found that a general nutrition limitation, not a specific methionine deficiency, triggered predatory activity. Nevertheless, by means of genome-wide transcriptome analysis we observed dramatic responses to methionine deprivation, which were alleviated when S. cerevisiae was available as prey, and therefore postulate that S. schoenii acquired methionine from its prey cells. During predation, both proteomic and transcriptomic analyses revealed that S. schoenii highly upregulated and translated aspartic protease genes, probably used to break down prey cell walls. With these fundamental insights into the predatory behavior of S. schoenii, we open up for further exploitation of this yeast as a biocontrol yeast and/or source for novel antifungal agents. [ABSTRACT FROM AUTHOR]

Formalin-fixed paraffin-embedded (FFPE) tissues are among the most widely available clinical specimens. Their potential utility as a source of RNA for transcriptome studies would greatly enhance population-based cancer studies. Although preliminary studies suggest FFPE tissue may be used for RNA sequencing, the effect of storage time on these specimens needs to be determined. We conducted this study to determine whether RNA in archived FFPE high-grade ovarian serous adenocarcinomas from Surveillance, Epidemiology and End Results (SEER) registries was present in sufficient quantity and quality for RNA-Seq analysis. FFPE tissues, stored from 7 to 32 years, were obtained from three SEER sites. RNA was extracted, quantified, quality assessed, and subjected to RNA-Seq (a whole transcriptome sequencing technology). FFPE specimens stored for longer periods of time had poorer RNA sample quality as indicated by negative correlations between specimen storage time and fragment distribution values (DV). In addition, sample contamination was a common issue among the RNA, with 41 of 67 samples having 5% to 48% bacterial contamination. However, regardless of specimen storage time and bacterial contamination, 60% of the samples yielded data that enabled gene expression quantification, identifying more than 10,000 genes, with the correlations among most biological replicates above 0.7. This study demonstrates that FFPE high-grade ovarian serous adenocarcinomas specimens stored in repositories for up to 32 years and under varying storage conditions are a promising source of RNA for RNA-Seq. We also describe certain caveats to be considered when designing RNA-Seq studies using archived FFPE tissues. [ABSTRACT FROM AUTHOR]