Sequence Comparison Research Articles

The soybean Rpp3 gene encodes a TIR-NBS-LRR protein that confers resistance to Phakopsora pachyrhizi.

Soybean rust is an economically significant disease caused by the fungus Phakopsora pachyrhizi that negatively impacts soybean (Glycine max (L.) Merr.) production throughout the world. Susceptible plants infected by P. pachyrhizi develop tan-colored lesions on the leaf surface that give rise to funnel-shaped uredinia as the disease progresses. While most soybean germplasm is susceptible, seven genetic loci (Rpp1 to Rpp7) that provide race-specific resistance to P. pachyrhizi (Rpp) have been identified. Rpp3 was first discovered and characterized in the soybean accession PI 462312 (Ankur), and it was also determined to be one of two Rpp genes present in PI 506764 (Hyuuga). Genetic crosses with PI 506764 were later used to fine-map the Rpp3 locus to a 371 kb region on chromosome 6. The corresponding region in the susceptible Williams 82 (Wm82) reference genome contains several homologous nucleotide binding site-leucine rich repeat (NBS-LRR) genes. To identify Rpp3, we designed oligonucleotide primers to amplify Rpp3 candidate (Rpp3C) NBS-LRR genes at this locus from PI 462312, PI 506764, and Wm82 using polymerase chain reaction (PCR). Five Rpp3C genes were identified in both Rpp3-resistant soybean lines, and co-silencing these genes compromised resistance to P. pachyrhizi. Gene expression analysis and sequence comparisons of the Rpp3C genes in PI 462312 and PI 506764 suggest that a single candidate gene, Rpp3C3, is responsible for Rpp3-mediated resistance.

Characterization of Plasmids Isolated from Bacterial Flora of CRC Patients – Proof of Concept

Abstract Introduction Colorectal cancer (CRC) is one of the leading oncogenic disorders, both in terms of incidence and mortality. The etiology of the disease is certainly multifactorial. Various risk factors like alcohol consumption, smoking, CRC family history, inflammatory bowel disease, hormone therapy, aspirin/nonsteroidal anti-inflammatory drugs use, higher body mass index, consumption of red and/or processed meat, insufficient physical activity, and decreased intake of fruit and vegetables have been pointed out; however, there is not enough support evidence for a single particular causative mechanism. Recently, gut bacterial microbiota has been shown to influence significantly the pathogenesis of CRC. However, little attention is paid to the putative impact of plasmids in gut flora. Material and methods We have designed and tested the workflow for semi-selective isolation and amplification of random circular sequences. The exploitation of rolling circle amplification (RCA) with a random hexamers protocol is crucial for the outcome. Results Our results suggest that it is possible to isolate and amplify plasmid DNA from gut flora and further process, sequence, and identify them. Discussion Little is known about the interactions between bacterial plasmids and human cells. The collection of plasmid sequencing data and the comparison of CRC patients and healthy control sequences can be the first step to elucidating this phenomenon.

Abundance and diversity of fungal endophytes isolated from monk fruit (Siraitia grosvenorii) grown in a Canadian research greenhouse.

Monk fruit (Siraitia grosvenorii) is an herbaceous perennial vine of the Cucurbitaceae family cultivated commercially mainly in southern China. There is very little information available about the fungal endophytes in monk fruit. In this study, monk fruit plants were grown from seeds in a research greenhouse at Kwantlen Polytechnic University in British Columbia, Canada to explore the abundance and diversity of their fungal endophytes. Fungal endophytes were isolated from seeds, seedlings, mature monk fruit plants, and fruits, and cultured on potato dextrose agar and water agar media. Isolates were identified by microscopic examination and BLAST comparison of ITS sequences to published sequences in GenBank. At least 150 species of fungal endophytes representing 60 genera and 20 orders were recovered from monk fruit tissues. Non-metric multidimensional scaling (NMDS) was carried out to explore the similarity of fungal communities among roots, stems, leaves, flowers, fruits, and seeds based on fungal orders. Our study showed that monk fruit plants are a rich source of fungal endophytes with the greatest abundance and diversity in leaves. This work has deepened our understanding of the intricate interactions between plants and fungi that sustain ecosystems and underpin plant health and resilience.

Porphyromonas levii – A Boon for Periodontists?

