Sequence Reads Research Articles

Novel variants and genotype-phenotype correlation in a multicentre cohort of GNE myopathy in China

BackgroundGlcNAc2-epimerase (GNE) myopathy is a rare autosomal recessive disorder caused by pathogenic variants in the GNE gene, which is essential for the sialic acid biosynthesis pathway.ObjectiveThis multi-centre study aimed to...

Contrasting and combining transcriptome complexity captured by short and long RNA sequencing reads.

Mapping transcriptomic variations using either short- or long-reads RNA sequencing is a staple of genomic research. Long reads are able to capture entire isoforms and overcome repetitive regions, while short reads still provide improved coverage and error rates. Yet how to quantitatively compare the technologies, can we combine those, and what may be the benefit of such a combined view remain open questions. We tackle these questions by first creating a pipeline to assess matched long and short reads data using a variety of transcriptome statistics. We find that across datasets, algorithms, and technologies, matched short reads data detects roughly 30% more splice junctions such that 10-30% of the splice junctions included at 20% or more by short reads are missed by long reads. In contrast, long reads detect many more intron retention events and can detect full isoforms, pointing to the benefit of combining the technologies. We introduce MAJIQ-L, an extension of the MAJIQ software to enable a unified view of transcriptome variations from both technologies and demonstrate its benefits. Our software can be used to assess any future long-read technology or algorithm, and combine it with short reads data for improved transcriptome analysis.

FamLink2 – A comprehensive tool for likelihood computations in pedigrees analyses involving linked DNA markers accounting for genotype uncertainties

