Ultra-deep sequencing enables high-fidelity recovery of biodiversity for bulk arthropod samples without PCR amplification

BACKGROUND

Next-generation-sequencing (NGS) technologies combined with a classic DNA barcoding approach have enabled fast and credible measurement for biodiversity of mixed environmental samples. However, the PCR amplification involved in nearly all existing NGS protocols inevitably introduces taxonomic biases. In the present study, we developed new Illumina pipelines without PCR amplifications to analyze terrestrial arthropod communities.

RESULTS

Mitochondrial enrichment directly followed by Illumina shotgun sequencing, at an ultra-high sequence volume, enabled the recovery of Cytochrome c Oxidase subunit 1 (COI) barcode sequences, which allowed for the estimation of species composition at high fidelity for a terrestrial insect community. With 15.5 Gbp Illumina data, approximately 97% and 92% were detected out of the 37 input Operational Taxonomic Units (OTUs), whether the reference barcode library was used or not, respectively, while only 1 novel OTU was found for the latter. Additionally, relatively strong correlation between the sequencing volume and the total biomass was observed for species from the bulk sample, suggesting a potential solution to reveal relative abundance.

CONCLUSIONS

The ability of the new Illumina PCR-free pipeline for DNA metabarcoding to detect small arthropod specimens and its tendency to avoid most, if not all, false positives suggests its great potential in biodiversity-related surveillance, such as in biomonitoring programs. However, further improvement for mitochondrial enrichment is likely needed for the application of the new pipeline in analyzing arthropod communities at higher diversity.

Application of genomic tools in plant breeding

Plant breeding has been very successful in developing improved varieties using conventional tools and methodologies. Nowadays, the availability of genomic tools and resources is leading to a new revolution of plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. Next Generation Sequencing (NGS) technologies are allowing the mass sequencing of genomes and transcriptomes, which is producing a vast array of genomic information. The analysis of NGS data by means of bioinformatics developments allows discovering new genes and regulatory sequences and their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Genomic approaches include TILLING and EcoTILLING, which make possible to screen mutant and germplasm collections for allelic variants in target genes. Re-sequencing of genomes is very useful for the genome-wide discovery of markers amenable for high-throughput genotyping platforms, like SSRs and SNPs, or the construction of high density genetic maps. All these tools and resources facilitate studying the genetic diversity, which is important for germplasm management, enhancement and use. Also, they allow the identification of markers linked to genes and QTLs, using a diversity of techniques like bulked segregant analysis (BSA), fine genetic mapping, or association mapping. These new markers are used for marker assisted selection, including marker assisted backcross selection, ‘breeding by design’, or new strategies, like genomic selection. In conclusion, advances in genomics are providing breeders with new tools and methodologies that allow a great leap forward in plant breeding, including the ‘superdomestication’ of crops and the genetic dissection and breeding for complex traits.

Evaluation of methodologies for microRNA biomarker detection by next generation sequencing

In recent years, microRNAs (miRNAs) in tissues and biofluids have emerged as a new class of promising biomarkers for numerous diseases. Blood-based biomarkers are particularly desirable since serum or plasma is easily accessible and can be sampled repeatedly. To comprehensively explore the biomarker potential of miRNAs, sensitive, accurate and cost-efficient miRNA profiling techniques are required. Next generation sequencing (NGS) is emerging as the preferred method for miRNA profiling; offering high sensitivity, single-nucleotide resolution and the possibility to profile a considerable number of samples in parallel. Despite the excitement about miRNA biomarkers, challenges associated with insufficient characterization of the sequencing library preparation efficacy, precision and method-related quantification bias have not been addressed in detail and are generally underappreciated in the wider research community. Here, we have tested in parallel four commercially available small RNA sequencing kits against a cohort of samples comprised of human plasma, human serum, murine brain tissue and a reference library containing ~ 950 synthetic miRNAs. We discuss the advantages and limits of these methodologies for massive parallel microRNAs profiling. This work can serve as guideline for choosing an adequate library preparation method, based on sensitivity, specificity and accuracy of miRNA quantification, workflow convenience and potential for automation.

Let-7a-5p, miR-100-5p, miR-101-3p, and miR-199a-3p Hyperexpression as Potential Predictive Biomarkers in Early Breast Cancer Patients

BACKGROUND

