Search

Spontaneous preterm birth (sPTB) is a leading cause of maternal and neonatal morbidity and mortality, yet its prevention and early risk stratification are limited. Previous investigations have suggested that vaginal microbes and metabolites may be implicated in sPTB. Here we performed untargeted metabolomics on 232 second-trimester vaginal samples, 80 from pregnancies ending preterm. We find multiple associations between vaginal metabolites and subsequent preterm birth, and propose that several of these metabolites, including diethanolamine and ethyl glucoside, are exogenous. We observe associations between the metabolome and microbiome profiles previously obtained using 16S ribosomal RNA amplicon sequencing, including correlations between bacteria considered suboptimal, such as Gardnerella vaginalis, and metabolites enriched in term pregnancies, such as tyramine. We investigate these associations using metabolic models. We use machine learning models to predict sPTB risk from metabolite levels, weeks to months before birth, with good accuracy (area under receiver operating characteristic curve of 0.78). These models, which we validate using two external cohorts, are more accurate than microbiome-based and maternal covariates-based models (area under receiver operating characteristic curve of 0.55–0.59). Our results demonstrate the potential of vaginal metabolites as early biomarkers of sPTB and highlight exogenous exposures as potential risk factors for prematurity.

Cancers assume a variety of distinct histologies and may originate from a myriad of sites including solid organs, hematopoietic cells, and connective tissue. Clinical decision making based on consensus guidelines such as NCCN is often predicated on a specific histologic and anatomic diagnosis, supported by clinical features and pathologist interpretation of morphology and immunohistochemical (IHC) staining patterns. However, in patients with nonspecific morphologic and IHC findings—in addition to ambiguous clinical presentations such as recurrence versus new primary—a definitive diagnosis may not be possible, resulting in the patient being categorized as having a cancer of unknown primary (CUP). Therapeutic options and clinical outcomes are poor for CUP patients, with a median survival of 8-11 months. Here we describe and validate the Tempus Tumor Origin (Tempus TO) assay, an RNA-seq-based machine learning classifier capable of discriminating between 68 clinically relevant cancer subtypes. We show that the Tempus TO model is 91% accurate when assessed on retrospectively and prospectively held out cohorts of containing 9,210 samples with known diagnoses. When evaluated on a cohort of CUPs, the model recapitulated established associations between genomic alterations and cancer subtype. Combining diagnostic prediction tests (e.g., Tempus TO) with sequencing-based variant reporting (e.g., Tempus xT) may expand therapeutic options for patients with cancers of unknown primary or uncertain histology.

The oral–gut microbiota axis in health is a fundamentally important and clinically relevant topic; however, our current understanding of it remains biased and incomplete. By characterizing the bacterial and fungal microbiomes in paired oral rinse and stool samples from a large cohort of healthy Chinese adults, here we provided new evidence that oral–gut microbiota transmission is limited in non-Western population and across biological domains.

Species and subspecies within the Salmonella genus have been defined for public health purposes by biochemical properties; however, reference laboratories have increasingly adopted sequence-based, and especially whole genome sequence (WGS), methods for surveillance and routine identification. This leads to potential disparities in subspecies definitions, routine typing, and the ability to detect novel subspecies. A large-scale analysis of WGS data from the routine sequencing of clinical isolates was employed to define and characterise Salmonella subspecies population structure, demonstrating that the Salmonella species and subspecies were genetically distinct, including those previously identified through phylogenetic approaches, namely: S. enterica subspecies londinensis (VII), subspecies brasiliensis (VIII), subspecies hibernicus (IX) and subspecies essexiensis (X). The analysis also identified an additional novel subspecies, reptilium (XI). Further, these analyses indicated that S. enterica subspecies arizonae (IIIa) isolates were divergent from the other S. enterica subspecies, which clustered together and, on the basis of ANI analysis, subspecies IIIa was sufficiently distinct to be classified as a separate species, S. arizonae. Multiple phylogenetic and statistical approaches generated congruent results, suggesting that the proposed species and subspecies structure was sufficiently biologically robust for routine application. Biochemical analyses demonstrated that not all subspecies were distinguishable by these means and that biochemical approaches did not capture the genomic diversity of the genus. We recommend the adoption of standardised genomic definitions of species and subspecies and a genome sequence-based approach to routine typing for the identification and definition of novel subspecies., Highlights • A large-scale analysis of genomic data demonstrate Salmonella species and subspecies are genetically distinct. • Biochemical analysis does not capture the genomic diversity of the Salmonella genus but routine species and subspecies identification can be achieved with rMLST • Average Nucleotide Identify (ANI) with a 95% criteria was suitable to distinguish species and 98% to distinguish subspecies. • Five novel S. enteric subspecies (VII-XI) type strains are defined. • Reclassification of S. arizonae as a separate species is recommended.

Genetic diversity is a fundamental component of biological diversity, and its conservation is considered key to ensure the long-term survival of natural populations and species. National and international legislation increasingly mandates a monitoring of genetic diversity. Examples are the United Nation's Convention on Biological Diversity (CBD) Aichi target 13 and the current post-2020 negotiations to specify a new target for maintaining genetic diversity. To date, only a few pilot projects have been launched that systematically monitor genetic diversity over time in natural populations of a broad variety of wild species. The Swiss Federal Office for the Environment mandated a feasibility study in 2019 for implementing a national monitoring of genetic diversity in natural populations. To obtain information on whether stakeholders are interested in such a systematic monitoring, what they would expect from such a monitoring and where they see respective caveats, we conducted an online survey, which 138 (42% of those surveyed) Swiss stakeholders answered. We find that Swiss stakeholders are generally aware of the lacking evidence regarding the status of genetic diversity in wild populations and species. Accordingly, most stakeholders are interested in a monitoring of genetic diversity and see opportunities for the application of its results in their work. Nevertheless, stakeholders also expressed concerns regarding financial resources and that the results of a genetic diversity monitoring program would not benefit conservation practice. Our findings highlight the importance of stakeholder engagement and demonstrate the value of a detailed stakeholder analysis prior to developing and implementing a genetic diversity monitoring program. A powerful tool for examining the constellation and interactions of the different stakeholders are social network analyses (SNAs). Finally, it is particularly important to communicate transparently about the possibilities and limitations of a genetic diversity monitoring program as well as to closely involve stakeholders from the beginning to increase the acceptance of genetic diversity monitoring and facilitate its implementation., Conservation Genetics, 22 (5), ISSN:1572-9737, ISSN:1566-0621