There is an increasing demand for software that can handle an arbitrary number of linked markers in forensic genetics; primarily with application to inference of relationships and direct matching but also in applications such as ancestry inference and mixture interpretation. With the emergence of sequencing technologies, denser sets of SNP markers are generated and analyzed. Additionally, sequence data of low quality and quantity DNA generate uncertainty about the underlying true genotype. We provide an efficient implementation of a general model for pedigree likelihood computations with genetic marker data using a three-layered approach. The top and first layer is the population model where allele frequencies and population substructure are accounted for. The second layer is the inheritance model which efficiently handles linked markers using an IBD model. The third and bottom layer is the observational level where we model the likelihood of the true genotype given underlying reads as well as parameters for errors. We exemplify the utility of our implementation as well as provide validation according to guidelines established by the ISFG using a combination of two published SNP panels. We demonstrate that computations are feasible for panels encompassing 10,000 markers and we argue that, due to the properties of the underlying algorithm, extending the number of markers will result in a linear increase in computation time. In addition we study the impact of parameters used in our model and suggest some guidelines pertaining to their values. The results demonstrate that a probabilistic model for low coverage sequence read data is needed instead of relying on an a threshold based genotype and applying our general model for inference of relationships on a real case can be superior, i.e. higher information content, to other methods relying on either fixed genotypes with low quality sequence data or simple pair wise relationship tests. In summary, the implementation, FamLink2 (freely available at https://famlink.se), can jointly handle genetic linkage, genotype uncertainty and population substructure for an arbitrary pedigree with data for any number of individuals. Whereas the current study will focus on calculations disregarding mutations, FamLink2 has the ability to model mutations for certain built-in pedigrees.

Scalable and unsupervised discovery from raw sequencing reads using SPLASH2.

We introduce SPLASH2, a fast, scalable implementation of SPLASH based on an efficient k-mer counting approach for regulated sequence variation detection in massive datasets from a wide range of sequencing technologies and biological contexts. We demonstrate biological discovery by SPLASH2 in single-cell RNA sequencing (RNA-seq) data and in bulk RNA-seq data from the Cancer Cell Line Encyclopedia, including unannotated alternative splicing in cancer transcriptomes and sensitive detection of circular RNA.

Unveiling viral pathogens in acute respiratory disease: Insights from viral metagenomics in patients from the State of Alagoas, Brazil.

Respiratory illness affects individuals across all age demographics on a global scale, often precipitated by viral infections. The symptomatic manifestations of these diseases bear clinical resemblance, complicating the accurate determination of their etiological origins. Furthermore, the diagnostic panels for respiratory pathogens used within local medical practices, may not encompass the full spectrum of viral agents responsible for such ailments. Consequently, a significant number of clinically important viral pathogens may remain undetected. In the light of this, we conducted a metagenomic examination of 66 nasopharyngeal swab specimens, obtained from patients presenting with acute respiratory conditions yet tested negative by the standard diagnostic panels available locally. These specimens were obtained from the Public Health Laboratory, Maceio, State of Alagoas. Our findings indicate a predominant diagnostic escape of rhinoviruses and notably enterovirus D68. Moreover, our study identified a substantial quantity of sequence reads attributed to human respirovirus 3 (human parainfluenza 3) along with various herpresviruses including human herpesvirus-1, Epstein-Barr virus (Human herpesvirus-4), Human herpesviruses 6 and 7 and human parvovirus B19 (B19V). Notably, the metagenomic analysis uncovered a widespread presence of the emerging human vientovirus FB in most of sample pools, though its clinical importance remains to be elucidated. The obtained results in this study underscore the invaluable role of viral metagenomics in the identification of underrecognized viruses bearing clinical relevance. Furthermore, it offers insights into the dissemination of these pathogens within the studied area, thereby informing public health strategies aimed at enhancing diagnostic accuracy and improving patient care.

3t-seq: automatic gene expression analysis of single-copy genes, transposable elements, and tRNAs from RNA-seq data.

RNA sequencing is the gold-standard method to quantify transcriptomic changes between two conditions. The overwhelming majority of data analysis methods available are focused on polyadenylated RNA transcribed from single-copy genes and overlook transcripts from repeated sequences such as transposable elements (TEs). These self-autonomous genetic elements are increasingly studied, and specialized tools designed to handle multimapping sequencing reads are available. Transfer RNAs are transcribed by RNA polymerase III and are essential for protein translation. There is a need for integrated software that is able to analyze multiple types of RNA. Here, we present 3t-seq, a Snakemake pipeline for integrated differential expression analysis of transcripts from single-copy genes, TEs, and tRNA. 3t-seq produces an accessible report and easy-to-use results for downstream analysis starting from raw sequencing data and performing quality control, genome mapping, gene expression quantification, and statistical testing. It implements three methods to quantify TEs expression and one for tRNA genes. It provides an easy-to-configure method to manage software dependencies that lets the user focus on results. 3t-seq is released under MIT license and is available at https://github.com/boulardlab/3t-seq.

Root Endophytic Microorganisms Contribute to the Attribute of Full-Year Shooting in Woody Bamboo Cephalostachyum pingbianense.

Cephalostachyum pingbianense (Hsueh & Y.M. Yang ex Yi et al.) D.Z. Li & H.Q. Yang is unique among bamboo species for its ability to produce bamboo shoots in all seasons under natural conditions. Apart from the physiological mechanism, information regarding the effects of endophytic microorganisms on this full-year shooting characteristic is limited. We hypothesize that root endophytic microorganisms will have a positive impact on the full-year bamboo shooting characteristic of C. pingbianense by increasing the availability or supply of nutrients. To identify the seasonal variations in the root endophytic bacterial and fungal communities of C. pingbianense, and to assess their correlation with bamboo shoot productivity, the roots of C. pingbianense were selected as research materials, and the 16S rRNA and ITS rDNA genes of root endophytic microorganisms were sequenced using the Illumina platform. Following this sequencing, raw sequencing reads were processed, and OTUs were annotated. Alpha and beta diversity, microbial composition, and functional predictions were analyzed, with correlations to bamboo shoot numbers assessed. The results showed that seasonal changes significantly affected the community diversity and structure of root endophytic microbes of C. pingbianense. Bacterial communities in root samples from all seasons contained more nitrogen-fixing microorganisms, with members of the Burkholderiales and Rhizobiales predominating. The relative abundances of ectomycorrhizal and arbuscular mycorrhizal fungi in the autumn sample were significantly higher than in other seasons. Correlation analysis revealed that the bamboo shoot productivity was significantly and positively correlated with bacterial functions of nitrogen fixation, arsenate detoxification, and ureolysis, as well as with symbiotrophic fungi, ectomycorrhizal fungi, and arbuscular mycorrhizal fungi. At the genus level, the bacterial genus Herbaspirillum and the fungal genera Russula, unclassified_f_Acaulosporaceae, and unclassified_f_Glomeraceae were found to have a significant positive correlation with bamboo shoot number. Our study provides an ecological perspective for understanding the highly productive attribute of C. pingbianense and offers new insights into the forest management of woody bamboos.

Abstract C102: Polyethnic-1000: Advancing cancer genomics by studying ethnically diverse, underserved patient populations in New York

Abstract Prompted by the considerable biases toward European ancestry cases in large cancer genomics repositories, the New York Genome Center initiated an effort to diversify the public datasets that are forming the foundations of current cancer research, in order to optimize the unbiased applications of modern genomically-driven discoveries to a larger population. Our effort was focused on 8 cancer types (breast, prostate, bladder, lung, pancreatic, colon and endometrial cancer, as well as multiple myeloma) and generated Illumina tumor-normal whole-genome sequencing and tumor transcriptome sequencing. The majority of the samples were retrospectively obtained in clinics and hospitals of the New York region and the analyses are conducted with collaborators in these institutions. Our dataset is currently composed of 851 tumor-normal WGS and 471 RNASeq, and is progressively being distributed to the scientific community in partnership with ICGC-ARGO. Overall, 75% of the cohort is predominantly of non-European ancestry, as estimated by genetic ancestry analysis using the 1000Genomes continental references.We analyzed the samples with NYGC’s somatic variant calling pipeline and will present results from the 8 cancer types related to germline and somatic variants (single-nucleotide variants, indels, copy-number alterations, complex structural variants) as well as mutational signatures, neoantigen prediction, fusion genes and expression profiles, in relation to genetic ancestry estimation and clinical features of the patients. We anticipate that the dataset will significantly diversify the genetic ancestry of tumor profiles in public databases, and serve to form hypotheses regarding the contribution of ancestry on cancer disparities. The experience of the Polyethnic-1000 project informs our contribution to the Cancer Grand Challenge SAMBAI, which focuses on breast, prostate and pancreatic cancers in the African Diaspora. SAMBAI will further investigate the causes of cancer disparity by combining genomic data with social determinants of health and exposomics information, deeper immune profiling and new genomics technologies such as long reads sequencing and the use of a pangenome reference. Citation Format: Nicolas Robine, Lara Winterkorn, Tim Chu, Zoe Goldstein, Will Hooper, Ali Oku, Heather M. Geiger, Melissa B. Davis, Onyinye Balogun. Polyethnic-1000: Advancing cancer genomics by studying ethnically diverse, underserved patient populations in New York [abstract]. In: Proceedings of the 17th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2024 Sep 21-24; Los Angeles, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2024;33(9 Suppl):Abstract nr C102.

Evaluating the quantitative performance of environmental DNA metabarcoding for freshwater zooplankton community: a case study in Lake Biwa, Japan.

Zooplankton monitoring is important for understanding their population dynamics and life history, ecosystem health, and environmental changes. Compared with traditional morphological identification, environmental DNA (eDNA) analysis allows for more sensitive and efficient monitoring of zooplankton diversity. Previous eDNA studies have primarily used metabarcoding approaches to reveal their richness and composition, whereas its performance in predicting zooplankton abundance remains understudied. We conducted water and bulk sampling in Lake Biwa, Japan, showing that the number of sequence reads by metabarcoding moderately correlated with eDNA concentrations estimated by quantitative real-time PCR (qPCR). In addition, the eDNA read number was significantly related to cladoceran and copepod abundance estimated by microscopy sorting, although there remained too much uncertainty in the read-abundance relationship. Moreover, there was a significant difference in species composition between eDNA metabarcoding and sorting. Although our results indicated the potential applicability of eDNA metabarcoding for quantifying multiple zooplankton abundance, several methodological validations in eDNA metabarcoding would also be required to optimize its performance in the future.

Characterization of the complete mitogenome of Tiarella polyphylla, commonly known as Asian foamflower: insights into the multi-chromosomes structure and DNA transfers

BackgroundTiarella polyphylla D. Don has been traditionally used to cure asthma and skin eruptions. However, the sequence and the structure of the mitogenome of T. polyphylla remained elusive, limiting the genomic and evolution analysis based on the mitogenome.ResultsUsing a combination of Illumina and Nanopore sequencing reads, we de novo assembled the complete mitogenome of T. polyphylla. In addition to unveiling the major configuration of the T. polyphylla mitogenome was three circular chromosomes with lengths of 430,435 bp, 126,943 bp, and 55,269 bp, we revealed five (R01-R05) and one (R06) repetitive sequence could mediate the intra- and inter-chromosomal recombination, respectively. Furthermore, we identified 208 short and 25 long tandem segments, seven cp-derived mtDNAs, 106 segments of mtDNAs transferred to the nuclear genome, and 653 predicted RNA editing sites. Based on the sequence of the mitogenomes, we obtained the resolved phylogeny of the seven Saxifragales species.ConclusionsThese results presented the mitogenome features and expanded its potential applications in phylogenetics, species identification, and cytoplasmic male sterility (CMS) in the future.

Whole genome re-sequencing reveals high altitude adaptation signatures and admixture in Ladakhi cattle

Ladakhi cattle, native to the high-altitude region of Ladakh in northern India (ranging from 3,000 to 5,000 m above sea level), have evolved unique genetic adaptations to thrive in harsh environmental conditions, such as hypoxia, extreme cold, and low humidity. This study explored the genome of Ladakhi cattle to investigate genetic structure, selection signatures, and adaptive mechanisms. Whole genome sequencing reads, generated on Illumina NovaSeq 6000 platform, were aligned to the Bos taurus reference genome with BWA-MEM. SNPs were identified and filtered using GATK and bcftools, and functionally annotated with SnpEff. For population genomic analysis, PCA and admixture modeling assessed genetic structure, while Neighbor-Joining trees, LD decay, nucleotide diversity (π), and FST evaluated phylogenetic relationships and genetic variation. Selective sweeps were detected using RAiSD, and gene-set enrichment and protein–protein interaction analyses were conducted to explore functional pathways related to adaptation. The study revealed 3,759,279 unique SNPs and demonstrated that Ladakhi cattle form a distinct genetic cluster with an estimated admixture of 68 % Bos indicus and 32 % Bos taurus ancestry. Key findings include rapid linkage disequilibrium decay, low inbreeding level, and the identification of selection signatures and genes associated with hypoxia response, energy metabolism, and cold adaptation. Mean nucleotide diversity (π, 0.0037) and FST values indicated moderate genetic differentiation from other breeds. The analysis highlighted selection signatures for genes like HIF1A, ENO4, ANGPT1, EPO, NOS3, MAPK3, HMOX1, BCL2,CAMK2D, MTOR, AKT2,PIK3CB, and MAP2K1, among others, including various keratin and heat shock proteins. The interaction between genes associated with hypoxia signaling (HIF-1) and other enriched pathways such as PI3K, mTOR, NFκB, ERK, and ER stress, reveals a complex mechanism for managing hypoxic stress in Ladakhi cattle. These findings offer valuable insights for breeding programs aimed at enhancing livestock resilience in extreme environments and enhance understanding of mammalian adaptation to high-altitude conditions.

Chromosome-level genome assembly of Chinese water Scorpion Ranatra chinensis (Heteroptera: Nepidae)

Heteroptera (the true bugs), one of the most diverse lineages of insects, diversified in feeding strategies and living habitats, and thus become an ideal lineage for studies on adaptive evolution. Chinese water scorpion Ranatra chinensis (Heteroptera: Nepidae) is a predaceous bug living in lentic water systems, representing an ideal model for studying habitat transition and adaptation to water environment. However, genetic studies on this water bug remain limited. Here, we obtained a chromosome-level genome of R. chinensis using PacBio HiFi long reads and Hi-C sequencing reads. The total assembly size of genome is 867.89 Mb, with a scaffold N50 length of 26.48 Mb and the GC content of 39.50%. All contigs were assembled into 23 pseudo-chromosomes (N = 19 A + X1X2X3X4), and we predicted 18,424 protein-coding genes in this genome. This study will provide valuable genomic resources for future studies on the biology, water adaptation, and genome evolution of water bugs.

Deciphering transcript architectural complexity in bacteria and archaea.

RNA transcripts are potential therapeutic targets, yet bacterial transcripts have uncharacterized biodiversity. We developed an algorithm for transcript prediction called tp.py using it to predict transcripts (mRNA and other RNAs) in Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria), Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria). From >5 million E. coli K12 and >3 million E. coli E2348/69 newly generated Oxford Nanopore Technologies direct RNA sequencing reads, 2,487 K12 mRNAs and 1,844 E2348/69 mRNAs were predicted, with the K12 mRNAs containing more than half of the predicted E. coli K12 proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used, across all strains examined, the predicted average size of the mRNAs was 1.6-1.7 kbp, while the median size of the 5'- and 3'-untranslated regions (UTRs) were 30-90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions were of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100-200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, and reproducible method will facilitate the presentation of transcripts, and UTR predictions alongside coding sequences and protein predictions in bacterial genome annotation as important resources for the research community.IMPORTANCEOur understanding of bacterial and archaeal genes and genomes is largely focused on proteins since there have only been limited efforts to describe bacterial/archaeal RNA diversity. This contrasts with studies on the human genome, where transcripts were sequenced prior to the release of the human genome over two decades ago. We developed software for the quick, easy, and reproducible prediction of bacterial and archaeal transcripts from Oxford Nanopore Technologies direct RNA sequencing data. These predictions are urgently needed for more accurate studies examining bacterial/archaeal gene regulation, including regulation of virulence factors, and for the development of novel RNA-based therapeutics and diagnostics to combat bacterial pathogens, like those with extreme antimicrobial resistance.