The aim of this study is to identify miRNAs able to predict the outcomes in breast cancer patients after neoadjuvant chemotherapy (NAC).

PATIENTS AND METHODS

We retrospectively analyzed 24 patients receiving NAC and not reaching pathologic complete response (pCR). miRNAs were analyzed using an Illumina Next-Generation-Sequencing (NGS) system.

RESULTS

Event-free survival (EFS) and overall survival (OS) were significantly higher in patients with up-regulation of let-7a-5p (EFS p = 0.006; OS p = 0.0001), mirR-100-5p (EFS s p = 0.01; OS p = 0.03), miR-101-3p (EFS p = 0.05; OS p = 0.01), and miR-199a-3p (EFS p = 0.02; OS p = 0.01) in post-NAC samples, independently from breast cancer subtypes. At multivariate analysis, only let-7a-5p was significantly associated with EFS (p = 0.009) and OS (p = 0.0008).

CONCLUSION

Up-regulation of the above miRNAs could represent biomarkers in breast cancer.

Non-Invasive Screening Tools for Down's Syndrome: A Review

Down's syndrome (DS) is the most common genetic cause of developmental delay with an incidence of 1 in 800 live births, and is the predominant reason why women choose to undergo invasive prenatal diagnosis. However, as invasive tests are associated with around a 1% risk of miscarriage new non-invasive tests have been long sought after. Recently, the most promising approach for non-invasive prenatal diagnosis (NIPD) has been provided by the introduction of next generation sequencing (NGS) technologies. The clinical application of NIPD for DS detection is not yet applicable, as large scale validation studies in low-risk pregnancies need to be completed. Currently, prenatal screening is still the first line test for the detection of fetal aneuploidy. Screening cannot diagnose DS, but developing a more advanced screening program can help to improve detection rates, and therefore reduce the number of women offered invasive tests. This article describes how the prenatal screening program has developed since the introduction of maternal age as the original "screening" test, and subsequently discusses recent advances in detecting new screening markers with reference to both proteomic and bioinformatic techniques.

Complete chloroplast genome of the orchid Cattleya crispata (Orchidaceae:Laeliinae), a Neotropical rupiculous species

A partial genome dataset was sequenced for the orchid Cattleya crispata using both Illumina and 454 technologies. The chloroplast genome was assembled using iterative runs of MIRA software that yielded a circular molecule with 148,343 bp in length and deposited in GenBank database (Accession Number KP168671). The plastid genome conserved the quadripartite structure present in most Orchidaceae chloroplasts and was composed by 79 protein-coding genes, 39 tRNAs and 8 rRNAs. Genome structure, gene order and orientation were similar to previously described chloroplasts for Cymbidium orchids, differing in gene order for petN and psbM genes. Data described here contain the first report of a complete chloroplast for the Neotropical subtribe Laeliinae and may contribute to improve the phylogenetic resolution and allow the development of new molecular markers for population genetic studies of orchids.

Evaluating the Effect of QIIME Balanced Default Parameters on Metataxonomic Analysis Workflows With a Mock Community

Metataxonomic analysis represents a fast and cost-effective approach for acquiring informative insight into the composition of the microbiome of samples with variable diversity, such as wine samples. Nevertheless, it comprises a vast amount of laboratory procedures and bioinformatic frameworks each one associated with an inherent variability of protocols and algorithms, respectively. As a solution to the bioinformatic maze, QIIME bioinformatic framework has incorporated benchmarked, and balanced parameters as default parameters. In the current study, metataxonomic analysis of two types of mock community standards with the same microbial composition has been performed for evaluating the effectivess of QIIME balanced default parameters on a variety of aspects related to different laboratory and bioinformatic workflows. These aspects concern NGS platforms, PCR protocols, bioinformatic pipelines, and taxonomic classification algorithms. Several qualitative performance expectations have been the outcome of the analysis, rendering the mock community a useful evaluation tool.

The complete mitochondrial genome of the southern purple-spotted gudgeon Mogurnda adspersa (Perciformes: Eleotridae) through pyrosequencing

