Long-read Sequencing Research Articles

CGC1, a new reference genome for Caenorhabditis elegans.

The original 100.3 Mb reference genome for Caenorhabditis elegans, generated from the wild-type laboratory strain N2, has been crucial for analysis of C. elegans since 1998 and has been considered complete since 2005. Unexpectedly, this long-standing reference was shown to be incomplete in 2019 by a genome assembly from the N2-derived strain VC2010. Moreover, genetically divergent versions of N2 have arisen over decades of research and hindered reproducibility of C. elegans genetics and genomics. Here we provide a 106.4 Mb gap-free, telomere-to-telomere genome assembly of C. elegans, generated from CGC1, an isogenic derivative of the N2 strain. We used improved long-read sequencing and manual assembly of 43 recalcitrant genomic regions to overcome deficiencies of prior N2 and VC2010 assemblies, and to assemble tandem repeat loci including a 772-kb sequence for the 45S rRNA genes. While many differences from earlier assemblies came from repeat regions, unique additions to the genome were also found. Of 19,972 protein-coding genes in the N2 assembly, 19,790 (99.1%) encode products that are unchanged in the CGC1 assembly. The CGC1 assembly also may encode 183 new protein-coding and 163 new ncRNA genes. CGC1 thus provides both a completely defined reference genome and corresponding isogenic wild-type strain for C. elegans, allowing unique opportunities for model and systems biology.

Metagenomic Analysis of the Buccal Microbiome by Nanopore Sequencing Reveals Structural Differences in the Microbiome of a Patient with Molar Incisor Hypomineralization (MIH) Compared to a Healthy Child-Case Study.

Molar incisor hypomineralization (MIH) is a qualitative developmental defect that affects the enamel tissue of permanent molars and can also occur in permanent incisors. Enamel affected by MIH has reduced hardness, increased porosity, and a higher organic content than unaffected enamel. These characteristics predispose the enamel to accumulation of bacteria and a higher prevalence of caries lesions. Through a groundbreaking metagenomic analysis of the buccal mucosal sample from a patient with MIH, we explored the intricacies of its microbiome compared to a healthy control using state-of-the-art nanopore long-read sequencing. Out of the 210 bacterial taxa identified in the MIH microbiome, we found Streptococcus and Haemophilus to be the most abundant genera. The bacteria with the highest read counts in the patient with MIH included Streptococcus mitis, Haemophilus parainfluenzae, Streptococcus pneumoniae, Rothia dentocariosa, and Gemella haemolysans. Our results revealed a striking contrast between healthy and MIH affected children, with a higher dominance and number of pathogenic species (S. pneumoniae, H. influenzae, and N. meningitidis) and reduced diversity in the MIH-affected patient. This distinct microbial profile not only sheds light on MIH-affected patients, but paves the way for future research, inspiring deeper understanding and larger scale studies.

Gut microbiome and metabolome library construction based on age group using short-read and long-read sequencing techniques in Korean traditional canine species Sapsaree.

This study investigated age-related changes in the gut microbiota and metabolome of Sapsaree dogs through metagenomic and metabolomic analyses. Using Illumina (short-read) and Nanopore (long-read) sequencing technologies, we identified both common and unique bacterial genera in the dogs across different age groups. In metagenomic analysis, Firmicutes were predominant at the family level. At the genus level, Lactobacillus, Streptococcus, Romboutsia, and Clostridium XI were the most abundant, and the bacterial genera typically considered beneficial were less prevalent in senior dogs, whereas the genera associated with pathogenicity were more abundant. These findings suggest age-related shifts in gut microbiota composition. Metabolomic analysis showed distinct clustering of metabolites based on the age group, with changes in metabolite profiles correlating with metagenomic findings. Although Illumina and Nanopore methods provided distinctive results, the genera detected by both methods exhibited similar trends across all age groups in Sapsaree dogs. These findings highlight the relationship between ages, metabolite profiles and gut microbiota composition in dogs, suggesting the need for further research to explore this relation in greater depth.

Multi-omics analysis detail a submicroscopic inv(15)(q14q15) generating fusion transcripts and MEIS2 and NUSAP1 haploinsufficiency

