A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds.

BACKGROUND: The ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map. RESULTS: Each of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features. CONCLUSIONS: The improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Limited Parallelism in Genetic Adaptation to Brackish Water Bodies in European Sprat and Atlantic Herring.

The European sprat is a small plankton-feeding clupeid present in the northeastern Atlantic Ocean, in the Mediterranean Sea, and in the brackish Baltic Sea and Black Sea. This species is the target of a major fishery and, therefore, an accurate characterization of its genetic population structure is crucial to delineate proper stock assessments that aid ensuring the fishery's sustainability. Here, we present (i) a draft genome assembly, (ii) pooled whole genome sequencing of 19 population samples covering most of the species' distribution range, and (iii) the design and test of a single nucleotide polymorphism (SNP)-chip resource and use this to validate the population structure inferred from pooled sequencing. These approaches revealed, using the populations sampled here, three major groups of European sprat: Oceanic, Coastal, and Brackish with limited differentiation within groups even over wide geographical stretches. Genetic structure is largely driven by six large putative inversions that differentiate Oceanic and Brackish sprats, while Coastal populations display intermediate frequencies of haplotypes at each locus. Interestingly, populations from the Baltic and the Black Seas share similar frequencies of haplotypes at these putative inversions despite their distant geographic location. The closely related clupeids European sprat and Atlantic herring both show genetic adaptation to the brackish Baltic Sea, providing an opportunity to explore the extent of genetic parallelism. This analysis revealed limited parallelism because out of 125 independent loci detected in the Atlantic herring, three showed sharp signals of selection that overlapped between the two species and contained single genes such as PRLRA, which encodes the receptor for prolactin, a freshwater-adapting hormone in euryhaline species, and THRB, a receptor for thyroid hormones, important both for metabolic regulation and the development of red cone photoreceptors.

Genomic, transcriptomic and epigenomic sequencing data of the B-cell leukemia cell line REH.

OBJECTIVES: The aim of this data paper is to describe a collection of 33 genomic, transcriptomic and epigenomic sequencing datasets of the B-cell acute lymphoblastic leukemia (ALL) cell line REH. REH is one of the most frequently used cell lines for functional studies of pediatric ALL, and these data provide a multi-faceted characterization of its molecular features. The datasets described herein, generated with short- and long-read sequencing technologies, can both provide insights into the complex aberrant karyotype of REH, and be used as reference datasets for sequencing data quality assessment or for methods development. DATA DESCRIPTION: This paper describes 33 datasets corresponding to 867 gigabases of raw sequencing data generated from the REH cell line. These datasets include five different approaches for whole genome sequencing (WGS) on four sequencing platforms, two RNA sequencing (RNA-seq) techniques on two different sequencing platforms, DNA methylation sequencing, and single-cell ATAC-sequencing.

The Carniolan Honeybee from Slovenia-A Complete and Annotated Mitochondrial Genome with Comparisons to Closely Related Apis mellifera Subspecies.

The complete mitochondrial genome of the Carniolan honeybee (Apis mellifera carnica) from Slovenia, a homeland of this subspecies, was acquired in two contigs from WGS data and annotated. The newly obtained mitochondrial genome is a circular closed loop of 16,447 bp. It comprises 37 genes (13 protein coding genes, 22 tRNA genes, and 2 rRNA genes) and an AT-rich control region. The order of the tRNA genes resembles the order characteristic of A. mellifera. The mitogenomic sequence of A. m. carnica from Slovenia contains 44 uniquely coded sites in comparison to the closely related subspecies A. m. ligustica and to A. m. carnica from Austria. Furthermore, 24 differences were recognised in comparison between A. m. carnica and A. m. ligustica subspecies. Among them, there are three SNPs that affect translation in the nd2, nd4, and cox2 genes, respectively. The phylogenetic placement of A. m. carnica from Slovenia within C lineage deviates from the expected position and changes the perspective on relationship between C and O lineages. The results of this study represent a valuable addition to the information available in the phylogenomic studies of A. mellifera-a pollinator species of worldwide importance. Such genomic information is essential for this local subspecies' conservation and preservation as well as its breeding and selection.

