Search

The unexplained association between infection and autoimmune disease is strongest for hepatitis C virus-induced cryoglobulinemic vasculitis (HCV-cryovas). To analyze its origins, we traced the evolution of pathogenic rheumatoid factor (RF) autoantibodies in four HCV-cryovas patients by deep single-cell multi-omic analysis, revealing three sources of B cell somatic mutation converged to drive the accumulation of a large disease-causing clone. A method for quantifying low-affinity binding revealed recurring antibody variable domain combinations created by V(D)J recombination that bound self-immunoglobulin G (IgG) but not viral E2 antigen. Whole-genome sequencing revealed thousands of somatic mutations, numerically comparable to chronic lymphocytic leukemia and normal memory B cells, but with 1-2 corresponding to driver mutations found recurrently in B cell leukemia and lymphoma. V(D)J hypermutation created autoantibodies with compromised solubility in complex with self-IgG. In this virus-induced autoimmune disease, infection promotes a catastrophic confluence of somatic mutagenesis in the descendants of a single B cell., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2025 The Authors. Published by Elsevier Inc. All rights reserved.)

It is estimated that two billion people are latently infected with Mycobacterium tuberculosis ( Mtb ), the causative agent of tuberculosis (TB). Latent Mtb infection (LTBI) can occur in multiple organs, including the lymphatics. The risk of LTBI reactivation increases in immunocompromised conditions, such as coinfection with human immunodeficiency virus (HIV), and during treatment of autoimmune diseases and organ transplantation. The immunological correlates of protection against TB, including against reactivation of LTBI, remain largely elusive. Here, we used a mouse model of latent lymphatic Mtb infection to dissect the immunological mechanisms underlying LTBI containment versus reactivation. We show that immunosuppression-mediated reactivation of lymphatic LTBI and the subsequent spread to non-lymphatic organs can be prevented by vaccination with multiple recombinant BCG (rBCG) strains despite the deficiency of CD4 + T cells. Using spatial transcriptomics, multi-parameter imaging, network analysis and bioinformatic integration of histopathological images, we reveal that immunosuppression is associated with a distinct repositioning of non-CD4 immune cells at the edge of TB lesions within the infection-draining cervical lymph nodes. While B cells increased in numbers, they are dispensable for the containment of LTBI. Lymphatic Mtb infection in different immune cell-deficient mouse strains, antibody-mediated cell depletion and adoptive transfer experiments into highly susceptible mice unequivocally show that vaccination-mediated prevention of LTBI reactivation is critically dependent on CD8 + T cells. These findings have profound implications for our understanding of immunity to TB and the management of LTBI.

Background: The targeted application of cutting-edge high-throughput molecular data technologies provides an enormous opportunity to address key health, economic and environmental issues in the tropics within the One Health framework. The Earth's tropical regions are projected to contain > 50% of the world's population by 2050 coupled with 80% of its biodiversity however these regions are relatively less developed economically, with agricultural productivity substantially lower than temperate zones, a large percentage of its population having limited health care options and much of its biodiversity understudied and undescribed. The generation of high-throughput molecular data and bespoke bioinformatics capability to address these unique challenges offers an enormous opportunity for people living in the tropics. MAIN: In this review we discuss in depth solutions to challenges to populations living in tropical zones across three critical One Health areas: human health, biodiversity and food production. This review will examine how some of the challenges in the tropics can be addressed through the targeted application of advanced omics and bioinformatics and will discuss how local populations can embrace these technologies through strategic outreach and education ensuring the benefits of the One Health approach is fully realised through local engagement., Conclusion: Within the context of the One Health framework, we will demonstrate how genomic technologies can be utilised to improve the overall quality of life for half the world's population., Competing Interests: Declarations Ethics approval and consent to participate Not applicable. Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Whole-exome sequencing of two unrelated kindreds with systemic autoimmune disease featuring antinuclear antibodies with IgG4 elevation uncovered an identical ultrarare heterozygous TNIP1 Q333P variant segregating with disease. Mice with the orthologous Q346P variant developed antinuclear autoantibodies, salivary gland inflammation, elevated IgG2c, spontaneous germinal centers and expansion of age-associated B cells, plasma cells and follicular and extrafollicular helper T cells. B cell phenotypes were cell-autonomous and rescued by ablation of Toll-like receptor 7 (TLR7) or MyD88. The variant increased interferon-β without altering nuclear factor kappa-light-chain-enhancer of activated B cells signaling, and impaired MyD88 and IRAK1 recruitment to autophagosomes. Additionally, the Q333P variant impaired TNIP1 localization to damaged mitochondria and mitophagosome formation. Damaged mitochondria were abundant in the salivary epithelial cells of Tnip1 Q346P mice. These findings suggest that TNIP1-mediated autoimmunity may be a consequence of increased TLR7 signaling due to impaired recruitment of downstream signaling molecules and damaged mitochondria to autophagosomes and may thus respond to TLR7-targeted therapeutics., (© 2024. The Author(s).)

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

Genetic management is a critical component of threatened species conservation. Understanding spatial patterns of genetic diversity is essential for evaluating the resilience of fragmented populations to accelerating anthropogenic threats. Nowhere is this more relevant than on the Australian continent, which is experiencing an ongoing loss of biodiversity that exceeds any other developed nation. Using a proprietary genome complexity reduction-based method (DArTSeq), we generated a data set of 3239 high quality Single Nucleotide Polymorphisms (SNPs) to investigate spatial patterns and indices of genetic diversity in the koala ( Phascolarctos cinereus ), a highly specialised folivorous marsupial that is experiencing rapid and widespread population declines across much of its former range. Our findings demonstrate that current management divisions across the state of New South Wales (NSW) do not fully represent the distribution of genetic diversity among extant koala populations, and that care must be taken to ensure that translocation paradigms based on these frameworks do not inadvertently restrict gene flow between populations and regions that were historically interconnected. We also recommend that koala populations should be prioritised for conservation action based on the scale and severity of the threatening processes that they are currently faced with, rather than placing too much emphasis on their perceived value (e.g., as reservoirs of potentially adaptive alleles), as our data indicate that existing genetic variation in koalas is primarily partitioned among individual animals. As such, the extirpation of koalas from any part of their range represents a potentially critical reduction of genetic diversity for this iconic Australian species., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have influenced the work reported in this manuscript., (© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.)

Background: Past studies suggest that there are changes in peripheral blood cell gene expression in response to ischaemic stroke; however, the specific changes which occur during the acute phase are poorly characterised. The current study aimed to identify peripheral blood cell genes specifically associated with the early response to ischaemic stroke using whole blood samples collected from participants diagnosed with ischaemic stroke (n = 29) or stroke mimics (n = 27) following emergency presentation to hospital. Long non-coding RNA (lncRNA), mRNA and micro-RNA (miRNA) abundance was measured by RNA-seq, and the consensusDE package was used to identify genes which were differentially expressed between groups. A sensitivity analysis excluding two participants with metastatic disease was also conducted., Results: The mean time from symptom onset to blood collection was 2.6 h. Most strokes were mild (median NIH stroke scale score 2.0). Ten mRNAs (all down-regulated in samples provided by patients experiencing ischaemic stroke) and 30 miRNAs (14 over-expressed and 16 under-expressed in participants with ischaemic stroke) were significantly different between groups in the whole cohort and sensitivity analyses. No significant over-representation of gene ontology categories by the differentially expressed genes was observed. Random forest analysis suggested a panel of differentially expressed genes (ADGRG7 and miRNAs 96, 532, 6766, 6798 and 6804) as potential ischaemic stroke biomarkers, although modelling analyses demonstrated that these genes had poor diagnostic performance., Conclusions: This study provides evidence suggesting that the early response to minor ischaemic stroke is predominantly reflected by changes in the expression of miRNAs in peripheral blood cells. Further work in independent cohorts particularly in patients with more severe stroke is needed to validate these findings and investigate their clinical relevance., (© 2023. The Author(s).)

