Whole genome sequencing

RanBALL: An Ensemble Random Projection Model for Identifying Subtypes of B-cell Acute Lymphoblastic Leukemia

As the most common pediatric malignancy, B-cell acute lymphoblastic leukemia (B-ALL) has multiple distinct subtypes characterized by recurrent and sporadic somatic and germline genetic alterations. Identification of B-ALL subtypes can facilitate risk stratification and enable tailored therapeutic approaches. Existing methods for B-ALL subtyping primarily depend on immunophenotypic, cytogenetic and genomic analyses, which would be costly, complicated, and laborious in clinical practice applications. To overcome these challenges, we present RanBALL (an Ensemble Random Projection-Based Model for Identifying B-Cell Acute Lymphoblastic Leukemia Subtypes), an accurate and cost-effective model for B-ALL subtype identification based on transcriptomic profiling only. RanBALL leverages random projection (RP) to construct an ensemble of dimension-reduced multi-class support vector machine (SVM) classifiers for B-ALL subtyping. Results based on 100 times 5-fold cross validation tests for >1700 B-ALL patients demonstrated that the proposed model achieved an accuracy of 93.35%, indicating promising prediction capabilities of RanBALL for B-ALL subtyping. The high accuracies of RanBALL suggested that our model could effectively capture underlying patterns of transcriptomic profiling for accurate B-ALL subtype identification. We believe RanBALL will facilitate the discovery of B-ALL subtype-specific marker genes and therapeutic targets, and eventually have consequential positive impacts on downstream risk stratification and tailored treatment design.

Molecular typing methods to characterize Brucella spp. from animals: A review

Brucellosis is an infectious disease of animals that can infect humans. The disease causes significant economic losses and threatens human health. A timely and accurate disease diagnosis plays a vital role in the identification of brucellosis. In addition to traditional diagnostic methods, molecular methods allow diagnosis and typing of the causative agent of brucellosis. This review will discuss various methods, such as Bruce-ladder, Suiladder, high-resolution melt analysis, restriction fragment length polymorphism, multilocus sequence typing, multilocus variable-number tandem repeat analysis, and whole-genome sequencing single-nucleotide polymorphism, for the molecular typing of Brucella and discuss their advantages and disadvantages.

A haplotype-resolved genome assembly of Malus domestica 'Red Fuji'

The 'Red Fuji' apple (Malus domestica), is one of the most important and popular economic crops worldwide in the fruit industry. Using PacBio HiFi long reads and Hi-C reads, we assembled a high-quality haplotype-resolved genome of 'Red Fuji', with sizes of 668.7 and 668.8 Mb, and N50 sizes of 34.1 and 31.4 Mb. About 97.2% of sequences were anchored in 34 chromosomes. We annotated both haploid genomes, identifying a total of 95,439 protein-coding genes in the two haplotype genomes, with 98% functional annotation. The haplotype-resolved genome of 'Red Fuji' apple stands as a precise benchmark for an array of analyses, such as comparative genomics, transcriptomics, and allelic expression studies. This comprehensive resource is paramount in unraveling variations in allelic expression, advancing quality improvements, and refining breeding efforts.

Resolving tumor evolution: a phylogenetic approach

The evolutionary dynamics of cancer, characterized by its profound heterogeneity, demand sophisticated tools for a holistic understanding. This review delves into tumor phylogenetics, an essential approach bridging evolutionary biology with oncology, offering unparalleled insights into cancer's evolutionary trajectory. We provide an overview of the workflow, encompassing study design, data acquisition, and phylogeny reconstruction. Notably, the integration of diverse data sets emerges as a transformative step, enhancing the depth and breadth of evolutionary insights. With this integrated perspective, tumor phylogenetics stands poised to redefine our understanding of cancer evolution and influence therapeutic strategies.

Single Cell Atlas: a single-cell multi-omics human cell encyclopedia

Single-cell sequencing datasets are key in biology and medicine for unraveling insights into heterogeneous cell populations with unprecedented resolution. Here, we construct a single-cell multi-omics map of human tissues through in-depth characterizations of datasets from five single-cell omics, spatial transcriptomics, and two bulk omics across 125 healthy adult and fetal tissues. We construct its complement web-based platform, the Single Cell Atlas (SCA, www.singlecellatlas.org ), to enable vast interactive data exploration of deep multi-omics signatures across human fetal and adult tissues. The atlas resources and database queries aspire to serve as a one-stop, comprehensive, and time-effective resource for various omics studies.

Multimodal Omics Approaches to Aging and Age-Related Diseases

Aging is associated with a progressive decline in physiological capacities and an increased risk of aging-associated disorders. An increasing body of experimental evidence shows that aging is a complex biological process coordinately regulated by multiple factors at different molecular layers. Thus, it is difficult to delineate the overall systematic aging changes based on single-layer data. Instead, multimodal omics approaches, in which data are acquired and analyzed using complementary omics technologies, such as genomics, transcriptomics, and epigenomics, are needed for gaining insights into the precise molecular regulatory mechanisms that trigger aging. In recent years, multimodal omics sequencing technologies that can reveal complex regulatory networks and specific phenotypic changes have been developed and widely applied to decode aging and age-related diseases. This review summarizes the classification and progress of multimodal omics approaches, as well as the rapidly growing number of articles reporting on their application in the field of aging research, and outlines new developments in the clinical treatment of age-related diseases based on omics technologies.

Bilateral inflammatory recurrence of HER-2 positive breast cancer: a unique case report and literature review

Inflammatory breast cancer (IBC) is an aggressive and rare form of breast cancer with a poor prognosis. The occurrence of bilateral IBC in a short period of time is extremely rare. In this case report, a 54-year-old woman diagnosed with invasive ductal carcinoma of the left breast underwent lumpectomy, lymph node dissection, chemotherapy, and radiotherapy but opted against trastuzumab treatment. Four years later, she experienced bilateral breast inflammation, skin changes, edema, and heat (calor). Biopsies confirmed breast cancer metastasis to both breasts. Whole-Exome Sequencing revealed genetic mutations, including PIK3CA and C4orf54, in both primary and recurrent tumors, with significant downregulation in the recurrent tumors. KEGG analysis suggested potential enrichment of axon guidance signal pathways in both tumors. The patient showed a partial response after treatment with liposome paclitaxel, along with targeted therapy using trastuzumab and pertuzumab. This case report sheds light on the rare occurrence of bilateral inflammatory breast cancer post-HER-2 treatment and highlights the importance of genetic profiling in understanding the disease. Further research on clinical targets for breast cancer management is warranted.

Seroprevalence and Molecular Characterization of Brucella abortus from the Himalayan Marmot in Qinghai, China

Objective

Brucellosis is a serious public health issue in Qinghai (QH), China. Surveying the seroprevalence and isolation of B. abortus strains from marmots is key to understanding the role of wildlife in the maintenance and spread of brucellosis.

Methods

In this study, a set of methods, including a serology survey, bacteriology, antibiotic susceptibility, molecular genotyping (MLST and MLVA), and genome sequencing, were employed to characterize the two B. abortus strains.

Results

The seroprevalence of brucellosis in marmots was 7.0% (80/1146) by serum tube agglutination test (SAT); one Brucella strain was recovered from these positive samples, and another Brucella strain from a human. Two strains were identified as B. abortus bv. 1 and were susceptible to all eight drugs examined. The distribution patterns of the accessory genes, virulence associated genes, and resistance genes of the two strains were consistent, and there was excellent collinearity between the two strains on chromosome I, but they had significant SVs in chromosome II, including inversions and translocations. MLST genotyping identified two B. abortus strains as ST2, and MLVA-16 analysis showed that the two strains clustered with strains from northern China. WGS-SNP phylogenetic analysis showed that the strains were genetically homogeneous with strains from the northern region, implying that strains from a common lineage were spread continuously in different regions and hosts.

