Identification of LINE retrotransposons and long non-coding RNAs expressed in the octopus brain.

BACKGROUND: Transposable elements (TEs) widely contribute to the evolution of genomes allowing genomic innovations, generating germinal and somatic heterogeneity, and giving birth to long non-coding RNAs (lncRNAs). These features have been associated to the evolution, functioning, and complexity of the nervous system at such a level that somatic retrotransposition of long interspersed element (LINE) L1 has been proposed to be associated to human cognition. Among invertebrates, octopuses are fascinating animals whose nervous system reaches a high level of complexity achieving sophisticated cognitive abilities. The sequencing of the genome of the Octopus bimaculoides revealed a striking expansion of TEs which were proposed to have contributed to the evolution of its complex nervous system. We recently found a similar expansion also in the genome of Octopus vulgaris. However, a specific search for the existence and the transcription of full-length transpositionally competent TEs has not been performed in this genus. RESULTS: Here, we report the identification of LINE elements competent for retrotransposition in Octopus vulgaris and Octopus bimaculoides and show evidence suggesting that they might be transcribed and determine germline and somatic polymorphisms especially in the brain. Transcription and translation measured for one of these elements resulted in specific signals in neurons belonging to areas associated with behavioral plasticity. We also report the transcription of thousands of lncRNAs and the pervasive inclusion of TE fragments in the transcriptomes of both Octopus species, further testifying the crucial activity of TEs in the evolution of the octopus genomes. CONCLUSIONS: The neural transcriptome of the octopus shows the transcription of thousands of putative lncRNAs and of a full-length LINE element belonging to the RTE class. We speculate that a convergent evolutionary process involving retrotransposons activity in the brain has been important for the evolution of sophisticated cognitive abilities in this genus.

Therapeutic homology-independent targeted integration in retina and liver.

Challenges to the widespread application of gene therapy with adeno-associated viral (AAV) vectors include dominant conditions due to gain-of-function mutations which require allele-specific knockout, as well as long-term transgene expression from proliferating tissues, which is hampered by AAV DNA episomal status. To overcome these challenges, we used CRISPR/Cas9-mediated homology-independent targeted integration (HITI) in retina and liver as paradigmatic target tissues. We show that AAV-HITI targets photoreceptors of both mouse and pig retina, and this results in significant improvements to retinal morphology and function in mice with autosomal dominant retinitis pigmentosa. In addition, we show that neonatal systemic AAV-HITI delivery achieves stable liver transgene expression and phenotypic improvement in a mouse model of a severe lysosomal storage disease. We also show that HITI applications predominantly result in on-target editing. These results lay the groundwork for the application of AAV-HITI for the treatment of diseases affecting various organs.

Cantù syndrome: Report of a patient with a novel variant in KCNJ8 and revision of literature.

Cantù syndrome (CS) is a rare multisystemic disorder, characterized by congenital hypertrichosis, macrocephaly, facial dysmorphisms, cardiomegaly, vascular, and skeletal anomalies. From the cognitive point of view, most of the patients show a mild speech delay and a few of them present intellectual disability and learning difficulties. To date, most CS-reported cases are caused by heterozygous ABCC9 gene mutations. Only three patients with CS and heterozygous KCNJ8 gene variants have been reported. The authors here present the fourth case of CS with a variant in KCNJ8 in a 6-month-old baby. Diagnosis was reached through Trio-Whole Exome analysis that revealed a de novo missense variant in KCNJ8.

VarGenius-HZD Allows Accurate Detection of Rare Homozygous or Hemizygous Deletions in Targeted Sequencing Leveraging Breadth of Coverage.