Most genes in photosynthetic organisms remain functionally uncharacterized. Here, using a barcoded mutant library of the model eukaryotic alga Chlamydomonas reinhardtii, we determined the phenotypes of more than 58,000 mutants under more than 121 different environmental growth conditions and chemical treatments. A total of 59% of genes are represented by at least one mutant that showed a phenotype, providing clues to the functions of thousands of genes. Mutant phenotypic profiles place uncharacterized genes into functional pathways such as DNA repair, photosynthesis, the CO2-concentrating mechanism and ciliogenesis. We illustrate the value of this resource by validating phenotypes and gene functions, including three new components of an actin cytoskeleton defense pathway. The data also inform phenotype discovery in land plants; mutants in Arabidopsis thaliana genes exhibit phenotypes similar to those we observed in their Chlamydomonas homologs. We anticipate that this resource will guide the functional characterization of genes across the tree of life.

Minimal absent words (MAWs) are minimal-length oligomers absent from a genome or proteome. Although some artificially synthesized MAWs have deleterious effects, there is still a lack of a strategy for the classification of non-occurring sequences as potentially malicious or benign. In this work, by using Markovian models with multiple-testing correction, we reveal significant absent oligomers, which are statistically expected to exist. This suggests that their absence is due to negative selection. We survey genomes and proteomes covering the diversity of life and find thousands of significant absent sequences. Common significant MAWs are often mono- or dinucleotide tracts, or palindromic. Significant viral MAWs are often restriction sites and may indicate unknown restriction motifs. Surprisingly, significant mammal genome MAWs are often present, but rare, in other mammals, suggesting that they are suppressed but not completely forbidden. Significant human MAWs are frequently present in prokaryotes, suggesting immune function, but rarely present in human viruses, indicating viral mimicry of the host. More than one-fourth of human proteins are one substitution away from containing a significant MAW, with the majority of replacements being predicted harmful. We provide a web-based, interactive database of significant MAWs across genomes and proteomes.

Numerous studies have examined the composition of and factors shaping the oral bacterial microbiota in healthy adults; however, similar studies on the less dominant yet ecologically and clinically important fungal microbiota are scarce. In this study, we characterized simultaneously the oral bacterial and fungal microbiomes in a large cohort of systemically healthy Chinese adults by sequencing the bacterial 16S rRNA gene and fungal internal transcribed spacer. We showed that different factors shaped the oral bacterial and fungal microbiomes in healthy adults. Sex and age were associated with the alpha diversity of the healthy oral bacterial microbiome but not that of the fungal microbiome. Age was also a major factor affecting the beta diversity of the oral bacterial microbiome; however, it only exerted a small effect on the oral fungal microbiome when compared with other variables. After controlling for age and sex, the bacterial microbiota structure was most affected by marital status, recent oral conditions and oral hygiene-related factors, whereas the fungal microbiota structure was most affected by education level, fruits and vegetables, and bleeding gums. Bacterial-fungal interactions were limited in the healthy oral microbiota, with the strongest association existing between Pseudomonas sp. and

Opsins, the protein moieties of animal visual photo-pigments, have emerged as moonlighting proteins with diverse, light-dependent and -independent physiological functions. This raises the need to revise some basic assumptions concerning opsin expression, structure, classification, and evolution.

Whole-exome sequencing (WES) can detect not only single-nucleotide variants in causal genes, but also pathogenic copy-number variations using several methods. However, there may be overlooked pathogenic variations in the out of target genome regions of WES analysis (e.g., promoters), leaving many patients undiagnosed. Whole-genome sequencing (WGS) can potentially analyze such regions. We applied long-read nanopore WGS and our recently developed analysis pipeline "dnarrange" to a patient who was undiagnosed by trio-based WES analysis, and identified a heterozygous 97-kb deletion partially involving 5'-untranslated exons of MBD5, which was outside the WES target regions. The phenotype of the patient, a 32-year-old male, was consistent with haploinsufficiency of MBD5. The transcript level of MBD5 in the patient's lymphoblastoid cells was reduced. We therefore concluded that the partial MBD5 deletion is the culprit for this patient. Furthermore, we found other rare structural variations (SVs) in this patient, i.e., a large inversion and a retrotransposon insertion, which were not seen in 33 controls. Although we considered that they are benign SVs, this finding suggests that our pipeline using long-read WGS is useful for investigating various types of potentially pathogenic SVs. In conclusion, we identified a 97-kb deletion, which causes haploinsufficiency of MBD5 in a patient with neurodevelopmental disorder, demonstrating that long-read WGS is a powerful technique to discover pathogenic SVs.

Yeast and cancer cells are metabolically similar as they use fermentation of glucose as a primary means of generating energy. Reliance on glucose fermentation makes both of these cell types highly sensitive to the toxic glucose analog, 2-deoxyglucose. Here we review the cellular and metabolic pathways that play a role in 2-deoxyglucose sensitivity and discuss how the modifications to these pathways result in acquisition of 2-deoxyglucose resistance. Insights gained from genetic and proteomic studies in yeast provide new ideas for the design of combinatorial therapies for cancer treatment.