Conclusion

Seroprevalence and molecular clues demonstrated frequent direct or indirect contact between sheep/goats, cattle, and marmots, implying that wildlife plays a vital role in the maintenance and spread of B. abortus in the Qinghai-Tibet Plateau.

Proteome-wide structural analysis quantifies structural conservation across distant species

Traditional evolutionary biology research mainly relies on sequence information to infer evolutionary relationships between genes or proteins. In contrast, protein structural information has long been overlooked, although structures are more conserved and closely linked to the functions than the sequences. To address this gap, we conducted a proteome-wide structural analysis using experimental and computed protein structures for organisms from the three distinct domains, including Homo sapiens (eukarya), Escherichia coli (bacteria), and Methanocaldococcus jannaschii (archaea). We reveal the distribution of structural similarity and sequence identity at the genomic level and characterize the twilight zone, where signals obtained from sequence alignment are blurred and evolutionary relationships cannot be inferred unambiguously. We find that structurally similar homologous protein pairs in the twilight zone account for ∼0.004%-0.021% of all possible protein pair combinations, which translates to ∼8%-32% of the protein-coding genes, depending on the species under comparison. In addition, by comparing the structural homologs, we show that human proteins involved in the energy supply are more similar to their E. coli homologs, whereas proteins relating to the central dogma are more similar to their M. jannaschii homologs. We also identify a bacterial GPCR homolog in the E. coli proteome that displays distinctive domain architecture. Our results shed light on the characteristics of the twilight zone and the origin of different pathways from a protein structure perspective, highlighting an exciting new frontier in evolutionary biology.

The Development of Plant Genome Sequencing Technology and Its Conservation and Application in Endangered Gymnosperms

Genome sequencing is widely recognized as a fundamental pillar in genetic research and legal studies of biological phenomena, providing essential insights for genetic investigations and legal analyses of biological events. The field of genome sequencing has experienced significant progress due to rapid improvements in scientific and technological developments. These advancements encompass not only significant improvements in the speed and quality of sequencing but also provide an unparalleled opportunity to explore the subtle complexities of genomes, particularly in the context of rare species. Such a wide range of possibilities has successfully supported the validation of plant gene functions and the refinement of precision breeding methodologies. This expanded scope now includes a comprehensive exploration of the current state and conservation efforts of gymnosperm gene sequencing, offering invaluable insights into their genomic landscapes. This comprehensive review elucidates the trajectory of development and the diverse applications of genome sequencing. It encompasses various domains, including crop breeding, responses to abiotic stress, species evolutionary dynamics, biodiversity, and the unique challenges faced in the conservation and utilization of gymnosperms. It highlights both ongoing challenges and the unveiling of forthcoming developmental trajectories.

Comparative transcriptome analysis reveals deep molecular landscapes in stony coral Montipora clade

Introduction: Coral reefs, among the most invaluable ecosystems in the world, face escalating threats from climate change and anthropogenic activities. To decipher the genetic underpinnings of coral adaptation and resilience, we undertook comprehensive transcriptome profiling of two emblematic coral species, Montipora foliosa and Montipora capricornis, leveraging PacBio Iso-Seq technology. These species were strategically selected for their ecological significance and their taxonomic proximity within the Anthozoa class. Methods: Our study encompassed the generation of pristine transcriptomes, followed by thorough functional annotation via diverse databases. Subsequently, we quantified transcript abundance and scrutinized gene expression patterns, revealing notable distinctions between the two species. Results: Intriguingly, shared orthologous genes were identified across a spectrum of coral species, highlighting a substantial genetic conservation within scleractinian corals. Importantly, a subset of genes, integral to biomineralization processes, emerged as exclusive to scleractinian corals, shedding light on their intricate evolutionary history. Furthermore, we discerned pronounced upregulation of genes linked to immunity, stress response, and oxidative-reduction processes in M. foliosa relative to M. capricornis. These findings hint at the presence of more robust mechanisms in M. foliosa for maintaining internal equilibrium and effectively navigating external challenges, underpinning its potential ecological advantage. Beyond elucidating genetic adaptation in corals, our research underscores the urgency of preserving genetic diversity within coral populations. Discussion: These insights hold promise for informed conservation strategies aimed at safeguarding these imperiled ecosystems, bearing ecological and economic significance. In synthesis, our study seamlessly integrates genomic inquiry with ecological relevance, bridging the gap between molecular insights and the imperative to conserve coral reefs in the face of mounting threats.

Structure-based pathogenicity relationship identifier for predicting effects of single missense variants and discovery of higher-order cancer susceptibility clusters of mutations

We report the structure-based pathogenicity relationship identifier (SPRI), a novel computational tool for accurate evaluation of pathological effects of missense single mutations and prediction of higher-order spatially organized units of mutational clusters. SPRI can effectively extract properties determining pathogenicity encoded in protein structures, and can identify deleterious missense mutations of germ line origin associated with Mendelian diseases, as well as mutations of somatic origin associated with cancer drivers. It compares favorably to other methods in predicting deleterious mutations. Furthermore, SPRI can discover spatially organized pathogenic higher-order spatial clusters (patHOS) of deleterious mutations, including those of low recurrence, and can be used for discovery of candidate cancer driver genes and driver mutations. We further demonstrate that SPRI can take advantage of AlphaFold2 predicted structures and can be deployed for saturation mutation analysis of the whole human proteome.

Predicting Pathology of Missense Mutations through Protein-Specific Evolutionary Pattern

Missense mutations, which are single base pair genetic alternation resulting in a different amino acid, are among the most common occurring variants in exon regions of the human genome and may lead to diseases. Thus to assess the effects of missense mutations, it is essential to investigate the evolutionary history of the protein under selection pressures. In this study, we employ a continuous-time Markov model to investigate the evolutionary patterns in protein sequences and a Bayesian Markov chain Monte Carlo method to estimate the substitution rates for protein of interest, from which we obtain scoring matrices. Specifically, we examined the evolutionary patterns of protein sequences containing missense mutations using a species tree to define the phylogeny of the protein of interest. We thoroughly studied the evolutionary pattern of human muscle glycogen phosphorylase containing 127 known missense mutations, and identified characteristic evolutionary patterns in 63 proteins with 2,238 missense mutations, including both deleterious and neutral effects. Our results show that the estimated protein-specific evolutionary pattern-based scoring matrices (PSM) lead to higher sensitivity in detecting the pathological effects of missense mutations, compared to the general evolutionary pattern-based scoring matrix of Blosum62 (BL62) matrix. By incorporating PSM, the performance of a recently released structure-based model SPRI for evaluating missense mutations is further improved.

An Optimal Stacked ResNet-BiLSTM-Based Accurate Detection and Classification of Genetic Disorders

Gene is located inside the nuclease and the genetic data is contained in deoxyribonucleic acid (DNA). A person's gene count ranges from 20,000 to 30,000. Even a minor alteration to the DNA sequence can be harmful if it affects the cell's fundamental functions. As a result, the gene begins to act abnormally. The sorts of genetic abnormalities brought on by mutation include chromosomal disorders, complex disorders, and single-gene disorders. Therefore, a detailed diagnosis method is required. Thus, we proposed an Elephant Herd Optimization-Whale Optimization Algorithm (EHO-WOA) optimized Stacked ResNet-Bidirectional Long Term Short Memory (ResNet-BiLSTM) model for detecting genetic disorders. Here, a hybrid EHO-WOA algorithm is presented to assess the Stacked ResNet-BiLSTM architecture's fitness. The ResNet-BiLSTM design uses the genotype and gene expression phenotype as input data. Furthermore, the proposed method identifies rare genetic disorders such as Angelman Syndrome, Rett Syndrome, and Prader-Willi Syndrome. It demonstrates the effectiveness of the developed model with greater accuracy, recall, specificity, precision, and f1-score. Thus, a wide range of DNA deficiencies including Prader-Willi syndrome, Marfan syndrome, Early Onset Morbid Obesity, Rett syndrome, and Angelman syndrome are predicted accurately.