Homozygous deletions (HDs) may be the cause of rare diseases and cancer, and their discovery in targeted sequencing is a challenging task. Different tools have been developed to disentangle HD discovery but a sensitive caller is still lacking. We present VarGenius-HZD, a sensitive and scalable algorithm that leverages breadth-of-coverage for the detection of rare homozygous and hemizygous single-exon deletions (HDs). To assess its effectiveness, we detected both real and synthetic rare HDs in fifty exomes from the 1000 Genomes Project obtaining higher sensitivity in comparison with state-of-the-art algorithms that each missed at least one event. We then applied our tool on targeted sequencing data from patients with Inherited Retinal Dystrophies and solved five cases that still lacked a genetic diagnosis. We provide VarGenius-HZD either stand-alone or integrated within our recently developed software, enabling the automated selection of samples using the internal database. Hence, it could be extremely useful for both diagnostic and research purposes.

Linked-Read Whole Genome Sequencing Solves a Double DMD Gene Rearrangement.

Next generation sequencing (NGS) has changed our approach to diagnosis of genetic disorders. Nowadays, the most comprehensive application of NGS is whole genome sequencing (WGS) that is able to detect virtually all DNA variations. However, even after accurate WGS, many genetic conditions remain unsolved. This may be due to the current NGS protocols, based on DNA fragmentation and short reads. To overcome these limitations, we applied a linked-read sequencing technology that combines single-molecule barcoding with short-read WGS. We were able to assemble haplotypes and distinguish between alleles along the genome. As an exemplary case, we studied the case of a female carrier of X-linked muscular dystrophy with an unsolved genetic status. A deletion of exons 16-29 in DMD gene was responsible for the disease in her family, but she showed a normal dosage of these exons by Multiplex Ligation-dependent Probe Amplification (MLPA) and array CGH. This situation is usually considered compatible with a "non-carrier" status. Unexpectedly, the girl also showed an increased dosage of flanking exons 1-15 and 30-34. Using linked-read WGS, we were able to distinguish between the two X chromosomes. In the first allele, we found the 16-29 deletion, while the second allele showed a 1-34 duplication: in both cases, linked-read WGS correctly mapped the borders at single-nucleotide resolution. This duplication in trans apparently restored the normal dosage of exons 16-29 seen by quantitative assays. This had a dramatic impact in genetic counselling, by converting a non-carrier into a double carrier status prediction. We conclude that linked-read WGS should be considered as a valuable option to improve our understanding of unsolved genetic conditions.

Spectrum of Disease Severity in Patients With X-Linked Retinitis Pigmentosa Due to RPGR Mutations.

Purpose: The purpose of this study was to perform a detailed longitudinal phenotyping of X-linked retinitis pigmentosa (RP) caused by mutations in the RPGR gene during a long follow-up period. Methods: An Italian cohort of 48 male patients (from 31 unrelated families) with RPGR-associated RP was clinically assessed at a single center (mean follow-up = 6.5 years), including measurements of best-corrected visual acuity (BCVA), Goldmann visual field (GVF), optical coherence tomography (OCT), fundus autofluorescence (FAF), microperimetry, and full-field electroretinography (ERG). Results: Patients (29.6 ± 15.2 years) showed a mean BCVA of 0.6 ± 0.7 logMAR, mostly with myopic refraction (79.2%). Thirty patients (62.5%) presented a typical RP fundus, while the remaining sine pigmento RP. Over the follow-up, BCVA significantly declined at a mean rate of 0.025 logMAR/year. Typical RP and high myopia were associated with a significantly faster decline of BCVA. Blindness was driven primarily by GVF loss. ERG responses with a rod-cone pattern of dysfunction were detectable in patients (50%) that were significantly younger and more frequently presented sine pigmento RP. Thirteen patients (27.1%) had macular abnormalities without cystoid macular edema. Patients (50%) with a perimacular hyper-FAF ring were significantly younger, had a higher BCVA and a better-preserved ellipsoid zone band than those with markedly decreased FAF. Patients harboring pathogenic variants in exons 1 to 14 showed a milder phenotype compared to those with ORF15 mutations. Conclusions: Our monocentric, longitudinal retrospective study revealed a spectrum disease progression in male patients with RPGR-associated RP. Slow disease progression correlated with sine pigmento RP, absence of high myopia, and mutations in RPGR exons 1 to 14.