Antimicrobial resistance in Neisseria gonorrhoeae is a global health concern. Strains from two internationally circulating sequence types, ST-7363 and ST-1901, have acquired resistance to treatment with third-generation cephalosporins mainly due to the emergence of mosaic penA alleles. These two STs were first detected in Japan; however, when and how the mosaic penA alleles emerged and spread to other countries remains unknown. Here, we addressed the evolution of penA alleles by obtaining complete genomes from three Japanese ST-1901 clinical isolates harboring mosaic penA allele 34 (penA-34) dating from 2005 and generating a phylogenetic representation of 1,075 strains sampled from 37 countries. We also sequenced the genomes of 103 Japanese ST-7363 N. gonorrhoeae isolates from 1996-2005 and reconstructed a phylogeny including 88 previously sequenced genomes. Based on an estimate of the time of emergence of ST-1901 harboring mosaic penA-34 and ST-7363 harboring mosaic penA-10, and >300 additional genome sequences of Japanese strains representing multiple STs isolated in 1996-2015, we suggest that penA-34 in ST-1901 was generated from penA-10 via recombination with another Neisseria species, followed by a second recombination event with a gonococcal strain harboring wildtype penA-1. Following the acquisition of penA-10 in ST-7363, a dominant sub-lineage rapidly acquired fluoroquinolone resistance mutations at GyrA 95 and ParC 87-88, possibly due to independent mutations rather than horizontal gene transfer. Literature data suggest the emergence of these resistance determinants may reflect selection from the standard treatment regimens in Japan at that time. Our findings highlight how recombination and antibiotic use across and within Neisseria species intersect in driving the emergence and spread of drug-resistant gonorrhea.Author summaryAntimicrobial resistance is recognized as one of the greatest threats to human health, and Neisseria gonorrhoeae resistance is classified as one of the most urgent. The two major internationally spreading lineages resistant. to first line drugs likely originated in Japan, but when and how their genetic resistance determinants emerged remain unknown. In this study, we conducted an evolutionary analysis using clinical N. gonorrhoeae isolates from 37 countries, including a historical collection of Japanese isolates, to investigate the emergence of resistance in each of the two major lineages. We showed that the penA allele responsible for resistance to cephalosporins, the first-line treatment for gonorrhea, was possibly generated by two recombination events, one from another Neisseria species and one from another N. gonorrhoeae lineage. We also showed that mutations responsible for resistance to a previously widely used antibiotic treatment occurred twice independently in one of the two major lineages. The emergence of the genetic resistance determinants potentially reflects selection from the standard treatment regimen at that time. Our findings highlight how recombination (horizontal gene transfer) and antibiotic use across and within a bacterial species intersect in driving the emergence and spread of antimicrobial resistance genes and mutations.

BACKGROUND We previously determined that polyplexes formed by linear H2K peptides were more effective in transfecting tumors in vivo than polyplexes formed by branched H2K4b-20 peptides. Based on trypsin digest and salt displacement studies, the linear H2K polyplexes were less stable than the branched H2K4b-20 polyplexes. Because binding and release of the polymer and DNA from the H2K4b-20 polyplex may account for the ineffectiveness, we investigated whether four-branched histidine-lysine (HK) peptides with varying numbers of amino acids in their branches would be more effective in their ability to increase gene expression in tumors in vivo. METHODS Linear and branched peptides with multiple -KHHK- motifs were synthesized by solid-phase synthesis. The branched H2K4b-20, -18, -14 and 12 peptides had 20, 18, 14 and 12 amino acids in their branches, respectively. These peptides were examined for their ability to carry luciferase-expressing plasmids to human breast cancer xenografts in a mouse model. With gel retardation and in vivo transfection, the incorporation of a targeting ligand and an endosomal lysis peptide into these polyplexes was also examined. A blocking antibody was pre-injected prior to the polyplexes to determine the role of neuropilin 1 in the uptake of these polyplexes by the tumor. The size of the polyplexes was measured by dynamic light scattering. RESULTS Of the four negative surface-charge polyplexes formed by the branched carriers, the H2K4b-14 polyplex was determined to be the most effective plasmid delivery platform to tumors. The incorporation of a targeting ligand and an endosomal lysis peptide into H2K4b-14 polyplexes further enhanced their ability to transfect tumors in vivo. Furthermore, after pre-injecting tumor-bearing mice with a blocking antibody to the neuropilin-1 receptor (NRP-1), there was a marked reduction of tumor gene expression with the modified H2K4b-14 polyplexes, suggesting that NRP-1 mediated their transport into the tumor. CONCLUSIONS The present study established that branched peptides intermediate in length were very efficient in delivering plasmids to tumors in vivo.

Background Many genetic/genomic disorders are caused by genomic rearrangements. Standard methods can often characterize these variations only partly, e.g., copy number changes or breakpoints. It is important to fully understand the order and orientation of rearranged fragments, with precise breakpoints, to know the pathogenicity of the rearrangements. Methods We performed whole-genome-coverage nanopore sequencing of long DNA reads from four patients with chromosomal translocations. We identified rearrangements relative to a reference human genome, subtracted rearrangements shared by any of 33 control individuals, and determined the order and orientation of rearranged fragments, with our newly developed analysis pipeline. Results We describe the full characterization of complex chromosomal rearrangements, by filtering out genomic rearrangements seen in controls without the same disease, reducing the number of loci per patient from a few thousand to a few dozen. Breakpoint detection was very accurate; we usually see ~ 0 ± 1 base difference from Sanger sequencing-confirmed breakpoints. For one patient with two reciprocal chromosomal translocations, we find that the translocation points have complex rearrangements of multiple DNA fragments involving 5 chromosomes, which we could order and orient by an automatic algorithm, thereby fully reconstructing the rearrangement. A rearrangement is more than the sum of its parts: some properties, such as sequence loss, can be inferred only after reconstructing the whole rearrangement. In this patient, the rearrangements were evidently caused by shattering of the chromosomes into multiple fragments, which rejoined in a different order and orientation with loss of some fragments. Conclusions We developed an effective analytic pipeline to find chromosomal aberration in congenital diseases by filtering benign changes, only from long read sequencing. Our algorithm for reconstruction of complex rearrangements is useful to interpret rearrangements with many breakpoints, e.g., chromothripsis. Our approach promises to fully characterize many congenital germline rearrangements, provided they do not involve poorly understood loci such as centromeric repeats.