The Genome of a Pigeonpea Compatible Rhizobial Strain '10ap3' Appears to Lack Common Nodulation Genes

The symbiotic fixation of atmospheric nitrogen (N) in root nodules of tropical legumes such as pigeonpea (Cajanus cajan) is a complex process, which is regulated by multiple genetic factors at the host plant genotype microsymbiont interface. The process involves multiple genes with various modes of action and is accomplished only when both organisms are compatible. Therefore, it is necessary to develop tools for the genetic manipulation of the host or bacterium towards improving N fixation. In this study, we sequenced the genome of a robust rhizobial strain, Rhizobium tropici '10ap3' that was compatible with pigeonpea, and we determined its genome size. The genome consisted of a large circular chromosome (6,297,373 bp) and contained 6013 genes of which 99.13% were coding sequences. However only 5833 of the genes were associated with proteins that could be assigned to specific functions. The genes for nitrogen, phosphorus and iron metabolism, stress response and the adenosine monophosphate nucleoside for purine conversion were present in the genome. However, the genome contained no common nod genes, suggesting that an alternative pathway involving a purine derivative was involved in the symbiotic association with pigeonpea.

Retinitis Pigmentosa: Current Clinical Management and Emerging Therapies

Retinitis pigmentosa (RP) comprises a group of inherited retinal dystrophies characterized by the degeneration of rod photoreceptors, followed by the degeneration of cone photoreceptors. As a result of photoreceptor degeneration, affected individuals experience gradual loss of visual function, with primary symptoms of progressive nyctalopia, constricted visual fields and, ultimately, central vision loss. The onset, severity and clinical course of RP shows great variability and unpredictability, with most patients already experiencing some degree of visual disability in childhood. While RP is currently untreatable for the majority of patients, significant efforts have been made in the development of genetic therapies, which offer new hope for treatment for patients affected by inherited retinal dystrophies. In this exciting era of emerging gene therapies, it remains imperative to continue supporting patients with RP using all available options to manage their condition. Patients with RP experience a wide variety of physical, mental and social-emotional difficulties during their lifetime, of which some require timely intervention. This review aims to familiarize readers with clinical management options that are currently available for patients with RP.

The study of selection signature and its applications on identification of candidate genes using whole genome sequencing data in chicken - a review

Chicken is a major source of protein for the increasing human population and is useful for research purposes. There are almost 1,600 distinct regional breeds of chicken across the globe, among which a large body of genetic and phenotypic variations has been accumulated due to extensive natural and artificial selection. Moreover, natural selection is a crucial force for animal domestication. Several approaches have been adopted to detect selection signatures in different breeds of chicken using whole genome sequencing (WGS) data including integrated haplotype score (iHS), cross-populated extend haplotype homozygosity test (XP-EHH), fixation index (F_ST), cross-population composite likelihood ratio (XP-CLR), nucleotide diversity (Pi), and others. In addition, gene enrichment analyses are utilized to determine KEGG pathways and gene ontology (GO) terms related to traits of interest in chicken. Herein, we review different studies that have adopted diverse approaches to detect selection signatures in different breeds of chicken. This review systematically summarizes different findings on selection signatures and related candidate genes in chickens. Future studies could combine different selection signatures approaches to strengthen the quality of the results thereby providing more affirmative inference. This would further aid in deciphering the importance of selection in chicken conservation for the increasing human population.

De novo assembly of a chromosome-level reference genome of the ornamental butterfly Sericinus montelus based on nanopore sequencing and Hi-C analysis

Sericinus montelus (Lepidoptera, Papilionidae, Parnassiinae) is a high-value ornamental swallowtail butterfly species widely distributed in Northern and Central China, Japan, Korea, and Russia. The larval stage of this species feeds exclusively on Aristolochia plants. The Aristolochia species is well known for its high levels of aristolochic acids (AAs), which have been found to be carcinogenic for numerous animals. The swallowtail butterfly is among the few that can feed on these toxic host plants. However, the genetic adaptation of S. montelus to confer new abilities for AA tolerance has not yet been well explored, largely due to the limited genomic resources of this species. This study aimed to present a chromosome-level reference genome for S. montelus using the Oxford Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C technology. The final assembly was composed of 581.44 Mb with an expected genome size of 619.27 Mb. Further, 99.98% of the bases could be anchored onto 30 chromosomes. The N50 of contigs and scaffolds was 5.74 and 19.12 Mb, respectively. Approximately 48.86% of the assembled genome was suggested to be repeat elements, and 13,720 protein-coding genes were predicted in the current assembly. The phylogenetic analysis indicated that S. montelus diverged from the common ancestor of swallowtails about 58.57-80.46 million years ago. Compared with related species, S. montelus showed a significant expansion of P450 gene family members, and positive selections on eloa, heatr1, and aph1a resulted in the AA tolerance for S. montelus larva. The de novo assembly of a high-quality reference genome for S. montelus provided a fundamental genomic tool for future research on evolution, genome genetics, and toxicology of the swallowtail butterflies.

The consequences of recurrent genetic and epigenetic variants in human pluripotent stem cells

It is well established that human pluripotent stem cells (hPSCs) can acquire genetic and epigenetic changes during culture in vitro. Given the increasing use of hPSCs in research and therapy and the vast expansion in the number of hPSC lines available for researchers, the International Society for Stem Cell Research has recognized the need to reassess quality control standards for ensuring the genetic integrity of hPSCs. Here, we summarize current knowledge of the nature of recurrent genetic and epigenetic variants in hPSC culture, the methods for their detection, and what is known concerning their effects on cell behavior in vitro or in vivo. We argue that the potential consequences of low-level contamination of cell therapy products with cells bearing oncogenic variants are essentially unknown at present. We highlight the key challenges facing the field with particular reference to safety assessment of hPSC-derived cellular therapeutics.

Three highly variable genome regions of the four dengue virus serotypes can accurately recapitulate the CDS phylogeny

The circulation of the four-dengue virus (DENV) serotypes has significantly increased in recent years, accompanied by an increase in viral genetic diversity. In order to conduct disease surveillance and understand DENV evolution and its effects on virus transmission and disease, efficient and accurate methods for phylogenetic classification are required. Phylogenetic analysis of different viral genes sequences is the most used method, the envelope gene (E) being the most frequently selected target. We explored the genetic variability of the four DENV serotypes throughout their complete coding sequence (CDS) of sequences available in GenBank and used genomic regions of different variability rate to recapitulate the phylogeny obtained with the DENV CDS. Our results indicate that the use of high or low variable regions accurately recapitulate the phylogeny obtained with CDS of sequences from different DENV genotypes. However, when analyzing the phylogeny of a single genotype, highly variable regions performed better in recapitulating the distance branch length, topology, and support of the CDS phylogeny. The use of three concatenated highly variable regions was not statistically different in distance branch length and support to that obtained in CDS phylogeny.•This study demonstrated the ability of highly variable regions of the DENV genome to recapitulate the phylogeny obtained with the full coding sequence (CDS).•The use of genomic regions of high or low variability did not affect the performance in recapitulating the phylogeny obtained with CDS from different genotypes. However, when phylogeny was analyzed for sequences from a single genotype, highly variable regions performed better in recapitulating the distance branch length, topology, and support of the CDS phylogeny.•The use of concatenated highly variable genome regions represent a useful option for recapitulating genome-wide phylogenies in analyses of sequences belonging to the same DENV genotype.