Inversions are balanced structural variants that often remain undetected in genetic diagnostics. We present a female proband with a de novo Chromosome 15 paracentric inversion, disrupting MEIS2 and NUSAP1. The inversion was detected by short-read genome sequencing and confirmed with adaptive long-read sequencing. The breakpoint junction analysis revealed a 96 bp (bp) deletion and an 18 bp insertion in the two junctions, suggesting that the rearrangement arose through a replicative error. Transcriptome sequencing of cultured fibroblasts revealed normal MEIS2 levels and 0.61-fold decreased expression of NUSAP1. Furthermore, three fusion transcripts were detected and confirmed by Sanger sequencing. Heterozygous loss of MEIS2 (MIM# 600987) is associated with a cleft palate, heart malformations, and intellectual impairment, which overlap with the clinical symptoms observed in the proband. The observed fusion transcripts are likely non-functional, and MEIS2 haploinsufficiency is the likely disease causative mechanism. Altogether, this study’s findings illustrate the importance of including inversions in rare disease diagnostic testing and highlight the value of long read sequencing for the validation and characterization of such variants.

The ZFHX3 GGC Repeat Expansion Underlying Spinocerebellar Ataxia Type 4 has a Common Ancestral Founder.

The identification of a heterozygous exonic GGC repeat expansion in ZFHX3 underlying spinocerebellar ataxia type 4 (SCA4) has solved a 25-year diagnostic conundrum. We used adaptive long-read sequencing to decipher the pathogenic expansion in the index Utah family and an unrelated family from Iowa of Swedish ancestry. Contemporaneous to our discovery, other groups identified the same repeat expansion in affected individuals from Utah, Sweden, and Germany, highlighting the current pivotal time for detection of novel repeat expansion disorders. Given that the pathogenic repeat expansion is rare on a population level, we proposed a common ancestor across all families. Here, we employed targeted long-read sequencing through adaptive sampling, enriching for the chr16q22 region of interest. Using phased sequencing results from individuals from Utah, Iowa, and Southern Sweden, we confirmed a common ~2000-year-old ancestral haplotype harbouring the repeat expansion. This study provides further insight into the genetic architecture of SCA4. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Multispecies emergence of dual blaKPC/NDM carbapenemase-producing Enterobacterales recovered from invasive infections in Chile.

Carbapenemase-producing carbapenem-resistant Enterobacterales (CP-CRE) represent a significant global threat. The emergence of dual CP-CRE is particularly alarming, as they can potentially compromise the efficacy of newer antibiotics, further decreasing therapeutic alternatives. Herein, we report the emergence of multiple species of CP-CRE recovered from invasive infections in Chile that simultaneously harbor blaKPC and blaNDM and provide an in-depth genomic characterization of these worrisome pathogens. We collected carbapenem-resistant Enterobacterales (CRE) isolates from invasive infections over a 4-year period, across 11 healthcare centers in Chile. Bacterial species and the presence of carbapenemase genes were confirmed using MALDI-TOF and PCR assays, respectively. Antimicrobial susceptibility testing was conducted through disk diffusion and broth microdilution methods. Dual CP-CRE isolates were subjected to short- and long-read whole genome sequencing to perform a detailed genomic characterization of the isolates and of the mobile genetic elements harboring the enzymes. From a total of 1,335 CRE isolates, we observed an increase in the prevalence of CP-CRE, from 11% in 2019 to 38% in 2022. A total of 11 dual CP-CRE isolates were recovered, all of them harboring blaKPC and blaNDM. Species corresponded to Escherichia coli (n = 6), Klebsiella pneumoniae (n = 2), Klebsiella oxytoca (n = 2), and Citrobacter freundii (n = 1). Dual CP-CRE isolates exhibited resistance to all tested β-lactams except for cefiderocol. The blaKPC and blaNDM encoding genes were located on independent plasmids. Platforms harboring blaKPC were diverse and included IncN, IncF, and IncFIB plasmids. In contrast, blaNDM-7 was only found on fairly conserved IncX3 plasmids. We report that a rapid increase of CP-CRE in Chile, alongside with the emergence of multiple bacterial species of CP-CRE co-harboring blaKPC-2/3 and blaNDM-7, underscores a critical public health challenge. Our data suggest that the dissemination of blaNDM-7 was predominantly facilitated by IncX3 plasmids, whereas the spread of blaKPC involved multiple plasmid backbones. Active surveillance and genomic monitoring are critical to inform public policy and curtail the spread of these highly resistant pathogens.

Rearrangements of viral and human genomes at human papillomavirus integration events and their allele-specific impacts on cancer genome regulation.

Human papillomavirus (HPV) integration has been implicated in transforming HPV infection into cancer. To resolve genome dysregulation associated with HPV integration, we performed Oxford Nanopore long-read sequencing on 72 cervical cancer genomes from an Ugandan dataset that was previously characterized using short-read sequencing. We found recurrent structural rearrangement patterns at HPV integration events, which we categorized as: del(etion)-like, dup(lication)-like, translocation, multibreakpoint, or repeat region integrations. Integrations involving amplified HPV-human concatemers, particularly multibreakpoint events, frequently harbored heterogeneous forms and copy numbers of the viral genome. Transcriptionally active integrants were characterized by unmethylated regions in both the viral and human genomes downstream from the viral transcription start site, resulting in HPV-human fusion transcripts. In contrast, integrants without evidence of expression lacked consistent methylation patterns. Furthermore, whereas transcriptional dysregulation was limited to genes within 200 kilobases of an HPV integrant, dysregulation of the human epigenome in the form of allelic differentially methylated regions affected megabase expanses of the genome, irrespective of the integrant's transcriptional status. By elucidating the structural, epigenetic, and allele-specific impacts of HPV integration, we provide insight into the role of integrated HPV in cervical cancer.

Coupling of alternative splicing and alternative polyadenylation.

RNA splicing and 3'-cleavage and polyadenylation (CPA) are essential processes for the maturation of RNA. There have been extensive independent studies of these regulated processing events, including alternative splicing (AS) and alternative polyadenylation (APA). However, growing evidence suggests potential crosstalk between splicing and 3'-end processing in regulating AS or APA. Here, we first provide a brief overview of the molecular machines involved in splicing and 3'-end processing events, and then review recent studies on the functions and mechanisms of the crosstalk between the two processes. On one hand, 3'-end processing can affect splicing, as 3'-end processing factors and CPA-generated polyA tail promote the splicing of the last intron. Beyond that, 3'-end processing factors can also influence the splicing of internal and terminal exons. Those 3'-end processing factors can also interact with different RNA-binding proteins (RBPs) to exert their effects on AS. The length of 3' untranslated region (3' UTR) can affect the splicing of upstream exons. On the other hand, splicing and CPA may compete within introns in generating different products. Furthermore, splicing within the 3' UTR is a significant factor contributing to 3' UTR diversity. Splicing also influences 3'-end processing through the actions of certain splicing factors. Interestingly, some classical RBPs play dual roles in both splicing and 3'-end processing. Finally, we discuss how long-read sequencing technologies aid in understanding the coordination of AS-APA events and envision that these findings may potentially promote the development of new strategies for disease diagnosis and treatment.

Genetic Analysis and Reproductive Interventions for Two Rare Families Affected by Severe Haemophilia A.

Haemophilia A (HA) is a rare bleeding disorder caused by variants in F8. Although traditional mutational analyses have identified numerous pathogenic variants, the aetiology of HA in certain patients remains unclear. Furthermore, female patients with severe HA are rare. To investigate the molecular defects underlying severe HA in two patients and provide personalised reproductive interventions for their families. Two patients diagnosed with severe HA without other clinical phenotypes were enrolled in the study. A combination of whole-exome sequencing, real-time quantitative polymerase chain reaction and long-read sequencing (LR-sequencing) was performed to reveal the molecular defects of them, followed by the application of different reproductive intervention strategies. Proband 1, a 29-year-old man with FVIII activity of 0.8%, did not exhibit common F8 variants, including Inv1 or Inv22, in the coding region. However, he carried a rare maternal novel inversion on ChrX:154148973_154170321, spanning approximately 21.345Kbp, with breakpoints in introns 13 and 14 of F8. Finally, the couple of Proband 1 opted for assisted reproductive technology using preimplantation genetic testing and successfully conceived. Proband 2, a 20-year-old female with severe HA and FVIII activity of 0.6%, carried inv22 of F8. Further investigation combining whole exome sequencing (WES) and pedigree analysis revealed that she carried a maternal cross-deletion encompassing exons 1-22 of F8, FUNDC2, BRCC3 and CLIC2, along with a de novo missense variant c.5852T>C (p.Leu1951Ser) on her paternal X-chromosome. Chromosome X-inactivation (XCI) analysis demonstrated a highly skewed inactivation of the maternal X chromosome, with a ratio of 98:2. Subsequently, prenatal diagnosis confirmed that the third child in this family did not carry any of the F8 variants present in Proband 2. Our findings provide novel insights into the genetic aetiology of HA and emphasise the importance of a definitive diagnosis in guiding genetic counselling and personalised reproductive interventions for affected individuals and their families.

Exploring the Allelic Diversity of the Self-Incompatibility Gene Across Natural Populations in Petunia (Solanaceae).

Self-incompatibility (SI) is a genetic mechanism to prevent self-fertilization and thereby promote outcrossing in hermaphroditic plant species through discrimination of self and nonself-pollen by pistils. In many SI systems, recognition between pollen and pistils is controlled by a single multiallelic locus (called the S-locus), in which multiple alleles (called S-alleles) are segregating. Because of the extreme level of polymorphism of the S-locus, identification of S-alleles has been a major issue in many SI studies for decades. Here, we report an RNA-seq-based method to explore allelic diversity of the S-locus by employing the long-read sequencing technology of the Oxford Nanopore MinION and applied it for the gametophytic SI system of Petunia (Solanaceae), in which the female determinant is a secreted ribonuclease called S-RNase that inhibits the elongation of self-pollen tubes by degrading RNA. We developed a method to identify S-alleles by the search of S-RNase sequences, using the previously reported sequences as queries, and found in total 62 types of S-RNase including 45 novel types. We validated this method through Sanger sequencing and crossing experiments, confirming the sequencing accuracy and SI phenotypes corresponding to genotypes. Then, using the obtained sequence data together with polymerase chain reaction-based genotyping in a larger sample set of 187 plants, we investigated the diversity, frequency, and the level of shared polymorphism of S-alleles across populations and species. The method and the dataset obtained in Petunia will be an important basis for further studying the evolution of S-RNase-based gametophytic SI systems in natural populations.

Long-read sequencing sheds light on key bacteria contributing to deadwood decomposition processes

BackgroundDeadwood decomposition is an essential ecological process in forest ecosystems, playing a key role in nutrient cycling and carbon sequestration by enriching soils with organic matter. This process is driven by diverse microbial communities encompassing specialized functions in breaking down organic matter, but the specific roles of individual microorganisms in this process are still not fully understood.ResultsHere, we characterized the deadwood microbiome in a natural mixed temperate forest in Central Europe using PacBio HiFi long-read sequencing and a genome-resolved transcriptomics approach in order to uncover key microbial contributors to wood decomposition. We obtained high quality assemblies, which allowed attribution of complex microbial functions such as nitrogen fixation to individual microbial taxa and enabled the recovery of metagenome-assembled genomes (MAGs) from both abundant and rare deadwood bacteria. We successfully assembled 69 MAGs (including 14 high-quality and 7 single-contig genomes) from 4 samples, representing most of the abundant bacterial phyla in deadwood. The MAGs exhibited a rich diversity of carbohydrate-active enzymes (CAZymes), with Myxococcota encoding the highest number of CAZymes and the full complement of enzymes required for cellulose decomposition. For the first time we observed active nitrogen fixation by Steroidobacteraceae, as well as hemicellulose degradation and chitin recycling by Patescibacteria. Furthermore, PacBio HiFi sequencing identified over 1000 biosynthetic gene clusters, highlighting a vast potential for secondary metabolite production in deadwood, particularly in Pseudomonadota and Myxococcota.ConclusionsPacBio HiFi long-read sequencing offers comprehensive insights into deadwood decomposition processes by advancing the identification of functional features involving multiple genes. It represents a robust tool for unraveling novel microbial genomes in complex ecosystems and allows the identification of key microorganisms contributing to deadwood decomposition.

Whole-Genome Resequencing and Evolutionary Analysis of Wild Morel Mushroom Morchella sp.

Morels (Morchella sp.) are important edible fungi cultivated mainly in China. Although the relevant culture technology for Morchella is now fundamentally mature, it is limited to the Elata and Rufobrunnea clades, and the artificial culture technology for the Esculenta clade, which also has economic value, has not been extensively studied. In this study, we selected a wild morel belonging to the Esculenta clade as the research material and performed de novo sequencing and assembly of the Morchella sp. (Mosp) genome using second- and third-generation sequencing. The whole-genome size of Mosp was 55.17 Mb with a contig N50 of 1.89 Mb, and the GC content was 47.49%. A total of 10,896 protein-coding genes were identified. The non-coding RNA prediction results showed that there were 329 tRNAs, 65 rRNAs, and 37 snRNAs in the Mosp genome. The functional annotation of the Mosp genes showed that most of the genes were related to the reproductive and metabolic processes of the cells and participated in nutrient digestion, absorption, utilization, and catabolism in morels. There was a high degree of repetition (21.58%) in the Mosp genome, and the sizes of the DNA transposons and the long terminal repeats were 0.55 Mb and 5.85 Mb, respectively. The phylogeny analysis showed that Mosp clusters together with four other Morchella species: Morchella importuna, Morchella conica, Morchella sextelata, and Morchella snyderi. Molecular dating indicated that the differentiation of Mosp and the black morels group occurred about 147.0 million years ago (MYA). In addition, the evolutionary analysis showed that 296 gene families were contracted and 96 gene families were expanded in Mosp versus the related morel species. The results of this study provide new insights into the genome evolution of Mosp and lay the foundation for future in-depth research into the molecular biology and breeding of the genus Morchella.

Landscape of alternative splicing and polyadenylation during growth and development of muscles in pigs

Alternative polyadenylation (APA) is emerging as a post-transcriptional regulatory mechanism, similar as that of alternative splicing (AS), and plays a prominent role in regulating gene expression and increasing the complexity of the transcriptome and proteome. We use polyadenylation selected long-read isoform sequencing to obtain full-length transcript sequences in porcine muscles at five developmental stages. We identify numerous novel transcripts unannotated in the existing pig genome, including transcripts mapping to known and unknown gene loci, and widespread transcript diversity in porcine muscles. The top 100 most isoformic genes are mainly enriched in Gene Ontology terms related to muscle growth and development. It is revealed that intron retention/exon inclusion and the usage of distal polyadenylation site (PAS) are associated with ageing through analyzing changes of AS and PAS during muscle development. We also identify developmental changes in major transcripts and major PASs. Furthermore, genes/transcripts important for muscle development are identified. The results confirm the importance of AS and APA in pig muscles, substantially increasing transcriptional diversity and showing an important mechanism underlying gene regulation in muscles.

The Genome Architecture of the Copepod Eurytemora carolleeae -the Highly Invasive Atlantic Clade of the Eurytemoraaffinis Species Complex.

Copepods are among the most abundant organisms on the planet and play critical functions in aquatic ecosystems. Among copepods, populations of the Eurytemora affinis species complex are numerically dominant in many coastal habitats and serve as food sources for major fisheries. Intriguingly, certain populations possess the unusual capacity to invade novel salinities on rapid time scales. Despite their ecological importance, high-quality genomic resources have been absent for calanoid copepods, limiting our ability to comprehensively dissect the genome architecture underlying the highly invasive and adaptive capacity of certain populations. Here, we present the first chromosome-level genome of a calanoid copepod, from the Atlantic clade (Eurytemora carolleeae) of the E. affinis species complex. This genome was assembled using high-coverage PacBio long-read and Hi-C sequences of an inbred line, generated through 30 generations of full-sib mating. This genome, consisting of 529.3 Mb (contig N50 = 4.2 Mb, scaffold N50 = 140.6 Mb), was anchored onto four chromosomes. Genome annotation predicted 20,262 protein-coding genes, of which ion transport-related gene families were substantially expanded based on comparative analyses of 12 additional arthropod genomes. Also, we found genome-wide signatures of historical gene body methylation of the ion transport-related genes and the significant clustering of these genes on each chromosome. This genome represents one of the most contiguous copepod genomes to date and is among the highest quality marine invertebrate genomes. As such, this genome provides an invaluable resource to help yield fundamental insights into the ability of this copepod to adapt to rapidly changing environments.

High-resolution profiling reveals coupled transcriptional and translational regulation of transgenes.

Concentrations of RNAs and proteins provide important determinants of cell fate. Robust gene circuit design requires an understanding of how the combined actions of individual genetic components influence both mRNA and protein levels. Here, we simultaneously measure mRNA and protein levels in single cells using HCR Flow-FISH for a set of commonly used synthetic promoters. We find that promoters generate differences in both the mRNA abundance and the effective translation rate of these transcripts. Stronger promoters not only transcribe more RNA but also show higher effective translation rates. While the strength of the promoter is largely preserved upon genome integration with identical elements, the choice of polyadenylation signal and coding sequence can generate large differences in the profiles of the mRNAs and proteins. We used long-read direct RNA sequencing to characterize full-length mRNA isoforms and observe remarkable uniformity of mRNA isoforms from the transgenic units. Together, our high-resolution profiling of transgenic mRNAs and proteins offers insight into the impact of common synthetic genetic components on transcriptional and translational mechanisms. By developing a novel framework for quantifying expression profiles of transgenes, we have established a system for comparing native and synthetic gene regulation and for building more robust transgenic systems.

The whole chromosome-level genome provides resources and insights into the endangered fish Percocypris pingi evolution and conservation

BackgroundPercocypris pingi (Tchang) was classified as Endangered on the Red List of China′s Vertebrates in 2015 and is widely distributed in the Upper Yangtze River. Although breeding and release into wild habitats have been performed for this commercially important fish in recent years, low genetic diversity has been found in wild populations. Genomic resources are strongly recommended before formulating and carrying out conservation strategies for P. pingi. Thus, there is an urgent need to conserve germplasm resources and improve the population diversity of P. pingi. To date, the whole genome of P. pingi has not been reported.ResultsIn our study, we constructed the first chromosome-level genome of P. pingi by high-throughput chromosome conformation capture (Hi-C) technology and PacBio long-read sequencing. The assembled genome was 1.7 Gb in size, with an N50 of 17,692 bp and a GC content from circular consensus sequencing of 37.67%. The Hi-C results again demonstrated that P. pingi was tetraploid (n = 98), with the genome consisting of 24-type and 25-type chromosomes. Chr.19 of the 24-type chromosomes in P. pingi resulted from the fusion of chr.19 and chr.22 in zebrafish. The divergence times between 24-type and 25-type chromosomes was around 6.1 million years ago. A total of 25,198 and 25,291 protein-coding genes were obtained from the 24-type and 25-type chromosomes, respectively. The ploidy of P. pingi is an allotetraploid. A total of 8,741 genes of P. pingi were clustered into 4,378 gene families that were shared with 14 other species, and the P. pingi genome had 68 unique gene families. Phylogenetic analyses indicated that P. pingi was most closely related to Schizothorax oconnori, and the genes were clustered on one branch. We identified 166 significantly expanded gene families and 173 significantly contracted gene families in P. pingi. The most enriched positive protein-coding genes, such as Bmp-4, Etfdh, homeobox protein HB9, and ATG3, were screened.ConclusionOur study provides a high-quality chromosome-anchored reference genome for P. pingi and provides sufficient information on the chromosomes, which will lead to valuable resources for genetic, genomic, and biological studies of P. pingi and for improving the genetic diversity, population size, and scientific conservation of endangered fish and other key cyprinid species in aquaculture.

MetaCONNET: A metagenomic polishing tool for long-read assemblies.

Accurate and high coverage genome assemblies are the basis for downstream analysis of metagenomic studies. Long-read sequencing technology is an ideal tool to facilitate the assemblies of metagenome, except for the drawback of usually producing reads with high sequencing error rate. Many polishing tools were developed to correct the sequencing error, but most are designed on the ground of one or two species. Considering the complexity and uneven depth of metagenomic study, we present a novel deep-learning polishing tool named MetaCONNET for polishing metagenomic assemblies. We evaluate MetaCONNET against Medaka, CONNET and NextPolish in accuracy, coverage, contiguity and resource consumption. Our results demonstrate that MetaCONNET provides a valuable polishing tool and can be applied to many metagenomic studies.

Identification and characterization of human KALRN mRNA and Kalirin protein isoforms.

Kalirin is a multidomain protein with important roles in neurite outgrowth, and synaptic spine formation and remodeling. Genetic and pathophysiological links with various neuropsychiatric disorders associated with synaptic dysfunction and cognitive impairment have sparked interest in its potential as a pharmacological target. Multiple Kalirin proteoforms are detected in the adult human brain, yet we know little about the diversity of the transcripts that encode them or their tissue profiles. Here, we characterized full-length KALRN transcripts expressed in the adult human frontal lobe and hippocampus using rapid amplification of complementary DNA (cDNA) ends and nanopore long-read sequencing. For comparison with non-neural tissue, we also analyzed KALRN transcripts in the aorta. Multiple novel isoforms were identified and were largely similar between the two brain regions analyzed. Alternative splicing in the brain results in preferential inclusion of exon 37, which encodes 32 amino acids upstream of the second guanine nucleotide exchange factor (GEF) domain. Structural modeling predicts that a subset of these amino acids forms a conserved alpha helix. Although deletion of these amino acids had little effect on GEF activity, it did alter Kalirin-induced neurite outgrowth suggesting that this brain-enriched splicing event may be important for neural function. These data indicate that alternative splicing is potentially important for regulating Kalirin actions in the human brain.

The rate and spectrum of new mutations in mice inferred by long-read sequencing.

All forms of genetic variation originate from new mutations, making it crucial to understand their rates and mechanisms. Here, we use long-read sequencing from Pacific Biosciences (PacBio) to investigate de novo mutations that accumulated in 12 inbred mouse lines derived from three commonly used inbred strains (C3H, C57BL/6, and FVB) maintained for 8 to 15 generations in a mutation accumulation (MA) experiment. We built chromosome-level genome assemblies based on the MA line founders' genomes and then employed a combination of read and assembly-based methods to call the complete spectrum of new mutations. On average, there are about 45 mutations per haploid genome per generation, about half of which (54%) are insertions and deletions shorter than 50 bp (indels). The remainder are single-nucleotide mutations (SNMs; 44%) and large structural mutations (SMs; 2%). We found that the degree of DNA repetitiveness is positively correlated with SNM and indel rates and that a substantial fraction of SMs can be explained by homology-dependent mechanisms associated with repeat sequences. Most (90%) indels can be attributed to microsatellite contractions and expansions, and there is a marked bias toward 4 bp indels. Among the different types of SMs, tandem repeat mutations have the highest mutation rate, followed by insertions of transposable elements (TEs). We uncover a rich landscape of active TEs, notable differences in their spectrum among MA lines and strains, and a high rate of gene retroposition. Our study offers novel insights into mammalian genome evolution and highlights the importance of repetitive elements in shaping genomic diversity.

Adaptable and comprehensive approaches for long-read nanopore sequencing of polyadenylated and non-polyadenylated RNAs.

The advent of long-read (LR) sequencing technologies has provided a direct opportunity to determine the structure of transcripts with potential for end-to-end sequencing of full-length RNAs. LR methods that have been described to date include commercial offerings from Oxford Nanopore Technologies (ONT) and Pacific Biosciences. These kits are based on selection of polyadenylated (polyA+) RNAs and/or oligo-dT priming of reverse transcription. Thus, these approaches do not allow comprehensive interrogation of the transcriptome due to their exclusion of non-polyadenylated (polyA-) RNAs. In addition, polyA + specificity also results in 3'-biased measurements of PolyA+ RNAs especially when the RNA input is partially degraded. To address these limitations of current LR protocols, we modified rRNA depletion protocols that have been used in short-read sequencing: one approach representing a ligation-based method and the other a template-switch cDNA synthesis-based method to append ONT-specific adaptor sequences and by removing any deliberate fragmentation/shearing of RNA/cDNA. Here, we present comparisons with poly+ RNA-specific versions of the two approaches including the ONT PCR-cDNA Barcoding kit. The rRNA depletion protocols displayed higher proportions (30%-50%) of intronic content compared to that of the polyA-specific protocols (5%-8%). In addition, the rRNA depletion protocols enabled ∼20-50% higher detection of expressed genes. Other metrics that were favourable to the rRNA depletion protocols include better coverage of long transcripts, and higher accuracy and reproducibility of expression measurements. Overall, these results indicate that the rRNA depletion-based protocols described here allow the comprehensive characterization of polyadenylated and non-polyadenylated RNAs. While the resulting reads are long enough to help decipher transcript structures, future endeavors are warranted to improve the proportion of individual reads representing end-to-end spanning of transcripts.