Background Gonorrhea, caused by the bacterium Neisseria gonorrhoeae, is a globally prevalent sexually transmitted infection. The dynamics of gonococcal population biology have been poorly defined due to a lack of resolution in strain typing methods. Methods In this study, we assess how the core genome can be used to improve our understanding of gonococcal population structure compared with current typing schemes. Results A total of 1668 loci were identified as core to the gonococcal genome. These were organized into a core genome multilocus sequence typing scheme (N gonorrhoeae cgMLST v1.0). A clustering algorithm using a threshold of 400 allelic differences between isolates resolved gonococci into discrete and stable core genome groups, some of which persisted for multiple decades. These groups were associated with antimicrobial genotypes and non-overlapping NG-STAR and NG-MAST sequence types. The MLST-STs were more widely distributed among core genome groups. Conclusions Clustering with cgMLST identified globally distributed, persistent, gonococcal lineages improving understanding of the population biology of gonococci and revealing its population structure. These findings have implications for the emergence of antimicrobial resistance in gonococci and how this is associated with lineages, some of which are more predisposed to developing antimicrobial resistance than others.

Recently, a recessively inherited intronic repeat expansion in replication factor C1 (RFC1) was identified in cerebellar ataxia with neuropathy and bilateral vestibular areflexia syndrome (CANVAS). Here, we describe a Japanese case of genetically confirmed CANVAS with autonomic failure and auditory hallucination. The case showed impaired uptake of iodine-123-metaiodobenzylguanidine and 123I-ioflupane in the cardiac sympathetic nerve and dopaminergic neurons, respectively, by single-photon emission computed tomography. Long-read sequencing identified biallelic pathogenic (AAGGG)n nucleotide repeat expansion in RFC1 and heterozygous benign (TAAAA)n and (TAGAA)n expansions in brain expressed, associated with NEDD4 (BEAN1). Enrichment of the repeat regions in RFC1 and BEAN1 using a Cas9-mediated system clearly distinguished between pathogenic and benign repeat expansions. The haplotype around RFC1 indicated that the (AAGGG)n expansion in our case was on the same ancestral allele as that of European cases. Thus, long-read sequencing facilitates precise genetic diagnosis of diseases with complex repeat structures and various expansions.

Invasive meningococcal disease (IMD) due to serogroup Y Neisseria meningitidis emerged in Europe during the 2000s. Draft genomes of serogroup Y isolates in Sweden revealed that although the population structure of these isolates was similar to other serogroup Y isolates internationally, a distinct strain (YI) and more specifically a sublineage (1) of this strain was responsible for the increase of serogroup Y IMD in Sweden. We performed single molecule real-time (SMRT) sequencing on eight serogroup Y isolates from different sublineages to unravel the genetic and epigenetic factors delineating them, in order to understand the serogroup Y emergence. Extensive comparisons between the serogroup Y sublineages of all coding sequences, complex genomic regions, intergenic regions, and methylation motifs revealed small point mutations in genes mainly encoding hypothetical and metabolic proteins, and non-synonymous variants in genes involved in adhesion, iron acquisition, and endotoxin production. The methylation motif CACNNNNNTAC was only found in isolates of sublineage 2. Only seven genes were putatively differentially expressed, and another two genes encoding hypothetical proteins were only present in sublineage 2. These data suggest that the serogroup Y IMD increase in Sweden was most probably due to small changes in genes important for colonization and transmission.

Many invasive bacterial diseases are caused by organisms that are ordinarily harmless components of the human microbiome. Effective interventions against these microbes require an understanding of the processes whereby symbiotic or commensal relationships transition into pathology. Here, we describe bacterial genome-wide association studies (GWAS) of Neisseria meningitidis, a common commensal of the human respiratory tract that is nevertheless a leading cause of meningitis and sepsis. An initial GWAS discovered bacterial genetic variants, including single nucleotide polymorphisms (SNPs), associated with invasive meningococcal disease (IMD) versus carriage in several loci across the meningococcal genome, encoding antigens and other extracellular components, confirming the polygenic nature of the invasive phenotype. In particular, there was a significant peak of association around the fHbp locus, encoding factor H binding protein (fHbp), which promotes bacterial immune evasion of human complement by recruiting complement factor H (CFH) to the meningococcal surface. The association around fHbp with IMD was confirmed by a validation GWAS, and we found that the SNPs identified in the validation affected the 5’ region of fHbp mRNA, altering secondary RNA structures, thereby increasing fHbp expression and enhancing bacterial escape from complement-mediated killing. This finding is consistent with the known link between complement deficiencies and CFH variation with human susceptibility to IMD. These observations demonstrate the importance of human and bacterial genetic variation across the fHbp:CFH interface in determining IMD susceptibility, the transition from carriage to disease., Author summary The human bacterial pathogen Neisseria meningitidis (the meningococcus) causes sepsis and meningitis worldwide. Paradoxically, it is carried harmlessly in the nose and throat far more commonly than it causes invasive disease, raising the question of whether the bacteria causing disease are genetically the same as those that do not. We looked for systematic differences in the DNA of the bacteria from invasive disease and those that were carried harmlessly. Whilst controlling for population structure as a confounder, we identified multiple genes apparently contributing to invasion and identified a signal around the gene encoding an important component of some meningococcal vaccines, called factor H binding protein (fHbp). This protein helps the meningococcus evade killing by the human immune response by binding to a human defence protein called complement factor H. We validated this result in an independent population of meningococci and identified that DNA variants that control the amount of fHbp produced by the bacteria were important, with high levels of fHbp expression associated with invasion. Given that human DNA variation in complement factor H increases the risk of meningococcal invasion, this highlights the important connection between human and bacterial genetic variation and helps explain why some people infected with meningococci suffer meningococcal disease whilst others do not.