ABSTRACT The majority of species previously categorized as Bacteroides have been reassigned into new genera. Bacteroides levii (Holdeman, Cato, and Mooretaxonomic)’s status has remained uncertain. This species shares a high degree of similarity with members of the genus Porphyromonas based on biochemical, chemical, and comparative 16s rRNA sequence analysis. As a result, Bacteroides levii (Holdeman, Cato, and Moore) was reclassified as Porphyromonas levii comb. now under the genus Porphyromonas.

Comparative single-cell RNA sequencing analysis of immune response to inactivated vaccine and natural SARS-CoV-2 infection.

Uncovering the immune response to an inactivated SARS-CoV-2 vaccine (In-Vac) and natural infection is crucial for comprehending COVID-19 immunology. Here we conducted an integrated analysis of single-cell RNA sequencing (scRNA-seq) data from serial peripheral blood mononuclear cell (PBMC) samples derived from 12 individuals receiving In-Vac compared with those from COVID-19 patients. Our study reveals that In-Vac induces subtle immunological changes in PBMC, including cell proportions and transcriptomes, compared with profound changes for natural infection. In-Vac modestly upregulates IFN-α but downregulates NF-κB pathways, while natural infection triggers hyperactive IFN-α and NF-κB pathways. Both In-Vac and natural infection alter T/B cell receptor repertoires, but COVID-19 has more significant change in preferential VJ gene, indicating a vigorous immune response. Our study reveals distinct patterns of cellular communications, including a selective activation of IL-15RA/IL-15 receptor pathway after In-Vac boost, suggesting its potential role in enhancing In-Vac-induced immunity. Collectively, our study illuminates multifaceted immune responses to In-Vac and natural infection, providing insights for optimizing SARS-CoV-2 vaccine efficacy.

A hypervirulent Acinetobacter baumannii strain has robust anti-phagocytosis ability

BackgroundAcinetobacter baumannii (A. baumannii) is associated with both hospital-acquired infections (HAP) and community-acquired pneumonia (CAP). In this study, we present a novel CAP-associated A. baumannii (CAP-AB) strain causing severe pneumonia in an afore healthy male patient without underlying conditions. Subsequently, we investigated the pathogenicity and immunogenicity of this CAP-AB strain using a mice pneumonia model.ResultsA 58-year-old male patient with no underlying conditions experienced worsening symptoms of a productive cough, sputum, and fever that developed acutely, in just 24 h. The diagnosis was severe community-acquired pneumonia (CAP) and type-1 respiratory failure. An A. baumannii strain was isolated from his sputum and blood cultures. To gain a deeper understanding of the rapid progression of its pathology, we utilized the CAP-associated A. baumannii strain YC128, a previously obtained hospital-acquired pneumonia A. baumannii (HAP-AB) strain YC156, and a highly virulent A. baumannii control strain LAC-4 to construct a mouse pneumonia model, and subsequently compared the mortality rate of the three groups. Following inoculation with 107 CFU of A. baumannii, the mortality rate for the YC128, LAC-4, and YC156 groups was 60% (6/10), 30% (3/10), and 0%, respectively. The bacterial burden within the pulmonary, liver, and spleen tissues of mice in the YC128 group was significantly higher than that of the YC156 group, and slightly higher than that of the LAC-4 group. Pathological analysis of lung tissue using HE-staining revealed that the inflammatory pathological changes in mice from the YC128 group were significantly more severe than those in the YC156 group. Additionally, CT scan images displayed more pronounced inflammation in the lungs of mice from the YC128 group compared to the YC156 group. Local levels of cytokines/chemokines such as IL-1β, IL-6, TNF-α, and CXCL1 were assessed via RT-qPCR in lung tissues. In comparison with the YC156 strain, the highly virulent YC128 strain induced the expression of proinflammatory cytokines more rapidly and severely. Furthermore, we examined the in vitro anti-phagocytosis ability of YC128 and YC156 strains against mice peritoneal macrophages, revealing that the highly virulent YC128 isolate displayed greater resistance to macrophage uptake in contrast to YC156. Results from Whole Genome Sequencing (WGS) indicated that YC128 harbored a complete type VI secretion system (T6SS) gene cluster, while YC156 lacked the majority of genes within the T6SS gene cluster. The other virulence-related genes exhibited minimal differences between YC128 and YC156. Drawing from previous studies, we postulated that the T6SS is linked to the hypervirulence and robust anti-phagocytic ability of YC128.ConclusionsThis article reports on the isolation of a novel hypervirulent CAP-AB strain, YC128, from a severe CAP patient. The results demonstrate that this CAP-AB strain, YC128, is capable of inducing fatal pneumonia and extrapulmonary dissemination in a mouse pneumonia model. Moreover, this highly virulent CAP-AB strain exhibits significantly stronger anti-phagocytic abilities compared to the HAP-AB YC156 strain. Genome sequencing comparisons reveal that the heightened hypervirulence and enhanced anti-phagocytosis abilities observed in YC128 may be attributed to the presence of the T6SS.

