Comparison of qPCR and metagenomic sequencing methods for quantifying antibiotic resistance genes in wastewater

Surveillance methods of circulating antibiotic resistance genes (ARGs) are of utmost importance in order to tackle what has been described as one of the greatest threats to humanity in the 21st century. In order to be effective, these methods have to be accurate, quickly deployable, and scalable. In this study, we compare metagenomic shotgun sequencing (TruSeq DNA sequencing) of wastewater samples with a state-of-the-art PCR-based method (Resistomap HT-qPCR) on four wastewater samples that were taken from hospital, industrial, urban and rural areas. ARGs that confer resistance to 11 antibiotic classes have been identified in these wastewater samples using both methods, with the most abundant observed classes of ARGs conferring resistance to aminoglycoside, multidrug-resistance (MDR), macrolide-lincosamide-streptogramin B (MLSB), tetracycline and beta-lactams. In comparing the methods, we observed a strong correlation of relative abundance of ARGs obtained by the two tested methods for the majority of antibiotic classes. Finally, we investigated the source of discrepancies in the results obtained by the two methods. This analysis revealed that false negatives were more likely to occur in qPCR due to mutated primer target sites, whereas ARGs with incomplete or low coverage were not detected by the sequencing method due to the parameters set in the bioinformatics pipeline. Indeed, despite the good correlation between the methods, each has its advantages and disadvantages which are also discussed here. By using both methods together, a more robust ARG surveillance program can be established. Overall, the work described here can aid wastewater treatment plants that plan on implementing an ARG surveillance program.

Roflumilast inhibits tumor growth and migration in STK11/LKB1 deficient pancreatic cancer

Pancreatic cancer is a malignant tumor of the digestive system. It is highly aggressive, easily metastasizes, and extremely difficult to treat. This study aimed to analyze the genes that might regulate pancreatic cancer migration to provide an essential basis for the prognostic assessment of pancreatic cancer and individualized treatment. A CRISPR knockout library directed against 915 murine genes was transfected into TB 32047 cell line to screen which gene loss promoted cell migration. Next-generation sequencing and PinAPL.py- analysis was performed to identify candidate genes. We then assessed the effect of serine/threonine kinase 11 (STK11) knockout on pancreatic cancer by wound-healing assay, chick agnosia (CAM) assay, and orthotopic mouse pancreatic cancer model. We performed RNA sequence and Western blotting for mechanistic studies to identify and verify the pathways. After accelerated Transwell migration screening, STK11 was identified as one of the top candidate genes. Further experiments showed that targeted knockout of STK11 promoted the cell migration and increased liver metastasis in mice. Mechanistic analyses revealed that STK11 knockout influences blood vessel morphogenesis and is closely associated with the enhanced expression of phosphodiesterases (PDEs), especially PDE4D, PDE4B, and PDE10A. PDE4 inhibitor Roflumilast inhibited STK11-KO cell migration and tumor size, further demonstrating that PDEs are essential for STK11-deficient cell migration. Our findings support the adoption of therapeutic strategies, including Roflumilast, for patients with STK11-mutated pancreatic cancer in order to improve treatment efficacy and ultimately prolong survival.

Critical assessment of on-premise approaches to scalable genome analysis

BACKGROUND

Plummeting DNA sequencing cost in recent years has enabled genome sequencing projects to scale up by several orders of magnitude, which is transforming genomics into a highly data-intensive field of research. This development provides the much needed statistical power required for genotype-phenotype predictions in complex diseases.

METHODS

In order to efficiently leverage the wealth of information, we here assessed several genomic data science tools. The rationale to focus on on-premise installations is to cope with situations where data confidentiality and compliance regulations etc. rule out cloud based solutions. We established a comprehensive qualitative and quantitative comparison between BCFtools, SnpSift, Hail, GEMINI, and OpenCGA. The tools were compared in terms of data storage technology, query speed, scalability, annotation, data manipulation, visualization, data output representation, and availability.

RESULTS

Tools that leverage sophisticated data structures are noted as the most suitable for large-scale projects in varying degrees of scalability in comparison to flat-file manipulation (e.g., BCFtools, and SnpSift). Remarkably, for small to mid-size projects, even lightweight relational database.

CONCLUSION

The assessment criteria provide insights into the typical questions posed in scalable genomics and serve as guidance for the development of scalable computational infrastructure in genomics.

Self-oxygenation of engineered living tissues orchestrates osteogenic commitment of mesenchymal stem cells