ACE2 gene variants may underlie interindividual variability and susceptibility to COVID-19 in the Italian population.

In December 2019, an initial cluster of interstitial bilateral pneumonia emerged in Wuhan, China. A human-to-human transmission was assumed and a previously unrecognized entity, termed coronavirus disease-19 (COVID-19) due to a novel coronavirus (SARS-CoV-2) was described. The infection has rapidly spread out all over the world and Italy has been the first European country experiencing the endemic wave with unexpected clinical severity in comparison with Asian countries. It has been shown that SARS-CoV-2 utilizes angiotensin converting enzyme 2 (ACE2) as host receptor and host proteases for cell surface binding and internalization. Thus, a predisposing genetic background can give reason for interindividual disease susceptibility and/or severity. Taking advantage of the Network of Italian Genomes (NIG), here we mined whole-exome sequencing data of 6930 Italian control individuals from five different centers looking for ACE2 variants. A number of variants with a potential impact on protein stability were identified. Among these, three more common missense changes, p.(Asn720Asp), p.(Lys26Arg), and p.(Gly211Arg) were predicted to interfere with protein structure and stabilization. Rare variants likely interfering with the internalization process, namely p.(Leu351Val) and p.(Pro389His), predicted to interfere with SARS-CoV-2 spike protein binding, were also observed. Comparison of ACE2 WES data between a cohort of 131 patients and 258 controls allowed identifying a statistically significant (P value < 0.029) higher allelic variability in controls compared with patients. These findings suggest that a predisposing genetic background may contribute to the observed interindividual clinical variability associated with COVID-19, allowing an evidence-based risk assessment leading to personalized preventive measures and therapeutic options.

Sinus pericranii, skull defects, and structural brain anomalies in TRAF7-related disorder.

BACKGROUND: Several somatic mutations in TRAF7 have been reported in cancers, whereas a few germline heterozygous mutations have been recently linked to a neurodevelopmental disorder, characterized by craniofacial dysmorphisms, congenital heart defects, and digital anomalies. CASES: We report two subjects harboring de novo heterozygous missense variants in TRAF7, namely the recurrent 1964G>A(p.Arg655Gln) and the novel missense c.1204C>G(p.Leu402Val) variants. In addition to the typical hallmarks of the TRAF7-related disorder, both subjects presented with a recognizable "pear-shaped" skull due to multiple craniosynostosis, sinus pericranii, skull base/cranio-cervical junction anomalies, dysgyria, and inferior cerebellar vermis hypoplasia. CONCLUSIONS: Hence, we expand the genotypic and phenotypic spectrum of this neurodevelopmental disorder, discussing possible implications for clinical management of subjects with germline TRAF7 mutations.

Congenital posterior cervical spine malformation due to biallelic c.240-4T>G RIPPLY2 variant: A discrete entity.

The clinical and radiological spectrum of spondylocostal dysostosis syndromes encompasses distinctive costo-vertebral anomalies. RIPPLY2 biallelic pathogenic variants were described in two distinct cervical spine malformation syndromes: Klippel-Feil syndrome and posterior cervical spine malformation. RIPPLY2 is involved in the determination of rostro-caudal polarity and somite patterning during development. To date, only four cases have been reported. The current report aims at further delineating the posterior malformation in three new patients. Three patients from two unrelated families underwent clinical and radiological examination through X-ray, 3D computed tomography and brain magnetic resonance imaging. After informed consent was obtained, family-based whole exome sequencing (WES) was performed. Complex vertebral segmentation defects in the cervico-thoracic spine were observed in all patients. WES led to the identification of the homozygous splicing variant c.240-4T>G in all subjects. This variant is predicted to result in aberrant splicing of Exon 4. The current report highlights a subtype of cervical spine malformation with major atlo-axoidal malformation compromising spinal cord integrity. This distinctive mutation-specific pattern of malformation differs from Klippel-Feil syndrome and broadens the current classification, defining a sub-type of RIPPLY2-related skeletal disorder. Of note, the phenotype of one patient overlaps with oculo-auriculo-vertebral spectrum disorder.