First Report of Powdery Mildew Caused by Golovinomyces adenophorae on Adenophora triphylla var. japonica in Korea.

Adenophora triphylla var. japonica (Campanulaceae), known as Japanese lady bell, is native to East Asia. It has been used as a medicinal plant but is widely cultivated in Korea as an indigenous vegetable (Park et al. 2011). In the summer of 2020, about 100 plants in an experimental plot at the National Institute of Forest Science, Seoul, Korea, showed powdery mildew symptoms with a 100% disease incidence. Signs first appeared as white colonies, subsequently expanding over the leaves, stems, and inflorescences. Infected young shoots were elongated and became slender. Chasmothecia were found in late October. Voucher specimens were deposited in the Korea University Herbarium (KUS-F). Conidiophores arising from the lateral part of the hyphae were upright, 100 to 220 × 10 to 12 µm, and produced 2 to 5 immature conidia in chains with sinuate edge lines. Basal parts of foot-cells in conidiophores were curved. Conidia were barrel-shaped to ellipsoid, 26 to 40 × 14 to 20 µm, and produced germ tubes on the perihilar position of the conidia. Chasmothecia with short mycelioid appendages were gregarious, 144 to176 µm in diam., and contained 8 to 22 asci. Asci were clavate-saccate with short stalks, 60 to 82 × 28 to 42 µm, and contained two spores. Ascospores were broadly ellipsoid, cytoplasm-dense without vacuoles, colorless, and 22 to 28 × 12 to 18 µm. The structures and measurements were consistent with those of Golovinomyces adenophorae (R.Y. Zheng & G.Q. Chen) Heluta (Braun & Cook, 2012). To confirm the morphology-based identification, two herbarium specimens (KUS-F29252 and F31898) were sequenced for the internal transcribed spacer (ITS) and large subunit (LSU) regions with PM10/ITS4 and PM3/TW14 primers, respectively (Bradshaw and Tobin, 2020). A Blastn search revealed high similarities in the ITS and LSU sequences, with 99.81% (538/539 bp) and 99.86% (697/698 bp) to G. adenophorae sequences (AB077633 and AB077632), respectively. All resulting sequences were deposited in GenBank under accession numbers OR841069-70 for ITS and OR841071 for LSU. A pathogenicity test was performed through inoculation by gently dusting the conidia from a detached symptomatic leaf onto the leaves of five healthy plants. Five non-inoculated plants served as controls. Following inoculation, plants were covered with plastic film and maintained in a greenhouse (24 to 32°C) until symptoms developed. Powdery mildew colonies developed on the inoculated plants after twelve days, whereas the control plants remained symptomless. The inoculated pathogen was confirmed morphologically and molecularly by the sequence comparison aforementioned, fulfilling Koch's postulates. Based on morphological characteristics and the sequencing data, the powdery mildew was identified as G. adenophorae. The association of G. adenophorae and Adenophora spp. has been known in China, Japan, Kazakhstan, and the Far East of Russia (Farr and Rossman, 2023). This is the first report of powdery mildew caused by G. adenophorae on A. triphylla var. japonica in Korea. Since the commercial cultivation of this plant aims to harvest young shoots as one of the most popular vegetables in Korea, appropriate control measures for the powdery mildew should be considered.

Viral metagenomic analysis reveals diverse viruses and a novel bocaparvovirus in the enteric virome of snow leopard (Panthera uncia)