Glucose is central to many biological processes, serving as an energy source and a building block for biosynthesis. After glucose enters the cell, hexokinases convert it to glucose-6-phosphate (Glc-6P) for use in anaerobic fermentation, aerobic oxidative phosphorylation, and the pentose-phosphate pathway. We here describe a genetic screen in Saccharomyces cerevisiae that generated a novel spontaneous mutation in hexokinase-2, hxk2G238V, that confers resistance to the toxic glucose analog 2-deoxyglucose (2DG). Wild-type hexokinases convert 2DG to 2-deoxyglucose-6-phosphate (2DG-6P), but 2DG-6P cannot support downstream glycolysis, resulting in a cellular starvation-like response. Curiously, though the hxk2G238V mutation encodes a loss-of-function allele, the affected amino acid does not interact directly with bound glucose, 2DG, or ATP. Molecular dynamics simulations suggest that Hxk2G238V impedes sugar binding by altering the protein dynamics of the glucose-binding cleft, as well as the large-scale domain-closure motions required for catalysis. These findings shed new light on Hxk2 dynamics and highlight how allosteric changes can influence catalysis, providing new structural insights into this critical regulator of carbohydrate metabolism. Given that hexokinases are upregulated in some cancers and that 2DG and its derivatives have been studied in anti-cancer trials, the present work also provides insights that may apply to cancer biology and drug resistance.

Digestive Chagas disease (DCD) is an enteric neuropathy caused by Trypanosoma cruzi infection. The mechanism of pathogenesis is poorly understood and the lack of a robust, predictive animal model has held back research. We screened a series of mouse models using gastrointestinal tracer assays and in vivo infection imaging systems to discover a subset exhibiting chronic digestive transit dysfunction and significant retention of faeces in both sated and fasted conditions. The colon was a specific site of both tissue parasite persistence, delayed transit and dramatic loss of myenteric neurons as revealed by whole-mount immunofluorescence analysis. DCD mice therefore recapitulated key clinical manifestations of human disease. We also exploited dual reporter transgenic parasites to home in on locations of rare chronic infection foci in the colon by ex vivo bioluminescence imaging and then used fluorescence imaging in tissue microdomains to reveal co-localisation of infection and enteric nervous system lesions. This indicates that long-term T. cruzi-host interactions in the colon drive DCD pathogenesis, suggesting that the efficacy of anti-parasitic chemotherapy against chronic disease progression warrants further pre-clinical investigation., Author summary Chagas disease (American trypanosomiasis) is caused by the protozoan parasite Trypanosoma cruzi. Chagas disease has two types, the cardiac form and the digestive form; some patients have symptoms of both. How the parasite causes digestive disease is poorly understood. It is known that damage to the gut’s nervous system is an important factor, but it has been unclear exactly where and when this damage occurs during the course of an infection and also why only a subset of infected people suffer from this outcome. We studied infections in mice and found certain combinations of strains of parasites and mice that exhibited symptoms similar to human digestive Chagas patients, including a problem with peristalsis that localised specifically to the colon. Using parasites that were genetically engineered to emit both bioluminescent and fluorescent light, we tracked infections over time and were able to analyse rare infected cells deep within the muscle tissue of the wall of the colon. We found evidence of damaged neurons in the same location as these infection foci over 6 months after initial infection. Our results show that digestive Chagas disease probably develops as a result of chronic infection and inflammation, which potentially changes approaches to treatment.

Pseudogenes have been classified as functionless and their annotation is an ongoing problem. The Adh6-ps1—a mouse pseudogene belonging to the alcohol dehydrogenase gene complex (Adh) was analyzed to review the conservation, homology, expression and interactions, and identify any role it plays in disease phenotypes using bioinformatics databases. Results showed that Adh6-ps1 have 2 transcripts (processed and unprocessed) which may have emerged from a transposition and duplication event, respectively, and that induced inversions (Uox gene, In(3)11Rk) involving gene complexes associated with Adh6-ps1 have been implicated in a diverse range of diseases. Adh6-ps1 is highly conserved in vertebrates particularly rodents, and expressed in the liver. The top 5 MirRNA targets identified are Mir455, Mir511, Mir1903, Mir361, and Mir669o markers. While much is unknown about Mir1903 and Mir669o, the silencing of Mir455 and Mir511 is linked with hepatocellular carcinoma (HCC), and Mir361 is implicated in endometrial cancers. Given the identified MirRNA interactions with Adh6-ps1 and their expression in cancer phenotypes, and Adh6-ps1 associations with induced inversions, it may well have a role in tumorigenesis and disease phenotypes. Nonetheless, further studies are required to establish these facts to add to the growing efforts to understand pseudogenes and their potential involvement in disease conditions.

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole‐genome pooled population sequencing to study genome‐wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over‐representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species‐specific.