Oxygenating biomaterials can alleviate anoxic stress, stimulate vascularization, and improve engraftment of cellularized implants. However, the effects of oxygen-generating materials on tissue formation have remained largely unknown. Here, we investigate the impact of calcium peroxide (CPO)-based oxygen-generating microparticles (OMPs) on the osteogenic fate of human mesenchymal stem cells (hMSCs) under a severely oxygen deficient microenvironment. To this end, CPO is microencapsulated in polycaprolactone to generate OMPs with prolonged oxygen release. Gelatin methacryloyl (GelMA) hydrogels containing osteogenesis-inducing silicate nanoparticles (SNP hydrogels), OMPs (OMP hydrogels), or both SNP and OMP (SNP/OMP hydrogels) are engineered to comparatively study their effect on the osteogenic fate of hMSCs. OMP hydrogels associate with improved osteogenic differentiation under both normoxic and anoxic conditions. Bulk mRNAseq analyses suggest that OMP hydrogels under anoxia regulate osteogenic differentiation pathways more strongly than SNP/OMP or SNP hydrogels under either anoxia or normoxia. Subcutaneous implantations reveal a stronger host cell invasion in SNP hydrogels, resulting in increased vasculogenesis. Furthermore, time-dependent expression of different osteogenic factors reveals progressive differentiation of hMSCs in OMP, SNP, and SNP/OMP hydrogels. Our work demonstrates that endowing hydrogels with OMPs can induce, improve, and steer the formation of functional engineered living tissues, which holds potential for numerous biomedical applications, including tissue regeneration and organ replacement therapy.

Allelic Variants Within the ABO Blood Group Phenotype Confer Protection Against Critical COVID-19 Hospital Presentation

Introduction: Coronavirus disease 2019 (COVID-19) disease severity differs widely due to numerous factors including ABO gene-derived susceptibility or resistance. The objective of this study was to investigate the association of the ABO blood group and genetic variations of the ABO gene with COVID-19 severity in a heterogeneous hospital population sample from the United Arab Emirates, with the use of an epidemiological and candidate gene approach from a genome-wide association study (GWAS). Methods: In this cross-sectional study, a total of 646 participants who tested positive for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) were recruited from multiple hospitals and population-based (quarantine camps) recruitment sites from March 2020 to February 2021. The participants were divided into two groups based on the severity of COVID-19: noncritical (n = 453) and critical [intensive care unit (ICU) patients] (n = 193), as per the COVID-19 Reporting and Data System (CO-RADS) classification. The multivariate logistic regression analysis demonstrated the association of ABO blood type as well as circulating anti-A antibodies and anti-B antibodies as well as A and B antigens, in association with critical COVID-19 hospital presentation. A candidate gene analysis approach was conducted from a GWAS where we examined 240 single nucleotide polymorphisms (SNPs) (position in chr9: 136125788-136150617) in the ABO gene, in association with critical COVID-19 hospital presentation. Results: Patients with blood group O [odds ratio (OR): 0.51 (0.33, 0.79); p = 0.003] were less likely to develop critical COVID-19 symptoms. Eight alleles have been identified to be associated with a protective effect of blood group O in ABO 3'untranslated region (UTR): rs199969472 (p = 0.0052), rs34266669 (p = 0.0052), rs76700116 (p = 0.0052), rs7849280 (p = 0.0052), rs34039247 (p = 0.0104), rs10901251 (p = 0.0165), rs9411475 (p = 0.0377), and rs13291798 (p = 0.0377). Conclusion: Our findings suggest that there are novel allelic variants that link genetic variants of the ABO gene and ABO blood groups contributing to the reduced risk of critical COVID-19 disease. This study is the first study to combine genetic and serological evidence of the involvement of the ABO blood groups and the ABO gene allelic associations with COVID-19 severity within the Middle Eastern population.

A Population-Specific Major Allele Reference Genome From The United Arab Emirates Population

The ethnic composition of the population of a country contributes to the uniqueness of each national DNA sequencing project and, ideally, individual reference genomes are required to reduce the confounding nature of ethnic bias. This work represents a representative Whole Genome Sequencing effort of an understudied population. Specifically, high coverage consensus sequences from 120 whole genomes and 33 whole exomes were used to construct the first ever population specific major allele reference genome for the United Arab Emirates (UAE). When this was applied and compared to the archetype hg19 reference, assembly of local Emirati genomes was reduced by ∼19% (i.e., some 1 million fewer calls). In compiling the United Arab Emirates Reference Genome (UAERG), sets of annotated 23,038,090 short (novel: 1,790,171) and 137,713 structural (novel: 8,462) variants; their allele frequencies (AFs) and distribution across the genome were identified. Population-specific genetic characteristics including loss-of-function variants, admixture, and ancestral haplogroup distribution were identified and reported here. We also detect a strong correlation between F ST and admixture components in the UAE. This baseline study was conceived to establish a high-quality reference genome and a genetic variations resource to enable the development of regional population specific initiatives and thus inform the application of population studies and precision medicine in the UAE.