The enteric virome, comprising a complex community of viruses inhabiting the gastrointestinal tract, plays a significant role in health and disease dynamics. In this study, the fecal sample of a wild snow leopard was subjected to viral metagenomic analysis using a double barcode Illumina MiSeq platform. The resulting reads were de novo assembled into contigs with SOAPdenovo2 version r240. Additional bioinformatic analysis of the assembled genome and genome annotation was done using the Geneious prime software (version 2022.0.2). Following viral metagenomic analysis and bioinformatic analysis, a total of 7 viral families and a novel specie of bocaparvovirus tentatively named Panthera uncia bocaparvovirus (PuBOV) with GenBank accession number OQ627713 were identified. The complete genome of PuBOV was predicted to contain 3 open reading frames (ORFs), contains 5433 nucleotides and has a G + C content of 47.40 %. BLASTx analysis and pairwise sequence comparison indicated the novel virus genome was a new species in the genus Bocaparvovirus based on the species demarcation criteria of the International Committee on the Taxonomy of Viruses. This study provides valuable insights into the diversity and composition of the enteric virome in wild endangered snow leopards. The identification and characterization of viruses in wildlife is crucial for developing effective strategies to manage and mitigate potential zoonotic and other viral disease threats to human and animal health.

A conserved NR5A1-responsive enhancer regulates SRY in testis-determination

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5’ regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

Genetic diversity of mitochondrial DNA D-loop sequence in Bang Troi chicken breed

Bang Troi chicken (GBT) is a local chicken that is free-gardening raised in Quang Ninh province. A 1,268 bp of mitochondrial DNA (mtDNA) D-loop of 17 Bang Troi chicken individuals was sequenced and analyzed. Comparative multiple sequence alignment revealed 17 base substitutions within 17 Bang Troi chicken individuals. Bang Troi chicken had a relatively high level of genetic diversity, expressed in the number of 13 haplotypes, haplotype diversity of 0.949. The Tajima’s test indicated that the analysed population in balancing selection with a D-value of 2.39. A phylogenetic tree based on 455 bp of hypervariable mtD-loop sequence classified 17 Bang Troi chicken individuals into three clades A, B, and E, which are common groups of Vietnamese and Asia indigenous chickens. This result provides the first genetic information and maternal lineages of Bang Troi chickens and it can support for registration of local chicken breeds.

Molecular characterization of emerging distinct watermelon bud necrosis virus (Orthotospovirus citrullonecrosis) variant causing necrosis disease on tomato in India

Watermelon bud necrosis virus (WBNV) cause serious economic impact on the cultivation of vegetable crops in the Indian subcontinent. In this study, an Orthotospovirus that infect tomato plants with typical severe necrotic symptoms had been identified as WBNV based on N gene sequences. The complete genome sequence of WBNV isolate (22ToNV) was determined and analyzed. The nucleotide sequences of the viral RNA segments in 22ToNV exhibit <91% similarity compared to the published WBNV genome sequence. Whereas, the amino acid sequences encoded by these segments, shared identity of about 96–98.54% among the other WBNV isolates. There is a notable divergence observed specifically within the intergenic regions (IGRs) of the S and M RNAs, with a lower similarity range of 78–84%. Nucleotide sequence comparisons and phylogenetic analysis showed that 22ToNV is a unique WBNV isolate having distinct genetic differences with other known WBNV isolates reported in different geographic regions of the country.

Persistent homology reveals strong phylogenetic signal in 3D protein structures.

Changes that occur in proteins over time provide a phylogenetic signal that can be used to decipher their evolutionary history and the relationships between organisms. Sequence comparison is the most common way to access this phylogenetic signal, while those based on 3D structure comparisons are still in their infancy. In this study, we propose an effective approach based on Persistent Homology Theory (PH) to extract the phylogenetic information contained in protein structures. PH provides efficient and robust algorithms for extracting and comparing geometric features from noisy datasets at different spatial resolutions. PH has a growing number of applications in the life sciences, including the study of proteins (e.g. classification, folding). However, it has never been used to study the phylogenetic signal they may contain. Here, using 518 protein families, representing 22,940 protein sequences and structures, from 10 major taxonomic groups, we show that distances calculated with PH from protein structures correlate strongly with phylogenetic distances calculated from protein sequences, at both small and large evolutionary scales. We test several methods for calculating PH distances and propose some refinements to improve their relevance for addressing evolutionary questions. This work opens up new perspectives in evolutionary biology by proposing an efficient way to access the phylogenetic signal contained in protein structures, as well as future developments of topological analysis in the life sciences.

Comparison of polymerase chain reaction and next-generation sequencing with conventional urine culture for the diagnosis of urinary tract infections: A meta-analysis.