Cancer is a heterogeneous disease with a strong genetic component making it suitable for precision medicine approaches aimed at identifying the underlying molecular drivers within a tumour. Large scale population-level cancer sequencing consortia have identified many actionable mutations common across both cancer types and sub-types, resulting in an increasing number of successful precision medicine programs. Nonetheless, such approaches fail to consider the effects of mutations unique to an individual patient and may miss rare driver mutations, necessitating personalised approaches to driver-gene prioritisation. One approach is to quantify the functional importance of individual mutations in a single tumour based on how they affect the expression of genes in a gene interaction network (GIN). These GIN-based approaches can be broadly divided into those that utilise an existing reference GIN and those that construct de novo patient-specific GINs. These single-tumour approaches have several limitations that likely influence their results, such as use of reference cohort data, network choice, and approaches to mathematical approximation, and more research is required to evaluate the in vitro and in vivo applicability of their predictions. This review examines the current state of the art methods that identify driver genes in single tumours with a focus on GIN-based driver prioritisation., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Rhys Gillman reports financial support was provided by Tour de Cure Ltd. Lionel Hebbard reports financial support was provided by Townsville Hospital and Health Service. Lionel Hebbard reports financial support was provided by Tropical Australian Academic Health Centre., (© 2023 The Authors.)

The reduced prevalence of insulin resistance and type 2 diabetes in countries with endemic parasitic worm infections suggests a protective role for worms against metabolic disorders, however clinical evidence has been non-existent. This 2-year randomised, double-blinded clinical trial in Australia of hookworm infection in 40 male and female adults at risk of type 2 diabetes assessed the safety and potential metabolic benefits of treatment with either 20 (n = 14) or 40 (n = 13) Necator americanus larvae (L3) or Placebo (n = 13) (Registration ACTRN12617000818336). Primary outcome was safety defined by adverse events and completion rate. Homoeostatic model assessment of insulin resistance, fasting blood glucose and body mass were key secondary outcomes. Adverse events were more frequent in hookworm-treated participants, where 44% experienced expected gastrointestinal symptoms, but completion rates were comparable to Placebo. Fasting glucose and insulin resistance were lowered in both hookworm-treated groups at 1 year, and body mass was reduced after L3-20 treatment at 2 years. This study suggests hookworm infection is safe in people at risk of type 2 diabetes and associated with improved insulin resistance, warranting further exploration of the benefits of hookworms on metabolic health., (© 2023. The Author(s).)

Background: One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long-read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed., Findings: We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on chromosomes 11, 16, 25, and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and 9 previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mitochondrial DNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified 2 differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphologic data, comprising geometric morphometric assessment of cranial morphology, place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue shows she had a larger cranial capacity than a similar-sized domestic dog., Conclusions: These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphologic characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney., (© The Author(s) 2023. Published by Oxford University Press GigaScience.)

Non-tuberculous mycobacterial pulmonary disease (NTM-PD) is a chronic, progressive, and growing worldwide health burden associated with mounting morbidity, mortality, and economic costs. Improvements in NTM-PD management are urgently needed, which requires a better understanding of fundamental immunopathology. Here, we examine temporal dynamics of the immune compartment during NTM-PD caused by Mycobacterium avium complex (MAC) and Mycobactereoides abscessus complex (MABS). We show that active MAC infection is characterized by elevated T cell immunoglobulin and mucin-domain containing-3 expression across multiple T cell subsets. In contrast, active MABS infection was characterized by increased expression of cytotoxic T-lymphocyte-associated protein 4. Patients who failed therapy closely mirrored the healthy individual immune phenotype, with circulating immune network appearing to 'ignore' infection in the lung. Interestingly, immune biosignatures were identified that could inform disease stage and infecting species with high accuracy. Additionally, programmed cell death protein 1 blockade rescued antigen-specific IFN-γ secretion in all disease stages except persistent infection, suggesting the potential to redeploy checkpoint blockade inhibitors for NTM-PD. Collectively, our results provide new insight into species-specific 'immune chatter' occurring during NTM-PD and provide new targets, processes and pathways for diagnostics, prognostics, and treatments needed for this emerging and difficult to treat disease., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Ratnatunga, Tungatt, Proietti, Halstrom, Holt, Lutzky, Price, Doolan, Bell, Field, Kupz, Thomson and Miles.)

Background and Aims: The nacht domain, leucine-rich repeat, and pyrin domain-containing protein 3 (NLRP3) inflammasome is upregulated in human abdominal aortic aneurysm (AAA), but its pathogenic role is unclear. The aims of this study were firstly to examine whether the inflammasome was upregulated in a mouse model of AAA and secondly to test whether the inflammasome inhibitor colchicine limited AAA growth., Methods: AAA was induced in eight-week-old male C57BL6/J mice with topical application of elastase to the infrarenal aorta and oral 3-aminopropionitrile (E-BAPN). For aim one, inflammasome activation, abdominal aortic diameter, and rupture were compared between mice with AAA and sham controls. For aim two, 3 weeks after AAA induction, mice were randomly allocated to receive colchicine ( n = 28, 0.2 mg/kg/d) or vehicle control ( n = 29). The primary outcome was the rate of maximum aortic diameter increase measured by ultrasound over 13 weeks., Results: There was upregulation of NLRP3 markers interleukin- (IL-) 1 β (median, IQR; 15.67, 7.11-22.60 pg/mg protein versus 6.87, 4.54-11.60 pg/mg protein, p = .048) and caspase-1 (109, 83-155 relative luminosity units (RLU) versus 45, 38-65 RLU, p < .001) in AAA samples compared to controls. Aortic diameter increase over 80 days (mean difference, MD, 4.3 mm, 95% CI 3.3, 5.3, p < .001) was significantly greater in mice in which aneurysms were induced compared to sham controls. Colchicine did not significantly limit aortic diameter increase over 80 days (MD -0.1 mm, 95% CI -1.1, 0.86, p = .922)., Conclusions: The inflammasome was activated in this mouse model of AAA; however, daily oral administration of colchicine did not limit AAA growth., Competing Interests: The authors declare that they have no conflicts of interest., (Copyright © 2022 James Phie et al.)

A decline in the prevalence of parasites such as hookworms appears to be correlated with the rise in non-communicable inflammatory conditions in people from high- and middle-income countries. This correlation has led to studies that have identified proteins produced by hookworms that can suppress inflammatory bowel disease (IBD) and asthma in animal models. Hookworms secrete a family of abundant netrin-domain containing proteins referred to as AIPs (Anti-Inflammatory Proteins), but there is no information on the structure-function relationships. Here we have applied a downsizing approach to the hookworm AIPs to derive peptides of 20 residues or less, some of which display anti-inflammatory effects when co-cultured with human peripheral blood mononuclear cells and oral therapeutic activity in a chemically induced mouse model of acute colitis. Our results indicate that a conserved helical region is responsible, at least in part, for the anti-inflammatory effects. This helical region has potential in the design of improved leads for treating IBD and possibly other inflammatory conditions., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Cobos, Bansal, Wilson, Jones, Zhao, Field, Eichenberger, Pickering, Ryan, Ratnatunga, Miles, Ruscher, Giacomin, Navarro, Loukas and Daly.)

Parasitic helminth infections, while a major cause of neglected tropical disease burden, negatively correlate with the incidence of immune-mediated inflammatory diseases such as inflammatory bowel diseases (IBD). To evade expulsion, helminths have developed sophisticated mechanisms to regulate their host's immune responses. Controlled experimental human helminth infections have been assessed clinically for treating inflammatory conditions; however, such a radical therapeutic modality has challenges. An alternative approach is to harness the immunomodulatory properties within the worm's excretory-secretory (ES) complement, its secretome. Here, we report a biologics discovery and validation pipeline to generate and screen in vivo a recombinant cell-free secretome library of helminth-derived immunomodulatory proteins. We successfully expressed 78 recombinant ES proteins from gastrointestinal hookworms and screened the crude in vitro translation reactions for anti-IBD properties in a mouse model of acute colitis. After statistical filtering and ranking, 20 proteins conferred significant protection against various parameters of colitis. Lead candidates from distinct protein families, including annexins, transthyretins, nematode-specific retinol-binding proteins, and SCP/TAPS were identified. Representative proteins were produced in mammalian cells and further validated, including ex vivo suppression of inflammatory cytokine secretion by T cells from IBD patient colon biopsies. Proteins identified herein offer promise as novel, safe, and mechanistically differentiated biologics for treating the globally increasing burden of inflammatory diseases.