Automated filtering of genome-wide large deletions through an ensemble deep learning framework

Computational methods based on whole genome linked-reads and short-reads have been successful in genome assembly and detection of structural variants (SVs). Numerous variant callers that rely on linked-reads and short reads can detect genetic variations, including SVs. A shortcoming of existing tools is a propensity for overestimating SVs, especially for deletions. Optimizing the advantages of linked-read and short-read sequencing technologies would thus benefit from an additional step to effectively identify and eliminate false positive large deletions. Here, we introduce a novel tool, AquilaDeepFilter, aiming to automatically filter genome-wide false positive large deletions. Our approach relies on transforming sequencing data into an image and then relying on convolutional neural networks to improve classification of candidate deletions as such. Input data take into account multiple alignment signals including read depth, split reads and discordant read pairs. We tested the performance of AquilaDeepFilter on five linked-reads and short-read libraries sequenced from the well-studied NA24385 sample, validated against the Genome in a Bottle benchmark. To demonstrate the filtering ability of AquilaDeepFilter, we utilized the SV calls from three upstream SV detection tools including Aquila, Aquila_stLFR and Delly as the baseline. We showed that AquilaDeepFilter increased precision while preserving the recall rate of all three tools. The overall F1-score improved by an average 20% on linked-reads and by an average of 15% on short-read data. AquilaDeepFilter also compared favorably to existing deep learning based methods for SV filtering, such as DeepSVFilter. AquilaDeepFilter is thus an effective SV refinement framework that can improve SV calling for both linked-reads and short-read data.

Sample-multiplexing approaches for single-cell sequencing

Single-cell sequencing is widely used in biological and medical studies. However, its application with multiple samples is hindered by inefficient sample processing, high experimental costs, ambiguous identification of true single cells, and technical batch effects. Here, we introduce sample-multiplexing approaches for single-cell sequencing in transcriptomics, epigenomics, genomics, and multiomics. In single-cell transcriptomics, sample multiplexing uses variants of native or artificial features as sample markers, enabling sample pooling and decoding. Such features include: (1) natural genetic variation, (2) nucleotide-barcode anchoring on cellular or nuclear membranes, (3) nucleotide-barcode internalization to the cytoplasm or nucleus, (4) vector-based barcode expression in cells, and (5) nucleotide-barcode incorporation during library construction. Other single-cell omics methods are based on similar concepts, particularly single-cell combinatorial indexing. These methods overcome current challenges, while enabling super-loading of single cells. Finally, selection guidelines are presented that can accelerate technological application.

The genetic basis for adult onset glaucoma: Recent advances and future directions

Glaucoma, a diverse group of eye disorders that results in the degeneration of retinal ganglion cells, is the world's leading cause of irreversible blindness. Apart from age and ancestry, the major risk factor for glaucoma is increased intraocular pressure (IOP). In primary open-angle glaucoma (POAG), the anterior chamber angle is open but there is resistance to aqueous outflow. In primary angle-closure glaucoma (PACG), crowding of the anterior chamber angle due to anatomical alterations impede aqueous drainage through the angle. In exfoliation syndrome and exfoliation glaucoma, deposition of white flaky material throughout the anterior chamber directly interfere with aqueous outflow. Observational studies have established that there is a strong hereditable component for glaucoma onset and progression. Indeed, a succession of genome wide association studies (GWAS) that were centered upon single nucleotide polymorphisms (SNP) have yielded more than a hundred genetic markers associated with glaucoma risk. However, a shortcoming of GWAS studies is the difficulty in identifying the actual effector genes responsible for disease pathogenesis. Building on the foundation laid by GWAS studies, research groups have recently begun to perform whole exome-sequencing to evaluate the contribution of protein-changing, coding sequence genetic variants to glaucoma risk. The adoption of this technology in both large population-based studies as well as family studies are revealing the presence of novel, protein-changing genetic variants that could enrich our understanding of the pathogenesis of glaucoma. This review will cover recent advances in the genetics of primary open-angle glaucoma, primary angle-closure glaucoma and exfoliation glaucoma, which collectively make up the vast majority of all glaucoma cases in the world today. We will discuss how recent advances in research methodology have uncovered new risk genes, and how follow up biological investigations could be undertaken in order to define how the risk encoded by a genetic sequence variant comes into play in patients. We will also hypothesise how data arising from characterising these genetic variants could be utilized to predict glaucoma risk and the manner in which new therapeutic strategies might be informed.

Breeding of African sheep reared under low-input/output smallholder production systems for trypanotolerance

Trypanosomiasis is a disease caused by unicellular protozoan parasites. Small ruminants succumb to trypanosomiasis in areas of high tsetse fly challenge, resulting in serious economic loss often to farmers in low-input smallholder systems. At present, trypanosomiasis is treated with trypanocidal drugs, but access to these can be limited, and increasing parasite resistance raises questions about their efficacy. The development of trypanotolerance in small ruminant flocks through targeted breeding strategies is considered a sustainable and economical option for controlling African trypanosomiasis. Recently, quantitative trait loci (QTLs) associated with trypanotolerance traits in sheep have been reported. The results of these studies form the basis for more studies to identify QTLs associated with trypanosomiasis resistance, particularly in African livestock species. For example, signatures of positive selection for trypanotolerance have been identified using genome-wide single-nucleotide polymorphism data. However, there are several challenges in performing genetic analyses using data from low-input smallholder systems, including a lack of recorded pedigree and production records and the need for large sample sizes when flock sizes are often fewer than 50 animals. Breeding strategies to improve trypanotolerance should also preserve existing genetic diversity as well as minimize excessive genetic introgression by trypanosusceptible breeds. This review discusses the possibilities of breeding for trypanosome tolerance/resistance in low-input/low-output small ruminant production systems. Potential challenges are outlined, and potential available genetic resources are described as a foundation for future work.

Genetics and sports performance: the present and future in the identification of talent for sports based on DNA testing

The impact of genetics on physiology and sports performance is one of the most debated research aspects in sports sciences. Nearly 200 genetic polymorphisms have been found to influence sports performance traits, and over 20 polymorphisms may condition the status of the elite athlete. However, with the current evidence, it is certainly too early a stage to determine how to use genotyping as a tool for predicting exercise/sports performance or improving current methods of training. Research on this topic presents methodological limitations such as the lack of measurement of valid exercise performance phenotypes that make the study results difficult to interpret. Additionally, many studies present an insufficient cohort of athletes, or their classification as elite is dubious, which may introduce expectancy effects. Finally, the assessment of a progressively higher number of polymorphisms in the studies and the introduction of new analysis tools, such as the total genotype score (TGS) and genome-wide association studies (GWAS), have produced a considerable advance in the power of the analyses and a change from the study of single variants to determine pathways and systems associated with performance. The purpose of the present study was to comprehensively review evidence on the impact of genetics on endurance- and power-based exercise performance to clearly determine the potential utility of genotyping for detecting sports talent, enhancing training, or preventing exercise-related injuries, and to present an overview of recent research that has attempted to correct the methodological issues found in previous investigations.

A comprehensive WGS-based pipeline for the identification of new candidate genes in inherited retinal dystrophies