A draft genome sequence of the elusive giant squid, Architeuthis dux.

BACKGROUND: The giant squid (Architeuthis dux; Steenstrup, 1857) is an enigmatic giant mollusc with a circumglobal distribution in the deep ocean, except in the high Arctic and Antarctic waters. The elusiveness of the species makes it difficult to study. Thus, having a genome assembled for this deep-sea-dwelling species will allow several pending evolutionary questions to be unlocked. FINDINGS: We present a draft genome assembly that includes 200 Gb of Illumina reads, 4 Gb of Moleculo synthetic long reads, and 108 Gb of Chicago libraries, with a final size matching the estimated genome size of 2.7 Gb, and a scaffold N50 of 4.8 Mb. We also present an alternative assembly including 27 Gb raw reads generated using the Pacific Biosciences platform. In addition, we sequenced the proteome of the same individual and RNA from 3 different tissue types from 3 other species of squid (Onychoteuthis banksii, Dosidicus gigas, and Sthenoteuthis oualaniensis) to assist genome annotation. We annotated 33,406 protein-coding genes supported by evidence, and the genome completeness estimated by BUSCO reached 92%. Repetitive regions cover 49.17% of the genome. CONCLUSIONS: This annotated draft genome of A. dux provides a critical resource to investigate the unique traits of this species, including its gigantism and key adaptations to deep-sea environments.

Clinical and Genetic Analysis of a European Cohort with Pericentral Retinitis Pigmentosa.

Retinitis pigmentosa (RP) is a clinically heterogenous disease that comprises a wide range of phenotypic and genetic subtypes. Pericentral RP is an atypical form of RP characterized by bone-spicule pigmentation and/or atrophy confined in the near mid-periphery of the retina. In contrast to classic RP, the far periphery is better preserved in pericentral RP. The aim of this study was to perform the first detailed clinical and genetic analysis of a cohort of European subjects with pericentral RP to determine the phenotypic features and the genetic bases of the disease. A total of 54 subjects from 48 independent families with pericentral RP, non-syndromic and syndromic, were evaluated through a full ophthalmological examination and underwent clinical exome or retinopathy gene panel sequencing. Disease-causative variants were identified in 22 of the 35 families (63%) in 10 different genes, four of which are also responsible for syndromic RP. Thirteen of the 34 likely pathogenic variants were novel. Intra-familiar variability was also observed. The current study confirms the mild phenotype of pericentral RP and extends the spectrum of genes associated with this condition.

Intraspecific Diversity in the Cold Stress Response of Transposable Elements in the Diatom Leptocylindrus aporus.

Transposable elements (TEs), activated as a response to unfavorable conditions, have been proposed to contribute to the generation of genetic and phenotypic diversity in diatoms. Here we explore the transcriptome of three warm water strains of the diatom Leptocylindrus aporus, and the possible involvement of TEs in their response to changing temperature conditions. At low temperature (13 °C) several stress response proteins were overexpressed, confirming low temperature to be unfavorable for L. aporus, while TE-related transcripts of the LTR retrotransposon superfamily were the most enriched transcripts. Their expression levels, as well as most of the stress-related proteins, were found to vary significantly among strains, and even within the same strains analysed at different times. The lack of overexpression after many months of culturing suggests a possible role of physiological plasticity in response to growth under controlled laboratory conditions. While further investigation on the possible central role of TEs in the diatom stress response is warranted, the strain-specific responses and possible role of in-culture evolution draw attention to the interplay between the high intraspecific variability and the physiological plasticity of diatoms, which can both contribute to the adaptation of a species to a wide range of conditions in the marine environment.

Novel CNS malformations and skeletal anomalies in a patient with Beaulieu-boycott-Innes syndrome.