A Haplotype-resolved, Chromosome-scale Genome for Malus domestica Borkh. 'WA 38'.

Genome sequencing for agriculturally important Rosaceous crops has made rapid progress both in completeness and annotation quality. Whole genome sequence and annotation gives breeders, researchers, and growers information about cultivar specific traits such as fruit quality and disease resistance, and informs strategies to enhance postharvest storage. Here we present a haplotype-phased, chromosomal level genome of Malus domestica, 'WA 38', a new apple cultivar released to market in 2017 as Cosmic Crisp®. Using both short and long read sequencing data with a k-mer based approach, chromosomes originating from each parent were assembled and segregated. This is the first pome fruit genome fully phased into parental haplotypes in which chromosomes from each parent are identified and separated into their unique, respective haplomes. The two haplome assemblies, 'Honeycrisp' originated HapA and 'Enterprise' originated HapB, are about 650 Megabases each, and both have a BUSCO score of 98.7% complete. A total of 53,028 and 54,235 genes were annotated from HapA and HapB, respectively. Additionally, we provide genome-scale comparisons to 'Gala', 'Honeycrisp', and other relevant cultivars highlighting major differences in genome structure and gene family circumscription. This assembly and annotation was done in collaboration with the American Campus Tree Genomes project that includes 'WA 38' (Washington State University), 'd'Anjou' pear (Auburn University), and many more. To ensure transparency, reproducibility, and applicability for any genome project, our genome assembly and annotation workflow is recorded in detail and shared under a public GitLab repository. All software is containerized, offering a simple implementation of the workflow.