The limitations of conventional urine culture methods can be avoided by using culture-independent approaches like polymerase chain reaction (PCR) and next-generation sequencing (NGS). However, the efficacy of these approaches in this setting is still subject to contention. PRISMA-compliant searches were performed on MEDLINE/PubMed, EMBASE, Web of Sciences, and the Cochrane Database until March 2023. The included articles compared PCR or NGS to conventional urine culture for the detection of urinary tract infections (UTIs). RevMan performed meta-analysis, and the Cochrane Risk of Bias Assessment Tool assessed study quality. A total of 10 selected studies that involved 1,291 individuals were included in this meta-analysis. The study found that PCR has a 99% sensitivity and a 94% specificity for diagnosing UTIs. Furthermore, NGS was shown to have a sensitivity of 90% for identifying UTIs and a specificity of 86%. The odds ratio (OR) for PCR to detect Gram-positive bacteria is 0.50 (95% confidence interval [CI] 0.41-0.61), while the OR for NGS to detect Gram-negative bacteria is 0.23 [95% CI 0.09-0.59]. UTIs are typically caused by Gram-negative bacteria like Escherichia coli and Gram-positive bacteria like Staphylococci and Streptococci. PCR and NGS are reliable, culture-free molecular diagnostic methodsthat, despite being expensive, are essential for UTI diagnosis and prevention due to their high sensitivity and specificity.

Multi-omics analysis provides new insights into the changes of important nutrients and fructose metabolism in loquat bud sport mutant.

Bud sport is a common and stable somatic variation in perennial fruit trees, and often leads to significant modification of fruit traits and affects the breeding value. To investigate the impact of bud sport on the main metabolites in the fruit of white-fleshed loquat, we conducted a multi-omics analysis of loquat fruits at different developmental stages of a white-fleshed bud sport mutant of Dongting loquat (TBW) and its wild type (TBY). The findings from the detection of main fruit quality indices and metabolites suggested that bud sport resulted in a reduction in the accumulation of carotenoids, fructose, titratable acid and terpenoids at the mature stage of TBW, while leading to the accumulation of flavonoids, phenolic acids, amino acids and lipids. The comparably low content of titratable acid further enhances the balanced and pleasent taste profile of TBW. Expression patterns of differentially expressed genes involved in fructose metabolism exhibited a significant increase in the expression level of S6PDH (EVM0006243, EVM0044405) prior to fruit maturation. The comparison of protein sequences and promoter region of S6PDH between TBY and TBW revealed no structural variations that would impact gene function or expression, indicating that transcription factors may be responsible for the rapid up-regulation of S6PDH before maturation. Furthermore, correlation analysis helped to construct a comprehensive regulatory network of fructose metabolism in loquat, including 23 transcription factors, six structural genes, and nine saccharides. Based on the regulatory network and existing studies, it could be inferred that transcription factors such as ERF, NAC, MYB, GRAS, and bZIP may promote fructose accumulation in loquat flesh by positively regulating S6PDH. These findings improve our understanding of the nutritional value and breeding potential of white-fleshed loquat bud sport mutant, as well as serve as a foundation for exploring the genes and transcription factors that regulate fructose metabolism in loquat.

Plasticity of repetitive sequences demonstrated by the complete mitochondrial genome of Eucalyptus camaldulensis.

The tree Eucalyptus camaldulensis is a ubiquitous member of the Eucalyptus genus, which includes several hundred species. Despite the extensive sequencing and assembly of nuclear genomes from various eucalypts, the genus has only one fully annotated and complete mitochondrial genome (mitogenome). Plant mitochondria are characterized by dynamic genomic rearrangements, facilitated by repeat content, a feature that has hindered the assembly of plant mitogenomes. This complexity is evident in the paucity of available mitogenomes. This study, to the best of our knowledge, presents the first E. camaldulensis mitogenome. Our findings suggest the presence of multiple isomeric forms of the E. camaldulensis mitogenome and provide novel insights into minor rearrangements triggered by nested repeat sequences. A comparative sequence analysis of the E. camaldulensis and E. grandis mitogenomes unveils evolutionary changes between the two genomes. A significant divergence is the evolution of a large repeat sequence, which may have contributed to the differences observed between the two genomes. The largest repeat sequences in the E. camaldulensis mitogenome align well with significant yet unexplained structural variations in the E. grandis mitogenome, highlighting the adaptability of repeat sequences in plant mitogenomes.