THO/TREX (transcription/export) is a conserved eukaryotic complex that plays a crucial role in gene expression and prevents DNA damage during mitosis and meiosis. In mammals, TREX is essential during embryogenesis, determining stem cell fate specification by regulating posttranscriptional self-renewal and differentiation in several tissues. It is composed of a core called THO, consisting of THOC1, 2, 5, 6, 7, and additional proteins. Bi-allelic mutations in THOC6 have been associated to Beaulieu-Boycott-Innes syndrome (BBIS), a syndromic form of intellectual disability (ID). To date, nine patients harbouring homozygous or compound heterozygous mutations in THOC6 have been reported. Despite the clinical heterogenity and subtle dysmorphic features in some individuals, distinctive facial features are tall forehead, short and upslanting palpebral fissures, deep set eyes, flat philtrum, and malocclusion. Nonlife threatening congenital anomalies are common, including cardiac and renal malformations, anteriorly displaced anus, cryptorchidism in males, submucous cleft palate, and corpus callosum dysgenesis. Affected patients usually have short stature, mild microcephaly, and mild to moderate ID. Here, we describe an Italian patient with BBIS, carrying two compound heterozygous loss-of-function (LoF) variants in THOC6 (c.577C > T, p.R193* and c.792_793delCA, p.V264Vfs*48). In addition to the common phenotype, she displays cerebellar hypoplasia with severe vermian dysgenesis and hydrocephalus due to aqueductal stenosis, multiple skeletal anomalies and hypergonadotropic hypogonadism. Thus, we review the previous cases and discuss the phenotypic spectrum of BBIS, providing further evidence regarding the pivotal role of TREX complex in human development.

De novo assembly of a transcriptome from the eggs and early embryos of Astropecten aranciacus.

Starfish have been instrumental in many fields of biological and ecological research. Oocytes of Astropecten aranciacus, a common species native to the Mediterranean Sea and the East Atlantic, have long been used as an experimental model to study meiotic maturation, fertilization, intracellular Ca2+ signaling, and cell cycle controls. However, investigation of the underlying molecular mechanisms has often been hampered by the overall lack of DNA or protein sequences for the species. In this study, we have assembled a transcriptome for this species from the oocytes, eggs, zygotes, and early embryos, which are known to have the highest RNA sequence complexity. Annotation of the transcriptome identified over 32,000 transcripts including the ones that encode 13 distinct cyclins and as many cyclin-dependent kinases (CDK), as well as the expected components of intracellular Ca2+ signaling toolkit. Although the mRNAs of cyclin and CDK families did not undergo significant abundance changes through the stages from oocyte to early embryo, as judged by real-time PCR, the transcript encoding Mos, a negative regulator of mitotic cell cycle, was drastically reduced during the period of rapid cleavages. Molecular phylogenetic analysis using the homologous amino acid sequences of cytochrome oxidase subunit I from A. aranciacus and 30 other starfish species indicated that Paxillosida, to which A. aranciacus belongs, is not likely to be the most basal order in Asteroidea. Taken together, the first transcriptome we assembled in this species is expected to enable us to perform comparative studies and to design gene-specific molecular tools with which to tackle long-standing biological questions.

Genomewide transcriptional reprogramming in the seagrass Cymodocea nodosa under experimental ocean acidification.