The complete mitochondrial genome of the Critically Endangered southern purple-spotted gudgeon, Mogurnda adspersa, was determined for the first time using 1/8 of a 454 pyrosequencing plate. The mitogenome was assembled using the bioinformatic software MIRA. The M. adspersa genome organization was very similar to most vertebrates, being 16,523 bp in length. It contained 13 protein-coding genes, 22 transfer RNA, 2 ribosomal RNA genes, and 1 non-coding control region. The current methodology was effective to assemble the whole mitogenome using 2008 mitochondrial reads representing 0.23% of total genomic reads produced and providing average mitogenome coverage of 17.4 reads per site. The whole mitogenome sequence provided here is deposited in NCBI with the accession number KJ130031 and may benefit systematics studies and conservation programs of M. adspersa.

Variable In Vivo Hepatitis D Virus (HDV) RNA Editing Rates According to the HDV Genotype

Human hepatitis delta virus (HDV) is a small defective RNA satellite virus that requires hepatitis B virus (HBV) envelope proteins to form its own virions. The HDV genome possesses a single coding open reading frame (ORF), located on a replicative intermediate, the antigenome, encoding the small (s) and the large (L) isoforms of the delta antigen (s-HDAg and L-HDAg). The latter is produced following an editing process, changing the amber/stop codon on the s-HDAg-ORF into a tryptophan codon, allowing L-HDAg synthesis by the addition of 19 (or 20) C-terminal amino acids. The two delta proteins play different roles in the viral cell cycle: s-HDAg activates genome replication, while L-HDAg blocks replication and favors virion morphogenesis and propagation. L-HDAg has also been involved in HDV pathogenicity. Understanding the kinetics of viral editing rates in vivo is key to unravel the biology of the virus and understand its spread and natural history. We developed and validated a new assay based on next-generation sequencing and aimed at quantifying HDV RNA editing in plasma. We analyzed plasma samples from 219 patients infected with different HDV genotypes and showed that HDV editing capacity strongly depends on the genotype of the strain.

Complete mitochondrial genome of the versicoloured emerald hummingbird Amazilia versicolor, a polymorphic species

The genome of the versicoloured emerald hummingbird (Amazilia versicolor) was partially sequenced in one-sixth of an Illumina HiSeq lane. The mitochondrial genome was assembled using MIRA and MITObim software, yielding a circular molecule of 16,861 bp in length and deposited in GenBank under the accession number KF624601. The mitogenome contained 13 protein-coding genes, 22 transfer tRNAs, 2 ribosomal RNAs and 1 non-coding control region. The molecule was assembled using 21,927 sequencing reads of 100 bp each, resulting in ∼130 × coverage of uniformly distributed reads along the genome. This is the forth mitochondrial genome described for this highly diverse family of birds and may benefit further phylogenetic, phylogeographic, population genetic and species delimitation studies of hummingbirds.

Association between Dysbiosis in the Gut Microbiota of Primary Osteoporosis Patients and Bone Loss

In recent decades, gut microbiome research has experienced significant growth, driven by technological advances that enable quantifying bacterial taxa with greater precision. Age, diet, and living environment have emerged as three key factors influencing gut microbes. Dysbiosis, resulting from alterations in these factors, may lead to changes in bacterial metabolites that regulate pro- and anti-inflammatory processes and consequently impact bone health. Restoration of a healthy microbiome signature could mitigate inflammation and potentially reduce bone loss associated with osteoporosis or experienced by astronauts during spaceflight. However, current research is hindered by contradictory findings, insufficient sample sizes, and inconsistency in experimental conditions and controls. Despite progress in sequencing technology, defining a healthy gut microbiome across global populations remains elusive. Challenges persist in identifying accurate gut bacterial metabolics, specific taxa, and their effects on host physiology. We suggest greater attention be directed towards this issue in Western countries as the cost of treating osteoporosis in the United States reaches billions of dollars annually, with expenses projected to continue rising.

Comprehensive Analysis of Hepatitis Delta Virus Assembly Determinants According to Genotypes: Lessons From a Study of 526 Hepatitis Delta Virus Clinical Strains