Background: Aedes albopictus is a highly invasive species and an important vector of dengue and chikungunya viruses. Indigenous to Southeast Asia, Ae. albopictus has successfully invaded every inhabited continent, except Antarctica, in the past 80 years. Vector surveillance and control at points of entry (PoE) is the most critical front line of defence against the introduction of Ae. albopictus to new areas. Identifying the pathways by which Ae. albopictus are introduced is the key to implementing effective vector surveillance to rapidly detect introductions and to eliminate them., Methods: A literature review was conducted to identify studies and data sources reporting the known and suspected dispersal pathways of human-mediated Ae. albopictus dispersal between 1940-2020. Studies and data sources reporting the first introduction of Ae. albopictus in a new country were selected for data extraction and analyses., Results: Between 1940-2020, Ae. albopictus was reported via various dispersal pathways into 86 new countries. Two main dispersal pathways were identified: (1) at global and continental spatial scales, maritime sea transport was the main dispersal pathway for Ae. albopictus into new countries in the middle to late 20th Century, with ships carrying used tyres of particular importance during the 1980s and 1990s, and (2) at continental and national spatial scales, the passive transportation of Ae. albopictus in ground vehicles and to a lesser extent the trade of used tyres and maritime sea transport appear to be the major drivers of Ae. albopictus dispersal into new countries, especially in Europe. Finally, the dispersal pathways for the introduction and spread of Ae. albopictus in numerous countries remains unknown, especially from the 1990s onwards., Conclusions: This review identified the main known and suspected dispersal pathways of human-mediated Ae. albopictus dispersal leading to the first introduction of Ae. albopictus into new countries and highlighted gaps in our understanding of Ae. albopictus dispersal pathways. Relevant advances in vector surveillance and genomic tracking techniques are presented and discussed in the context of improving vector surveillance., (© 2022. The Author(s).)

Although circumstantial evidence supports enhanced Toll-like receptor 7 (TLR7) signalling as a mechanism of human systemic autoimmune disease 1-7 , evidence of lupus-causing TLR7 gene variants is lacking. Here we describe human systemic lupus erythematosus caused by a TLR7 gain-of-function variant. TLR7 is a sensor of viral RNA 8 , 9 and binds to guanosine 10 - 12 . We identified a de novo, previously undescribed missense TLR7 Y264H variant in a child with severe lupus and additional variants in other patients with lupus. The TLR7 Y264H variant selectively increased sensing of guanosine and 2',3'-cGMP 10-12 , and was sufficient to cause lupus when introduced into mice. We show that enhanced TLR7 signalling drives aberrant survival of B cell receptor (BCR)-activated B cells, and in a cell-intrinsic manner, accumulation of CD11c + age-associated B cells and germinal centre B cells. Follicular and extrafollicular helper T cells were also increased but these phenotypes were cell-extrinsic. Deficiency of MyD88 (an adaptor protein downstream of TLR7) rescued autoimmunity, aberrant B cell survival, and all cellular and serological phenotypes. Despite prominent spontaneous germinal-centre formation in Tlr7 Y264H mice, autoimmunity was not ameliorated by germinal-centre deficiency, suggesting an extrafollicular origin of pathogenic B cells. We establish the importance of TLR7 and guanosine-containing self-ligands for human lupus pathogenesis, which paves the way for therapeutic TLR7 or MyD88 inhibition., (© 2022. The Author(s).)

Dogs are uniquely associated with human dispersal and bring transformational insight into the domestication process. Dingoes represent an intriguing case within canine evolution being geographically isolated for thousands of years. Here, we present a high-quality de novo assembly of a pure dingo (CanFam&#95;DDS). We identified large chromosomal differences relative to the current dog reference (CanFam3.1) and confirmed no expanded pancreatic amylase gene as found in breed dogs. Phylogenetic analyses using variant pairwise matrices show that the dingo is distinct from five breed dogs with 100% bootstrap support when using Greenland wolf as the outgroup. Functionally, we observe differences in methylation patterns between the dingo and German shepherd dog genomes and differences in serum biochemistry and microbiome makeup. Our results suggest that distinct demographic and environmental conditions have shaped the dingo genome. In contrast, artificial human selection has likely shaped the genomes of domestic breed dogs after divergence from the dingo.

Precision medicine programs to identify clinically relevant genetic variation have been revolutionized by access to increasingly affordable high-throughput sequencing technologies. A decade of continual drops in per-base sequencing costs means it is now feasible to sequence an individual patient genome and interrogate all classes of genetic variation for < $1,000 USD. However, while advances in these technologies have greatly simplified the ability to obtain patient sequence information, the timely analysis and interpretation of variant information remains a challenge for the rollout of large-scale precision medicine programs. This review will examine the challenges and potential solutions that exist in identifying predictive genetic biomarkers and pharmacogenetic variants in a patient and discuss the larger bioinformatic challenges likely to emerge in the future. It will examine how both software and hardware development are aiming to overcome issues in short read mapping, variant detection and variant interpretation. It will discuss the current state of the art for genetic disease and the remaining challenges to overcome for complex disease. Success across all types of disease will require novel statistical models and software in order to ensure precision medicine programs realize their full potential now and into the future., Competing Interests: The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Field.)

About half of the world's population and 80% of the world's biodiversity can be found in the tropics. Many diseases are specific to the tropics, with at least 41 diseases caused by endemic bacteria, viruses, parasites, and fungi. Such diseases are of increasing concern, as the geographic range of tropical diseases is expanding due to climate change, urbanization, change in agricultural practices, deforestation, and loss of biodiversity. While traditional medicines have been used for centuries in the treatment of tropical diseases, the active natural compounds within these medicines remain largely unknown. In this review, we describe infectious diseases specific to the tropics, including their causative pathogens, modes of transmission, recent major outbreaks, and geographic locations. We further review current treatments for these tropical diseases, carefully consider the biodiscovery potential of the tropical biome, and discuss a range of technologies being used for drug development from natural resources. We provide a list of natural products with antimicrobial activity, detailing the source organisms and their effectiveness as treatment. We discuss how technological advancements, such as next-generation sequencing, are driving high-throughput natural product screening pipelines to identify compounds with therapeutic properties. This review demonstrates the impact natural products from the vast tropical biome have in the treatment of tropical infectious diseases and how high-throughput technical capacity will accelerate this discovery process.

Tuberculosis (TB) is the leading infectious cause of death globally. The only licensed TB vaccine, Bacille Calmette-Guérin (BCG), has low efficacy against TB in adults and is not recommended in people with impaired immunity. The incorporation of the Mycobacterium tuberculosis (Mtb) secretion system ESX-1 into BCG improves immunogenicity and protection against TB in animal models, which is associated with the secretion of the ESX-1-dependent protein ESAT-6. However, the resulting strain, BCG::ESX1 Mtb , has been deemed unsafe as a human vaccine, due to prolonged persistence and increased virulence in immunocompromised mice. In this study, we describe a new recombinant BCG strain that uncouples the beneficial aspects of ESAT-6 secretion from the detrimental ESX-1effects on virulence and persistence. The strain was constructed by fusing the ESAT-6-encoding gene esxA to the general secretion signal for the mycobacterial type VII secretion pathway protein PE25. This new strain, BCG::ESAT6-PE25SS, secretes full-length ESAT-6 via the ESX-5 secretion system, which in contrast to ESX-1 is also present in BCG. In vivo testing revealed that ESX-5-targeted ESAT-6 export, induces cytosolic contact, generates ESAT-6-specific T cells and enhances the protective efficacy against TB disease, but is associated with low virulence and reduced persistence in immunocompetent and immunocompromised mice. Additionally, compared to BCG::ESX1 Mtb and parental BCG, mucosal administration of BCG::ESAT6-PE25SS is associated with more rapid clearance from the lung. These results warrant further studies to evaluate BCG::ESAT6-PE25SS as a potential live attenuated vaccine candidate for TB., Competing Interests: Declaration of Competing Interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Andreas Kupz reports financial support was provided by Australian National Health and Medical Research Council. Aliabbas Husain reports financial support was provided by Indian Government. Andreas Kupz reports a relationship with Collaboration for Tuberculosis Vaccine Discovery (Bill and Melinda Gates Foundation) that includes: board membership and travel reimbursement. Andreas Kupz has patent filing # number 2021900320 pending to James Cook University., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