Performance of DNA metabarcoding, standard barcoding and morphological approaches in the identification of insect biodiversity.

For two decades, DNA barcoding and, more recently, DNA metabarcoding have been used for molecular species identification and estimating biodiversity. Despite their growing use, few studies have systematically evaluated these methods. This study aims to evaluate the efficacy of barcoding methods in identifying species and estimating biodiversity, by assessing their consistency with traditional morphological identification and evaluating how assignment consistency is influenced by taxonomic group, sequence similarity thresholds and geographic distance. We first analysed 951 insect specimens across three taxonomic groups: butterflies, bumblebees and parasitic wasps, using both morphological taxonomy and single-specimen COI DNA barcoding. An additional 25,047 butterfly specimens were identified by COI DNA metabarcoding. Finally, we performed a systematic review of 99 studies to assess average consistency between insect species identity assigned via morphology and COI barcoding and to examine the distribution of research effort. Species assignment consistency was influenced by taxonomic group, sequence similarity thresholds and geographic distance. An average assignment consistency of 49% was found across taxonomic groups, with parasitic wasps displaying lower consistency due to taxonomic impediment. The number of missing matches doubled with a 100% sequence similarity threshold and COI intraspecific variation increased with geographic distance. Metabarcoding results aligned well with morphological biodiversity estimates and a strong positive correlation between sequence reads and species abundance was found. The systematic review revealed an 89% average consistency and also indicated taxonomic and geographic biases in research effort. Together, our findings demonstrate that while problems persist, barcoding approaches offer robust alternatives to traditional taxonomy for biodiversity assessment.

