Cycles of heat exposure elevate metabolic enzyme genes and alters digestion in mussels

The intertidal sea mussel Mytilus californianus inhabits the Pacific coastline of North America. As a sessile organism it must cope with daily fluctuations of the marine and terrestrial environments. Organisms in stressful environments are commonly faced with energetic trade-offs between somatic and reproductive growth and stress management. Although, this energetic theory is generally accepted for mussels as well, the spectrum of mechanisms underlying this framework have not been widely investigated. In the current study we hypothesized that mussels acclimated to a cyclical moderately warm aerial environment would display enhanced transcript abundance of genes related to metabolism and exhibit resilient digestive enzyme activity (energy acquisition). Following acclimation to simulated tidal regimes in the laboratory we observed higher gene-expression of citrate synthase (CS), citrate lyase (ACLY), and mammalian target of rapamycin (MTOR) in heat stressed mussels. The expression of CS and MTOR was not elevated under acute thermal stress, suggestive that repeated stress is required for robust expression of these genes given that all other environmental variables are constant. We also observed reduced activity of the digestive enzyme, amylase in heat-shocked acclimated mussels (a proxy for energy acquisition). Our results suggest that mussels that settle high on shore not only face the challenge of thermal stress repair and limited access to food but may also be compromised by reduced digestive performance. Mussels may have adapted to cyclical energetic stress by overexpressing particular energy-related genes that can mitigate the disturbance to energy balance once the abundant transcripts are translated into functional proteins

Time elapsed after contrast injection is crucial to determine infarct transmurality and myocardial functional recovery after an acute myocardial infarction

BACKGROUND: In acute myocardial infarction (MI), late Gadolinium enhancement (LGE) has been proposed to include the infarcted myocardium and area at risk. However, little information is available on the optimal timing after contrast injection to differentiate these 2 areas. Our aim was to determine in acute and chronic MI whether imaging time after contrast injection influences the LGE size that better predicts infarct size and functional recovery. METHODS: Subjects were evaluated by cardiovascular magnetic resonance (CMR) the first week (n = 60) and 3 months (n = 47) after a percutaneously revascularized STEMI. Inversion-recovery single-shot (ss-IR) imaging was acquired at multiple time points following contrast administration and compared to segmented inversion-recovery (seg-IR) sequences. Inversion time was properly adjusted and images were blinded, randomized and measured for LGE volumes. RESULTS: In acute MI, LGE volume decreased over several minutes (p = 0.005) with the greatest volume occurring at 3 minutes and the smallest at 25 minutes post-contrast injection; however, LGE volume remained constant over time in chronic MI (p = 0.886). Depending on the imaging time, in acute phase, a change in the transmurality index was also observed. A transmural infarction (>75%) at 25 minutes better predicted the absence of improvement in the wall motion score index (WMSI), a higher increase in left ventricular volumes and a lower ejection fraction compared to 10 minutes. CONCLUSIONS: A change was observed in LGE volume in the minutes following contrast administration in acute but not in chronic MI. Infarct transmurality 25 minutes post-contrast injection better predicted infarct size and functional recovery at follow-up

HapSolo: an optimization approach for removing secondary haplotigs during diploid genome assembly and scaffolding.

BackgroundDespite marked recent improvements in long-read sequencing technology, the assembly of diploid genomes remains a difficult task. A major obstacle is distinguishing between alternative contigs that represent highly heterozygous regions. If primary and secondary contigs are not properly identified, the primary assembly will overrepresent both the size and complexity of the genome, which complicates downstream analysis such as scaffolding.ResultsHere we illustrate a new method, which we call HapSolo, that identifies secondary contigs and defines a primary assembly based on multiple pairwise contig alignment metrics. HapSolo evaluates candidate primary assemblies using BUSCO scores and then distinguishes among candidate assemblies using a cost function. The cost function can be defined by the user but by default considers the number of missing, duplicated and single BUSCO genes within the assembly. HapSolo performs hill climbing to minimize cost over thousands of candidate assemblies. We illustrate the performance of HapSolo on genome data from three species: the Chardonnay grape (Vitis vinifera), with a genome of 490 Mb, a mosquito (Anopheles funestus; 200 Mb) and the Thorny Skate (Amblyraja radiata; 2650 Mb).ConclusionsHapSolo rapidly identified candidate assemblies that yield improvements in assembly metrics, including decreased genome size and improved N50 scores. Contig N50 scores improved by 35%, 9% and 9% for Chardonnay, mosquito and the thorny skate, respectively, relative to unreduced primary assemblies. The benefits of HapSolo were amplified by down-stream analyses, which we illustrated by scaffolding with Hi-C data. We found, for example, that prior to the application of HapSolo, only 52% of the Chardonnay genome was captured in the largest 19 scaffolds, corresponding to the number of chromosomes. After the application of HapSolo, this value increased to ~ 84%. The improvements for the mosquito's largest three scaffolds, representing the number of chromosomes, were from 61 to 86%, and the improvement was even more pronounced for thorny skate. We compared the scaffolding results to assemblies that were based on PurgeDups for identifying secondary contigs, with generally superior results for HapSolo

HapSolo: an optimization approach for removing secondary haplotigs during diploid genome assembly and scaffolding.