Most genes associated with neurodevelopmental disorders (NDDs) were identified with an excess of de novo mutations (DNMs) but the significance in case–control mutation burden analysis is unestablished. Here, we sequence 63 genes in 16,294 NDD cases and an additional 62 genes in 6,211 NDD cases. By combining these with published data, we assess a total of 125 genes in over 16,000 NDD cases and compare the mutation burden to nonpsychiatric controls from ExAC. We identify 48 genes (25 newly reported) showing significant burden of ultra-rare (MAF, For many neurodevelopmental disorder (NDD) risk genes, the significance for mutational burden is unestablished. Here, the authors sequence 125 candidate NDD genes in over 16,000 NDD cases; case-control mutational burden analysis identifies 48 genes with a significant burden of severe ultra-rare mutations.

This human campylobacteriosis sentinel surveillance project was based at two sites in Oxfordshire and North East England chosen (i) to be representative of the English population on the Office for National Statistics urban-rural classification and (ii) to provide continuity with genetic surveillance started in Oxfordshire in October 2003. Between October 2015 and September 2018 epidemiological questionnaires and genome sequencing of isolates from human cases was accompanied by sampling and genome sequencing of isolates from possible food animal sources. The principal aim was to estimate the contributions of the main sources of human infection and to identify any changes over time. An extension to the project focussed on antimicrobial resistance in study isolates and older archived isolates. These older isolates were from earlier years at the Oxfordshire site and the earliest available coherent set of isolates from the national archive at Public Health England (1997/8). The aim of this additional work was to analyse the emergence of the antimicrobial resistance that is now present among human isolates and to describe and compare antimicrobial resistance in recent food animal isolates. Having identified the presence of bias in population genetic attribution, and that this was not addressed in the published literature, this study developed an approach to adjust for bias in population genetic attribution, and an alternative approach to attribution using sentinel types. Using these approaches the study estimated that approximately 70% of Campylobacter jejuni and just under 50% of C. coli infection in our sample was linked to the chicken source and that this was relatively stable over time. Ruminants were identified as the second most common source for C. jejuni and the most common for C. coli where there was also some evidence for pig as a source although less common than ruminant or chicken. These genomic attributions of themselves make no inference on routes of transmission. However, those infected with isolates genetically typical of chicken origin were substantially more likely to have eaten chicken than those infected with ruminant types. Consumption of lamb’s liver was very strongly associated with infection by a strain genetically typical of a ruminant source. These findings support consumption of these foods as being important in the transmission of these infections and highlight a potentially important role for lamb’s liver consumption as a source of Campylobacter infection. Antimicrobial resistance was predicted from genomic data using a pipeline validated by Public Health England and using BIGSdb software. In C. jejuni this showed a nine-fold increase in resistance to fluoroquinolones from 1997 to 2018. Tetracycline resistance was also common, with higher initial resistance (1997) and less substantial change over time. Resistance to aminoglycosides or macrolides remained low in human cases across all time periods. Among C. jejuni food animal isolates, fluoroquinolone resistance was common among isolates from chicken and substantially less common among ruminants, ducks or pigs. Tetracycline resistance was common across chicken, duck and pig but lower among ruminant origin isolates. In C. coli resistance to all four antimicrobial classes rose from low levels in 1997. The fluoroquinolone rise appears to have levelled off earlier and among animals, levels are high in duck as well as chicken isolates, although based on small sample sizes, macrolide and aminoglycoside resistance, was substantially higher than for C. jejuni among humans and highest among pig origin isolates. Tetracycline resistance is high in isolates from pigs and the very small sample from ducks. Antibiotic use following diagnosis was relatively high (43.4%) among respondents in the human surveillance study. Moreover, it varied substantially across sites and was highest among non-elderly adults compared to older adults or children suggesting opportunities for improved antimicrobial stewardship. The study also found evidence for stable lineages over time across human and source animal species as well as some tighter genomic clusters that may represent outbreaks. The genomic dataset will allow extensive further work beyond the specific goals of the study. This has been made accessible on the web, with access supported by data visualisation tools.

Many copy number variants (CNVs) confer risk for the same range of neurodevelopmental symptoms and psychiatric conditions including autism and schizophrenia. Yet, to date neuroimaging studies have typically been carried out one mutation at a time, showing that CNVs have large effects on brain anatomy. Here, we aimed to characterize and quantify the distinct brain morphometry effects and latent dimensions across 8 neuropsychiatric CNVs. We analyzed T1-weighted MRI data from clinically and non-clinically ascertained CNV carriers (deletion/duplication) at the 1q21.1 (n = 39/28), 16p11.2 (n = 87/78), 22q11.2 (n = 75/30), and 15q11.2 (n = 72/76) loci as well as 1296 non-carriers (controls). Case-control contrasts of all examined genomic loci demonstrated effects on brain anatomy, with deletions and duplications showing mirror effects at the global and regional levels. Although CNVs mainly showed distinct brain patterns, principal component analysis (PCA) loaded subsets of CNVs on two latent brain dimensions, which explained 32 and 29% of the variance of the 8 Cohen’s d maps. The cingulate gyrus, insula, supplementary motor cortex, and cerebellum were identified by PCA and multi-view pattern learning as top regions contributing to latent dimension shared across subsets of CNVs. The large proportion of distinct CNV effects on brain morphology may explain the small neuroimaging effect sizes reported in polygenic psychiatric conditions. Nevertheless, latent gene brain morphology dimensions will help subgroup the rapidly expanding landscape of neuropsychiatric variants and dissect the heterogeneity of idiopathic conditions.

The recent Chandos House meeting of the Alport Variant Collaborative extended the indications for screening for pathogenic variants in the COL4A5, COL4A3 and COL4A4 genes beyond the classical Alport phenotype (haematuria, renal failure; family history of haematuria or renal failure) to include persistent proteinuria, steroid-resistant nephrotic syndrome, focal and segmental glomerulosclerosis (FSGS), familial IgA glomerulonephritis and end-stage kidney failure without an obvious cause. The meeting refined the ACMG criteria for variant assessment for the Alport genes (COL4A3–5). It identified ‘mutational hotspots’ (PM1) in the collagen IV α5, α3 and α4 chains including position 1 Glycine residues in the Gly-X-Y repeats in the intermediate collagenous domains; and Cysteine residues in the carboxy non-collagenous domain (PP3). It considered that ‘well-established’ functional assays (PS3, BS3) were still mainly research tools but sequencing and minigene assays were commonly used to confirm splicing variants. It was not possible to define the Minor Allele Frequency (MAF) threshold above which variants were considered Benign (BA1, BS1), because of the different modes of inheritances of Alport syndrome, and the occurrence of hypomorphic variants (often Glycine adjacent to a non-collagenous interruption) and local founder effects. Heterozygous COL4A3 and COL4A4 variants were common ‘incidental’ findings also present in normal reference databases. The recognition and interpretation of hypomorphic variants in the COL4A3–COL4A5 genes remains a challenge.

Viral diseases pose major threats to humans and other animals, including the billions of chickens that are an important food source as well as a public health concern due to zoonotic pathogens. Unlike humans and other typical mammals, the major histocompatibility complex (MHC) of chickens can confer decisive resistance or susceptibility to many viral diseases. An iconic example is Marek’s disease, caused by an oncogenic herpesvirus with over 100 genes. Classical MHC class I and class II molecules present antigenic peptides to T lymphocytes, and it has been hard to understand how such MHC molecules could be involved in susceptibility to Marek’s disease, given the potential number of peptides from over 100 genes. We used a new in vitro infection system and immunopeptidomics to determine peptide motifs for the 2 class II molecules expressed by the MHC haplotype B2, which is known to confer resistance to Marek’s disease. Surprisingly, we found that the vast majority of viral peptide epitopes presented by chicken class II molecules arise from only 4 viral genes, nearly all having the peptide motif for BL2*02, the dominantly expressed class II molecule in chickens. We expressed BL2*02 linked to several Marek’s disease virus (MDV) peptides and determined one X-ray crystal structure, showing how a single small amino acid in the binding site causes a crinkle in the peptide, leading to a core binding peptide of 10 amino acids, compared to the 9 amino acids in all other reported class II molecules. The limited number of potential T cell epitopes from such a complex virus can explain the differential MHC-determined resistance to MDV, but raises questions of mechanism and opportunities for vaccine targets in this important food species, as well as providing a basis for understanding class II molecules in other species including humans., This study shows that chicken MHC class II molecules present peptides from only a handful of the more than 100 genes of the oncogenic herpesvirus Marek’s disease virus, explaining the strong genetic association of chicken MHC with resistance and susceptibility to this and other economically-important pathogens.

Here, we report the application of a portable sequencer, MinION, for genotyping the malaria parasite Plasmodium falciparum. In the present study, an amplicon mixture of nine representative genes causing resistance to anti-malaria drugs is diagnosed. First, we developed the procedure for four laboratory strains (3D7, Dd2, 7G8, and K1), and then applied the developed procedure to ten clinical samples. We sequenced and re-sequenced the samples using the obsolete flow cell R7.3 and the most recent flow cell R9.4. Although the average base-call accuracy of the MinION sequencer was 74.3%, performing >50 reads at a given position improves the accuracy of the SNP call, yielding a precision and recall rate of 0.92 and 0.8, respectively, with flow cell R7.3. These numbers increased significantly with flow cell R9.4, in which the precision and recall are 1 and 0.97, respectively. Based on the SNP information, the drug resistance status in ten clinical samples was inferred. We also analyzed K13 gene mutations from 54 additional clinical samples as a proof of concept. We found that a novel amino-acid changing variation is dominant in this area. In addition, we performed a small population-based analysis using 3 and 5 cases (K13) and 10 and 5 cases (PfCRT) from Thailand and Vietnam, respectively. We identified distinct genotypes from the respective regions. This approach will change the standard methodology for the sequencing diagnosis of malaria parasites, especially in developing countries.

Photosynthetic organisms provide food and energy for nearly all life on Earth, yet half of their protein-coding genes remain uncharacterized1,2. Characterization of these genes could be greatly accelerated by new genetic resources for unicellular organisms. Here, we generated a genome-wide, indexed library of mapped insertion mutants for the unicellular alga Chlamydomonas reinhardtii. The 62,389 mutants in the library, covering 83% of nuclear, protein-coding genes, are available to the community. Each mutant contains unique DNA barcodes, allowing the collection to be screened as a pool. We performed a genome-wide survey of genes required for photosynthesis, which identified 303 candidate genes. Characterization of one of these genes, the conserved predicted phosphatase-encoding gene CPL3, showed it is important for accumulation of multiple photosynthetic protein complexes. Notably, 21 of the 43 highest-confidence genes are novel, opening new opportunities for advances in our understanding of this biogeochemically fundamental process. This library will accelerate the characterization of thousands of genes in algae, plants and animals., Editorial summary: Generation of a library of 62,389 mapped insertion mutants for the unicellular alga Chlamydomonas reinhardtii enables screening for genes required for photosynthesis and the identification of 303 candidate genes.

The phosphorylation status of a protein is highly regulated and is determined by the opposing activities of protein kinases and protein phosphatases within the cell. While much is known about the protein kinases found in Saccharomyces cerevisiae, the protein phosphatases are much less characterized. Of the 127 protein kinases in yeast, over 90% are in the same evolutionary lineage. In contrast, protein phosphatases are fewer in number (only 43 have been identified in yeast) and comprise multiple, distinct evolutionary lineages. Here we review the protein phosphatase families of yeast with regard to structure, catalytic mechanism, regulation, and signal transduction participation.

The propensity of cancer cells to transition between epithelial and mesenchymal phenotypic states via the epithelial–mesenchymal transition (EMT) program can regulate metastatic processes, cancer progression, and treatment resistance. Transcriptional investigations using reversible models of EMT, revealed the mesenchymal-to-epithelial reverting transition (MErT) to be enriched in clinical samples of metastatic castrate resistant prostate cancer (mCRPC). From this enrichment, a metastasis-derived gene signature was identified that predicted more rapid cancer relapse and reduced survival across multiple human carcinoma types. Additionally, the transcriptional profile of MErT is not a simple mirror image of EMT as tumour cells retain a transcriptional “memory” following a reversible EMT. This memory was also enriched in mCRPC samples. Cumulatively, our studies reveal the transcriptional profile of epithelial–mesenchymal plasticity and highlight the unique transcriptional properties of MErT. Furthermore, our findings provide evidence to support the association of epithelial plasticity with poor clinical outcomes in multiple human carcinoma types.

Background This study examined current physical activity levels and preferences for exercise settings and activities among adult survivors of childhood cancers as a strategy to inform the feasibility and design of such programs. Methods A mixed-methods design was used to investigate current activity levels as well as barriers to and preferences for physical activity among 20 adult survivors of pediatric cancer. Results One-half of participants reported engaging in regular physical activity, although the frequency, intensity, and duration varied. Overall, 17 of the 20 participants (85%) stated they would be interested in participating in a structured exercise intervention, and they expressed a strong interest in walking (76%), bicycling (53%), and weight training (53%). Common barriers to participation in a potential structured exercise program were insufficient time, current health issues, and program location/distance. Nearly all participants agreed that information on nutrition and diet should be included as part of an exercise intervention. Conclusions These findings will help inform the design and implementation of future exercise programs to enhance physical activity among this high-risk group of cancer survivors.

Significance Chloroplasts are of vital importance in photosynthetic eukaryotic organisms. Like mitochondria, they contain their own genomes. Nevertheless, most chloroplast proteins are encoded by nuclear genes, translated in the cytosol, and must cross the chloroplast envelope membranes to reach their proper destinations inside the organelle. Despite its fundamental role, our knowledge of the machinery catalyzing chloroplast protein import is incomplete and controversial. Here, we address the evolutionary conservation, composition, and function of the chloroplast protein import machinery using the green alga Chlamydomonas reinhardtii as a model system. Our findings help clarify the current debate regarding the composition of the chloroplast protein import machinery, provide evidence for cross-compartmental coordination of its biogenesis, and open promising avenues for its structural characterization., In photosynthetic eukaryotes, thousands of proteins are translated in the cytosol and imported into the chloroplast through the concerted action of two translocons—termed TOC and TIC—located in the outer and inner membranes of the chloroplast envelope, respectively. The degree to which the molecular composition of the TOC and TIC complexes is conserved over phylogenetic distances has remained controversial. Here, we combine transcriptomic, biochemical, and genetic tools in the green alga Chlamydomonas (Chlamydomonas reinhardtii) to demonstrate that, despite a lack of evident sequence conservation for some of its components, the algal TIC complex mirrors the molecular composition of a TIC complex from Arabidopsis thaliana. The Chlamydomonas TIC complex contains three nuclear-encoded subunits, Tic20, Tic56, and Tic100, and one chloroplast-encoded subunit, Tic214, and interacts with the TOC complex, as well as with several uncharacterized proteins to form a stable supercomplex (TIC-TOC), indicating that protein import across both envelope membranes is mechanistically coupled. Expression of the nuclear and chloroplast genes encoding both known and uncharacterized TIC-TOC components is highly coordinated, suggesting that a mechanism for regulating its biogenesis across compartmental boundaries must exist. Conditional repression of Tic214, the only chloroplast-encoded subunit in the TIC-TOC complex, impairs the import of chloroplast proteins with essential roles in chloroplast ribosome biogenesis and protein folding and induces a pleiotropic stress response, including several proteins involved in the chloroplast unfolded protein response. These findings underscore the functional importance of the TIC-TOC supercomplex in maintaining chloroplast proteostasis.

Photosynthetic organisms are essential for human life, yet most of their genes remain functionally uncharacterized. Single-celled photosynthetic model systems have the potential to accelerate our ability to connect genes to functions. Here, using a barcoded mutant library of the model eukaryotic alga Chlamydomonas reinhardtii, we determined the phenotypes of more than 58,000 mutants under more than 121 different environmental growth conditions and chemical treatments. 78% of genes are represented by at least one mutant that showed a phenotype, providing clues to the functions of thousands of genes. Mutant phenotypic profiles allow us to place known and previously uncharacterized genes into functional pathways such as DNA repair, photosynthesis, the CO2-concentrating mechanism, and ciliogenesis. We illustrate the value of this resource by validating novel phenotypes and gene functions, including the discovery of three novel components of a defense pathway that counteracts actin cytoskeleton inhibitors released by other organisms. The data also inform phenotype discovery in land plants: mutants in Arabidopsis thaliana genes exhibit similar phenotypes to those we observed in their Chlamydomonas homologs. We anticipate that this resource will guide the functional characterization of genes across the tree of life.

Viral diseases pose major threats to humans and other animals, including the billions of chickens that are an important food source as well as a public health concern due to zoonotic pathogens. Unlike humans and other typical mammals, the major histocompatibility complex (MHC) of chickens can confer decisive resistance or susceptibility to many viral diseases. An iconic example is Marek’s disease, caused by an oncogenic herpesvirus with over 100 genes. Classical MHC class I and class II molecules present antigenic peptides to T lymphocytes, and it has been hard to understand how such MHC molecules could be involved in susceptibility to Marek’s disease, given the potential number of peptides from over 100 genes. We used a new in vitro infection system and immunopeptidomics to determine peptide motifs for the two class II molecules expressed by the MHC haplotype B2, which is known to confer resistance to Marek’s disease. Surprisingly, we found that the vast majority of viral peptide epitopes presented by chicken class II molecules arise from only four viral genes, nearly all having the peptide motif for BL2, the dominantly-expressed class II molecule in chickens. We expressed BL2 linked to several MDV peptides, and determined one X-ray crystal structure, showing how a single small amino acid in the binding site causes a crinkle in the peptide, leading to core binding peptide of ten amino acids, compared to the nine amino acids in all other reported class II molecules. The limited number of potential T cell epitopes from such a complex virus can explain the differential MHC-determined resistance to MDV, but raises questions of mechanism and opportunities for vaccine targets in this important food species, as well as providing a basis for understanding class II molecules in other species including humans.