Explore the vast range of titles available.

In non-mammalian vertebrates, the molecular mechanisms involved in the transformation of haploid germ cells (HGCs) into spermatozoa (spermiogenesis) are largely unknown. Here, we investigated this process in the marine teleost gilthead seabream (Sparus aurata) through the examination of the changes in the transcriptome between cell-sorted HGCs and ejaculated sperm (SPZEJ). Samples were collected under strict quality controls employing immunofluorescence microscopy as well as by determining the sperm motion kinematic parameters by computer-assisted sperm analysis. Deep sequencing by RNA-seq identified a total of 7286 differentially expressed genes (DEGs) (p-value < 0.01) between both cell types, of which nearly half were upregulated in SPZEJ compared to HCGs. In addition, approximately 9000 long non-coding RNAs (lncRNAs) were found, of which 56% were accumulated or emerged de novo in SPZEJ. The upregulated transcripts are involved in transcriptional and translational regulation, chromatin and cytoskeleton organization, metabolic processes such as glycolysis and oxidative phosphorylation, and also include a number of ion and water channels, exchangers, transporters and receptors. Pathway analysis conducted on DEGs identified 37 different signaling pathways enriched in SPZEJ, including 13 receptor pathways, from which the most predominant correspond to the chemokine and cytokine, gonadotropin-releasing hormone receptor and platelet derived growth factor signaling pathways. Our data provide new insight into the mRNA and lncRNA cargos of teleost spermatozoa and uncover the possible involvement of novel endocrine mechanisms during the differentiation and maturation of spermatozoa.

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research. Accurate analysis of mitochondrial DNA is important for mitochondrial disease clinical research and diagnostics. Here, authors present a method using Cas9 cleavage, nanopore sequencing and a custom pipeline to identify pathogenic variants, deletions and accurately quantify heteroplasmy to below 1%.

The molecular classification of glioblastoma (GBM) based on gene expression might better explain outcome and response to treatment than clinical factors. Whole transcriptome sequencing using next-generation sequencing platforms is rapidly becoming accepted as a tool for measuring gene expression for both research and clinical use. Fresh frozen (FF) tissue specimens of GBM are difficult to obtain since tumor tissue obtained at surgery is often scarce and necrotic and diagnosis is prioritized over freezing. After diagnosis, leftover tissue is usually stored as formalin-fixed paraffin-embedded (FFPE) tissue. However, RNA from FFPE tissues is usually degraded, which could hamper gene expression analysis. We compared RNA-Seq data obtained from matched pairs of FF and FFPE GBM specimens. Only three FFPE out of eleven FFPE-FF matched samples yielded informative results. Several quality-control measurements showed that RNA from FFPE samples was highly degraded but maintained transcriptomic similarities to RNA from FF samples. Certain issues regarding mutation analysis and subtype prediction were detected. Nevertheless, our results suggest that RNA-Seq of FFPE GBM specimens provides reliable gene expression data that can be used in molecular studies of GBM if the RNA is sufficiently preserved.

Background Long-read sequencing technologies enable resolution of structural variants (SV) and long-range genome assembly, but require high molecular weight (HMW) DNA of both high quantity and quality to produce optimal sequencing results. New DNA extraction methods have been developed but these have not been assessed for use in routine testing. The interlaboratory study described here tested four commonly used methods: Fire Monkey, Nanobind, Puregene and Genomic-tip with a reference cell line containing known chromosomal alterations. Samples were assessed with commonly applied approaches for evaluating DNA purity and integrity as well as a method based on linkage using digital PCR. Sequencing performance was evaluated and the impact of extraction method on structural variant calling investigated. Results All methods generally produced samples of acceptable purity although yield varied considerably between laboratories. Library preparation and sequencing were successful for all four methods, with Fire Monkey extracts achieving the highest N50 values, Genomic Tip giving the highest sequencing yields and Nanobind, the highest proportion of ultra-long reads (> 100 kb). The dPCR assay with duplexes at 100 kb and 150 kb distances was predictive of ultra-long reads and provides a more quantitative read-out (% linkage) than pulse-field gel electrophoresis (PFGE) which varied in performance between instruments and gel dyes. Neither PFGE nor dPCR were predictive of the proportion of short reads (< 10 kb). Coverage was a key factor in the success of SV calling, but this was dependent on SV caller. Megabase scale SVs were challenging to analyse with SV callers and required confirmation based on coverage plots and mapping of junction sequences, and the findings of earlier studies were only partially confirmed. Conclusions This study highlights some of the challenges of HMW DNA extraction as well as the need for robust sample QC metrics to ensure optimal sequencing yield and read length which in turn influence the success of SV analysis. dPCR approaches for DNA integrity showed potential but require further development. As long-read methods are increasingly applied in routine settings such as clinical testing laboratories, cellular reference samples with well-characterised SVs are recommended as controls for the full long-read sequencing workflow.