BackgroundDespite marked recent improvements in long-read sequencing technology, the assembly of diploid genomes remains a difficult task. A major obstacle is distinguishing between alternative contigs that represent highly heterozygous regions. If primary and secondary contigs are not properly identified, the primary assembly will overrepresent both the size and complexity of the genome, which complicates downstream analysis such as scaffolding.ResultsHere we illustrate a new method, which we call HapSolo, that identifies secondary contigs and defines a primary assembly based on multiple pairwise contig alignment metrics. HapSolo evaluates candidate primary assemblies using BUSCO scores and then distinguishes among candidate assemblies using a cost function. The cost function can be defined by the user but by default considers the number of missing, duplicated and single BUSCO genes within the assembly. HapSolo performs hill climbing to minimize cost over thousands of candidate assemblies. We illustrate the performance of HapSolo on genome data from three species: the Chardonnay grape (Vitis vinifera), with a genome of 490 Mb, a mosquito (Anopheles funestus; 200 Mb) and the Thorny Skate (Amblyraja radiata; 2650 Mb).ConclusionsHapSolo rapidly identified candidate assemblies that yield improvements in assembly metrics, including decreased genome size and improved N50 scores. Contig N50 scores improved by 35%, 9% and 9% for Chardonnay, mosquito and the thorny skate, respectively, relative to unreduced primary assemblies. The benefits of HapSolo were amplified by down-stream analyses, which we illustrated by scaffolding with Hi-C data. We found, for example, that prior to the application of HapSolo, only 52% of the Chardonnay genome was captured in the largest 19 scaffolds, corresponding to the number of chromosomes. After the application of HapSolo, this value increased to ~ 84%. The improvements for the mosquito's largest three scaffolds, representing the number of chromosomes, were from 61 to 86%, and the improvement was even more pronounced for thorny skate. We compared the scaffolding results to assemblies that were based on PurgeDups for identifying secondary contigs, with generally superior results for HapSolo

SNP errors in maize

Whole genome assemblies of the NAM maize inbred lines were each aligned to the v5 reference genome using Anchorwave.  These were compared to SNPs identified from short-read data published with the original NAM assemblies. Short read SNPs were identified as erroneous if they overlapped with a deletion identified via whole genome assembly. The histogram shows the distribution of error rates across the NAM assemblies.</p

Rapid Low-Cost Assembly of the Drosophila melanogaster Reference Genome Using Low-Coverage, Long-Read Sequencing.

Accurate and comprehensive characterization of genetic variation is essential for deciphering the genetic basis of diseases and other phenotypes. A vast amount of genetic variation stems from large-scale sequence changes arising from the duplication, deletion, inversion, and translocation of sequences. In the past 10 years, high-throughput short reads have greatly expanded our ability to assay sequence variation due to single nucleotide polymorphisms. However, a recent de novo assembly of a second Drosophila melanogaster reference genome has revealed that short read genotyping methods miss hundreds of structural variants, including those affecting phenotypes. While genomes assembled using high-coverage long reads can achieve high levels of contiguity and completeness, concerns about cost, errors, and low yield have limited widespread adoption of such sequencing approaches. Here we resequenced the reference strain of D. melanogaster (ISO1) on a single Oxford Nanopore MinION flow cell run for 24 hr. Using only reads longer than 1 kb or with at least 30x coverage, we assembled a highly contiguous de novo genome. The addition of inexpensive paired reads and subsequent scaffolding using an optical map technology achieved an assembly with completeness and contiguity comparable to the D. melanogaster reference assembly. Comparison of our assembly to the reference assembly of ISO1 uncovered a number of structural variants (SVs), including novel LTR transposable element insertions and duplications affecting genes with developmental, behavioral, and metabolic functions. Collectively, these SVs provide a snapshot of the dynamics of genome evolution. Furthermore, our assembly and comparison to the D. melanogaster reference genome demonstrates that high-quality de novo assembly of reference genomes and comprehensive variant discovery using such assemblies are now possible by a single lab for under $1,000 (USD)

Supplemental Material for Solares et al., 2018

Supplemental figures and tables for "Rapid low-cost assembly of the <i>Drosophila melanogaster</i> reference genome using low-coverage, long-read sequencing"

Rapid Functional And Sequence Differentiation of a Tandemly-Repeated Species-Specific Multigene Family in Drosophila

Gene clusters of recently duplicated genes are hotbeds for evolutionary change. However, our understanding of how mutational mechanisms and evolutionary forces shape the structural and functional evolution of these clusters is hindered by the high sequence identity among the copies, which typically results in their inaccurate representation in genome assemblies. The presumed testis-specific, chimeric gene Sdic originated and tandemly expanded in Drosophila melanogaster, contributing to increased male-male competition. Using various types of massively parallel sequencing data, we studied the organization, sequence evolution, and functional attributes of the different Sdic copies. By leveraging long-read sequencing data, we uncovered both copy number and order differences from the currently accepted annotation for the Sdic region. Despite evidence for pervasive gene conversion affecting the Sdic copies, we also detected signatures of two episodes of diversifying selection, which have contributed to the evolution of a variety of C-termini and miRNA binding site compositions. Expression analyses involving RNA-seq datasets from 59 different biological conditions revealed distinctive expression breadths among the copies, with three copies being transcribed in females, opening the possibility to a sexually antagonistic effect. Phenotypic assays using Sdic knock-out strains indicated that should this antagonistic effect exist, it does not compromise female fertility. Our results strongly suggest that the genome consolidation of the Sdic gene cluster is more the result of a quick exploration of different paths of molecular tinkering by different copies than a mere dosage increase, which could be a recurrent evolutionary outcome in the presence of persistent sexual selection