Human hepatitis Delta virus (HDV) infection is associated to the most severe viral hepatic disease, including severe acute liver decompensation and progression to cirrhosis, and hepatocellular carcinoma. HDV is a satellite of hepatitis B virus (HBV) that requires the HBV envelope proteins for assembly of HDV virions. HDV and HBV exhibit a large genetic diversity that extends, respectively to eight (HDV-1 to -8) and to ten (HBV/A to/J) genotypes. Molecular determinants of HDV virion assembly consist of a C-terminal Proline-rich domain in the large Hepatitis Delta Antigen (HDAg) protein, also known as the Delta packaging domain (DPD) and of a Tryptophan-rich domain, the HDV matrix domain (HMD) in the C-terminal region of the HBV envelope proteins. In this study, we performed a systematic genotyping of HBV and HDV in a cohort 1,590 HDV-RNA-positive serum samples collected between 2001 to 2014, from patients originated from diverse parts of the world, thus reflecting a large genetic diversity. Among these samples, 526 HBV (HBV/A, B, C, D, E, and G) and HDV (HDV-1, 2, 3, and 5 to -8) genotype couples could be obtained. We provide results of a comprehensive analysis of the amino-acid sequence conservation within the HMD and structural and functional features of the DPD that may account for the yet optimal interactions between HDV and its helper HBV.

Diagnostic Implications of NGS-Based Molecular Profiling in Mature B-Cell Lymphomas with Potential Bone Marrow Involvement

Background: Methods such as cytogenetics and immunocytology/immunohistology provide essential diagnostic insights but may be limited in ambiguous cases of mature B-cell lymphoma. Next-generation sequencing (NGS) has emerged as a potential tool to improve diagnostics. Methods: We validated the analytical performance of a lymphoid customized NGS panel. Clinical validation was conducted in 226 patients with suspected mature B-cell lymphoma with potential bone marrow involvement across multiple clinically relevant scenarios. Results: NGS (1) achieved 100% sensitivity and specificity with high reproducibility (r = 0.995), confirming its analytical performance. (2) It reliably detected WHO-classified markers, including BRAF mutations in all hairy cell leukemia cases, MYD88/CXCR4 mutations in lymphoplasmacytic lymphoma, and absence of BRAF mutations in splenic B-cell lymphoma with prominent nucleoli. (3) In lymphoma exclusion diagnostics, NGS identified mutations in previously undiagnosed cases, including a BCORL1 mutation leading to reclassification as marginal zone lymphoma. (4) Among 105 confirmed lymphomas, 65% harbored mutations, with detection rates highest in HCL and LPL (100%) and CLL (62%), while follicular lymphoma showed no detectable mutations. (5) In cases with non-interpretable cytogenetics, NGS detected pathogenic variants in 61% of patients, compensating for inconclusive findings. (6) In cases with limited morphological assessment, NGS identified relevant mutations in 70%, outperforming cytogenetics (30%; p = 0.0256, OR = 5.44). Conclusions: NGS enhances the diagnostic accuracy of mature B-cell lymphomas by complementing traditional methods, refining WHO-classified subtypes, and improving detection in cases with inconclusive cytogenetics or morphology. NGS may reduce the need for unnecessary bone marrow re-punctures by providing additional information in ambiguous cases.

Concurrent Detection of a Papillomatous Lesion and Sequence Reads Corresponding to a Member of the Family Adintoviridae in a Bell's Hinge-Back Tortoise (Kinixys belliana)

An adult male Bell's hinge-back tortoise (Kinixys belliana) was admitted to a veterinary clinic due to a swelling in the oral cavity. Physical examination revealed an approximately 2.5 × 1.5 cm sized, irregularly shaped tissue mass with villiform projections extending from its surface located in the oropharyngeal cavity. An initial biopsy was performed, and the lesion was diagnosed as squamous papilloma. Swabs taken for virological examination tested negative with specific PCRs for papillomavirus and herpesvirus. Further analysis of the oropharyngeal mass via metagenomic sequencing revealed sequence reads corresponding to a member of the family Adintoviridae. The tissue mass was removed one week after the initial examination. The oral cavity remained unsuspicious in follow-up examinations performed after one, five and twenty weeks. However, a regrowth of the tissue was determined 23 months after the initial presentation. The resampled biopsy tested negative for sequence reads of Adintoviridae. Conclusively, this report presents the diagnostic testing and therapy of an oral cavity lesion of unknown origin. The significance of concurrent metagenomic determination of adintovirus sequence reads within the tissue lesion is discussed.

An updated database of virus circular RNAs provides new insights into the biogenesis mechanism of the molecule