Comparison of metagenomic next-generation sequencing and conventional culture for the diagnostic performance in febrile patients with suspected infections

BackgroundTimely and accurate identification of pathogens is crucial for appropriate treatment and prognosis of infectious diseases. As an increasingly popular pathogen detection method, the performance of metagenomic next-generation sequencing (mNGS) in detecting pathogens in febrile patients with suspected infection requires further exploration.MethodsThis study included 368 febrile patients with suspected infections who were admitted to the Infectious Disease Department of Qilu Hospital, Shandong University between January 5, 2021 and April 14, 2023. Both mNGS testing and conventional culture were performed in all patients. Clinical data of enrolled patients were collected, and the diagnostic performances of mNGS and culture were compared.ResultsOf the 368 enrolled patients, 231 were finally diagnosed with infection and 137 were with diseases other than infection. The sensitivity (58.01% vs. 21.65%, p < 0.001) and negative predictive value (54.67% vs. 42.9%) of mNGS were superior to those of culture. In contrast, the culture exhibited higher specificity (99.27% vs. 85.40%, p < 0.001) and positive predictive value (98.84% vs. 87.01%) than mNGS. Among infected patients with positive mNGS results, 64 received adjusted antibiotic therapy including treatment transitions, antibiotic downgrading, and combination therapy. Among them, 9 had additional antifungal drugs and 21 patients had a treatment turning point based on the mNGS results and these patients recovered and discharged due to timely antibiotic adjustment. Both positive rates of puncture fluid mNGS and tissue mNGS were higher than those of culture in the patients who had prior antibiotic use, and this difference was statistically significant (p = 0.000).ConclusionmNGS is more sensitive and accurate than traditional culture, making it ideal for identifying pathogens and screening infectious diseases, especially for those with uncultivated or difficult-to-cultivate species. Early diagnosis allows for prompt treatment with targeted antibiotics, and mNGS is recommended when samples are limited.

COSAP: Comparative Sequencing Analysis Platform

BackgroundRecent improvements in sequencing technologies enabled detailed profiling of genomic features. These technologies mostly rely on short reads which are merged and compared to reference genome for variant identification. These operations should be done with computers due to the size and complexity of the data. The need for analysis software resulted in many programs for mapping, variant calling and annotation steps. Currently, most programs are either expensive enterprise software with proprietary code which makes access and verification very difficult or open-access programs that are mostly based on command-line operations without user interfaces and extensive documentation. Moreover, a high level of disagreement is observed among popular mapping and variant calling algorithms in multiple studies, which makes relying on a single algorithm unreliable. User-friendly open-source software tools that offer comparative analysis are an important need considering the growth of sequencing technologies.ResultsHere, we propose Comparative Sequencing Analysis Platform (COSAP), an open-source platform that provides popular sequencing algorithms for SNV, indel, structural variant calling, copy number variation, microsatellite instability and fusion analysis and their annotations. COSAP is packed with a fully functional user-friendly web interface and a backend server which allows full independent deployment for both individual and institutional scales. COSAP is developed as a workflow management system and designed to enhance cooperation among scientists with different backgrounds. It is publicly available at https://cosap.bio and https://github.com/MBaysanLab/cosap/. The source code of the frontend and backend services can be found at https://github.com/MBaysanLab/cosap-webapi/ and https://github.com/MBaysanLab/cosap_frontend/ respectively. All services are packed as Docker containers as well. Pipelines that combine algorithms can be customized and new algorithms can be added with minimal coding through modular structure.ConclusionsCOSAP simplifies and speeds up the process of DNA sequencing analyses providing commonly used algorithms for SNV, indel, structural variant calling, copy number variation, microsatellite instability and fusion analysis as well as their annotations. COSAP is packed with a fully functional user-friendly web interface and a backend server which allows full independent deployment for both individual and institutional scales. Standardized implementations of popular algorithms in a modular platform make comparisons much easier to assess the impact of alternative pipelines which is crucial in establishing reproducibility of sequencing analyses.

Using a Handful of Transcriptomes to Detect Sex-Linked Markers and Develop Molecular Sexing Assays in a Species with Homomorphic Sex Chromosomes.