Here, we report the first use of massive-scale RNA-sequencing to explore seagrass response to CO2 -driven ocean acidification (OA). Large-scale gene expression changes in the seagrass Cymodocea nodosa occurred at CO2 levels projected by the end of the century. C. nodosa transcriptome was obtained using Illumina RNA-Seq technology and de novo assembly, and differential gene expression was explored in plants exposed to short-term high CO2 /low pH conditions. At high pCO2 , there was a significant increased expression of transcripts associated with photosynthesis, including light reaction functions and CO2 fixation, and also to respiratory pathways, specifically for enzymes involved in glycolysis, in the tricarboxylic acid cycle and in the energy metabolism of the mitochondrial electron transport. The upregulation of respiratory metabolism is probably supported by the increased availability of photosynthates and increased energy demand for biosynthesis and stress-related processes under elevated CO2 and low pH. The upregulation of several chaperones resembling heat stress-induced changes in gene expression highlighted the positive role these proteins play in tolerance to intracellular acid stress in seagrasses. OA further modifies C. nodosa secondary metabolism inducing the transcription of enzymes related to biosynthesis of carbon-based secondary compounds, in particular the synthesis of polyphenols and isoprenoid compounds that have a variety of biological functions including plant defence. By demonstrating which physiological processes are most sensitive to OA, this research provides a major advance in the understanding of seagrass metabolism in the context of altered seawater chemistry from global climate change.

Marine diatoms change their gene expression profile when exposed to microscale turbulence under nutrient replete conditions.

Diatoms are a fundamental microalgal phylum that thrives in turbulent environments. Despite several experimental and numerical studies, if and how diatoms may profit from turbulence is still an open question. One of the leading arguments is that turbulence favours nutrient uptake. Morphological features, such as the absence of flagella, the presence of a rigid exoskeleton and the micrometre size would support the possible passive but beneficial role of turbulence on diatoms. We demonstrate that in fact diatoms actively respond to turbulence in non-limiting nutrient conditions. TURBOGEN, a prototypic instrument to generate natural levels of microscale turbulence, was used to expose diatoms to the mechanical stimulus. Differential expression analyses, coupled with microscopy inspections, enabled us to study the morphological and transcriptional response of Chaetoceros decipiens to turbulence. Our target species responds to turbulence by activating energy storage pathways like fatty acid biosynthesis and by modifying its cell chain spectrum. Two other ecologically important species were examined and the occurrence of a morphological response was confirmed. These results challenge the view of phytoplankton as unsophisticated passive organisms.

Transcriptome sequencing of three Pseudo-nitzschia species reveals comparable gene sets and the presence of Nitric Oxide Synthase genes in diatoms.

Diatoms are among the most diverse eukaryotic microorganisms on Earth, they are responsible for a large fraction of primary production in the oceans and can be found in different habitats. Pseudo-nitzschia are marine planktonic diatoms responsible for blooms in coastal and oceanic waters. We analyzed the transcriptome of three species, Pseudo-nitzschia arenysensis, Pseudo-nitzschia delicatissima and Pseudo-nitzschia multistriata, with different levels of genetic relatedness. These species have a worldwide distribution and the last one produces the neurotoxin domoic acid. We were able to annotate about 80% of the sequences in each transcriptome and the analysis of the relative functional annotations allowed comparison of the main metabolic pathways, pathways involved in the biosynthesis of isoprenoids (MAV and MEP pathways), and pathways putatively involved in domoic acid synthesis. The search for homologous transcripts among the target species and other congeneric species resulted in the discovery of a sequence annotated as Nitric Oxide Synthase (NOS), found uniquely in Pseudo-nitzschia multistriata. The predicted protein product contained all the domains of the canonical metazoan sequence. Putative NOS sequences were found in other available diatom datasets, supporting a role for nitric oxide as signaling molecule in this group of microalgae.

Annocript: a flexible pipeline for the annotation of transcriptomes able to identify putative long noncoding RNAs.

UNLABELLED: The eukaryotic transcriptome is composed of thousands of coding and long non-coding RNAs (lncRNAs). However, we lack a software platform to identify both RNA classes in a given transcriptome. Here we introduce Annocript, a pipeline that combines the annotation of protein coding transcripts with the prediction of putative lncRNAs in whole transcriptomes. It downloads and indexes the needed databases, runs the analysis and produces human readable and standard outputs together with summary statistics of the whole analysis. AVAILABILITY AND IMPLEMENTATION: Annocript is distributed under the GNU General Public License (version 3 or later) and is freely available at https://github.com/frankMusacchia/Annocript. CONTACT: remo.sanges@szn.it.