This study documents disparities in the mental health burden of young adults in Rhode Island during the COVID-19 pandemic as it pertains to essential worker status, sexual orientation, gender identity, and childhood trauma. A cross-sectional web-based survey of young adults aged 18 to 25 years conducted between May and October of 2020 assessed anxiety symptoms, depressive symptoms, and thoughts of suicide. In fully adjusted regression models, (n =528 young adults) being an essential worker, a sexual minority, having lower relative SES, and having more adverse childhood experiences were significantly associated with negative mental health outcomes. In models adjusted for individual ACEs, exposure to mental illness in the household, physical violence between adults in the household, emotional abuse, and sexual abuse were independently associated with all three outcomes. Young adults most vulnerable and stigmatized in the community are also those who are most severely affected in terms of mental health.

Background: Sensitive diagnostics are needed for effective management and surveillance of schistosomiasis so that current transmission interruption goals set by WHO can be achieved. We aimed to screen the Schistosoma haematobium secretome to find antibody biomarkers of schistosome infection, validate their diagnostic performance in samples from endemic populations, and evaluate their utility as point of care immunochromatographic tests (POC-ICTs) to diagnose urogenital schistosomiasis in the field., Methods: We did a biomarker identification study, in which we constructed a proteome array containing 992 validated and predicted proteins from S haematobium and screened it with serum and urine antibodies from endemic populations in Gabon, Tanzania, and Zimbabwe. Arrayed antigens that were IgG-reactive and a select group of antigens from the worm extracellular vesicle proteome, predicted to be diagnostically informative, were then evaluated by ELISA using the same samples used to probe arrays, and samples from individuals residing in a low-endemicity setting (ie, Pemba and Unguja islands, Zanzibar, Tanzania). The two most sensitive and specific antigens were incorporated into POC-ICTs to assess their ability to diagnose S haematobium infection from serum in a field-deployable format., Findings: From array probing, in individuals who were infected, 208 antigens were the targets of significantly elevated IgG responses in serum and 45 antigens were the targets of significantly elevated IgG responses in urine. Of the five proteins that were validated by ELISA, Sh -TSP-2 (area under the curve [AUC] serum =0·98 [95% CI 0·95-1·00]; AUC urine =0·96 [0·93-0·99]), and MS3&#95;01370 (AUC serum =0·93 [0·89-0·97]; AUC urine =0·81 [0·72-0·89]) displayed the highest overall diagnostic performance in each biofluid and exceeded that of S haematobium -soluble egg antigen in urine (AUC=0·79 [0·69-0·90]). When incorporated into separate POC-ICTs, Sh -TSP-2 showed absolute specificity and a sensitivity of 75% and MS3&#95;01370 showed absolute specificity and a sensitivity of 89%., Interpretation: We identified numerous biomarkers of urogenital schistosomiasis that could form the basis of novel antibody diagnostics for this disease. Two of these antigens, Sh -TSP-2 and MS3&#95;01370, could be used as sensitive, specific, and field-deployable diagnostics to support schistosomiasis control and elimination initiatives, with particular focus on post-elimination surveillance., Funding: Australian Trade and Investment Commission and Merck Global Health Institute., Competing Interests: MSP, AL, and JS are inventors on a provisional patent filed by the funding body that captures the data presented in this study. All other authors declare no competing interests., (© 2021 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY-NC-ND 4.0 license.)

Prostate cancer (PrCa) is highly heritable, and although rare variants contribute significantly to PrCa risk, few have been identified to date. Herein, whole-genome sequencing was performed in a large PrCa family featuring multiple affected relatives spanning several generations. A rare, predicted splice site EZH2 variant, rs78589034 (G > A), was identified as segregating with disease in all but two individuals in the family, one of whom was affected with lymphoma and bowel cancer and a female relative. This variant was significantly associated with disease risk in combined familial and sporadic PrCa datasets (n = 1551; odds ratio [OR] = 3.55, P = 1.20 × 10 -5 ). Transcriptome analysis was performed on prostate tumour needle biopsies available for two rare variant carriers and two wild-type cases. Although no allele-dependent differences were detected in EZH2 transcripts, a distinct differential gene expression signature was observed when comparing prostate tissue from the rare variant carriers with the wild-type samples. The gene expression signature comprised known downstream targets of EZH2 and included the top-ranked genes, DUSP1, FOS, JUNB and EGR1, which were subsequently validated by qPCR. These data provide evidence that rs78589034 is associated with increased PrCa risk in Tasmanian men and further, that this variant may be associated with perturbed EZH2 function in prostate tissue. Disrupted EZH2 function is a driver of tumourigenesis in several cancers, including prostate, and is of significant interest as a therapeutic target., (© 2021 Union for International Cancer Control.)

Background: Pharmacogenetic variation is important to drug responses through diverse and complex mechanisms. Predictions of the functional impact of missense pharmacogenetic variants primarily rely on the degree of sequence conservation between species as a primary discriminator. However, idiosyncratic or off-target drug-variant interactions sometimes involve effects that are peripheral or accessory to the central systems in which a gene functions. Given the importance of sequence conservation to functional prediction tools-these idiosyncratic pharmacogenetic variants may violate the assumptions of predictive software commonly used to infer their effect., Methods: Here we exhaustively assess the effectiveness of eleven missense mutation functional inference tools on all known pharmacogenetic missense variants contained in the Pharmacogenomics Knowledgebase (PharmGKB) repository. We categorize PharmGKB entries into sub-classes to catalog likely off-target interactions, such that we may compare predictions across different variant annotations., Results: As previously demonstrated, functional inference tools perform variably across the complete set of PharmGKB variants, with large numbers of variants incorrectly classified as 'benign'. However, we find substantial differences amongst PharmGKB variant sub-classes, particularly in variants known to cause off-target, type B adverse drug reactions, that are largely unrelated to the main pharmacological action of the drug. Specifically, variants associated with off-target effects (hence referred to as off-target variants) were most often incorrectly classified as 'benign'. These results highlight the importance of understanding the underlying mechanism of pharmacogenetic variants and how variants associated with off-target effects will ultimately require new predictive algorithms., Conclusion: In this work we demonstrate that functional inference tools perform poorly on pharmacogenetic variants, particularly on subsets enriched for variants causing off-target, type B adverse drug reactions. We describe how to identify variants associated with off-target effects within PharmGKB in order to generate a training set of variants that is needed to develop new algorithms specifically for this class of variant. Development of such tools will lead to more accurate functional predictions and pave the way for the increased wide-spread adoption of pharmacogenetics in clinical practice., Competing Interests: The authors declare that they have no competing interests., (© 2021 McConnell et al.)