To enhance the use of Whole Genome Sequencing (WGS) in clinical practice, it is still necessary to standardize data analysis pipelines. Herein, we aimed to define a WGS-based algorithm for the accurate interpretation of variants in inherited retinal dystrophies (IRD). This study comprised 429 phenotyped individuals divided into three cohorts. A comparison of 14 pathogenicity predictors, and the re-definition of its cutoffs, were performed using panel-sequencing curated data from 209 genetically diagnosed individuals with IRD (training cohort). The optimal tool combinations, previously validated in 50 additional IRD individuals, were also tested in patients with hereditary cancer (n = 109), and with neurological diseases (n = 47) to evaluate the translational value of this approach (validation cohort). Then, our workflow was applied for the WGS-data analysis of 14 individuals from genetically undiagnosed IRD families (discovery cohort). The statistical analysis showed that the optimal filtering combination included CADDv1.6, MAPP, Grantham, and SIFT tools. Our pipeline allowed the identification of one homozygous variant in the candidate gene CFAP20 (c.337 C > T; p.Arg113Trp), a conserved ciliary gene, which was abundantly expressed in human retina and was located in the photoreceptors layer. Although further studies are needed, we propose CFAP20 as a candidate gene for autosomal recessive retinitis pigmentosa. Moreover, we offer a translational strategy for accurate WGS-data prioritization, which is essential for the advancement of personalized medicine.

5WBF: a low-cost and straightforward whole blood filtration method suitable for whole-genome sequencing of Plasmodium falciparum clinical isolates

Background

Whole-genome sequencing (WGS) is becoming increasingly helpful to assist malaria control programmes. A major drawback of this approach is the large amount of human DNA compared to parasite DNA extracted from unprocessed whole blood. As red blood cells (RBCs) have a diameter of about 7–8 µm and exhibit some deformability, it was hypothesized that cheap and commercially available 5 µm filters might retain leukocytes but much less of Plasmodium falciparum-infected RBCs. This study aimed to test the hypothesis that such a filtration method, named 5WBF (for 5 µm Whole Blood Filtration), may provide highly enriched parasite material suitable for P. falciparum WGS.

Methods

Whole blood was collected from five patients experiencing a P. falciparum malaria episode (ring-stage parasitaemia range: 0.04–5.5%) and from mock samples obtained by mixing synchronized, ring-stage cultured P. falciparum 3D7 parasites with uninfected human whole blood (final parasitaemia range: 0.02–1.1%). These whole blood samples (50 to 400 µL) were diluted in RPMI 1640 medium or PBS 1× buffer and filtered with a syringe connected to a 5 µm commercial filter. DNA was extracted from 5WBF-treated and unfiltered counterpart blood samples using a commercial kit. The 5WBF method was evaluated on the ratios of parasite:human DNA assessed by qPCR and by sequencing depth and percentages of coverage from WGS data (Illumina NextSeq 500). As a comparison, the popular selective whole-genome amplification (sWGA) method, which does not rely on blood filtration, was applied to the unfiltered counterpart blood samples.

Results

After applying 5WBF, qPCR indicated an average of twofold loss in the amount of parasite template DNA (Pf ARN18S gene) and from 4096- to 65,536-fold loss of human template DNA (human β actin gene). WGS analyses revealed that > 95% of the  parasite nuclear and organellar genomes were all covered at ≥ 10× depth for all samples tested. In sWGA counterparts, the organellar genomes were poorly covered and from 47.7 to 82.1% of the nuclear genome was covered at ≥ 10× depth depending on parasitaemia. Sequence reads were homogeneously distributed across gene sequences for 5WBF-treated samples (n = 5460 genes; mean coverage: 91×; median coverage: 93×; 5th percentile: 70×; 95th percentile: 103×), allowing the identification of gene copy number variations such as for gch1. This later analysis was not possible for sWGA-treated samples, as a much more heterogeneous distribution of reads across gene sequences was observed (mean coverage: 80×; median coverage: 51×; 5th percentile: 7×; 95th percentile: 245×).

Conclusions

The novel 5WBF leucodepletion method is simple to implement and based on commercially available, standardized 5 µm filters which cost from 1.0 to 1.7€ per unit depending on suppliers. 5WBF permits extensive genome-wide analysis of P. falciparum ring-stage isolates from minute amounts of whole blood even with parasitaemias as low as 0.02%.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04073-1.

Genomic and neuroimaging approaches to bipolar disorder

BACKGROUND

To date, besides genome-wide association studies, a variety of other genetic analyses (e.g. polygenic risk scores, whole-exome sequencing and whole-genome sequencing) have been conducted, and a large amount of data has been gathered for investigating the involvement of common, rare and very rare types of DNA sequence variants in bipolar disorder. Also, non-invasive neuroimaging methods can be used to quantify changes in brain structure and function in patients with bipolar disorder.

AIMS

To provide a comprehensive assessment of genetic findings associated with bipolar disorder, based on the evaluation of different genomic approaches and neuroimaging studies.

METHOD

We conducted a PubMed search of all relevant literatures from the beginning to the present, by querying related search strings.

RESULTS

ANK3, CACNA1C, SYNE1, ODZ4 and TRANK1 are five genes that have been replicated as key gene candidates in bipolar disorder pathophysiology, through the investigated studies. The percentage of phenotypic variance explained by the identified variants is small (approximately 4.7%). Bipolar disorder polygenic risk scores are associated with other psychiatric phenotypes. The ENIGMA-BD studies show a replicable pattern of lower cortical thickness, altered white matter integrity and smaller subcortical volumes in bipolar disorder.

CONCLUSIONS

The low amount of explained phenotypic variance highlights the need for further large-scale investigations, especially among non-European populations, to achieve a more complete understanding of the genetic architecture of bipolar disorder and the missing heritability. Combining neuroimaging data with genetic data in large-scale studies might help researchers acquire a better knowledge of the engaged brain regions in bipolar disorder.

FgSfl1 and Its Conserved PKA Phosphorylation Sites Are Important for Conidiation, Sexual Reproduction, and Pathogenesis in Fusarium graminearum

The fungal plant pathogen, Fusarium graminearum, contains two genes, FgCPK1 and FgCPK2, encoding the catalytic subunits of cAMP-dependent protein kinase A. FgCPK1 and FgCPK2 are responsible for most of the PKA activities and have overlapping functions in various cellular processes in F. graminearum. The cpk1 cpk2 double mutant was significantly reduced in growth, rarely produced conidia, and was non-pathogenic. In this study, we found that the cpk1 cpk2 double mutant was unstable and produced fast-growing spontaneous sectors that were defective in plant infection. All spontaneous suppressor strains had mutations in FgSFL1, a transcription factor gene orthologous to SFL1 in yeast. Thirteen suppressor strains had non-sense mutations at Q501, three suppressor strains had frameshift mutations at W198, and five suppressor strains had mutations in the HSF binding domain of FgSfl1. Only one suppressor strain had both a non-synonymous mutation at H225 and a non-sense mutation at R490. We generated the SFL1 deletion mutant and found that it produced less than 2% of conidia than that of the wild-type strain PH-1. The sfl1 mutant was significantly reduced in the number of perithecia on carrot agar plates at 7 days post-fertilization (dpf). When incubated for more than 12 days, ascospore cirrhi were observed on the sfl1 mutant perithecia. The infection ability of the sfl1 deletion mutant was also obviously defective. Furthermore, we found that in addition to the S223 and S559 phosphorylation sites, FgSFL1 had another predicted phosphorylation site: T452. Interestingly, the S223 phosphorylation site was responsible for sexual reproduction, and the T452 phosphorylation site was responsible for growth and sexual reproduction. Only the S559 phosphorylation site was found to play an important role in conidiation, sexual reproduction, and infection. Overall, our results indicate that FgSFL1 and its conserved PKA phosphorylation sites are important for vegetative growth, conidiation, sexual reproduction, and pathogenesis in F. graminearum.

JWES: a new pipeline for whole genome/exome sequence data processing, management, and gene-variant discovery, annotation, prediction, and genotyping