To understand the biology of a species, it is often crucial to be able to differentiate males and females. However, many species lack easily identifiable sexually dimorphic traits. In those that possess sex chromosomes, molecular sexing offers a good alternative, and molecular sexing assays can be developed through the comparison of male and female genomic sequences. However, in many nonmodel species, sex chromosomes are poorly differentiated, and identifying sex-linked sequences and developing sexing assays can be challenging. In this study, we highlight a simple transcriptome-based procedure for the detection of sex-linked markers suitable for the development of sexing assays that circumvents limitations of more commonly used approaches. We apply it to the spotted snow skink Carinascincus ocellatus, a viviparous lizard with homomorphic XY chromosomes that has environmentally induced sex reversal. With transcriptomes from three males and three females alone, we identify thousands of putative Y-linked sequences. We confirm linkage through alignment of assembled transcripts to a distantly related lizard genome and readily design multiple single locus polymerase chain reaction primers to sex C. ocellatus and related species. Our approach also facilitates valuable comparisons of sex determining systems on a broad taxonomic scale.

Comparison of MRI Sequences to Predict IDH Mutation Status in Gliomas Using Radiomics-Based Machine Learning.

Regarding the 2021 World Health Organization (WHO) classification of central nervous system (CNS) tumors, the isocitrate dehydrogenase (IDH) mutation status is one of the most important factors for CNS tumor classification. The aim of our study is to analyze which of the commonly used magnetic resonance imaging (MRI) sequences is best suited to obtain this information non-invasively using radiomics-based machine learning models. We developed machine learning models based on different MRI sequences and determined which of the MRI sequences analyzed yields the highest discriminatory power in predicting the IDH mutation status. In our retrospective IRB-approved study, we used the MRI images of 106 patients with histologically confirmed gliomas. The MRI images were acquired using the T1 sequence with and without administration of a contrast agent, the T2 sequence, and the Fluid-Attenuated Inversion Recovery (FLAIR) sequence. To objectively compare performance in predicting the IDH mutation status as a function of the MRI sequence used, we included only patients in our study cohort for whom MRI images of all four sequences were available. Seventy-one of the patients had an IDH mutation, and the remaining 35 patients did not have an IDH mutation (IDH wild-type). For each of the four MRI sequences used, 107 radiomic features were extracted from the corresponding MRI images by hand-delineated regions of interest. Data partitioning into training data and independent test data was repeated 100 times to avoid random effects associated with the data partitioning. Feature preselection and subsequent model development were performed using Random Forest, Lasso regression, LDA, and Naïve Bayes. The performance of all models was determined with independent test data. Among the different approaches we examined, the T1-weighted contrast-enhanced sequence was found to be the most suitable for predicting IDH mutations status using radiomics-based machine learning models. Using contrast-enhanced T1-weighted MRI images, our seven-feature model developed with Lasso regression achieved a mean area under the curve (AUC) of 0.846, a mean accuracy of 0.792, a mean sensitivity of 0.847, and a mean specificity of 0.681. The administration of contrast agents resulted in a significant increase in the achieved discriminatory power. Our analyses show that for the prediction of the IDH mutation status using radiomics-based machine learning models, among the MRI images acquired with the commonly used MRI sequences, the contrast-enhanced T1-weighted images are the most suitable.

Complete genome characterization of cacao leafroll virus, a newly described cacao-infecting polerovirus.

The complete genome sequence of cacao leafroll virus (CaLRV; family Solemoviridae, genus Polerovirus) was determined by high-throughput sequencing of total RNA isolated from symptomatic cacao Theobroma cacao L. plants (n = 4). The CaLRV genome sequences ranged from 5,976 to 5,997 nucleotides (nt) in length and contained seven open reading frames (ORFs). Nucleotide and amino acid (aa) sequence comparisons showed that, among selected well-characterized poleroviruses, the CaLRV genome shared the highest nt sequence identity of 62% with that of potato leafroll virus (PLRV, NC_076505). A comparison of the predicted aa sequence of the CaLRV coat protein indicated that cotton leafroll dwarf virus (CLRDV, NC_014545) and melon aphid-borne yellows virus (MABYV, NC_010809) were the closest relatives, sharing 57% aa sequence identity. Bayesian phylogenetic analysis based on complete genome sequences showed that CaLRV grouped with well-characterized poleroviruses that cause diseases of cereal and vegetable crops. During the course of publishing this work, the nearly complete genome sequence of a member of the same polerovirus species, referred to as "cacao polerovirus" (OR605721), with which CaLRV shares 99% nt sequence identity, was reported.