Background and Aims: Abdominal aortic aneurysm (AAA) is an important cause of death worldwide and has an estimated heritability between 70 and 77%. Genome-wide association studies (GWAS) are an established way to discover genetic risk variants. The aim of this study was to systematically review the findings and quality of previous AAA GWAS., Methods: The Medline, PubMed, Web of Science and relevant genetic databases were searched to identify previous AAA GWAS. A framework was developed to grade the methodological quality of the GWAS. Data from included studies were extracted to assess methods and findings., Results: Eight case-control studies were included. Thirty-three of the 38 total single nucleotide polymorphisms (SNPs) previously reported were associated with AAA diagnosis at genome-wide significance (p < 5.0 × 10 -8 ). The CDKN2B antisense RNA-1 gene had the most significant association with AAA diagnosis (p = 6.94 × 10 -29 and p = 1.54 × 10 -33 for rs4007642 and rs10757274 respectively). Age, sex and smoking history were not reported for the complete cohort in any of the included studies, although five of the eight studies adjusted or matched for at least two confounding variables. All included studies had important design limitations including lack of sample size estimation, inconsistent case and control ascertainment and limited phenotyping of the AAAs. AAA growth was assessed in one GWAS, however, no significant associations with the reported SNPs were found., Conclusions: This systematic review identified 33 SNPs associated with AAA diagnosis at genome-wide significance previously validated in multiple cohorts. The association between SNPs and AAA growth was not adequately examined. Previous GWAS have a number of design limitations., (Copyright © 2021. Published by Elsevier B.V.)

Tuberculosis (TB) is one of the deadliest infectious diseases in the world. The metabolic disease type 2 diabetes (T2D) significantly increases the risk of developing active TB. Effective new TB vaccine candidates and novel therapeutic interventions are required to meet the challenges of global TB eradication. Recent evidence suggests that the microbiota plays a significant role in how the host responds to infection, injury and neoplastic changes. Animal models that closely reflect human physiology are crucial in assessing new treatments and to decipher the underlying immunological defects responsible for increased TB susceptibility in comorbid patients. In this study, using a diet-induced murine T2D model that reflects the etiopathogenesis of clinical T2D and increased TB susceptibility, we investigated how the intestinal microbiota may impact the development of T2D, and how the gut microbial composition changes following a very low-dose aerosol infection with Mycobacterium tuberculosis ( Mtb ). Our data revealed a substantial intestinal microbiota dysbiosis in T2D mice compared to non-diabetic animals. The observed differences were comparable to previous clinical reports in TB patients, in which it was shown that Mtb infection causes rapid loss of microbial diversity. Furthermore, diversity index and principle component analyses demonstrated distinct clustering of Mtb -infected non-diabetic mice vs. Mtb -infected T2D mice. Our findings support a broad applicability of T2D mice as a tractable small animal model for studying distinct immune parameters, microbiota and the immune-metabolome of TB/T2D comorbidity. This model may also enable answers to be found to critical outstanding questions about targeted interventions of the gut microbiota and the gut-lung axis., Competing Interests: The authors declare no conflicts of interest., (© 2021 The Authors. Animal Models and Experimental Medicine published by John Wiley & Sons Australia, Ltd on behalf of The Chinese Association for Laboratory Animal Sciences.)