Whole genome and exome sequencing (WGS/WES) are the most popular next‐generation sequencing (NGS) methodologies and are at present often used to detect rare and common genetic variants of clinical significance. We emphasize that automated sequence data processing, management, and visualization should be an indispensable component of modern WGS and WES data analysis for sequence assembly, variant detection (SNPs, SVs), imputation, and resolution of haplotypes. In this manuscript, we present a newly developed findable, accessible, interoperable, and reusable (FAIR) bioinformatics‐genomics pipeline Java based Whole Genome/Exome Sequence Data Processing Pipeline (JWES) for efficient variant discovery and interpretation, and big data modeling and visualization. JWES is a cross‐platform, user‐friendly, product line application, that entails three modules: (a) data processing, (b) storage, and (c) visualization. The data processing module performs a series of different tasks for variant calling, the data storage module efficiently manages high‐volume gene‐variant data, and the data visualization module supports variant data interpretation with Circos graphs. The performance of JWES was tested and validated in‐house with different experiments, using Microsoft Windows, macOS Big Sur, and UNIX operating systems. JWES is an open‐source and freely available pipeline, allowing scientists to take full advantage of all the computing resources available, without requiring much computer science knowledge. We have successfully applied JWES for processing, management, and gene‐variant discovery, annotation, prediction, and genotyping of WGS and WES data to analyze variable complex disorders. In summary, we report the performance of JWES with some reproducible case studies, using open access and in‐house generated, high‐quality datasets.

Next generation sequencing approaches to evaluate water and wastewater quality

The emergence of next generation sequencing (NGS) is revolutionizing the potential to address complex microbiological challenges in the water industry. NGS technologies can provide holistic insight into microbial communities and their functional capacities in water and wastewater systems, thus eliminating the need to develop a new assay for each target organism or gene. However, several barriers have hampered wide-scale adoption of NGS by the water industry, including cost, need for specialized expertise and equipment, challenges with data analysis and interpretation, lack of standardized methods, and the rapid pace of development of new technologies. In this critical review, we provide an overview of the current state of the science of NGS technologies as they apply to water, wastewater, and recycled water. In addition, a systematic literature review was conducted in which we identified over 600 peer-reviewed journal articles on this topic and summarized their contributions to six key areas relevant to the water and wastewater fields: taxonomic classification and pathogen detection, functional and catabolic gene characterization, antimicrobial resistance (AMR) profiling, bacterial toxicity characterization, Cyanobacteria and harmful algal bloom identification, and virus characterization. For each application, we have presented key trends, noteworthy advancements, and proposed future directions. Finally, key needs to advance NGS technologies for broader application in water and wastewater fields are assessed.

Short and long-read genome sequencing methodologies for somatic variant detection; genomic analysis of a patient with diffuse large B-cell lymphoma

Recent advances in throughput and accuracy mean that the Oxford Nanopore Technologies PromethION platform is a now a viable solution for genome sequencing. Much of the validation of bioinformatic tools for this long-read data has focussed on calling germline variants (including structural variants). Somatic variants are outnumbered many-fold by germline variants and their detection is further complicated by the effects of tumour purity/subclonality. Here, we evaluate the extent to which Nanopore sequencing enables detection and analysis of somatic variation. We do this through sequencing tumour and germline genomes for a patient with diffuse B-cell lymphoma and comparing results with 150 bp short-read sequencing of the same samples. Calling germline single nucleotide variants (SNVs) from specific chromosomes of the long-read data achieved good specificity and sensitivity. However, results of somatic SNV calling highlight the need for the development of specialised joint calling algorithms. We find the comparative genome-wide performance of different tools varies significantly between structural variant types, and suggest long reads are especially advantageous for calling large somatic deletions and duplications. Finally, we highlight the utility of long reads for phasing clinically relevant variants, confirming that a somatic 1.6 Mb deletion and a p.(Arg249Met) mutation involving TP53 are oriented in trans.

Responsible, practical genomic data sharing that accelerates research

Data sharing anchors reproducible science, but expectations and best practices are often nebulous. Communities of funders, researchers and publishers continue to grapple with what should be required or encouraged. To illuminate the rationales for sharing data, the technical challenges and the social and cultural challenges, we consider the stakeholders in the scientific enterprise. In biomedical research, participants are key among those stakeholders. Ethical sharing requires considering both the value of research efforts and the privacy costs for participants. We discuss current best practices for various types of genomic data, as well as opportunities to promote ethical data sharing that accelerates science by aligning incentives.

Whole-Genome Sequencing of Yeast Cells

The budding yeast, Saccharomyces cerevisiae, has been widely used for genetic studies of fundamental cellular functions. The isolation and analysis of yeast mutants is a commonly used and powerful technique to identify the genes that are involved in a process of interest. Furthermore, natural genetic variation among wild yeast strains has been studied for analysis of polygenic traits by quantitative trait loci mapping. Whole-genome sequencing, often combined with bulk segregant analysis, is a powerful technique that helps determine the identity of mutations causing a phenotype. Here, we describe protocols for the construction of libraries for S. cerevisiae whole-genome sequencing. We also present a bioinformatic pipeline to determine the genetic variants in a yeast strain using whole-genome sequencing data. This pipeline can also be used for analyzing Schizosaccharomyces pombe mutants. © 2019 by John Wiley & Sons, Inc. Basic Protocol 1: Generation of haploid spores for bulk segregant analysis Basic Protocol 2: Extraction of genomic DNA from yeast cells Basic Protocol 3: Shearing of genomic DNA for library preparation Basic Protocol 4: Construction and amplification of DNA libraries Support Protocol 1: Annealing oligonucleotides for forming Y-adapters Support Protocol 2: Size selection and cleanup using SPRI beads Basic Protocol 5: Identification of genomic variants from sequencing data.

Genetics and Extracellular Vesicles of Pediatrics Sleep Disordered Breathing and Epilepsy

Sleep remains one of the least understood phenomena in biology, and sleep disturbances are one of the most common behavioral problems in childhood. The etiology of sleep disorders is complex and involves both genetic and environmental factors. Epilepsy is the most popular childhood neurological condition and is characterized by an enduring predisposition to generate epileptic seizures, and the neurobiological, cognitive, psychological, and social consequences of this condition. Sleep and epilepsy are interrelated, and the importance of sleep in epilepsy is less known. The state of sleep also influences whether a seizure will occur at a given time, and this differs considerably for various epilepsy syndromes. The development of epilepsy has been associated with single or multiple gene variants. The genetics of epilepsy is complex and disorders exhibit significant genetic heterogeneity and variability in the expressivity of seizures. Phenobarbital (PhB) is the most widely used antiepileptic drug. With its principal mechanism of action to prolong the opening time of the γ-aminobutyric acid (GABA)-A receptor-associated chloride channel, it enhances chloride anion influx into neurons, with subsequent hyperpolarization, thereby reducing excitability. Enzymes that metabolize pharmaceuticals including PhB are well known for having genetic polymorphisms that contribute to adverse drug–drug interactions. PhB metabolism is highly dependent upon the cytochrome P450 (CYP450) and genetic polymorphisms can lead to variability in active drug levels. The highly polymorphic CYP2C19 isozymes are responsible for metabolizing a large portion of routinely prescribed drugs and variants contribute significantly to adverse drug reactions and therapeutic failures. A limited number of CYP2C19 single nucleotide polymorphisms (SNPs) are involved in drug metabolism. Extracellular vesicles (EVs) are circular membrane fragments released from the endosomal compartment as exosomes are shed from the surfaces of the membranes of most cell types. Increasing evidence indicated that EVs play a pivotal role in cell-to-cell communication. Theses EVs may play an important role between sleep, epilepsy, and treatments. The discovery of exosomes provides potential strategies for the diagnosis and treatment of many diseases including neurocognitive deficit. The aim of this study is to better understand and provide further knowledge about the metabolism and interactions between phenobarbital and CYP2C19 polymorphisms in children with epilepsy, interplay between sleep, and EVs. Understanding this interplay between epilepsy and sleep is helpful in the optimal treatment of all patients with epileptic seizures. The use of genetics and extracellular vesicles as precision medicine for the diagnosis and treatment of children with sleep disorder will improve the prognosis and the quality of life in patients with epilepsy.