Virus circular RNAs (circRNA) have been reported to be extensively expressed and play important roles in viral infections. Previously we build the first database of virus circRNAs named VirusCircBase which has been widely used in the field. This study significantly improved the database on both the data quantity and database functionality: the number of virus circRNAs, virus species, host organisms was increased from 46440, 23, 9 to 60859, 43, 22, respectively, and 1902 full-length virus circRNAs were newly added; new functions were added such as visualization of the expression level of virus circRNAs and visualization of virus circRNAs in the Genome Browser. Analysis of the expression of virus circRNAs showed that they had low expression levels in most cells or tissues and showed strong expression heterogeneity. Analysis of the splicing of virus circRNAs showed that they used a much higher proportion of non-canonical back-splicing signals compared to those in animals and plants, and mainly used the A5SS (alternative 5' splice site) in alternative-splicing. Most virus circRNAs have no more than two isoforms. Finally, human genes associated with the virus circRNA production were investigated and more than 1000 human genes exhibited moderate correlations with the expression of virus circRNAs. Most of them showed negative correlations including 42 genes encoding RNA-binding proteins. They were significantly enriched in biological processes related to cell cycle and RNA processing. Overall, the study provides a valuable resource for further studies of virus circRNAs and also provides new insights into the biogenesis mechanisms of virus circRNAs.

Proteomic Characterization of Undifferentiated Small Round Cell Sarcomas With EWSR1 and CIC::DUX4 Translocations Reveals Diverging Tumor Biology and Distinct Diagnostic Markers

Undifferentiated small round cell sarcomas (USRS) of bone and soft tissue are a group of tumors with heterogenic genomic alterations sharing similar morphology. In the present study, we performed a comparative large-scale proteomic analysis of USRS (n = 42) with diverse genomic translocations including classic Ewing sarcomas with EWSR1::FLI1 fusions (n = 24) or EWSR1::ERG fusions (n = 4), sarcomas with an EWSR1 rearrangement (n = 2), CIC::DUX4 fusion (n = 8), as well as tumors classified as USRS with no genetic data available (n = 4). Proteins extracted from formalin-fixed, paraffin-embedded pretherapeutic biopsies were analyzed qualitatively and quantitatively using shotgun mass spectrometry (MS). More than 8000 protein groups could be quantified using data-independent acquisition. Unsupervised hierarchical cluster analysis based on proteomic data allowed stratification of the 42 cases into distinct groups reflecting the different molecular genotypes. Protein signatures that significantly correlated with the respective genomic translocations were identified and used to generate a heatmap of all 42 sarcomas with assignment of cases with unknown molecular genetic data to either the EWSR1- or CIC-rearranged groups. MS-based prediction of sarcoma subtypes was molecularly confirmed in 2 cases where next-generation sequencing was technically feasible. MS also detected proteins routinely used in the immunohistochemical approach for the differential diagnosis of USRS. BCL11B highly expressed in Ewing sarcomas, and BACH2 as well as ETS-1 highly expressed in CIC::DUX4-associated sarcomas, were among proteins identified by the present proteomic study, and were chosen for immunohistochemical confirmation of MS data in our study cohort. Differential expressions of these 3 markers in the 2 genetic groups were further validated in an independent cohort of n = 34 USRS. Finally, our proteomic results point toward diverging signaling pathways in the different USRS subgroups.

Case report: Complete paternal isodisomy on chromosome 18 induces methylation changes in PARD6G-AS1 promotor in a case with arthrogryposis

Uniparental disomy (UPD) is the inheritance of both alleles of a chromosome from only one parent. So far, the detection of UPDs in sequencing data is not well established and a known gap in next-generation sequencing (NGS) diagnostics. By developing a new tool for UPD detection, we re-evaluated an eight-year-old individual presenting with scoliosis, muscle weakness and global developmental delay. Previous panel analysis identified a homozygous likely pathogenic loss-of-function variant in the PIEZO2-gene associated with arthrogryposis (OMIM # 617146). Interestingly, during a re-evaluation process, we identified a region of homozygosity (ROH) covering over 95% of chromosome 18. Segregation and microsatellite analysis within the family revealed that only the father is a heterozygous carrier of the variant in PIEZO2 and confirmed paternal uniparental isodisomy (iUPD) on chromosome 18 in the individual. Further methylation analysis indicated demethylation of the promotor region of PARD6G-AS1, which is described to be maternally imprinted and could possibly influence the individuals' phenotype. Our report describes the first complete iUPD on chromosome 18 and highlights that UPDs can be a cause for homozygous pathogenic variants, which reduces the risk of reoccurrence in case of a new pregnancy in comparison to an autosomal recessive inheritance trait significantly.