Duplication in the microtubule-actin cross-linking factor 1 gene causes a novel neuromuscular condition.

Spectrins and plakins are important communicators linking cytoskeletal components to each other and to cellular junctions. Microtubule-actin cross-linking factor 1 (MACF1) belongs to the spectraplakin family and is involved in control of microtubule dynamics. Complete knock out of MACF1 in mice is associated with developmental retardation and embryonic lethality. Here we present a family with a novel neuromuscular condition. Genetic analyses show a heterozygous duplication resulting in reduced MACF1 gene product. The functional consequence is affected motility observed as periodic hypotonia, lax muscles and diminished motor skills, with heterogeneous presentation among the affected family members. To corroborate these findings we used RNA interference to knock down the VAB-10 locus containing the MACF1 homologue in C. elegans, and we could show that this also causes movement disturbances. These findings suggest that changes in the MACF1 gene is implicated in this neuromuscular condition, which is an important observation since MACF1 has not previously been associated with any human disease and thus presents a key to understanding the essential nature of this gene.

Reduced ceramide synthase 2 activity causes progressive myoclonic epilepsy.

OBJECTIVE: Ceramides are precursors of complex sphingolipids (SLs), which are important for normal functioning of both the developing and mature brain. Altered SL levels have been associated with many neurodegenerative disorders, including epilepsy, although few direct links have been identified between genes involved in SL metabolism and epilepsy. METHODS: We used quantitative real-time PCR, Western blotting, and enzymatic assays to determine the mRNA, protein, and activity levels of ceramide synthase 2 (CERS2) in fiibroblasts isolated from parental control subjects and from a patient diagnosed with progressive myoclonic epilepsy (PME). Mass spectrometry and fluorescence microscopy were used to examine the effects of reduced CERS2 activity on cellular lipid composition and plasma membrane functions. RESULTS: We identify a novel 27 kb heterozygous deletion including the CERS2 gene in a proband diagnosed with PME. Compared to parental controls, levels of CERS2 mRNA, protein, and activity were reduced by ˜50% in fibroblasts isolated from this proband, resulting in significantly reduced levels of ceramides and sphingomyelins containing the very long-chain fatty acids C24:0 and C26:0. The change in SL composition was also reflected in a reduction in cholera toxin B immunofluorescence, indicating that membrane composition and function are altered. INTERPRETATION: We propose that reduced levels of CERS2, and consequently diminished levels of ceramides and SLs containing very long-chain fatty acids, lead to development of PME.

Functional loss of two ceramide synthases elicits autophagy-dependent lifespan extension in C. elegans.

Ceramide and its metabolites constitute a diverse group of lipids, which play important roles as structural entities of biological membranes as well as regulators of cellular growth, differentiation, and development. The C. elegans genome comprises three ceramide synthase genes; hyl-1, hyl-2, and lagr-1. HYL-1 function is required for synthesis of ceramides and sphingolipids containing very long acyl-chains (≥C24), while HYL-2 is required for synthesis of ceramides and sphingolipids containing shorter acyl-chains (≤C22). Here we show that functional loss of HYL-2 decreases lifespan, while loss of HYL-1 or LAGR-1 does not affect lifespan. We show that loss of HYL-1 and LAGR-1 functions extend lifespan in an autophagy-dependent manner, as knock down of the autophagy-associated gene ATG-12 abolishes hyl-1;lagr-1 longevity. The transcription factors PHA-4/FOXA, DAF-16/FOXO, and SKN-1 are also required for the observed lifespan extension, as well as the increased number of autophagosomes in hyl-1;lagr-1 animals. Both autophagic events and the transcription factors PHA-4/FOXA, DAF-16, and SKN-1 have previously been associated with dietary restriction-induced longevity. Accordingly, we find that hyl-1;lagr-1 animals display reduced feeding, increased resistance to heat, and reduced reproduction. Collectively, our data suggest that specific sphingolipids produced by different ceramide synthases have opposing roles in determination of C. elegans lifespan. We propose that loss of HYL-1 and LAGR-1 result in dietary restriction-induced autophagy and consequently prolonged longevity.