Comparative 16S rDNA metagenomics study of two samples of cassava peel heap from Nigeria and India

The microbiology of many cassava products and the wastes generated during the processing have been reported; however, majority of these reports used culture-dependent methods. This has resulted in a dearth of information on the bacterial diversity of cassava peels and peel heaps. Large amounts of cassava peels generated during the processing of cassava root are usually discharged on land or water as wastes and are allowed to rot in the open, especially in some developing countries. Culture-independent methods such as PCR-based amplification and sequencing of 16S rRNA genes, among others have been used in recent times to study the diversity of microbes in different environmental samples. In this study, bacterial isolates were screened for cellulase and xylanase enzyme activities on minimal agar and genomic DNA was isolated from cassava peel samples; metagenomics was carried out using MiSeq 2 × 300 with primers specific for V3-V4 bacterial region. Samples collected from Nigeria (AAG) had more species compared with samples from India (JHA) with 793 and 525 observed OTUs (operational taxonomic units), respectively. Five bacterial isolates from cassava peel-heap samples obtained from Ogbomoso, Nigeria showed no ability to produce cellulase enzyme, seven isolates from the Nigeria samples and three from Jorhat samples were positive for xylanase production; the highest amylase activity was shown by isolate AG18 (10,055 U/mL), while the lowest was recorded for isolate JA2 (2333 U/mL) with a significant difference observed in the amylase activities of isolates (p ≤ 0.05). Comparing the most abundant taxonomy for each of the samples at different taxonomic levels, the most abundant for sample AAG were phylum Firmicutes (42.11%), class Bacilli (41.27%), order Lactobacillales (33.11%), family Acetobacteraceae (31.30%), genus Acetobacter (30.02%) and unclassified species of Acetobacter (29.88%), while sample JHA had Actinobacteria (47.47%) as the highest phylum and class, order Actinomycetales (47.47%), family Brevibacteriaceae (46.97%), genus Brevibacterium (46.97%) and unclassified species of Brevibacterium (46.89%). This study provides an insight into the vast diversity of the bacteria associated with cassava peel heaps and the ability of some of the bacteria to produce selected extracellular enzymes.

Effective identification of soybean candidate genes involved in resistance to soybean cyst nematode via direct whole genome re-sequencing of two segregating mutants

KEY MESSAGE

Three soybean candidate genes involved in resistance to soybean cyst nematode race 4 were identified via direct whole genome re-sequencing of two segregating mutants. The genes conferring resistance to soybean cyst nematode (SCN) race 4 (Hg type 1.2.3.5.7) in soybean (Glycine max L. Merr.) remains unknown. Next generation sequencing-based methods identify a wide range of targets, it is difficult to identify genes underlying traits. Use of the MutMap and QTL-seq methods to identify trait candidate genes needs backcrossing and is very time-consuming. Here we report a simple method to effectively identify candidate genes involved in resistance to SCN race 4. Two ethane methylsulfonate mutagenized mutants of soybean 'PI 437654', whose SCN race 4-infection phenotype altered, were selected. Six relevant whole genomes were re-sequenced, and then calling of genomic variants (SNPs and InDels) was conducted and compared to 'Williams 82'. The comparison eliminated many genomic variants from the mutant lines that overlapped two non-phenotypic but mutant progeny plants, wild-type PI 437654 and 'Zhonghuang 13'. Finally, only 27 mutations were found among 10 genes. Of these 10 genes, 3 genes, Glyma.09g054000, Glyma.16g065700 and Glyma.18g192200 were overlapped between two phenotypic mutant progeny plants. Therefore, the three genes may be the candidate genes involved in resistance of PI 437654 to soybean cyst nematode race 4. This method simplifies the effective identification of candidate genes.

Whole metagenomic shotgun sequencing of the subgingival microbiome of diabetics and non-diabetics with different periodontal conditions

OBJECTIVE

The aim of this study was to use high-resolution whole metagenomic shotgun sequencing to characterize the subgingival microbiome of patients with/without type 2 Diabetes Mellitus and with/without periodontitis.

DESIGN

Twelve subjects, falling into one of the four study groups based on the presence/absence of poorly controlled type 2 Diabetes Mellitus and moderate-severe periodontitis, were selected. For each eligible subject, subgingival plaque samples were collected at 4 sites, all representative of the periodontal condition of the individual (i.e., non-bleeding sulci in subjects without a history of periodontitis, bleeding pockets in patients with moderate-severe periodontitis). The subgingival microbiome was evaluated using high-resolution whole metagenomic shotgun sequencing.

RESULTS

The results showed that: (i) the presence of type 2 Diabetes Mellitus and/or periodontitis were associated with a tendency of the subgingival microbiome to decrease in richness and diversity; (ii) the presence of type 2 Diabetes Mellitus was not associated with significant differences in the relative abundance of one or more species in patients either with or without periodontitis; (iii) the presence of periodontitis was associated with a significantly higher relative abundance of Anaerolineaceae bacterium oral taxon 439 in type 2 Diabetes Mellitus patients.

CONCLUSIONS

Whole metagenomic shotgun sequencing of the subgingival microbiome was extremely effective in the detection of low-abundant taxon. Our results point out a significantly higher relative abundance of Anaerolineaceae bacterium oral taxon 439 in patients with moderate to severe periodontitis vs patients without history of periodontitis, which was maintained when the comparison was restricted to type 2 diabetics.

Development of a Novel Mule Deer Genomic Assembly and Species-Diagnostic SNP Panel for Assessing Introgression in Mule Deer, White-Tailed Deer, and Their Interspecific Hybrids

Mule deer (Odocoileus hemionus) are endemic to a wide variety of habitats in western North America, many of which are shared in sympatry with their closely related sister-species white-tailed deer (Odocoileus virginianus), whom they hybridize with in wild populations. Although mule deer meet many ideal conditions for a molecular ecological research species, such as high abundance, ecological importance, and broad dispersal and gene flow, conservation genetic studies have been limited by a relative lack of existing genomic resources and inherent difficulties caused by introgression with white-tailed deer. Many molecular tools currently available for the study of cervids were designed using reference assemblies of divergent model species, specifically cattle (Bos taurus). Bovidae and Cervidae diverged approximately 28 million years ago, therefore, we sought to ameliorate the available resources by contributing the first mule deer whole genome sequence draft assembly with an average genome-wide read depth of 25X, using the white-tailed genome assembly (Ovir.te_1.0) as a reference. Comparing the two assemblies, we identified ∼33 million single nucleotide polymorphisms (SNPs) and insertion/deletion variants. We then verified fixed SNP differences between the two species and developed a 40-loci SNP assay capable of identifying pure mule deer, white-tailed deer, and interspecific hybrids. Assignment capacity of the panel, which was tested on simulated datasets, is reliable up to and including the third backcross hybrid generation. Identification of post-F1 hybrids will be necessary for hybrid zone population studies going forward, and the new mule deer assembly will be a valuable resource for genetic and comparative genomics studies.

Mass Spectrometry-Based Biomarkers in Drug Development

Advances in mass spectrometry, proteomics, protein bioanalytical approaches, and biochemistry have led to a rapid evolution and expansion in the area of mass spectrometry-based biomarker discovery and development. The last decade has also seen significant progress in establishing accepted definitions, guidelines, and criteria for the analytical validation, acceptance and qualification of biomarkers. These advances have coincided with a decreased return on investment for pharmaceutical research and development and an increasing need for better early decision making tools. Empowering development teams with tools to measure a therapeutic interventions impact on disease state and progression, measure target engagement and to confirm predicted pharmacodynamic effects is critical to efficient data-driven decision making. Appropriate implementation of a biomarker or a combination of biomarkers can enhance understanding of a drugs mechanism, facilitate effective translation from the preclinical to clinical space, enable early proof of concept and dose selection, and increases the efficiency of drug development. Here we will provide descriptions of the different classes of biomarkers that have utility in the drug development process as well as review specific, protein-centric, mass spectrometry-based approaches for the discovery of biomarkers and development of targeted assays to measure these markers in a selective and analytically precise manner.