We apply a known algorithm for computing exactly inequalities between Beta distributions to assess whether a given position in a genome is differentially methylated across samples. We discuss the advantages brought by the adoption of this solution with respect to two approximations (Fisher's test and Z score). The same formalism presented here can be applied in a similar way to variant calling.

Geminiviruses are plant ssDNA viruses that replicate through dsDNA intermediates and form minichromosomes which carry the same epigenetic marks as the host chromatin. During the infection, geminiviruses are targets of the post-transcriptional and transcriptional gene silencing machinery. To obtain insights into the connection between virus-derived small RNAs (vsRNAs), viral genome methylation and gene expression, we obtained the transcriptome, sRNAome and methylome from the geminivirus Tomato yellow leaf curl virus -infected tomato plants. The results showed accumulation of transcripts just at the viral ORFs, while vsRNAs spanned the entire genome, showing a prevalent accumulation at regions where the viral ORFs overlapped. The viral genome was not homogenously methylated showing two highly methylated regions located in the C1 ORF and around the intergenic region (IR). The compilation of those results showed a partial correlation between vsRNA accumulation, gene expression and DNA methylation. We could distinguish different epigenetic scenarios along the viral genome, suggesting that in addition to its function as a plant defence mechanism, DNA methylation could have a role in viral gene regulation. To our knowledge, this is the first report that shows integrative single-nucleotide maps of DNA methylation, vsRNA accumulation and gene expression from a plant virus.

Understanding the adaptive capacity to current climate change of drought-sensitive tree species is mandatory, given their limited prospect of migration and adaptation as long-lived, sessile organisms. Knowledge about the molecular and eco-physiological mechanisms that control drought resilience is thus key, since water shortage appears as one of the main abiotic factors threatening forests ecosystems. However, our current background is scarce, especially in conifers, due to their huge and complex genomes. Here we investigated the eco-physiological and transcriptomic basis of drought response of the climate change-threatened conifer . We studied seedlings from two locations with contrasting drought conditions to investigate a local adaptation. Seedlings were subjected to experimental drought conditions, and were monitored at immediate (24 hours) and extended (20 days) times. In addition, post-drought recovery was investigated, depicting two contrasting responses in both locations (drought resilient and non-resilient). Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of drought resilience and investigate a rapid local adaptation of . transcriptome assembly was performed for the first time in this species, providing differences in gene expression between the immediate and extended treatments, as well as among the post-drought recovery phenotypes. Weighted gene co-expression network analysis showed a regulation of stomatal closing and photosynthetic activity during the immediate drought, consistent with an isohydric dynamic. During the extended drought, growth and flavonoid biosynthesis inhibition mechanisms prevailed, probably to increase root-to-shoot ratio and to limit the energy-intensive biosynthesis of secondary metabolites. Drought sensitive individuals failed in metabolism and photosynthesis regulation under drought stress, and in limiting secondary metabolite production. Moreover, genomic differences (SNPs) were found between drought resilient and sensitive seedlings, and between the two studied locations, which were mostly related to transposable elements. This work provides novel insights into the transcriptomic basis of drought response of , a set of candidate genes mechanistically involved in its drought sensitivity and evidence of a rapid local adaptation. Our results may help guide conservation programs for this threatened conifer, contribute to advance drought-resilience research and shed light on trees' adaptive potential to current climate change.

RNA-Sequencing (RNA-Seq) can identify gene fusions in tumors, but not all these fusions have functional consequences. Using multiple data bases, we have performed an in silico analysis of fusions detected by RNA-Seq in tumor samples from 139 newly diagnosed glioblastoma patients to identify in-frame fusions with predictable oncogenic potential. Among 61 samples with fusions, there were 103 different fusions, involving 167 different genes, including 20 known oncogenes or tumor suppressor genes (TSGs), 16 associated with cancer but not oncogenes or TSGs, and 32 not associated with cancer but previously shown to be involved in fusions in gliomas. After selecting in-frame fusions able to produce a protein product and running Oncofuse, we identified 30 fusions with predictable oncogenic potential and classified them into four non-overlapping categories: six previously described in cancer; six involving an oncogene or TSG; four predicted by Oncofuse to have oncogenic potential; and 14 other in-frame fusions. Only 24 patients harbored one or more of these 30 fusions, and only two fusions were present in more than one patient: FGFR3::TACC3 and EGFR::SEPTIN14 . This in silico study provides a good starting point for the identification of gene fusions with functional consequences in the pathogenesis or treatment of glioblastoma.

Rift Valley fever phlebovirus (RVFV) causes an emerging zoonotic disease and is mainly transmitted by Culex and Aedes mosquitoes. While Aedes aegypti -dengue virus (DENV) is the most studied model, less is known about the genes involved in infection-responses in other mosquito-arboviruses pairing. The main objective was to investigate the molecular responses of Cx . pipiens to RVFV exposure focusing mainly on genes implicated in innate immune responses. Mosquitoes were fed with blood spiked with RVFV. The fully-engorged females were pooled at 3 different time points: 2 hours post-exposure (hpe), 3- and 14-days post-exposure (dpe). Pools of mosquitoes fed with non-infected blood were also collected for comparisons. Total RNA from each mosquito pool was subjected to RNA-seq analysis and a de novo transcriptome was constructed. A total of 451 differentially expressed genes (DEG) were identified. Most of the transcriptomic alterations were found at an early infection stage after RVFV exposure. Forty-eight DEG related to immune infection-response were characterized. Most of them were related with the RNAi system, Toll and IMD pathways, ubiquitination pathway and apoptosis. Our findings provide for the first time a comprehensive view on Cx . pipiens -RVFV interactions at the molecular level. The early depletion of RNAi pathway genes at the onset of the RVFV infection would allow viral replication in mosquitoes. While genes from the Toll and IMD immune pathways were altered in response to RVFV none of the DEG were related to the JAK/STAT pathway. The fact that most of the DEG involved in the Ubiquitin-proteasome pathway (UPP) or apoptosis were found at an early stage of infection would suggest that apoptosis plays a regulatory role in infected Cx . pipiens midguts. This study provides a number of target genes that could be used to identify new molecular targets for vector control.

Recent advances in long-read sequencing technologies have allowed the generation and curation of more complete genome assemblies, enabling the analysis of traditionally neglected chromosomes, such as the human Y chromosome (chrY). Native DNA was sequenced on a MinION Oxford Nanopore Technologies sequencing device to generate genome assemblies for seven major chrY human haplogroups. We analyzed and compared the chrY enrichment of sequencing data obtained using two different selective sequencing approaches: adaptive sampling and flow cytometry chromosome sorting. We show that adaptive sampling can produce data to create assemblies comparable to chromosome sorting while being a less expensive and time-consuming technique. We also assessed haplogroup-specific structural variants, which would be otherwise difficult to study using short-read sequencing data only. Finally, we took advantage of this technology to detect and profile epigenetic modifications among the considered haplogroups. Altogether, we provide a framework to study complex genomic regions with a simple, fast, and affordable methodology that could be applied to larger population genomics datasets. This study generates genome assemblies for seven major haplogroups of the human Y chromosome at a unique resolution that overcomes the structurally complex regions within. These data enable the detection and profiling of epigenetic modifications among the considered haplogroups.