Epigenetic control and inheritance of rDNA arrays.

Ribosomal RNA (rRNA) genes exist in multiple copies arranged in tandem arrays known as ribosomal DNA (rDNA). The total number of gene copies is variable, and the mechanisms buffering this copy number variation remain unresolved. We surveyed the number, distribution, and activity of rDNA arrays at the level of individual chromosomes across multiple human and primate genomes. Each individual possessed a unique fingerprint of copy number distribution and activity of rDNA arrays. In some cases, entire rDNA arrays were transcriptionally silent. Silent rDNA arrays showed reduced association with the nucleolus and decreased interchromosomal interactions, indicating that the nucleolar organizer function of rDNA depends on transcriptional activity. Methyl-sequencing of flow-sorted chromosomes, combined with long read sequencing, showed epigenetic modification of rDNA promoter and coding region by DNA methylation. Silent arrays were in a closed chromatin state, as indicated by the accessibility profiles derived from Fiber-seq. Removing DNA methylation restored the transcriptional activity of silent arrays. Array activity status remained stable through the iPS cell re-programming. Family trio analysis demonstrated that the inactive rDNA haplotype can be traced to one of the parental genomes, suggesting that the epigenetic state of rDNA arrays may be heritable. We propose that the dosage of rRNA genes is epigenetically regulated by DNA methylation, and these methylation patterns specify nucleolar organizer function and can propagate transgenerationally.