AFLP-AFLP in silico-NGS approach reveals polymorphisms in repetitive elements in the malignant genome

The increasing interest in exploring the human genome and identifying genetic risk factors contributing to the susceptibility to and outcome of diseases has supported the rapid development of genome-wide techniques. However, the large amount of obtained data requires extensive bioinformatics analysis. In this work, we established an approach combining amplified fragment length polymorphism (AFLP), AFLP in silico and next generation sequencing (NGS) methods to map the malignant genome of patients with chronic myeloid leukemia. We compared the unique DNA fingerprints of patients generated by the AFLP technique approach with those of healthy donors to identify AFLP markers associated with the disease and/or the response to treatment with imatinib, a tyrosine kinase inhibitor. Among the statistically significant AFLP markers selected for NGS analysis and virtual fingerprinting, we identified the sequences of three fragments in the region of DNA repeat element OldhAT1, LINE L1M7, LTR MER90, and satellite ALR/Alpha among repetitive elements, which may indicate a role of these non-coding repetitive sequences in hematological malignancy. SNPs leading to the presence/absence of these fragments were confirmed by Sanger sequencing. When evaluating the results of AFLP analysis for some fragments, we faced the frequently discussed size homoplasy, resulting in co-migration of non-identical AFLP fragments that may originate from an insertion/deletion, SNP, somatic mutation anywhere in the genome, or combination thereof. The AFLP–AFLP in silico–NGS procedure represents a smart alternative to microarrays and relatively expensive and bioinformatically challenging whole-genome sequencing to detect the association of variable regions of the human genome with diseases.

Establishing a stable, repeatable platform for measuring changes in sperm DNA methylation

BACKGROUND

Several independent research groups have shown that alterations in human sperm methylation profiles correlate with decreased fecundity and an increased risk of poor embryo development. Moving these initial findings from the lab into a clinical setting where they can be used to measure male infertility though requires a platform that is stable and robust against batch effects that can occur between sample runs. Operating parameters must be established, performance characteristics determined, and guidelines set to ensure repeatability and accuracy. The standard for technical validation of a lab developed test (LDT) in the USA comes from the Clinical Laboratory Improvement Amendments (CLIA). However, CLIA was introduced in 1988, before the advent of genome-wide profiling and associated computational analysis. This, coupled with its intentionally general nature, makes its interpretation for epigenetic assays non-trivial.

RESULTS

Here, we present an interpretation of the CLIA technical validation requirements for profiling DNA methylation and calling aberrant methylation using the Illumina Infinium platform (e.g., the 450HM and MethylationEPIC). We describe an experimental design to meet these requirements, the experimental results obtained, and the operating parameters established.

CONCLUSIONS

The CLIA guidelines, although not intended for high-throughput assays, can be interpreted in a way that is consistent with modern epigenetic assays. Based on such an interoperation, Illumina's Infinium platform is quite amenable to usage in a clinical setting for diagnostic work.

The advances of genetics research on Hirschsprung's disease

Hirschsprung's disease (HSCR) is a rare and complex congenital disorder characterized by the absence of the enteric neurons in lower digestive tract with an incidence of 1/5 000. Affected infant usually suffer from severe constipation with megacolon and distended abdomen, and face long-term complications even after surgery. In the last 2 decades, great efforts and progresses have been made in understanding the genetics and molecular biological mechanisms that underlie HSCR. However, only a small fraction of the genetic risk can be explained by the identified mutations in the previously established genes. To search novel genetic alterations, new study designs with advanced technologies such as genome/exome-wide association studies (GWASs/EWASs) and next generation sequencing (NGS) on target genes or whole genome/exome, were applied to HSCR. In this review, we summaries the current development of the genetics researches on HSCR based on GWASs/EWASs and NGS, focusing on the newly discovered variants and genes, and their potential roles in HSCR pathogenesis.

Molecular footprints of inshore aquatic adaptation in Indo-Pacific humpback dolphin (Sousa chinensis)

The Indo-Pacific humpback dolphin, Sousa chinensis, being a member of cetaceans, had fully adapted to inshore waters. As a threatened marine mammal, little molecular information available for understanding the genetic basis of ecological adaptation. We firstly sequenced and obtained the draft genome map of S. chinensis. Phylogenetic analysis in this study, based on the single copy orthologous genes of the draft genome, is consistent with traditional phylogenetic classification. The comparative genomic analysis indicated that S. chinensis had 494 species-specific gene families, which involved immune, DNA repair and sensory systems associated with the potential adaption mechanism. We also identified the expansion and positive selection genes in S. chinensis lineage to investigate the potential adaptation mechanism. Our study provided the potential insight into the molecular bases of ecological adaptation in Indo-Pacific humpback dolphin and will be also valuable for future understanding the ecological adaptation and evolution of cetaceans at the genomic level.

Private payer coverage policies for exome sequencing (ES) in pediatric patients: trends over time and analysis of evidence cited

Purpose

Whole Exome Sequencing (WES) is being adopted for neurodevelopmental disorders in pediatric patients. However, little is known about current coverage policies or the evidence cited supporting these policies. Our study is the first in-depth review of private payer WES coverage policies for pediatric patients with neurodevelopmental disorders.

Methods

We reviewed private payer coverage policies and examined evidence cited in the policies of the 15 largest payers in 2017, and trends in coverage policies and evidence cited (2015 – 2017) for the five largest payers.

Results

There were four relevant policies (N=5 payers) in 2015 and 13 policies (N=15 payers) in 2017. In 2015, no payer covered WES, but by 2017, three payers from the original registry payers did. In 2017, eight of the 15 payers covered WES. We found variations in the number and types of evidence cited. Positive coverage policies tended to include a larger number and range of citations.

Conclusion

We conclude that more systematic assessment of evidence cited in coverage policies can provide a greater understanding of coverage policies and how evidence is used. Such assessments could facilitate the ability of researchers to provide the needed evidence, and the ability of clinicians to provide the most appropriate testing for patients.

Sepsis Diagnostics in the Era of "Omics" Technologies

Sepsis is a multifactorial clinical syndrome with an extremely dynamic clinical course and with high diverse clinical phenotype. Early diagnosis is crucial for the final clinical outcome. Previous studies have not identified a biomarker for the diagnosis of sepsis which would have sufficient sensitivity and specificity. Identification of the infectious agents or the use of molecular biology, next gene sequencing, has not brought significant benefit for the patient in terms of early diagnosis. Therefore, we are currently searching for biomarkers, through "omics" technologies with sufficient diagnostic specificity and sensitivity, able to predict the clinical course of the disease and the patient response to therapy. Current progress in the use of systems biology technologies brings us hope that by using big data from clinical trials such biomarkers will be found.

Bioinformatics-based tools in drug discovery: the cartography from single gene to integrative biological networks

Originally developed for the analysis of biological sequences, bioinformatics has advanced into one of the most widely recognized domains in the scientific community. Despite this technological evolution, there is still an urgent need for nontoxic and efficient drugs. The onus now falls on the 'omics domain to meet this need by implementing bioinformatics techniques that will allow for the introduction of pioneering approaches in the rational drug design process. Here, we categorize an updated list of informatics tools and explore the capabilities of integrative bioinformatics in disease control. We believe that our review will serve as a comprehensive guide toward bioinformatics-oriented disease and drug discovery research.