Whole Genome Sequencing of Four Representatives From the Admixed Population of the United Arab Emirates

Whole genome sequences (WGS) of four nationals of the United Arab Emirates (UAE) at an average coverage of 33X have been completed and described. The selection of suitable subpopulation representatives was informed by a preceding comprehensive population structure analysis. Representatives were chosen based on their central location within the subpopulation on a principal component analysis (PCA) and the degree to which they were admixed. Novel genomic variations among the different subgroups of the UAE population are reported here. Specifically, the WGS analysis identified 4,161,067-4,798,806 variants in the four individual samples, where approximately 80% were single nucleotide polymorphisms (SNPs) and 20% were insertions or deletions (indels). An average of 2.75% was found to be novel variants according to dbSNP (build 151). This is the first report of structural variants (SV) from WGS data from UAE nationals. There were 15,677-20,339 called SVs, of which around 13.5% were novel. The four samples shared 1,399,178 variants, each with distinct variants as follows: 1,085,524 (for the individual denoted as UAE S011), 1,228,559 (UAE S012), 791,072 (UAE S013), and 906,818 (UAE S014). These results show a previously unappreciated population diversity in the region. The synergy of WGS and genotype array data was demonstrated through variant annotation of the former using 2.3 million allele frequencies for the local population derived from the latter technology platform. This novel approach of combining breadth and depth of array and WGS technologies has guided the choice of population genetic representatives and provides complementary, regionalized allele frequency annotation to new genomes comprising millions of loci.

Genetic Diversity and Low Stratification of the Population of the United Arab Emirates

With high consanguinity rates on the Arabian Peninsula, it would not have been unexpected if the population of the United Arab Emirates (UAE) was shown to be relatively homogenous. However, this study of 1000 UAE nationals provided a contrasting perspective, one of a relatively heterogeneous population. Located at the apex of Europe, Asia, and Africa, the observed diversity could be explained by a plethora of migration patterns since the first Out-of-Africa movement. A strategy to explore the extent of genetic variation of the population of the UAE is presented. The first step involved a comprehensive population stratification study that was instructive for subsequent whole genome sequencing (WGS) of suitable representatives (which is described elsewhere). When these UAE data were compared to previous smaller studies from the region, the findings were consistent with a population that is a diverse and admixed group of people. However, rather than sharp and distinctive clusters, cluster analysis reveals low levels of stratification throughout the population. UAE emirates exhibit high within-Emirate-distance/among-Emirate distance ratios. Supervised admixture analysis showed a continuous gradient of ancestral populations, suggesting that admixture on the south eastern tip of the Arabian Peninsula occurred gradually. When visualized using a unique technique that combined admixture ratios and principal component analysis (PCA), unappreciated diversity was revealed while mitigating projection bias of conventional PCA. We observe low population stratification in the UAE in terms of homozygosity versus separation cluster coefficients. This holds for the UAE in a global context as well as for isolated cluster analysis of the Emirati birthplaces. However, the subtle clustering observed in the Emirates reflects geographic proximity and historic migration events. The analytical strategy used here highlights the complementary nature of data from genotype array and WGS for anthropological studies. Specifically, genotype array data were instructive to select representative subjects for WGS. Furthermore, from the 2.3 million allele frequencies obtained from genotype arrays, we identified 46,481 loci with allele frequencies that were significantly different with respect to other world populations. This comparison of allele frequencies facilitates variant prioritization in common diseases. In addition, these loci bear great potential as biomarkers in anthropological and forensic studies.

Introducing the first whole genomes of nationals from the United Arab Emirates

Whole Genome Sequencing (WGS) provides an in depth description of genome variation. In the era of large-scale population genome projects, the assembly of ethnic-specific genomes combined with mapping human reference genomes of underrepresented populations has improved the understanding of human diversity and disease associations. In this study, for the first time, whole genome sequences of two nationals of the United Arab Emirates (UAE) at >27X coverage are reported. The two Emirati individuals were predominantly of Central/South Asian ancestry. An in-house customized pipeline using BWA, Picard followed by the GATK tools to map the raw data from whole genome sequences of both individuals was used. A total of 3,994,521 variants (3,350,574 Single Nucleotide Polymorphisms (SNPs) and 643,947 indels) were identified for the first individual, the UAE S001 sample. A similar number of variants, 4,031,580 (3,373,501 SNPs and 658,079 indels), were identified for UAE S002. Variants that are associated with diabetes, hypertension, increased cholesterol levels, and obesity were also identified in these individuals. These Whole Genome Sequences has provided a starting point for constructing a UAE reference panel which will lead to improvements in the delivery of precision medicine, quality of life for affected individuals and a reduction in healthcare costs. The information compiled will likely lead to the identification of target genes that could potentially lead to the development of novel therapeutic modalities.