De novo assembly of sialotranscriptome of Hyalomma anatolicum and insights into expression dynamics in response to Theileria annulata infection.

Hyalomma anatolicum is a tick of significant one-health importance due to its role as a vector for various pathogens affecting humans, animals and the environment, such as Theileria annulata, which causes tropical theileriosis in cattle, leading to severe economic losses. When infected with pathogens like T. annulata, the salivary glands of H. anatolicum undergo gene expression changes, secrete modified proteins and activate immune responses, all of which facilitate pathogen survival and transmission by modulating the host immune response and optimizing conditions for pathogen development. Understanding these responses is crucial for developing control strategies for tick-borne diseases. To understand the interaction between H. anatolicum and T. annulata, we performed a differential gene expression analysis of H. anatolicum salivary glands. An average of approximately 25 million raw sequencing reads were generated in each replicate using Illumina Sequencing. The sequenced reads were de novo assembled and the assembled transcriptome yielded approximately 50,231 non-redundant transcripts after clustering with CD-HIT using a sequence identity of 95% and alignment coverage of 90%. The assembly quality was evaluated with BUSCO analysis and found to be 86% complete using the Arachnida dataset and then blasted against non-redundant protein sequence database from NCBI followed by counting of reads and differential expression analysis. Overall, around 2400 and 400 genes were found differentially expressed with logFC > 1 and logFC > 2 respectively at FDR < 0.05. Top up-regulated genes included Calpain, Papilin, Neprilysin, and Ankyrin repeat-containing protein. Top down-regulated genes included Scoloptoxin, and Selenoprotein S and other uncharacterized proteins. Many other up-regulated proteins with high significance were uncharacterized suggesting room for further H. anatolicum functional and structural characterization studies. To our best knowledge, this is the first study of H. anatolicum sialotranscriptome which greatly contributes to sialotranscriptome information not only as sequence database but also indicates the potential targets for development of vaccine against ticks and transmission-blocking vaccines against T. annulata.

A retrospective whole genome sequencing of SARS-CoV-2 isolate from respiratory sample of an individual from Ota – Ogun State, Nigeria

Severe acute respiratory syndrome coronavirus – 2 (SARS-CoV-2) which was responsible for the COVID-19 pandemic has now been considered an endemic virus, that will continually produce sporadic outbreaks in different communities around the globe. Although, there are several surveillance studies on SARS-CoV-2 globally, including Nigeria, there is still an important need to understand the uniqueness of the strains of the virus that was circulating immediately after the pandemic, to add to the global database of information that will aid vaccine production and future preparedness against another SARS - related CoV pandemic.Towards the end of the pandemic, between February – August 2022, SARS-CoV-2 was detected in a surveillance project in Ota – Ogun State Nigeria. We carried out a retrospective whole genome sequencing (WGS) analysis of SARS-CoV-2 in the previously reported positive sample, at Inqaba biotech in South Africa. RNA extraction was carried out using the Quick – RNA viral kit (Zymo) and library preparation was done by NEBNext ARTIC SARS-CoV-2 FS Library Prep kit (Illumina) according to the manufacturer's instruction. The WGS was carried out on the NextSeq500 platform by Illumina. The fastq file containing the unassembled raw sequence reads were submitted to NCBI SARS-CoV-2 resources, and a BioProject accession number PRJNA1076330 was issued. The DRAGEN Targeted Microbial, GISAID-CoVsurver mutation, BLAST and MAFFT applications were used for analysis.The spike (s) gene of the study sequence possessed seven mutations (G142D, A163V, V213G, D614G, H655Y, N679K, P681H) and shared 99.4 % identity with those of the Wuhan reference sequence WIV04. The non-structural proteins (NSP) – 7,8,9,10 and 16 shared 100 % identity with bat/Yunnan/RaTG13 (a SARS- related CoV found in bats) sequence on GISAID. Although, the study sequence was obtained from an asymptomatic individual, the S mutations observed are known to be related to virulence, antigenic shift, enhanced transmissibility and host change. Thus, upon successful exploitation of this asymptomatic variant, it may possibly be transformed to a potential candidate for SARS-CoV-2 vaccine in future.