Background: Basenjis are considered an ancient dog breed of central African origins that still live and hunt with tribesmen in the African Congo. Nicknamed the barkless dog, Basenjis possess unique phylogeny, geographical origins and traits, making their genome structure of great interest. The increasing number of available canid reference genomes allows us to examine the impact the choice of reference genome makes with regard to reference genome quality and breed relatedness., Results: Here, we report two high quality de novo Basenji genome assemblies: a female, China (CanFam&#95;Bas), and a male, Wags. We conduct pairwise comparisons and report structural variations between assembled genomes of three dog breeds: Basenji (CanFam&#95;Bas), Boxer (CanFam3.1) and German Shepherd Dog (GSD) (CanFam&#95;GSD). CanFam&#95;Bas is superior to CanFam3.1 in terms of genome contiguity and comparable overall to the high quality CanFam&#95;GSD assembly. By aligning short read data from 58 representative dog breeds to three reference genomes, we demonstrate how the choice of reference genome significantly impacts both read mapping and variant detection., Conclusions: The growing number of high-quality canid reference genomes means the choice of reference genome is an increasingly critical decision in subsequent canid variant analyses. The basal position of the Basenji makes it suitable for variant analysis for targeted applications of specific dog breeds. However, we believe more comprehensive analyses across the entire family of canids is more suited to a pangenome approach. Collectively this work highlights the importance the choice of reference genome makes in all variation studies.

Aims: In 2003, an Australian woman was convicted by a jury of smothering and killing her four children over a 10-year period. Each child died suddenly and unexpectedly during a sleep period, at ages ranging from 19 days to 18 months. In 2019 we were asked to investigate if a genetic cause could explain the children's deaths as part of an inquiry into the mother's convictions., Methods and Results: Whole genomes or exomes of the mother and her four children were sequenced. Functional analysis of a novel CALM2 variant was performed by measuring Ca2+-binding affinity, interaction with calcium channels and channel function. We found two children had a novel calmodulin variant (CALM2 G114R) that was inherited maternally. Three genes (CALM1-3) encode identical calmodulin proteins. A variant in the corresponding residue of CALM3 (G114W) was recently reported in a child who died suddenly at age 4 and a sibling who suffered a cardiac arrest at age 5. We show that CALM2 G114R impairs calmodulin's ability to bind calcium and regulate two pivotal calcium channels (CaV1.2 and RyR2) involved in cardiac excitation contraction coupling. The deleterious effects of G114R are similar to those produced by G114W and N98S, which are considered arrhythmogenic and cause sudden cardiac death in children., Conclusion: A novel functional calmodulin variant (G114R) predicted to cause idiopathic ventricular fibrillation, catecholaminergic polymorphic ventricular tachycardia, or mild long QT syndrome was present in two children. A fatal arrhythmic event may have been triggered by their intercurrent infections. Thus, calmodulinopathy emerges as a reasonable explanation for a natural cause of their deaths., (© The Author(s) 2020. Published by Oxford University Press on behalf of the European Society of Cardiology.)

Population genomic approaches can characterize dispersal across a single generation through to many generations in the past, bridging the gap between individual movement and intergenerational gene flow. These approaches are particularly useful when investigating dispersal in recently altered systems, where they provide a way of inferring long-distance dispersal between newly established populations and their interactions with existing populations. Human-mediated biological invasions represent such altered systems which can be investigated with appropriate study designs and analyses. Here we apply temporally restricted sampling and a range of population genomic approaches to investigate dispersal in a 2004 invasion of Aedes albopictus (the Asian tiger mosquito) in the Torres Strait Islands (TSI) of Australia. We sampled mosquitoes from 13 TSI villages simultaneously and genotyped 373 mosquitoes at genome-wide single nucleotide polymorphisms (SNPs): 331 from the TSI, 36 from Papua New Guinea (PNG) and four incursive mosquitoes detected in uninvaded regions. Within villages, spatial genetic structure varied substantially but overall displayed isolation by distance and a neighbourhood size of 232-577. Close kin dyads revealed recent movement between islands 31-203 km apart, and deep learning inferences showed incursive Ae. albopictus had travelled to uninvaded regions from both adjacent and nonadjacent islands. Private alleles and a co-ancestry matrix indicated direct gene flow from PNG into nearby islands. Outlier analyses also detected four linked alleles introgressed from PNG, with the alleles surrounding 12 resistance-associated cytochrome P450 genes. By treating dispersal as both an intergenerational process and a set of discrete events, we describe a highly interconnected invasive system., (© 2021 John Wiley & Sons Ltd.)

The distinct properties of allo-reactive T-cell repertoires are not well understood. To investigate whether auto-reactive and allo-reactive T-cell repertoires encoded distinct properties, we used dextramer enumeration, enrichment, single-cell T-cell receptor (TCR) sequencing and multiparameter analysis. We found auto-reactive and allo-reactive T-cells differed in mean ex vivo frequency which was antigen dependent. Allo-reactive T-cells showed clear differences in TCR architecture, with enriched usage of specific T-cell receptor variable ( TRBJ ) genes and broader use of T-cell receptor variable joining ( TRBJ ) genes. Auto-reactive T-cell repertoires exhibited complementary determining regions three (CDR3) lengths using a Gaussian distribution whereas allo-reactive T-cell repertoires exhibited distorted patterns in CDR3 length. CDR3 loops from allo-reactive T-cells showed distinct physical-chemical properties, tending to encode loops that were more acidic in charge. Allo-reactive T-cell repertoires differed in diversity metrics, tending to show increased overall diversity and increased homogeneity between repertoires. Motif analysis of CDR3 loops showed allo-reactive T-cell repertoires differed in motif preference which included broader motif use. Collectively, these data conclude that allo-reactive T-cell repertoires are indeed different to auto-reactive repertoires and provide tangible metrics for further investigations and validation. Given that the antigens studied here are overexpressed on multiple cancers and that allo-reactive TCRs often show increased ligand affinity, this new TCR bank also has translational potential for adoptive cell therapy, soluble TCR-based therapy and rational TCR design.

Sequencing the first human genome in 2003 took 15 years and cost $2.7 billion. Advances in sequencing technologies have since decreased costs to the point where it is now feasible to resequence a whole human genome for $1000 in a single day. These advances have allowed the generation of huge volumes of high-quality human sequence data used to construct increasingly large catalogs of both population-level and disease-causing variation. The existence of such databases, coupled with a high-quality human reference genome, means we are able to interrogate and annotate all types of genetic variation and identify pathogenic variants for many diseases. Increasingly, sequencing-based approaches are being used to elucidate the underlying genetic cause of autoimmune diseases, a group of roughly 80 polygenic diseases characterized by abnormal immune responses where healthy tissue is attacked. Although sequence data generation has become routine and affordable, significant challenges remain with no gold-standard methodology to identify pathogenic variants currently available. This review examines the latest methodologies used to identify pathogenic variants in autoimmune diseases and considers available sequencing options and subsequent bioinformatic methodologies and strategies. The development of reliable and robust sequencing and analytic workflows to detect pathogenic variants is critical to realize the potential of precision medicine programs where patient variant information is used to inform clinical practice., (© 2020 The Authors. Immunology & Cell Biology published by John Wiley & Sons Australia, Ltd on behalf of Australian and New Zealand Society for Immunology, Inc.)

Increasingly affordable sequencing technologies are revolutionizing the field of genomic medicine. It is now feasible to interrogate all major classes of variation in an individual across the entire genome for less than $1000 USD. While the generation of patient sequence information using these technologies has become routine, the analysis and interpretation of this data remains the greatest obstacle to widespread clinical implementation. This chapter summarizes the steps to identify, annotate, and prioritize variant information required for clinical report generation. We discuss methods to detect each variant class and describe strategies to increase the likelihood of detecting causal variant(s) in Mendelian disease. Lastly, we describe a sample workflow for synthesizing large amount of genetic information into concise clinical reports.

Genetic signatures caused by demographic and adaptive processes during past climatic shifts can inform predictions of species' responses to anthropogenic climate change. To identify these signatures in Acropora tenuis , a reef-building coral threatened by global warming, we first assembled the genome from long reads and then used shallow whole-genome resequencing of 150 colonies from the central inshore Great Barrier Reef to inform population genomic analyses. We identify population structure in the host that reflects a Pleistocene split, whereas photosymbiont differences between reefs most likely reflect contemporary (Holocene) conditions. Signatures of selection in the host were associated with genes linked to diverse processes including osmotic regulation, skeletal development, and the establishment and maintenance of symbiosis. Our results suggest that adaptation to post-glacial climate change in A. tenuis has involved selection on many genes, while differences in symbiont specificity between reefs appear to be unrelated to host population structure., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

Tuberculosis (TB) claims 1.5 million lives per year. This situation is largely due to the low efficacy of the only licensed TB vaccine, Bacillus Calmette-Guérin (BCG) against pulmonary TB. The metabolic disease type 2 diabetes (T2D) is a risk factor for TB and the mechanisms underlying increased TB susceptibility in T2D are not well understood. Furthermore, it is unknown if new TB vaccines will provide protection in the context of T2D. Here we used a diet-induced murine model of T2D to investigate the underlying mechanisms of TB/T2D comorbidity and to evaluate the protective capacity of two experimental TB vaccines in comparison to conventional BCG. Our data reveal a distinct immune dysfunction that is associated with diminished recognition of mycobacterial antigens in T2D. More importantly, we provide compelling evidence that mucosal delivery of recombinant BCG strains expressing the Mycobacterium tuberculosis ( Mtb ) ESX-1 secretion system (BCG::RD1 and BCG::RD1 ESAT-6 ∆92-95) are safe and confer superior immunity against aerosol Mtb infection in the context of T2D. Our findings suggest that the remarkable anti-TB immunity by these recombinant BCG strains is achieved via augmenting the numbers and functional capacity of antigen presenting cells in the lungs of diabetic mice., Competing Interests: The authors declare no competing interest.

Background: The German Shepherd Dog (GSD) is one of the most common breeds on earth and has been bred for its utility and intelligence. It is often first choice for police and military work, as well as protection, disability assistance, and search-and-rescue. Yet, GSDs are well known to be susceptible to a range of genetic diseases that can interfere with their training. Such diseases are of particular concern when they occur later in life, and fully trained animals are not able to continue their duties., Findings: Here, we provide the draft genome sequence of a healthy German Shepherd female as a reference for future disease and evolutionary studies. We generated this improved canid reference genome (CanFam&#95;GSD) utilizing a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. The GSD assembly is ∼80 times as contiguous as the current canid reference genome (20.9 vs 0.267 Mb contig N50), containing far fewer gaps (306 vs 23,876) and fewer scaffolds (429 vs 3,310) than the current canid reference genome CanFamv3.1. Two chromosomes (4 and 35) are assembled into single scaffolds with no gaps. BUSCO analyses of the genome assembly results show that 93.0% of the conserved single-copy genes are complete in the GSD assembly compared with 92.2% for CanFam v3.1. Homology-based gene annotation increases this value to ∼99%. Detailed examination of the evolutionarily important pancreatic amylase region reveals that there are most likely 7 copies of the gene, indicative of a duplication of 4 ancestral copies and the disruption of 1 copy., Conclusions: GSD genome assembly and annotation were produced with major improvement in completeness, continuity, and quality over the existing canid reference. This resource will enable further research related to canine diseases, the evolutionary relationships of canids, and other aspects of canid biology., (© The Author(s) 2020. Published by Oxford University Press.)

Pathogenic autoantibodies arise in many autoimmune diseases, but it is not understood how the cells making them evade immune checkpoints. Here, single-cell multi-omics analysis demonstrates a shared mechanism with lymphoid malignancy in the formation of public rheumatoid factor autoantibodies responsible for mixed cryoglobulinemic vasculitis. By combining single-cell DNA and RNA sequencing with serum antibody peptide sequencing and antibody synthesis, rare circulating B lymphocytes making pathogenic autoantibodies were found to comprise clonal trees accumulating mutations. Lymphoma driver mutations in genes regulating B cell proliferation and V(D)J mutation (CARD11, TNFAIP3, CCND3, ID3, BTG2, and KLHL6) were present in rogue B cells producing the pathogenic autoantibody. Antibody V(D)J mutations conferred pathogenicity by causing the antigen-bound autoantibodies to undergo phase transition to insoluble aggregates at lower temperatures. These results reveal a pre-neoplastic stage in human lymphomagenesis and a cascade of somatic mutations leading to an iconic pathogenic autoantibody., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Extensive evaluation of RNA-seq methods have demonstrated that no single algorithm consistently outperforms all others. Removal of unwanted variation (RUV) has also been proposed as a method for stabilizing differential expression (DE) results. Despite this, it remains a challenge to run multiple RNA-seq algorithms to identify significant differences common to multiple algorithms, whilst also integrating and assessing the impact of RUV into all algorithms. consensusDE was developed to automate the process of identifying significant DE by combining the results from multiple algorithms with minimal user input and with the option to automatically integrate RUV. consensusDE only requires a table describing the sample groups, a directory containing BAM files or preprocessed count tables and an optional transcript database for annotation. It supports merging of technical replicates, paired analyses and outputs a compendium of plots to guide the user in subsequent analyses. Herein, we assess the ability of RUV to improve DE stability when combined with multiple algorithms and between algorithms, through application to real and simulated data. We find that, although RUV increased fold change stability between algorithms, it demonstrated improved FDR in a setting of low replication for the intersect, the effect was algorithm specific and diminished with increased replication, reinforcing increased replication for recovery of true DE genes. We finish by offering some rules and considerations for the application of RUV in a consensus-based setting. consensusDE is freely available, implemented in R and available as a Bioconductor package, under the GPL-3 license, along with a comprehensive vignette describing functionality: http://bioconductor.org/packages/consensusDE/., Competing Interests: The authors declare there are no competing interests., (©2019 Waardenberg and Field.)

Background: Short-read resequencing of genomes produces abundant information of the genetic variation of individuals. Due to their numerous nature, these variants are rarely exhaustively validated. Furthermore, low levels of undetected variant miscalling will have a systematic and disproportionate impact on the interpretation of individual genome sequence information, especially should these also be carried through into in reference databases of genomic variation., Results: We find that sequence variation from short-read sequence data is subject to recurrent-yet-intermittent miscalling that occurs in a sequence intrinsic manner and is very sensitive to sequence read length. The miscalls arise from difficulties aligning short reads to redundant genomic regions, where the rate of sequencing error approaches the sequence diversity between redundant regions. We find the resultant miscalled variants to be sensitive to small sequence variations between genomes, and thereby are often intrinsic to an individual, pedigree, strain or human ethnic group. In human exome sequences, we identify 2-300 recurrent false positive variants per individual, almost all of which are present in public databases of human genomic variation. From the exomes of non-reference strains of inbred mice, we identify 3-5000 recurrent false positive variants per mouse - the number of which increasing with greater distance between an individual mouse strain and the reference C57BL6 mouse genome. We show that recurrently miscalled variants may be reproduced for a given genome from repeated simulation rounds of read resampling, realignment and recalling. As such, it is possible to identify more than two-thirds of false positive variation from only ten rounds of simulation., Conclusion: Identification and removal of recurrent false positive variants from specific individual variant sets will improve overall data quality. Variant miscalls arising are highly sequence intrinsic and are often specific to an individual, pedigree or ethnicity. Further, read length is a strong determinant of whether given false variants will be called for any given genome - which has profound significance for cohort studies that pool datasets collected and sequenced at different points in time.

Metagenomic sequencing is an increasingly common tool in environmental and biomedical sciences yet analysis workflows remain immature relative to other field such as DNASeq and RNASeq analysis pipelines.  While software for detailing the composition of microbial communities using 16S rRNA marker genes is constantly improving, increasingly researchers are interested in identifying changes exhibited within microbial communities under differing environmental conditions. In order to gain maximum value from metagenomic sequence data we must improve the existing analysis environment by providing accessible and scalable computational workflows able to generate reproducible results. Here we describe a complete end-to-end open-source metagenomics workflow running within Galaxy for 16S differential abundance analysis. The workflow accepts 454 or Illumina sequence data (either overlapping or non-overlapping paired end reads) and outputs lists of the operational taxonomic unit (OTUs) exhibiting the greatest change under differing conditions. A range of analysis steps and graphing options are available giving users a high-level of control over their data and analyses. Additionally, users are able to input complex sample-specific metadata information which can be incorporated into differential analysis and used for grouping / colouring within graphs.  Detailed tutorials containing sample data and existing workflows are available for three different input types: overlapping and non-overlapping read pairs as well as for pre-generated Biological Observation Matrix (BIOM) files. Using the Galaxy platform we developed MetaDEGalaxy, a complete metagenomics differential abundance analysis workflow. MetaDEGalaxy is designed for bench scientists working with 16S data who are interested in comparative metagenomics.  MetaDEGalaxy builds on momentum within the wider Galaxy metagenomics community with the hope that more tools will be added as existing methods mature., Competing Interests: No competing interests were disclosed., (Copyright: © 2019 Thang MWC et al.)

Systemic lupus erythematosus (SLE) is the prototypic systemic autoimmune disease. It is thought that many common variant gene loci of weak effect act additively to predispose to common autoimmune diseases, while the contribution of rare variants remains unclear. Here we describe that rare coding variants in lupus-risk genes are present in most SLE patients and healthy controls. We demonstrate the functional consequences of rare and low frequency missense variants in the interacting proteins BLK and BANK1, which are present alone, or in combination, in a substantial proportion of lupus patients. The rare variants found in patients, but not those found exclusively in controls, impair suppression of IRF5 and type-I IFN in human B cell lines and increase pathogenic lymphocytes in lupus-prone mice. Thus, rare gene variants are common in SLE and likely contribute to genetic risk.

Genetic mutations account for many devastating early onset immune deficiencies. In contrast, less severe and later onset immune diseases, including in patients with no prior family history, remain poorly understood. Whole exome sequencing in two cohorts of such patients identified a novel heterozygous de novo IKBKB missense mutation (c.607G>A) in two separate kindreds in whom probands presented with immune dysregulation, combined T and B cell deficiency, inflammation, and epithelial defects. IKBKB encodes IKK2, which activates NF-κB signaling. IKK2 V203I results in enhanced NF-κB signaling, as well as T and B cell functional defects. IKK2 V203 is a highly conserved residue, and to prove causation, we generated an accurate mouse model by introducing the precise orthologous codon change in Ikbkb using CRISPR/Cas9. Mice and humans carrying this missense mutation exhibit remarkably similar cellular and biochemical phenotypes. Accurate mouse models engineered by CRISPR/Cas9 can help characterize novel syndromes arising from de novo germline mutations and yield insight into pathogenesis., (© 2018 Crown copyright. The government of Australia, Canada, or the UK ("the Crown") owns the copyright interests of authors who are government employees. The Crown Copyright is not transferable.)

Gastrointestinal (GI) parasites, hookworms in particular, have evolved to cause minimal harm to their hosts, allowing them to establish chronic infections. This is mediated by creating an immunoregulatory environment. Indeed, hookworms are such potent suppressors of inflammation that they have been used in clinical trials to treat inflammatory bowel diseases (IBD) and celiac disease. Since the recent description of helminths (worms) secreting extracellular vesicles (EVs), exosome-like EVs from different helminths have been characterized and their salient roles in parasite-host interactions have been highlighted. Here, we analyze EVs from the rodent parasite Nippostrongylus brasiliensis , which has been used as a model for human hookworm infection. N. brasiliensis EVs ( Nb -EVs) are actively internalized by mouse gut organoids, indicating a role in driving parasitism. We used proteomics and RNA-Seq to profile the molecular composition of Nb -EVs. We identified 81 proteins, including proteins frequently present in exosomes (like tetraspanin, enolase, 14-3-3 protein, and heat shock proteins), and 27 sperm-coating protein-like extracellular proteins. RNA-Seq analysis revealed 52 miRNA species, many of which putatively map to mouse genes involved in regulation of inflammation. To determine whether GI nematode EVs had immunomodulatory properties, we assessed their potential to suppress GI inflammation in a mouse model of inducible chemical colitis. EVs from N. brasiliensis but not those from the whipworm Trichuris muris or control vesicles from grapes protected against colitic inflammation in the gut of mice that received a single intraperitoneal injection of EVs. Key cytokines associated with colitic pathology (IL-6, IL-1β, IFNγ, and IL-17a) were significantly suppressed in colon tissues from EV-treated mice. By contrast, high levels of the anti-inflammatory cytokine IL-10 were detected in Nb -EV-treated mice. Proteins and miRNAs contained within helminth EVs hold great potential application in development of drugs to treat helminth infections as well as chronic non-infectious diseases resulting from a dysregulated immune system, such as IBD.

Whipworms are parasitic nematodes that live in the gut of more than 500 million people worldwide. Owing to the difficulty in obtaining parasite material, the mouse whipworm Trichuris muris has been extensively used as a model to study human whipworm infections. These nematodes secrete a multitude of compounds that interact with host tissues where they orchestrate a parasitic existence. Herein we provide the first comprehensive characterization of the excretory/secretory products of T. muris . We identify 148 proteins secreted by T. muris and show for the first time that the mouse whipworm secretes exosome-like extracellular vesicles (EVs) that can interact with host cells. We use an Optiprep® gradient to purify the EVs, highlighting the suitability of this method for purifying EVs secreted by a parasitic nematode. We also characterize the proteomic and genomic content of the EVs, identifying >350 proteins, 56 miRNAs (22 novel) and 475 full-length mRNA transcripts mapping to T. muris gene models. Many of the miRNAs putatively mapped to mouse genes are involved in regulation of inflammation, implying a role in parasite-driven immunomodulation. In addition, for the first time to our knowledge, colonic organoids have been used to demonstrate the internalization of parasite EVs by host cells. Understanding how parasites interact with their host is crucial to develop new control measures. This first characterization of the proteins and EVs secreted by T. muris provides important information on whipworm-host communication and forms the basis for future studies., Competing Interests: No potential conflict of interest was reported by the authors.

The HOXB13 G84E variant is associated with risk of prostate cancer (PCa), however the role this variant plays in PCa development is unknown. This study examined 751 cases, 450 relatives and 355 controls to determine the contribution of this variant to PCa risk in Tasmania and investigated HOXB13 gene and protein expression in tumours from nine G84E heterozygote variant and 13 wild-type carriers. Quantitative PCR and immunohistochemistry showed that HOXB13 gene and protein expression did not differ between tumour samples from variant and wild-type carriers. Allele-specific transcription revealed that two of seven G84E carriers transcribed both the variant and wild-type allele, while five carriers transcribed the wild-type allele. Methylation of surrounding CpG sites was lower in the variant compared to the wild-type allele, however overall methylation across the region was very low. Notably, tumour characteristics were less aggressive in the two variant carriers that transcribed the variant allele compared to the five that did not. This study has shown that HOXB13 expression does not differ between tumour tissue of G84E variant carriers and non-carriers. Intriguingly, the G84E variant allele was rarely transcribed in carriers, suggesting that HOXB13 expression may be driven by the wild-type allele in the majority of carriers.

Melanoma of the skin is a common cancer only in Europeans, whereas it arises in internal body surfaces (mucosal sites) and on the hands and feet (acral sites) in people throughout the world. Here we report analysis of whole-genome sequences from cutaneous, acral and mucosal subtypes of melanoma. The heavily mutated landscape of coding and non-coding mutations in cutaneous melanoma resolved novel signatures of mutagenesis attributable to ultraviolet radiation. However, acral and mucosal melanomas were dominated by structural changes and mutation signatures of unknown aetiology, not previously identified in melanoma. The number of genes affected by recurrent mutations disrupting non-coding sequences was similar to that affected by recurrent mutations to coding sequences. Significantly mutated genes included BRAF, CDKN2A, NRAS and TP53 in cutaneous melanoma, BRAF, NRAS and NF1 in acral melanoma and SF3B1 in mucosal melanoma. Mutations affecting the TERT promoter were the most frequent of all; however, neither they nor ATRX mutations, which correlate with alternative telomere lengthening, were associated with greater telomere length. Most melanomas had potentially actionable mutations, most in components of the mitogen-activated protein kinase and phosphoinositol kinase pathways. The whole-genome mutation landscape of melanoma reveals diverse carcinogenic processes across its subtypes, some unrelated to sun exposure, and extends potential involvement of the non-coding genome in its pathogenesis.

Background. Massively parallel sequencing technology is being used to sequence highly diverse populations of DNA such as that derived from heterogeneous cell mixtures containing both wild-type and disease-related states. At the core of such molecule tagging techniques is the tagging and identification of sequence reads derived from individual input DNA molecules, which must be first computationally disambiguated to generate read groups sharing common sequence tags, with each read group representing a single input DNA molecule. This disambiguation typically generates huge numbers of reads groups, each of which requires additional variant detection analysis steps to be run specific to each read group, thus representing a significant computational challenge. While sequencing technologies for producing these data are approaching maturity, the lack of available computational tools for analysing such heterogeneous sequence data represents an obstacle to the widespread adoption of this technology. Results. Using synthetic data we successfully detect unique variants at dilution levels of 1 in a 1,000,000 molecules, and find DeeepSNVMiner obtains significantly lower false positive and false negative rates compared to popular variant callers GATK, SAMTools, FreeBayes and LoFreq, particularly as the variant concentration levels decrease. In a dilution series with genomic DNA from two cells lines, we find DeepSNVMiner identifies a known somatic variant when present at concentrations of only 1 in 1,000 molecules in the input material, the lowest concentration amongst all variant callers tested. Conclusions. Here we present DeepSNVMiner; a tool to disambiguate tagged sequence groups and robustly identify sequence variants specific to subsets of starting DNA molecules that may indicate the presence of a disease. DeepSNVMiner is an automated workflow of custom sequence analysis utilities and open source tools able to differentiate somatic DNA variants from artefactual sequence variants that likely arose during DNA amplification. The workflow remains flexible such that it may be customised to variants of the data production protocol used, and supports reproducible analysis through detailed logging and reporting of results. DeepSNVMiner is available for academic non-commercial research purposes at https://github.com/mattmattmattmatt/DeepSNVMiner.

Whole genome sequencing (WGS) of cancer patients' tumours offers the most comprehensive method of identifying both novel and known clinically-actionable genomic targets. However, the practicalities of performing WGS on clinical samples are poorly defined.This study was designed to test sample preparation, sequencing specifications and bioinformatic algorithms for their effect on accuracy and cost-efficiency in a large WGS analysis of human melanoma samples.WGS was performed on melanoma cell lines (n = 15) and melanoma fresh frozen tumours (n = 222). The appropriate level of coverage and the optimal mutation detection algorithm for the project pipeline were determined.An incremental increase in sequencing coverage from 36X to 132X in melanoma tissue samples and 30X to 103X for cell lines only resulted in a small increase (1-2%) in the number of mutations detected, and the quality scores of the additional mutations indicated a low probability that the mutations were real. The results suggest that 60X coverage for melanoma tissue and 40X for melanoma cell lines empower the detection of 98-99% of informative single nucleotide variants (SNVs), a sensitivity level at which clinical decision making or landscape research projects can be carried out with a high degree of confidence in the results. Likewise the bioinformatic mutation analysis methodology strongly influenced the number and quality of SNVs detected. Detecting mutations in the blood genomes separate to the tumour genomes generated 41% more SNVs than if the blood and melanoma tissue genomes were analysed simultaneously. Therefore, simultaneous analysis should be employed on matched melanoma tissue and blood genomes to reduce errors in mutation detection.This study provided valuable insights into the accuracy of SNV with WGS at various coverage levels in human clinical cancer specimens. Additionally, we investigated the accuracy of the publicly available mutation detection algorithms to detect cancer specific SNVs which will aid researchers and clinicians in study design and implementation of WGS for the identification of somatic mutations in other cancers.

A diversity of tools is available for identification of variants from genome sequence data. Given the current complexity of incorporating external software into a genome analysis infrastructure, a tendency exists to rely on the results from a single tool alone. The quality of the output variant calls is highly variable however, depending on factors such as sequence library quality as well as the choice of short-read aligner, variant caller, and variant caller filtering strategy. Here we present a two-part study first using the high quality 'genome in a bottle' reference set to demonstrate the significant impact the choice of aligner, variant caller, and variant caller filtering strategy has on overall variant call quality and further how certain variant callers outperform others with increased sample contamination, an important consideration when analyzing sequenced cancer samples. This analysis confirms previous work showing that combining variant calls of multiple tools results in the best quality resultant variant set, for either specificity or sensitivity, depending on whether the intersection or union, of all variant calls is used respectively. Second, we analyze a melanoma cell line derived from a control lymphocyte sample to determine whether software choices affect the detection of clinically important melanoma risk-factor variants finding that only one of the three such variants is unanimously detected under all conditions. Finally, we describe a cogent strategy for implementing a clinical variant detection pipeline; a strategy that requires careful software selection, variant caller filtering optimizing, and combined variant calls in order to effectively minimize false negative variants. While implementing such features represents an increase in complexity and computation the results offer indisputable improvements in data quality.