Cave Myotis (Myotis velifer) consume diverse prey items and provide important ecosystem services

Abstract Insectivorous bats play an essential role as predators in natural ecosystems and contribute to pest control in agricultural landscapes. However, characterizing diets of specific bat species is difficult using conventional methods that cannot capture detailed dietary information. In this study, we used metabarcoding of the cytochrome oxidase I mitochondrial gene to analyze fecal samples of Myotis velifer and provide insight into the seasonal variation of diet from a colony located in the Chihuahuan Desert region of Texas. After filtering sequence reads, we recovered and analyzed 706 molecular operational taxonomic units (MOTUs) in the diet of M. velifer. We found 484 taxa (species and genus level) belonging to 11 insect orders in 66 fecal samples collected from March to October 2021. The orders containing the most MOTUs were Diptera (n = 353), Lepidoptera (n = 160), and Blattodea (n = 59). In their diet, we identified important insect crop pests, non-native insects, and substantial consumption of mosquitoes. Dietary composition shifted throughout seasons. Bats captured in summer months consumed the highest diversity of arthropod orders, and those captured in the spring consumed more diet items but were less diverse in arthropod orders. Our results uncovered 2 additional orders consumed by M. velifer, Blattodea and Odonata, that have not been previously described in their diet.

Anaerobic bacteria dominate the cholesteatoma tissue of chronic suppurative otitis media patients

The aim of this study was to investigate both the microbial composition and absolute abundance of clinically relevant bacteria in tissue specimens from patients with chronic suppurative otitis media with cholesteatoma (CSOM with cholesteatoma). Mastoid mucosa and cholesteatoma tissue from eleven subjects with CSOM with cholesteatoma, and mastoid mucosa from ten controls were examined using standard hospital culture swabs, Gram staining, bacterial 16S rRNA gene sequencing, Droplet Digital PCR (ddPCR), and multiplex PCR. Positive results from culture swabs were reported in half the CSOM with cholesteatoma samples and 1 control sample. In contrast, ddPCR detected bacterial genes copies in all 11 mucosa and cholesteatoma of CSOM subjects and 3 control samples. The average bacterial gene copies in tissue samples with CSOM with cholesteaotoma (1.6 ± 0.7 log10) was significantly higher compared to healthy controls (0.3 ± 1.6). These results were corroborated with Gram-staining that identified the large presence of Gram-positive cocci cells in the cholesteatoma tissue of CSOM subjects which were not seen in the mucosa of controls. The most abundant genus detected by sequencing in the mucosa and cholesteatoma of CSOM samples was Anaerococcus (93.5% of all reads), and genus Meiothermus (0.9%) in the control sample. The 3 samples with the highest sequencing reads (>300) were further analysed using multiplex PCR to identify the dominant Anaerococcus species. Anaerococcus hydrogenalis was the dominant species identified in these samples. In contract, commonly named ear pathogens, genera Staphylococcus and Pseudomonas, were detected in low numbers (<0.001% of all sequencing reads) and low prevalence (2/16 samples) in the tissue samples of this study.The results show that culture severely underestimated the bacterial diversity in CSOM samples and investigating tissue rather than standard culture swabs might be advantageous to understanding the disease process. The high abundance of bacteria and the large presence of Gram-positive cells detected in the cholesteatoma tissue of CSOM compared to mucosa of CSOM or controls could be members from the genus Anaerococcus. Anaerococcus may well be a pathogen in CSOM with cholesteatoma, but their role in this condition requires further investigation.