ALBERT

Description: Motivation Accurate automatic annotation of protein function relies on both innovative models and robust data sets. Due to their importance in biological processes, the identification of DNA-binding proteins directly from protein sequence has been the focus of many studies. However, the data sets used to train and evaluate these methods have suffered from substantial flaws. We describe some of the weaknesses of the data sets used in previous DNA-binding protein literature and provide several new data sets addressing these problems. We suggest new evaluative benchmark tasks that more realistically assess real-world performance for protein annotation models. We propose a simple new model for the prediction of DNA-binding proteins and compare its performance on the improved data sets to two previously published models. Additionally, we provide extensive tests showing how the best models predict across taxonomies. Results Our new gradient boosting model, which uses features derived from a published protein language model, outperforms the earlier models. Perhaps surprisingly, so does a baseline nearest neighbor model using BLAST percent identity. We evaluate the sensitivity of these models to perturbations of DNA-binding regions and control regions of protein sequences. The successful data-driven models learn to focus on DNA-binding regions. When predicting across taxonomies, the best models are highly accurate across species in the same kingdom and can provide some information when predicting across kingdoms. Code and Data Availability The data and results for this paper can be found at https://doi.org/10.5281/zenodo.5153906. The code for this paper can be found at https://doi.org/10.5281/zenodo.5153683. The code, data and results can also be found at https://github.com/AZaitzeff/tools_for_dna_binding_proteins.

Description: Motivation The investigation of quantitative trait loci (QTL) is an essential component in our understanding of how organisms vary phenotypically. However, many important crop species are polyploid (carrying more than two copies of each chromosome), requiring specialized tools for such analyses. Moreover, deciphering meiotic processes at higher ploidy levels is not straightforward, but is necessary to understand the reproductive dynamics of these species, or uncover potential barriers to their genetic improvement. Results Here, we present polyqtlR, a novel software tool to facilitate such analyses in (auto)polyploid crops. It performs QTL interval mapping in F1 populations of outcrossing polyploids of any ploidy level using identity-by-descent probabilities. The allelic composition of discovered QTL can be explored, enabling favourable alleles to be identified and tracked in the population. Visualization tools within the package facilitate this process, and options to include genetic co-factors and experimental factors are included. Detailed information on polyploid meiosis including prediction of multivalent pairing structures, detection of preferential chromosomal pairing and location of double reduction events can be performed. Availabilityand implementation polyqtlR is freely available from the Comprehensive R Archive Network (CRAN) at http://cran.r-project.org/package=polyqtlR. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Glioblastoma (GB) is the most aggressive and common form of primary brain tumor characterized by fast proliferation, high invasion and resistance to current standard treatment. The average survival rate post-diagnosis is 14.6 months, despite the aggressive standard post-surgery radiotherapy concomitant with chemotherapy with temozolomide (TMZ). Currently, efforts are being endowed to develop new and more efficient therapeutic approaches capable to overcome chemoresistance, inhibit tumor progression and improve overall patient survival rate. Abnormal microRNA (miRNA) expression has been correlated with chemoresistance, proliferation and resistance to apoptosis, which result from their master regulatory role of gene expression. Altered cell metabolism, favoring glycolysis, was identified as an emerging cancer hallmark and has been described in GB, thus offering a new target for innovative GB therapies. In this work, we hypothesized that a gene therapy-based strategy consisting of the overexpression of a miRNA downregulated in GB and predicted to target crucial metabolic enzymes might promote a shift of GB cell metabolism, decreasing the glycolytic dependence of tumor cells and contributing to their sensitization to chemotherapy with TMZ. The increase of miR-200c levels in DBTRG cells resulted in downregulation of messenger RNA of enzymes involved in bioenergetics pathways and impaired cell metabolism and mobility. In addition, miR-200c overexpression prior to DBTRG cell exposure to TMZ resulted in cell cycle arrest. Overall, our results show that miR-200c overexpression could offer a way to overcome chemoresistance developed by GB cells in response to current standard chemotherapy, providing an improvement to current GB standard treatment, with benefit for patient outcome.

Description: Motivation The investigation of the structure of biological systems at the molecular level gives insights about their functions and dynamics. Shape and surface of biomolecules are fundamental to molecular recognition events. Characterizing their geometry can lead to more adequate predictions of their interactions. In the present work, we assess the performance of reference shape retrieval methods from the computer vision community on protein shapes. Results Shape retrieval methods are efficient in identifying orthologous proteins and tracking large conformational changes. This work illustrates the interest for the protein surface shape as a higher-level representation of the protein structure that (i) abstracts the underlying protein sequence, structure or fold, (ii) allows the use of shape retrieval methods to screen large databases of protein structures to identify surficial homologs and possible interacting partners and (iii) opens an extension of the protein structure–function paradigm toward a protein structure-surface(s)-function paradigm. Availabilityand implementation All data are available online at http://datasetmachat.drugdesign.fr. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Contactin 4 (CNTN4) is a crucial synaptic adhesion protein that belongs to the contactin superfamily. Evidence from both human genetics and mouse models suggests that synapse formation and structural deficits strongly correlate with neurodevelopmental disorders, including autism. In addition, several lines of evidence suggest that CNTN4 is associated with the risk of autism. However, the biological functions of CNTN4 in neural development and disease pathogenesis are poorly understood. In this study, we investigated whether and how CNTN4 is autonomously involved in the development of dendrites and dendritic spines in cortical neurons. Disruption of Cntn4 decreased the number of excitatory synapses, which led to a reduction in neural activity. Truncated proteins lacking the signal peptide, FnIII domains, or GPI domain lacked the ability to regulate dendritic spine formation, indicating that CNTN4 regulates dendritic spine density through a mechanism dependent on FnIII domains. Importantly, we revealed that autism-related variants lacked the ability to regulate spine density and neural activity. In conclusion, our study suggests that CNTN4 is essential for promoting dendrite growth and dendritic spine formation and that disruptive variants of CNTN4 interfere with abnormal synapse formation and may increase the risk of autism.

Description: Numerous genome-wide association studies (GWASs) have been conducted for the identification of genetic variants involved with human height. The vast majority of these studies, however, have been conducted in populations of European ancestry. Here, we report the first GWAS of adult height in the Taiwan Biobank using a discovery sample of 14 571 individuals and an independent replication sample of 20 506 individuals. From our analysis, we generalize to the Taiwanese population genome-wide significant associations with height and 18 previously identified genes in European and non-Taiwanese East Asian populations. We also identify and replicate, at the genome-wide significance level, associated variants for height in four novel genes at two loci that have not previously been reported: RASA2 on chromosome 3 and NABP2, RNF41 and SLC39A5 at 12q13.3 on chromosome 12. RASA2 and RNF41 are strong candidates for having a role in height with copy number and loss of function variants in RASA2 previously found to be associated with short stature disorders, and decreased expression of the RNF41 gene resulting in insulin resistance in skeletal muscle. The results from our analysis of the Taiwan Biobank underscore the potential for the identification of novel genetic discoveries in underrepresented worldwide populations, even for traits, such as height, that have been extensively investigated in large-scale studies of European ancestry populations.

Description: Motivation The mathematically optimal solution in computational protein folding simulations does not always correspond to the native structure, due to the imperfection of the energy force fields. There is therefore a need to search for more diverse suboptimal solutions in order to identify the states close to the native. We propose a novel multimodal optimization protocol to improve the conformation sampling efficiency and modeling accuracy of de novo protein structure folding simulations. Results A distance-assisted multimodal optimization sampling algorithm, MMpred, is proposed for de novo protein structure prediction. The protocol consists of three stages: The first is a modal exploration stage, in which a structural similarity evaluation model DMscore is designed to control the diversity of conformations, generating a population of diverse structures in different low-energy basins. The second is a modal maintaining stage, where an adaptive clustering algorithm MNDcluster is proposed to divide the populations and merge the modal by adjusting the annealing temperature to locate the promising basins. In the last stage of modal exploitation, a greedy search strategy is used to accelerate the convergence of the modal. Distance constraint information is used to construct the conformation scoring model to guide sampling. MMpred is tested on a large set of 320 non-redundant proteins, where MMpred obtains models with TM-score≥0.5 on 291 cases, which is 28% higher than that of Rosetta guided with the same set of distance constraints. In addition, on 320 benchmark proteins, the enhanced version of MMpred (E-MMpred) has 167 targets better than trRosetta when the best of five models are evaluated. The average TM-score of the best model of E-MMpred is 0.732, which is comparable to trRosetta (0.730). Availability and implementation The source code and executable are freely available at https://github.com/iobio-zjut/MMpred. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Asthenoteratospermia is a common cause of male infertility. Recent studies have revealed that CFAP65 mutations lead to severe asthenoteratospermia due to acrosome hypoplasia and flagellum malformations. However, the molecular mechanism underlying CFAP65-associated sperm malformation is largely unclear. Here, we initially examined the role of CFAP65 during spermiogenesis using Cfap65 knockout (Cfap65−/−) mice. The results showed that Cfap65−/− male mice exhibited severe asthenoteratospermia characterized by morphologically defective sperm heads and flagella. In Cfap65−/− mouse testes, hyper-constricted sperm heads were apparent in step 9 spermatids accompanied by abnormal manchette development, and acrosome biogenesis was abnormal in the maturation phase. Moreover, subsequent flagellar elongation was also severely affected and characterized by disrupted assembly of the mitochondrial sheath (MS) in Cfap65−/− male mice. Furthermore, the proteomic analysis revealed that the proteostatic system during acrosome formation, manchette organization and MS assembly was disrupted when CFAP65 was lost. Importantly, endogenous immunoprecipitation and immunostaining experiments revealed that CFAP65 may form a cytoplasmic protein network comprising MNS1, RSPH1, TPPP2, ZPBP1 and SPACA1. Overall, these findings provide insights into the complex molecular mechanisms of spermiogenesis by uncovering the essential roles of CFAP65 during sperm head shaping, acrosome biogenesis and MS assembly.

Description: Mitochondrial DNA (mtDNA) disorders are recognized as one of the most common causes of inherited metabolic disorders. The mitochondrial genome occurs in multiple copies resulting in both homoplasmic and heteroplasmic pathogenic mtDNA variants. A biochemical defect arises when the pathogenic variant level reaches a threshold, which differs between variants. Moreover, variants can segregate, clonally expand, or be lost from cellular populations resulting in a dynamic and tissue-specific mosaic pattern of oxidative deficiency. MtDNA is maternally inherited but transmission patterns of heteroplasmic pathogenic variants are complex. During oogenesis, a mitochondrial bottleneck results in offspring with widely differing variant levels to their mother, whilst highly deleterious variants, such as deletions, are not transmitted. Complemented by a complex interplay between mitochondrial and nuclear genomes, these peculiar genetics produce marked phenotypic variation, posing challenges to the diagnosis and clinical management of patients. Novel therapeutic compounds and several genetic therapies are currently under investigation, but proven disease-modifying therapies remain elusive. Women who carry pathogenic mtDNA variants require bespoke genetic counselling to determine their reproductive options. Recent advances in in vitro fertilization techniques, have greatly improved reproductive choices, but are not without their challenges. Since the first pathogenic mtDNA variants were identified over 30 years ago, there has been remarkable progress in our understanding of these diseases. However, many questions remain unanswered and future studies are required to investigate the mechanisms of disease progression and to identify new disease-specific therapeutic targets.

Description: Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP), and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the major histocompatibility complex region [Odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.66–0.83, lowest P-value = 4.53 × 10−8]. The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95% CI = 0.90–0.97, P-value = 3.54 × 10−4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the human leukocyte antigen genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease.

Description: Motivation Co-evolution analysis can be used to accurately predict residue–residue contacts from multiple sequence alignments. The introduction of machine-learning techniques has enabled substantial improvements in precision and a shift from predicting binary contacts to predict distances between pairs of residues. These developments have significantly improved the accuracy of de novo prediction of static protein structures. With AlphaFold2 lifting the accuracy of some predicted protein models close to experimental levels, structure prediction research will move on to other challenges. One of those areas is the prediction of more than one conformation of a protein. Here, we examine the potential of residue–residue distance predictions to be informative of protein flexibility rather than simply static structure. Results We used DMPfold to predict distance distributions for every residue pair in a set of proteins that showed both rigid and flexible behaviour. Residue pairs that were in contact in at least one reference structure were classified as rigid, flexible or neither. The predicted distance distribution of each residue pair was analysed for local maxima of probability indicating the most likely distance or distances between a pair of residues. We found that rigid residue pairs tended to have only a single local maximum in their predicted distance distributions while flexible residue pairs more often had multiple local maxima. These results suggest that the shape of predicted distance distributions contains information on the rigidity or flexibility of a protein and its constituent residues. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Antimicrobial peptides harboring S- and or O-linked glycans are known as glycocins. Glycocins were first discovered and best characterized in Firmicutes. S-glycosylation is an enzymatic process catalyzed by S-glycosyltransferases of the GT2 family. Using a heterologous expression system, here we describe an inverting S/O-HexNAc-transferase (SvGT), encoded by ORF AQF52_3101 of S. venezuelae ATCC 15439, along with its acceptor substrate (SvC), encoded by ORF AQF52_3099. Using in vitro and in vivo assays, we define the distinct donor specificity, acceptor specificity, regioselectivity, chemoselectivity, and Y(G/A/K/Q/E ≠ ΔG)(C/S/T ≠ Y/N)(G/A ≠ P/Q)G as the minimum acceptor sequon of SvGT. Although UDP-GlcNAc served as the donor in the cellular milieu, SvGT could also utilize UDP-Glc and UDP-GalNAc as donors in vitro. Using mass spectrometry and western blotting, we provide evidence that an anti-O-GlcNAc antibody (CTD110.6) cross-reacts with S-GlcNAc and may be used to detect S-GlcNAcylated glycoconjugates directly. With an understanding of enzyme specificities, we finally employed SvGT to generate two proof-of-concept neoglycocins against L. monocytogenes. In conclusion, this study provides the first experimental evidence for S-glycosylation in Actinobacteria and the application of its S/O-HexNAc-transferase in glycocin engineering.

Description: Truncated O-GalNAc glycosylation is an important feature of pancreatic ductal adenocarcinomas (PDAC) and expression of truncated O-GalNAc glycans is strongly associated with decreased survival and poor prognosis. It has been proven, that aberrant O-GalNAc glycosylation influence PDAC signaling to promote oncogenic properties, but elucidation of the influence of truncated O-GalNAc glycosylation on different signaling molecules has just been started. We herein elucidated the impact of aberrant O-GalNAc glycosylation on two important PDAC signaling pathways, namely AKT/mTOR and RAS/MAPK. In PDAC cells expressing truncated O-GalNAc glycans, we identified differentially expressed proteins associated with AKT/mTOR and RAS/MAPK pathways using quantitative proteomics. Since AKT, a key-signaling molecule in PDAC, was among the identified proteins, we analyzed AKT and found a strikingly enhanced S473 phosphorylation and identified a previously unknown O-GalNAc-modification. Consecutive analysis of COSMC knockdowns in PDAC revealed strong effects on AKT upstream and downstream effector molecules. Interestingly, truncated O-GalNAc glycans could facilitate an mTORC1 inhibitor resistance using AZD8055. In addition, as AKT/mTOR pathway has extensive cross talks with RAS/MAPK pathway we analyzed the pathways and found it negatively regulated. Finally, we found that the expression of epithelial-mesenchymal-transition markers, key features of aggressive PDACs cells, are enhanced and truncated O-GalNAc glycans enhance pancreatic cancer cell growth in a xenograft mouse model. Our study demonstrates that truncated O-GalNAc glycans have a strong impact on AKT/mTOR and RAS/MAPK signaling pathways, are modulated by EGF or IGF-1 signaling and should be considered for targeted therapy of these pathways in PDAC.

Description: Congenital disorders of glycosylation (CDG) include 150 disorders constituting in genetically and clinically heterogeneous diseases, showing significant glycoprotein hypoglycosylation that leads to pathological consequences on multiple organs and systems which underlying mechanisms are not yet understood. A few cellular and animal models have been used to study specific CDG characteristics although they have given limited information due to the few CDG mutations tested and the still missing comprehensive molecular and cellular basic research. Here we provide specific gene expression profiles, based on RNA microarray analysis, together with some biochemical and cellular characteristics of a total of 9 control EBV-transformed lymphoblastoid B cell lines (B-LCL) and 13 CDG B-LCL from patients carrying severe mutations in the PMM2 gene, strong serum protein hypoglycosylation and neurological symptoms. Significantly dysregulated genes in PMM2-CDG cells included those regulating stress responses, transcription factors, glycosylation, motility, cell junction and, importantly, those related to development and neuronal differentiation and synapse such as CA2 and ADAM23. PMM2-CDG associated biological consequences involved the unfolded protein response, RNA metabolism and the endoplasmic reticulum, Golgi apparatus and mitochondria components. Changes in transcriptional and CA2 protein levels are consistent with CDG physiopathology. These results demonstrate the global transcriptional impact in phosphomannomutase 2 deficient cells, reveal CA2 as a potential cellular biomarker and confirm B-LCL as an advantageous model for CDG studies.

Description: Motivation The well-known fact that protein structures are more conserved than their sequences forms the basis of several areas of computational structural biology. Methods based on the structure analysis provide more complete information on residue conservation in evolutionary processes. This is crucial for the determination of evolutionary relationships between proteins and for the identification of recurrent structural patterns present in biomolecules involved in similar functions. However, algorithmic structural alignment is much more difficult than multiple sequence alignment. This study is devoted to the development and applications of DAMA—a novel effective environment capable to compute and analyze multiple structure alignments. Results DAMA is based on local structural similarities, using local 3D structure descriptors and thus accounts for nearest-neighbor molecular environments of aligned residues. It is constrained neither by protein topology nor by its global structure. DAMA is an extension of our previous study (DEDAL) which demonstrated the applicability of local descriptors to pairwise alignment problems. Since the multiple alignment problem is NP-complete, an effective heuristic approach has been developed without imposing any artificial constraints. The alignment algorithm searches for the largest, consistent ensemble of similar descriptors. The new method is capable to capture most of the biologically significant similarities present in canonical test sets and is discriminatory enough to prevent the emergence of larger, but meaningless, solutions. Tests performed on the test sets, including protein kinases, demonstrate DAMA’s capability of identifying equivalent residues, which should be very useful in discovering the biological nature of proteins similarity. Performance profiles show the advantage of DAMA over other methods, in particular when using a strict similarity measure QC, which is the ratio of correctly aligned columns, and when applying the methods to more difficult cases. Availability and implementation DAMA is available online at http://dworkowa.imdik.pan.pl/EP/DAMA. Linux binaries of the software are available upon request. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Sperm fibrous sheath is closely related to sperm maturation, capacitation and motility, and A-kinase anchor protein 4 (AKAP4) is the most abundant protein in sperm fibrous sheath. Previous studies found incomplete sperm fibrous sheaths and abnormal flagella in Akap4 knockout (KO) mice. Meanwhile, it was reported that the partial deletion in AKAP4 is highly relevant to the dysplasia of the fibrous sheath in an infertile man, and so far, there is no report about male infertility caused by hemizygous AKAP4 variant. Furthermore, the specific mechanisms of how the variant is relevant to the phenotype remain elusive. In this study, we investigated three multiple morphological abnormalities of the sperm flagella (MMAF)-affected men from three independent families (including one consanguine family) carried hemizygous c.C1285T variant in AKAP4. The patients carried thisvariant showed dysplastic sperm fibrous sheath and the protein expression of AKAP4 was decreased in flagella which was further confirmed in HEK-293 T cells in vitro. In addition, the co-localization and interaction between AKAP4 and glutamine-rich protein 2 (QRICH2) on the molecular level were identified by immunofluorescence and Co-immunoprecipitation (CO-IP). The hemizygous c.1285C 〉 T variant in AKAP4 induced decreased protein expression of QRICH2 in spermatozoa. These results suggested that the normal expression of AKAP4 is required for maintaining the expression of QRICH2 and the decreased protein expression of AKAP4 and QRICH2，as well as the interaction between them induced by the hemizygous variant of AKAP4 caused dysplastic fibrous sheath, which eventually led to reduced sperm motility and male infertility.

Description: Sulfomucins are in some body-locations and species a normal occurrence, whereas in other situations a sign of pathology. Sulfomucin content on histological sections and isolated material is frequently analyzed with Alcian blue staining at pH 1.0. However, since the stain detects the charge, a high density of other charged molecules, such as sialic acids has potential to impede specificity. Here, we compared the outcome from four staining protocols with the level of sulfation determined by liquid chromatography–tandem mass spectrometry analysis (MS) on samples from various tissues with variable sulfation and sialylation levels. We found that a protocol we designed, including rinsing with MetOH and 0.5 M NaCl buffer at pH 1.0 eliminates false positive staining of tissues outperforming commonly recommended solutions. In tissues with low to moderately sulfated mucins (e.g. human stomach and salmonid epithelia) this method enables accurate relative quantification (e.g. sulfate scoring comparisons between healthy and diseased tissues) whereas the range of the method is not suitable for comparisons between tissues with high sulfomucin content (e.g. pig stomach and colon).

Description: Nephronophthisis (NPH) is the most prevalent monogenetic disorder leading to end-stage renal failure (ESRD) in childhood. Mutations in Nphp1, encoding a cilia-localized protein, account for the majority of NPH cases. Despite its identification many years ago, Nphp1 deletions targeting exon 4 or exon 20 have not reproduced the histological features of human NPH in murine models. In this study, we deleted exon 2–20 of Nphp1 by CRISPR/Cas9 gene editing to create a near-total knockout (KO) mouse model (Nphp1del2–20/del2–20). Nphp1del2–20/del2–20 mice faithfully reproduced the renal and extrarenal phenotypes associated with human NPH, including renal cyst development, tubular basement membrane thickening, retinal degeneration and abnormal spermatogenesis. Importantly, Nphp1 re-expression using an adenoviral-associated-virus-9 (AAV9) vector could partially rescue both renal and retinal phenotypes in Nphp1del2–20/del2–20 mice. Our results reported the first relevant Nphp1 mouse model with renal phenotypes for human disease. It will be a valuable model for future studies of Nphp1 function and to develop novel treatments for this common childhood disease.

Description: Recent studies show decreasing prostate-specific antigen utilization and increasing incidence of metastatic prostate cancer in the United States after national recommendations against screening in 2012. Yet whether the increasing incidence of metastatic prostate cancer is consistent in magnitude with the expected impact of decreased screening is unknown. We compared observed incidence of metastatic prostate cancer from the Surveillance, Epidemiology, and End Results program and published effects of continued historical screening and discontinued screening starting in 2013 projected by two models of disease natural history, screening, and diagnosis. The observed rate of new metastatic prostate cancer cases in 2017 was 44%-60% of the projected increase under discontinued screening relative to continued screening. Thus, the observed increase in incident metastatic prostate cancer is consistent with the expected impact of reduced screening. Although this comparison does not establish a causal relationship, it highlights the plausible role of decreased screening in the observed trend.

Description: Background Few studies have examined gene-specific associations with contralateral and/or second breast cancer (SBC). Methods The frequency of pathogenic/likely pathogenic (P/LP) variants in clinically-actionable genes (BRCA1, BRCA2, PTEN, TP53, CHEK2, CDH1, ATM, PALB2, NBN, and NF1) was compared between women with a primary breast cancer (PBC) and SBC who underwent multi-gene panel testing at a single diagnostic testing laboratory. Race/ethnicity-specific logistic regression burden tests adjusted for age at diagnosis of first breast cancer, histology, presence of first- or second-degree relatives with breast cancer, and prior testing for BRCA1/2 genes, were used to test for associations with SBC. All statistical tests were two-sided. Results The study was comprised of 75,550 women with PBC and 7,728 with SBC. Median time between breast cancers for SBC was 11 (IQR=6,17) years. Restricting to women tested for all actionable genes (n = 60,310), there were 4,231 (7.8%) carriers of P/LP variants in actionable genes among the controls (PBC) compared to 652 (11.1%) women with SBC (p〈 .001). Among Caucasians, exclusive of Ashkenazi Jewish women, those carrying a P/LP variant in a clinically actionable gene were 1.44 (95% CI: 1.30-1.60) times as likely to have SBC than non-carriers, after accounting for potential confounders. Among African American and Hispanic women, a P/LP variant in a clinically actionable gene was 1.88 (1.36-2.56) and 1.66 (1.02-2.58) times as likely to be associated with SBC respectively (p

Description: The Sahel/Savannah belt of Africa is a contact zone between two subsistence systems (nomadic pastoralism and sedentary farming) and of two groups of populations, namely Eurasians penetrating from northern Africa southwards and sub-Saharan Africans migrating northwards. Because pastoralism is characterised by a high degree of mobility, it leaves few significant archaeological traces. Demographic history seen through the lens of population genetic studies complements our historical and archaeological knowledge in this African region. In this review, we highlight recent advances in our understanding of demographic history in the Sahel/Savannah belt as revealed by genetic studies. We show the impact of food-producing subsistence strategies on population structure as well as the somewhat different migration patterns in the western and eastern part of the region. Genomic studies show that the gene pool of various groups of Sahelians consists in a complex mosaic of several ancestries. We also touch upon various signals of genetic adaptations such as lactase persistence, taste sensitivity, and malaria resistance, all of which have different distribution patterns among Sahelian populations. Overall, genetic studies contribute to gain a deeper understanding about the demographic and adaptive history of human populations in this specific African region and beyond.

Description: Talipes equinovarus (clubfoot, TEV) is a congenital rotational foot deformity occurring in 1 per 1000 births with increased prevalence in males compared to females. The genetic etiology of isolated clubfoot (iTEV) remains unclear. Using a genome-wide association study, we identified a locus within FSTL5, encoding follistatin-like 5, significantly associated with iTEV. FSTL5 is an uncharacterized gene whose potential role in embryonic and post-natal development were previously unstudied. Utilizing multiple model systems, we found that Fstl5 was expressed during later stages of embryonic hindlimb development, and, in mice, expression was restricted to the condensing cartilage anlage destined to form the limb skeleton. In the post-natal growth plate, Fstl5 was specifically expressed in pre-hypertrophic chondrocytes. As Fstl5 knock-out rats displayed no gross malformations, we engineered a conditional transgenic mouse line (Fstl5LSL) to over-express Fstl5 in skeletal osteochondroprogenitors. We observed that hindlimbs were slightly shorter and that bone mineral density was reduced in adult male, but not female, Prrx1-cre;Fstl5LSL mice compared to control. No overt clubfoot-like deformity was observed in Prrx1-cre;Fstl5LSL mice, suggesting FSTL5 may function in other cell types to contribute to iTEV pathogenesis. Interrogating published mouse embryonic single-cell expression data showed Fstl5 was expressed in cell lineage sub-clusters whose transcriptomes were associated with neural system development. Moreover, our results suggest lineage-specific expression of the Fstl genes correlates with their divergent roles as modulators of TGF-ß and BMP signaling. Results from this study associate FSTL5 with iTEV and suggest a potential sexually dimorphic role for Fstl5 in vivo.

Description: The in vivo Pig-a assay is being used in safety studies to evaluate the potential of chemicals to induce somatic cell gene mutations. Ongoing work is aimed at developing an Organisation for Economic Cooperation and Development (OECD) test guideline to support routine use for regulatory purposes (OECD project number 4.93). Among the details that will need to be articulated in an eventual guideline are recommended treatment and harvest schedules. With this in mind, experiments reported herein were performed with Wistar Han rats exposed to aristolochic acid I (AA), 1,3-propane sultone, chlorambucil, thiotepa or melphalan using each of two commonly used treatment schedules: 3 or 28 consecutive days. In the case of the 3-day studies, blood was collected for Pig-a analysis on days 15 or 16 and 29 or 30. For the 28-day studies blood was collected on day 29 or 30. The effect of treatment on mutant reticulocytes and mutant erythrocytes was evaluated with parametric pair-wise tests. While each of the five mutagens increased mutant phenotype cell frequencies irrespective of study design, statistical significance was consistently achieved at lower dose levels when the 28-day format was used (e.g. 2.75 vs 20 mg/kg/bw for AA). To more thoroughly investigate the dose–response relationships, benchmark dose (BMD) analyses were performed with PROAST software. These results corroborate the pair-wise testing results in that lower BMD values were obtained with the 28-day design. Finally, mutagenic potency, as measured by BMD analyses, most consistently correlated with the mutagens’ tumorigenic dose 50 values when the lengthier treatment schedule was used. Collectively, these results suggest that both 3- and 28-day treatment schedules have merit in hazard identification-type studies. That being said, for the purpose of regulatory safety assessments, there are clear advantages to study designs that utilise protracted exposures.

Description: The chemotherapeutic efficacy in colorectal cancer (CRC) is limited due to the inter-individual variability in drug response and the development of tumour resistance. ATP-binding cassette (ABC) transporters are crucial in the development of resistance by the efflux of anticancer agents from cancer cells. In this study, we identified 14 single nucleotide polymorphisms (SNPs) in 11 ABC transporter genes acting as an expression of quantitative trait loci (eQTLs), i.e. whose variation influence the expression of many downstream genes. These SNPs were genotyped in a case–control study comprising 1098 cases and 1442 healthy controls and analysed in relation to CRC development risk and patient survival. Considering a strict correction for multiple tests, we did not observe any significant association between SNPs and CRC risk. The rs3819720 polymorphism in the ABCB3/TAP2 gene was statistically significantly associated with shorter overall survival (OS) in the codominant, and dominant models [GA vs. GG, hazard ratio (HR) = 1.48; P = 0.002; AA vs. GG, HR = 1.70; P = 0.004 and GA + AA vs. GG, HR = 1.52; P = 0.0006]. Additionally, GA carriers of the same SNP displayed worse OS after receiving 5-FU based chemotherapy. The variant allele of rs3819720 polymorphism statistically significantly affected the expression of 36 downstream genes. Screening for eQTL polymorphisms in relevant genes such as ABC transporters that can regulate the expression of several other genes may help to identify the genetic background involved in the individual response to the treatment of CRC patients.

Description: Cancer deaths account for nearly 10 million deaths worldwide each year, with lung cancer (LCa) as the leading cause of cancer-related death. Smoking is one of the major LCa risk factors, and tobacco-related carcinogens are potent mutagens and epi-mutagens. In the present study, we aimed to analyse smoking-related epigenetic changes in lung tissues from LCa cases. The study cohort consisted of paired LCa and noncancerous lung tissues (NLT) from 104 patients, 90 of whom were smokers or ex-smokers (i.e. ever smokers) at the time of diagnosis. DNA methylation status of tumour suppressor genes DAPK1, MGMT, p16, RASSF1 and RARB was screened by means of methylation-specific PCR (MSP) and further analysed quantitatively by pyrosequencing. Methylation of at least one gene was detected in 59% (61 of 104) of LCa samples and in 39% (41 of 104) of NLT. DAPK1 and RASSF1 were more frequently methylated in LCa than in NLT (P = 0.022 and P = 0.041, respectively). The levels of DNA methylation were higher in LCa than NLT at most of the analysed CpG positions. More frequent methylation of at least one gene was observed in LCa samples of ever smokers (63%, 57 of 90) as compared with never smokers (36%, 5 of 14; P = 0.019). In the ever smokers group, methylation of the genes also occurred in NLT, but was rare or absent in the samples of never smokers. Among the current smokers, RASSF1 methylation in LCa showed association with the number of cigarettes smoked per day (P = 0.017), whereas in NLT it was positively associated with the duration of smoking (P = 0.039). Similarly, p16 methylation in LCa of current smokers correlated with the larger number of cigarettes smoked per day (P = 0.047). Overall, DNA methylation changes were present in both cancerous and noncancerous tissues of LCa patients and showed associations with smoking-related parameters.

Description: The consumption of fructose during pregnancy can cause hyperglycaemia and may stimulate production of reactive oxygen species; however, there are only a few studies reporting whether fructose consumption during pregnancy causes DNA damage. Therefore, the aim of this study was to evaluate the effects of fructose consumption on genetic and biochemical parameters in Swiss mice treated during pregnancy and lactation. For this, 15 couples of 60-day-old Swiss mice were divided into three groups of five couples: negative control (water) and two fructose groups (fructose dose of 10%/l and 20%/l). During this period, we evaluated food consumption, energy efficiency and body weight. Samples of blood were collected from the females before copulation, after the 15th day of conception and on the 21st day after the lactation period, for the glycaemic and lipid profiles as well as comet assay and micronucleus (MN) test. Comet assay and MN test evaluate DNA damage and clastogenicity, respectively. In the gestation and lactation period, the two fructose doses tested showed DNA damage as observed in the comet assay, which is associated with an increase in dietary intake, body weight, lipid profile and fasting glycaemia in females. Thus, it can be suggested that the high consumption of fructose during these periods is harmful for pregnancy and lactation.

Description: The antitumour activity of chrysin have been studied in several types of cancer cells. In urinary bladder cancer, its cytotoxic effects have already demonstrated; however, its mechanism of action is not completely understood and the role of tumour protein p53 (TP53) gene in these effects is unclear. In this study, we investigated the role of chrysin (10, 20, 40, 60 80 and 100 µM) in progression of bladder tumour cells with different status of the TP53 gene and different degrees of tumour (RT4, grade 1, TP53 wild type; 5637, grade 2, TP53 mutated and T24, grade 3, TP53 mutated). Results demonstrated that chrysin inhibited cell proliferation by increasing reactive oxygen species and DNA damage and inhibited cell migration in all cell lines. In TP53 wild-type cells, a sub-G1 apoptotic population was present. In mutated TP53 cells, chrysin caused arrest at the G2/M phase and morphological changes accompanied by downregulation of PLK1, SRC and HOXB3 genes. In addition, in Grade 2 cells, chrysin induced global DNA hypermethylation and, in the highest-grade cells, downregulated c-MYC, FGFR3 and mTOR gene expression. In conclusion, chrysin has antiproliferative and toxicogenetic activity in bladder tumour cells independently of TP53 status; however, the mechanisms of action are dependent on TP53 status.

Description: Tuberous sclerosis complex (TSC) is a rare autosomal dominant neurodevelopmental disorder characterized by variable expressivity. TSC results from inactivating variants within the TSC1 or TSC2 genes, leading to constitutive activation of mechanistic target of rapamycin complex 1 signaling. Using a mouse model of TSC (Tsc2-RG) in which the Tsc2 gene is deleted in radial glial precursors and their neuronal and glial descendants, we observed increased ornithine decarboxylase (ODC) enzymatic activity and concentration of its product, putrescine. To test if increased ODC activity and dysregulated polyamine metabolism contribute to the neurodevelopmental defects of Tsc2-RG mice, we used pharmacologic and genetic approaches to reduce ODC activity in Tsc2-RG mice, followed by histologic assessment of brain development. We observed that decreasing ODC activity and putrescine levels in Tsc2-RG mice worsened many of the neurodevelopmental phenotypes, including brain growth and neuronal migration defects, astrogliosis and oxidative stress. These data suggest a protective effect of increased ODC activity and elevated putrescine that modify the phenotype in this developmental Tsc2-RG model.

Description: Ageing is the major risk factor for Alzheimer’s disease (AD), a condition involving brain hypoxia. The majority of early-onset familial AD (EOfAD) cases involve dominant mutations in the gene PSEN1. PSEN1 null mutations do not cause EOfAD. We exploited putative hypomorphic and EOfAD-like mutations in the zebrafish psen1 gene to explore the effects of age and genotype on brain responses to acute hypoxia. Both mutations accelerate age-dependent changes in hypoxia-sensitive gene expression supporting that ageing is necessary, but insufficient, for AD occurrence. Curiously, the responses to acute hypoxia become inverted in extremely aged fish. This is associated with an apparent inability to upregulate glycolysis. Wild-type PSEN1 allele expression is reduced in post-mortem brains of human EOfAD mutation carriers (and extremely aged fish), possibly contributing to EOfAD pathogenesis. We also observed that age-dependent loss of HIF1 stabilization under hypoxia is a phenomenon conserved across vertebrate classes.

Description: Heterozygous de novo mutations in EEF1A2, encoding the tissue-specific translation elongation factor eEF1A2, have been shown to cause neurodevelopmental disorders including often severe epilepsy and intellectual disability. The mutational profile is unusual; ~50 different missense mutations have been identified but no obvious loss of function mutations, though large heterozygous deletions are known to be compatible with life. A key question is whether the heterozygous missense mutations operate through haploinsufficiency or a gain of function mechanism, an important prerequisite for design of therapeutic strategies. In order both to address this question and to provide a novel model for neurodevelopmental disorders resulting from mutations in EEF1A2, we created a new mouse model of the D252H mutation. This mutation causes the eEF1A2 protein to be expressed at lower levels in brain but higher in muscle in the mice. We compared both heterozygous and homozygous D252H and null mutant mice using behavioural and motor phenotyping alongside molecular modelling and analysis of binding partners. Although the proteomic analysis pointed to a loss of function for the D252H mutant protein, the D252H homozygous mice were more severely affected than null homozygotes on the same genetic background. Mice that are heterozygous for the missense mutation show no behavioural abnormalities but do have sex-specific deficits in body mass and motor function. The phenotyping of our novel mouse lines, together with analysis of molecular modelling and interacting proteins, suggest that the D252H mutation results in a gain of function.

Description: ADCA-DN and HSN-IE are rare neurodegenerative syndromes caused by dominant mutations in the replication foci targeting sequence (RFTS) of the DNA methyltransferase 1 (DNMT1) gene. Both phenotypes resemble mitochondrial disorders, and mitochondrial dysfunction was first observed in ADCA-DN. To explore mitochondrial involvement, we studied the effects of DNMT1 mutations in fibroblasts from four ADCA-DN and two HSN-IE patients. We documented impaired activity of purified DNMT1 mutant proteins, which in fibroblasts results in increased DNMT1 amount. We demonstrated that DNMT1 is not localized within mitochondria, but it is associated with the mitochondrial outer membrane. Concordantly, mitochondrial DNA failed to show meaningful CpG methylation. Strikingly, we found activated mitobiogenesis and OXPHOS with significant increase of H2O2, sharply contrasting with a reduced ATP content. Metabolomics profiling of mutant cells highlighted purine, arginine/urea cycle and glutamate metabolisms as the most consistently altered pathways, similar to primary mitochondrial diseases. The most severe mutations showed activation of energy shortage AMPK-dependent sensing, leading to mTORC1 inhibition. We propose that DNMT1 RFTS mutations deregulate metabolism lowering ATP levels, as a result of increased purine catabolism and urea cycle pathways. This is associated with a paradoxical mitochondrial hyper-function and increased oxidative stress, possibly resulting in neurodegeneration in non-dividing cells.

Description: Retinal degeneration is a common clinical feature of ciliopathies, a group of genetic diseases linked to ciliary dysfunction, and gene therapy is an attractive treatment option to prevent vision loss. Although the efficacy of retinal gene therapy is well established by multiple proof-of-concept preclinical studies, its long-term effect, particularly when treatments are given at advanced disease stages, is controversial. Incomplete treatment and intrinsic variability of gene delivery methods may contribute to the variable outcomes. Here, we used a genetic rescue approach to ‘optimally’ treat retinal degeneration at various disease stages and examined the long-term efficacy of gene therapy in a mouse model of ciliopathy. We used a Bardet–Biedl syndrome type 17 (BBS17) mouse model, in which the gene-trap that suppresses Bbs17 (also known as Lztfl1) expression can be removed by tamoxifen administration, restoring normal gene expression systemically. Our data indicate that therapeutic effects of retinal gene therapy decrease gradually as treatments are given at later stages. These results suggest the presence of limited time window for successful gene therapy in certain retinal degenerations. Our study also implies that the long-term efficacy of retinal gene therapy may depend on not only the timing of treatment but also other factors such as the function of mutated genes and residual activities of mutant alleles.

Description: Individuals with germline mutations in the gene encoding phosphatase and tensin homolog on chromosome ten (PTEN) are diagnosed with PTEN hamartoma tumor syndrome (PHTS) and are at high risk for developing breast, thyroid and other cancers and/or autoimmunity or neurodevelopmental issues including autism spectrum disorders. Although well recognized as a tumor suppressor, involvement of PTEN mutations in mediating such a diverse range of phenotypes indicates a more central involvement for PTEN in immunity than previously recognized. To address this, sequencing of the T-cell receptor variable-region β-chain was performed on peripheral blood from PHTS patients. Based on patient findings, we performed mechanistic studies in two Pten knock-in murine models, distinct from each other in cell compartment-specific predominance of Pten. We found that PTEN mutations in humans and mice are associated with a skewed T- and B-cell gene repertoire, characterized by increased prevalence of high-frequency clones. Immunological characterization showed that Pten mutants have increased B-cell proliferation and a proclivity towards increased T-cell reactivity upon Toll-like-receptor stimulation. Furthermore, decreases in nuclear but not cytoplasmic Pten levels associated with a reduction in expression of the autoimmune regulator (Aire), a critical mediator of central immune tolerance. Mechanistically, we show that nuclear PTEN most likely regulates Aire expression via its emerging role in splicing regulation. We conclude that germline disruption of PTEN, both in human and mouse, results in compromised central immune tolerance processes that may significantly impact individual stress responses and therefore predisposition to autoimmunity and cancer.

Description: Factor I (FI) is one of the main inhibitors of complement activity, and numerous rare coding variants have been reported in patients with age-related macular degeneration, atypical hemolytic uremic syndrome and C3 glomerulopathy. Since many of these variants are of unknown clinical significance, this study aimed to determine the effect of rare coding variants in the complement factor I (CFI) gene on FI expression. We measured FI levels in plasma samples of carriers of rare coding variants and in vitro in the supernatants of epithelial cells expressing recombinant FI. FI levels were measured in 177 plasma samples of 155 individuals, carrying 24 different rare coding variants in CFI. In carriers of the variants p.Gly119Arg, p.Leu131Arg, p.Gly188Ala and c.772G〉A (r.685_773del), significantly reduced FI plasma levels were detected. Furthermore, recombinant FI expression levels were determined for 126 rare coding variants. Of these variants 68 (54%) resulted in significantly reduced FI expression in supernatant compared to wildtype (WT). The recombinant protein expression levels correlated significantly with the FI level in plasma of carriers of CFI variants. In this study, we performed the most comprehensive FI expression level analysis of rare coding variants in CFI to date. More than half of CFI variants lead to reduced FI expression, which might impair complement regulation in vivo. Our study will aid the interpretation of rare coding CFI variants identified in clinical practice, which is in particular important in light of patient inclusion in ongoing clinical trials for CFI gene supplementation in AMD.

Description: Cyclin-dependent kinase-like 5 (CDKL5), a serine-threonine kinase encoded by an X-linked gene, is highly expressed in the mammalian forebrain. Mutations in this gene cause CDKL5 deficiency disorder, a neurodevelopmental encephalopathy characterized by early-onset seizures, motor dysfunction, and intellectual disability. We previously found that mice lacking CDKL5 exhibit hyperlocomotion and increased impulsivity, resembling the core symptoms in attention-deficit hyperactivity disorder (ADHD). Here, we report the potential neural mechanisms and treatment for hyperlocomotion induced by CDKL5 deficiency. Our results showed that loss of CDKL5 decreases the proportion of phosphorylated dopamine transporter (DAT) in the rostral striatum, leading to increased levels of extracellular dopamine and hyperlocomotion. Administration of methylphenidate (MPH), a DAT inhibitor clinically effective to improve symptoms in ADHD, significantly alleviated the hyperlocomotion phenotype in Cdkl5 null mice. In addition, the improved behavioral effects of MPH were accompanied by a region-specific restoration of phosphorylated dopamine- and cAMP-regulated phosphoprotein Mr 32 kDa, a key signaling protein for striatal motor output. Finally, mice carrying a Cdkl5 deletion selectively in DAT-expressing dopaminergic neurons, but not dopamine receptive neurons, recapitulated the hyperlocomotion phenotype found in Cdkl5 null mice. Our findings suggest that CDKL5 is essential to control locomotor behavior by regulating region-specific dopamine content and phosphorylation of dopamine signaling proteins in the striatum. The direct, as well as indirect, target proteins regulated by CDKL5 may play a key role in movement control and the therapeutic development for hyperactivity disorders.

Description: Individuals with pseudoexfoliation (PEX) syndrome exhibit various connective tissue pathologies associated with dysregulated extracellular matrix homeostasis. PEX glaucoma is a common, aggressive form of open-angle glaucoma resulting from the deposition of fibrillary material in the conventional outflow pathway. However, the molecular mechanisms that drive pathogenesis and genetic risk remain poorly understood. PEX glaucoma-associated single-nucleotide polymorphisms are located in and affect activity of the promoter of LOXL1-AS1, a long non-coding RNA (lncRNA). Nuclear and non-nuclear lncRNAs regulate a host of biological processes, and when dysregulated, contribute to disease. Here we report that LOXL1-AS1 localizes to the nucleus where it selectively binds to the mRNA processing protein, heterogeneous nuclear ribonucleoprotein-L (hnRNPL). Both components of this complex are critical for the regulation of global gene expression in ocular cells, making LOXL1-AS1 a prime target for investigation in PEX syndrome and glaucoma.

Description: Mutations in DEP domain containing 5 (DEPDC5) are increasingly appreciated as one of the most common causes of inherited focal epilepsy. Epilepsies due to DEPDC5 mutations are often associated with brain malformations, tend to be drug-resistant, and have been linked to an increased risk of sudden unexplained death in epilepsy (SUDEP). Generation of epilepsy models to define mechanisms of epileptogenesis remains vital for future therapies. Here, we describe a novel mouse model of Depdc5 deficiency with a severe epilepsy phenotype, generated by conditional deletion of Depdc5 in dorsal telencephalic neuroprogenitor cells. In contrast to control and heterozygous mice, Depdc5-Emx1-Cre conditional knockout (CKO) mice demonstrated macrocephaly, spontaneous seizures and premature death. Consistent with increased mTORC1 activation, targeted neurons were enlarged and both neurons and astrocytes demonstrated increased S6 phosphorylation. Electrophysiologic characterization of miniature inhibitory post-synaptic currents in excitatory neurons was consistent with impaired post-synaptic response to GABAergic input, suggesting a potential mechanism for neuronal hyperexcitability. mTORC1 inhibition with rapamycin significantly improved survival of CKO animals and prevented observed seizures, including for up to 40 days following rapamycin withdrawal. These data not only support a primary role for mTORC1 hyperactivation in epilepsy following homozygous loss of Depdc5, but also suggest a developmental window for treatment which may have a durable benefit for some time even after withdrawal.

Description: Mutations of the RNA granule component TDRD7 (OMIM: 611258) cause pediatric cataract. We applied an integrated approach to uncover the molecular pathology of cataract in Tdrd7−/− mice. Early postnatal Tdrd7−/− animals precipitously develop cataract suggesting a global-level breakdown/misregulation of key cellular processes. High-throughput RNA sequencing integrated with iSyTE-bioinformatics analysis identified the molecular chaperone and cytoskeletal modulator, HSPB1, among high-priority downregulated candidates in Tdrd7−/− lens. A protein fluorescence two-dimensional difference in-gel electrophoresis (2D-DIGE)-coupled mass spectrometry screen also identified HSPB1 downregulation, offering independent support for its importance to Tdrd7−/− cataractogenesis. Lens fiber cells normally undergo nuclear degradation for transparency, posing a challenge: how is their cell morphology, also critical for transparency, controlled post-nuclear degradation? HSPB1 functions in cytoskeletal maintenance, and its reduction in Tdrd7−/− lens precedes cataract, suggesting cytoskeletal defects may contribute to Tdrd7−/− cataract. In agreement, scanning electron microscopy (SEM) revealed abnormal fiber cell morphology in Tdrd7−/− lenses. Further, abnormal phalloidin and wheat germ agglutinin (WGA) staining of Tdrd7−/− fiber cells, particularly those exhibiting nuclear degradation, reveals distinct regulatory mechanisms control F-actin cytoskeletal and/or membrane maintenance in post-organelle degradation maturation stage fiber cells. Indeed, RNA immunoprecipitation identified Hspb1 mRNA in wild-type lens lysate TDRD7-pulldowns, and single-molecule RNA imaging showed co-localization of TDRD7 protein with cytoplasmic Hspb1 mRNA in differentiating fiber cells, suggesting that TDRD7–ribonucleoprotein complexes may be involved in optimal buildup of key factors. Finally, Hspb1 knockdown in Xenopus causes eye/lens defects. Together, these data uncover TDRD7’s novel upstream role in elevation of stress-responsive chaperones for cytoskeletal maintenance in post-nuclear degradation lens fiber cells, perturbation of which causes early-onset cataracts.

Description: In myotonic dystrophy type 1 (DM1), the CUG expansion (CUGexp) in the 3′ untranslated region of the dystrophia myotonica protein kinase messenger ribonucleic acid affects the homeostasis of ribonucleic acid-binding proteins, causing the multiple symptoms of DM1. We have previously reported that Staufen1 is increased in skeletal muscles from DM1 mice and patients and that sustained Staufen1 expression in mature mouse muscle causes a progressive myopathy. Here, we hypothesized that the elevated levels of Staufen1 contributes to the myopathic features of the disease. Interestingly, the classic DM1 mouse model human skeletal actin long repeat (HSALR) lacks overt atrophy while expressing CUGexp transcripts and elevated levels of endogenous Staufen1, suggesting a lower sensitivity to atrophic signaling in this model. We report that further overexpression of Staufen1 in the DM1 mouse model HSALR causes a myopathy via inhibition of protein kinase B signaling through an increase in phosphatase tensin homolog, leading to the expression of atrogenes. Interestingly, we also show that Staufen1 regulates the expression of muscleblind-like splicing regulator 1 and CUG-binding protein elav-like family member 1 in wild-type and DM1 skeletal muscle. Together, data obtained from these new DM1 mouse models provide evidence for the role of Staufen1 as an atrophy-associated gene that impacts progressive muscle wasting in DM1. Accordingly, our findings highlight the potential of Staufen1 as a therapeutic target and biomarker.

Description: The RNA exosome is an essential ribonuclease complex required for processing and/or degradation of both coding and non-coding RNAs. We identified five patients with biallelic variants in EXOSC5, which encodes a structural subunit of the RNA exosome. The clinical features of these patients include failure to thrive, short stature, feeding difficulties, developmental delays that affect motor skills, hypotonia and esotropia. Brain MRI revealed cerebellar hypoplasia and ventriculomegaly. While we ascertained five patients, three patients with distinct variants of EXOSC5 were studied in detail. The first patient had a deletion involving exons 5–6 of EXOSC5 and a missense variant, p.Thr114Ile, that were inherited in trans, the second patient was homozygous for p.Leu206His and the third patient had paternal isodisomy for chromosome 19 and was homozygous for p.Met148Thr. The additional two patients ascertained are siblings who had an early frameshift mutation in EXOSC5 and the p.Thr114Ile missense variant that were inherited in trans. We employed three complementary approaches to explore the requirement for EXOSC5 in brain development and assess consequences of pathogenic EXOSC5 variants. Loss of function for exosc5 in zebrafish results in shortened and curved tails/bodies, reduced eye/head size and edema. We modeled pathogenic EXOSC5 variants in both budding yeast and mammalian cells. Some of these variants cause defects in RNA exosome function as well as altered interactions with other RNA exosome subunits. These findings expand the number of genes encoding RNA exosome subunits linked to human disease while also suggesting that disease mechanism varies depending on the specific pathogenic variant.

Description: Olfactory dysfunction and altered neurogenesis are observed in several neurodegenerative disorders including Huntington disease (HD). These deficits occur early and correlate with a decline in global cognitive performance, depression and structural abnormalities of the olfactory system including the olfactory epithelium, bulb and cortices. However, the role of olfactory system dysfunction in the pathogenesis of HD remains poorly understood and the mechanisms underlying this dysfunction are unknown. We show that deficits in odour identification, discrimination and memory occur in HD individuals. Assessment of the olfactory system in an HD murine model demonstrates structural abnormalities in the olfactory bulb (OB) and piriform cortex, the primary cortical recipient of OB projections. Furthermore, a decrease in piriform neuronal counts and altered expression levels of neuronal nuclei and tyrosine hydroxylase in the OB are observed in the YAC128 HD model. Similar to the human HD condition, olfactory dysfunction is an early phenotype in the YAC128 mice and concurrent with caspase activation in the murine HD OB. These data provide a link between the structural olfactory brain region atrophy and olfactory dysfunction in HD and suggest that cell proliferation and cell death pathways are compromised and may contribute to the olfactory deficits in HD.

Description: Synaptojanin1 (synj1) is a phosphoinositide phosphatase with dual SAC1 and 5′-phosphatase enzymatic activities in regulating phospholipid signaling. The brain-enriched isoform has been shown to participate in synaptic vesicle (SV) recycling. More recently, recessive human mutations were identified in the two phosphatase domains of SYNJ1, including R258Q, R459P and R839C, which are linked to rare forms of early-onset Parkinsonism. We now demonstrate that Synj1 heterozygous deletion (Synj1+/−), which is associated with an impaired 5′-phosphatase activity, also leads to Parkinson’s disease (PD)-like pathologies in mice. We report that male Synj1+/− mice display age-dependent motor function abnormalities as well as alpha-synuclein accumulation, impaired autophagy and dopaminergic terminal degeneration. Synj1+/− mice contain elevated 5′-phosphatase substrate, PI(4,5)P2, particularly in the midbrain neurons. Moreover, pharmacological elevation of membrane PI(4,5)P2 in cultured neurons impairs SV endocytosis, specifically in midbrain neurons, and further exacerbates SV trafficking defects in Synj1+/− midbrain neurons. We demonstrate down-regulation of SYNJ1 transcript in a subset of sporadic PD brains, implicating a potential role of Synj1 deficiency in the decline of dopaminergic function during aging.

Description: We investigated the genetic origin of the phenotype displayed by three children from two unrelated Italian families, presenting with a previously unrecognized autosomal recessive disorder that included a severe form of spondylo-epiphyseal dysplasia, sensorineural hearing loss, intellectual disability and Leber congenital amaurosis (SHILCA), as well as some brain anomalies that were visible at the MRI. Autozygome-based analysis showed that these children shared a 4.76 Mb region of homozygosity on chromosome 1, with an identical haplotype. Nonetheless, whole-exome sequencing failed to identify any shared rare coding variants, in this region or elsewhere. We then determined the transcriptome of patients’ fibroblasts by RNA sequencing, followed by additional whole-genome sequencing experiments. Gene expression analysis revealed a 4-fold downregulation of the gene NMNAT1, residing indeed in the shared autozygous interval. Short- and long-read whole-genome sequencing highlighted a duplication involving 2 out of the 5 exons of NMNAT1 main isoform (NM_022787.3), leading to the production of aberrant mRNAs. Pathogenic variants in NMNAT1 have been previously shown to cause non-syndromic Leber congenital amaurosis (LCA). However, no patient with null biallelic mutations has ever been described, and murine Nmnat1 knockouts show embryonic lethality, indicating that complete absence of NMNAT1 activity is probably not compatible with life. The rearrangement found in our cases, presumably causing a strong but not complete reduction of enzymatic activity, may therefore result in an intermediate syndromic phenotype with respect to LCA and lethality.

Description: Statistical imputation applied to genome-wide array data is the most cost-effective approach to complete the catalog of genetic variation in a study population. However, imputed genotypes in underrepresented populations incur greater inaccuracies due to ascertainment bias and a lack of representation among reference individuals, further contributing to the obstacles to study these populations. Here we examined the consequences due to the lack of representation by genotyping in a large number of self-reported Native Hawaiians (N = 3693) a functionally important, Polynesian-specific variant in the CREBRF gene, rs373863828. We found the derived allele was significantly associated with several adiposity traits with large effects (e.g. ~ 1.28 kg/m2 per allele in body mass index as the most significant; P = 7.5 × 10−5), consistent with the original findings in Samoans. Due to the current absence of Polynesian representation in publicly accessible reference sequences, rs373863828 or its proxies could not be tested through imputation using these existing resources. Moreover, the association signals at the entire CREBRF locus could not be captured by alternative approaches, such as admixture mapping. In contrast, highly accurate imputation can be achieved even if a small number (

Description: Lysinuric protein intolerance (LPI) is an inborn error of cationic amino acid (arginine, lysine, ornithine) transport caused by biallelic pathogenic variants in SLC7A7, which encodes the light subunit of the y+LAT1 transporter. Treatments for the complications of LPI, including growth failure, renal disease, pulmonary alveolar proteinosis, autoimmune disorders and osteoporosis, are limited. Given the early lethality of the only published global Slc7a7 knockout mouse model, a viable animal model to investigate global SLC7A7 deficiency is needed. Hence, we generated two mouse models with global Slc7a7 deficiency (Slc7a7em1Lbu/em1Lbu; Slc7a7Lbu/Lbu and Slc7a7em1(IMPC)Bay/em1(IMPC)Bay; Slc7a7Bay/Bay) using CRISPR/Cas9 technology by introducing a deletion of exons 3 and 4. Perinatal lethality was observed in Slc7a7Lbu/Lbu and Slc7a7Bay/Bay mice on the C57BL/6 and C57BL/6NJ inbred genetic backgrounds, respectively. We noted improved survival of Slc7a7Lbu/Lbu mice on the 129 Sv/Ev × C57BL/6 F2 background, but postnatal growth failure occurred. Consistent with human LPI, these Slc7a7Lbu/Lbu mice exhibited reduced plasma and increased urinary concentrations of the cationic amino acids. Histopathological assessment revealed loss of brush border and lipid vacuolation in the renal cortex of Slc7a7Lbu/Lbu mice, which combined with aminoaciduria suggests proximal tubular dysfunction. Micro-computed tomography of L4 vertebrae and skeletal radiographs showed delayed skeletal development and suggested decreased mineralization in Slc7a7Lbu/Lbu mice, respectively. In addition to delayed skeletal development and delayed development in the kidneys, the lungs and liver were observed based on histopathological assessment. Overall, our Slc7a7Lbu/Lbu mouse model on the F2 mixed background recapitulates multiple human LPI phenotypes and may be useful for future studies of LPI pathology.

Description: Rare coding variants have been proven to be one of the significant factors contributing to spermatogenic failure in patients with non-obstructive azoospermia (NOA) and severe oligospermia (SO). To delineate the molecular characteristics of idiopathic NOA and SO, we performed whole-exome sequencing of 314 unrelated patients of Chinese Han origin and verified our findings by comparing to 400 fertile controls. We detected six pathogenic/likely pathogenic variants and four variants of unknown significance, in genes known to cause NOA/SO, and 9 of which had not been earlier reported. Additionally, we identified 20 novel NOA candidate genes affecting 25 patients. Among them, five (BRDT, CHD5, MCM9, MLH3 and ZFX) were considered as strong candidates based on the evidence obtained from murine functional studies and human single-cell (sc)RNA-sequencing data. These genetic findings provide insight into the aetiology of human NOA/SO and pave the way for further functional analysis and molecular diagnosis of male infertility.

Description: Dysfunction of the gonadotropin-releasing hormone (GnRH) axis causes a range of reproductive phenotypes resulting from defects in the specification, migration and/or function of GnRH neurons. To identify additional molecular components of this system, we initiated a systematic genetic interrogation of families with isolated GnRH deficiency (IGD). Here, we report 13 families (12 autosomal dominant and one autosomal recessive) with an anosmic form of IGD (Kallmann syndrome) with loss-of-function mutations in TCF12, a locus also known to cause syndromic and non-syndromic craniosynostosis. We show that loss of tcf12 in zebrafish larvae perturbs GnRH neuronal patterning with concomitant attenuation of the orthologous expression of tcf3a/b, encoding a binding partner of TCF12, and stub1, a gene that is both mutated in other syndromic forms of IGD and maps to a TCF12 affinity network. Finally, we report that restored STUB1 mRNA rescues loss of tcf12 in vivo. Our data extend the mutational landscape of IGD, highlight the genetic links between craniofacial patterning and GnRH dysfunction and begin to assemble the functional network that regulates the development of the GnRH axis.

Description: It has been reported that abnormal epigenetic modification is associated with the occurrence of Parkinson’s disease (PD). Here, we found that a ten-eleven translocation 2 (TET2), a staff of the DNA hydroxylases family, was increased in dopaminergic neurons in vitro and in vivo. Genome-wide mapping of DNA 5-hydroxymethylcytosine (5-hmC)-sequencing has revealed an aberrant epigenome 5-hmC landscape in 1-methyl-4-phenylpyridinium iodide (MPP+)-induced SH-SY5Y cells. The TET family of DNA hydroxylases could reverse DNA methylation by oxidization of 5-methylcytosine (5-mC) to 5-hmC. However, the relationship between modification of DNA hydroxymethylation and the pathogenesis of PD is not clear. According to the results of 5-hmC-sequencing studies, 5-hmC was associated with gene-rich regions in the genomes related to cell cycle, especially gene-cyclin-dependent kinase inhibitor 2A (Cdkn2A). Downregulation of TET2 expression could significantly rescue MPP+-stimulated SH-SY5Y cell damage and cell cycle arrest. Meanwhile, knockdown of Tet2 expression in the substantia nigra pars compacta of MPTP-induced PD mice resulted in attenuated MPTP-induced motor deficits and dopaminergic neuronal injury via p16 suppression. In this study, we demonstrated a critical function of TET2 in PD development via the CDKN2A activity-dependent epigenetic pathway, suggesting a potential new strategy for epigenetic therapy.

Description: Spinal muscular atrophy (SMA) occurs as a result of cell-ubiquitous depletion of the essential survival motor neuron (SMN) protein. Characteristic disease pathology is driven by a particular vulnerability of the ventral motor neurons of the spinal cord to decreased SMN. Perhaps not surprisingly, many other organ systems are also impacted by SMN depletion. The normal kidney expresses very high levels of SMN protein, equivalent to those found in the nervous system and liver, and levels are dramatically lowered by ~90–95% in mouse models of SMA. Taken together, these data suggest that renal pathology may be present in SMA. We have addressed this using an established mouse model of severe SMA. Nephron number, as assessed by gold standard stereological techniques, was significantly reduced. In addition, morphological assessment showed decreased renal vasculature, particularly of the glomerular capillary knot, dysregulation of nephrin and collagen IV, and ultrastructural changes in the trilaminar filtration layers of the nephron. To explore the molecular drivers underpinning this process, we correlated these findings with quantitative PCR measurements and protein analyses of glial cell-line-derived neurotrophic factor, a crucial factor in ureteric bud branching and subsequent nephron development. Glial cell-line-derived neurotrophic factor levels were significantly reduced at early stages of disease in SMA mice. Collectively, these findings reveal significant renal pathology in a mouse model of severe SMA, further reinforcing the need to develop and administer systemic therapies for this neuromuscular disease.

Description: Glycosuria is a condition where glucose is detected in urine at higher concentrations than normal (i.e. not detectable). Glycosuria at some point during pregnancy has an estimated prevalence of 50% and is associated with adverse outcomes in both mothers and offspring. Little is currently known about the genetic contribution to this trait or the extent to which it overlaps with other seemingly related traits, e.g. diabetes. We performed a genome-wide association study (GWAS) for self-reported glycosuria in pregnant mothers from the Avon Longitudinal Study of Parents and Children (cases/controls = 1249/5140). We identified two loci, one of which (lead SNP = rs13337037; chromosome 16; odds ratio of glycosuria per effect allele: 1.42; 95% CI: 1.30, 1.56; P = 1.97 × 10−13) was then validated using an obstetric measure of glycosuria measured in the same cohort (227/6639). We performed a secondary GWAS in the 1986 Northern Finland Birth Cohort (NFBC1986; 747/2991) using midwife-reported glycosuria and offspring genotype as a proxy for maternal genotype. The combined results revealed evidence for a consistent effect on glycosuria at the chromosome 16 locus. In follow-up analyses, we saw little evidence of shared genetic underpinnings with the exception of urinary albumin-to-creatinine ratio (Rg = 0.64; SE = 0.22; P = 0.0042), a biomarker of kidney disease. In conclusion, we identified a genetic association with self-reported glycosuria during pregnancy, with the lead SNP located 15kB upstream of SLC5A2, a target of antidiabetic drugs. The lack of strong genetic correlation with seemingly related traits such as type 2 diabetes suggests different genetic risk factors exist for glycosuria during pregnancy.

Description: Mutations in S-phase cyclin A-associated protein in the endoplasmic reticulum (SCAPER) cause a recessively inherited multisystemic disorder whose main features are retinal degeneration and intellectual disability. SCAPER, originally identified as a cell cycle regulator, was also suggested to be a ciliary protein. Because Scaper mutant males are sterile, we set up to characterize their phenotype. The testes of Scaper mutant mice are significantly smaller than those of WT mice. Histology revealed no signs of spermatogenesis, and seminiferous tubules contained mainly Sertoli cells with a few spermatogonia/spermatogonial stem cells (SSCs). In WT testes, SCAPER is expressed by SSCs and in the various stages of spermatogenesis, as well as in Sertoli cells. In WT spermatozoa SCAPER is not expressed in the flagellum but rather in the head compartment, where it is found both in the nucleus and in the perinuclear region. Scaper mutant females present reduced fertility, manifested by a significantly smaller litter size compared to WT females. Mutant ovaries are similar in size but comprised of significantly less primordial and antral follicles, compared to WT ovaries, while the number of atretic follicles is significantly higher. In WT ovarian follicles SCAPER is expressed in the somatic granulosa cells as well as in the oocyte. In conclusion, our data demonstrate that SCAPER is a crucial component in both male and female reproductive systems. We hypothesize that the reproductive phenotype observed in Scaper mutant mice is rooted in SCAPER’s interaction with cyclin A/Cdk2, which play an important role, however different, in male and female gonads.

Description: In mice, male sex determination depends on FGF9 signalling via FGFR2c in the bipotential gonads to maintain the expression of the key testis gene SOX9. In humans, however, while FGFR2 mutations have been linked to 46,XY disorders of sex development (DSD), the role of FGF9 is unresolved. The only reported pathogenic mutations in human FGF9, FGF9S99N and FGF9R62G, are dominant and result in craniosynostosis (fusion of cranial sutures) or multiple synostoses (fusion of limb joints). Whether these synostosis-causing FGF9 mutations impact upon gonadal development and DSD etiology has not been explored. We therefore examined embryonic gonads in the well-characterized Fgf9 missense mouse mutants, Fgf9S99N and Fgf9N143T, which phenocopy the skeletal defects of FGF9S99N and FGF9R62G variants, respectively. XY Fgf9S99N/S99N and XY Fgf9N143T/N143T fetal mouse gonads showed severely disorganized testis cords and partial XY sex reversal at 12.5 days post coitum (dpc), suggesting loss of FGF9 function. By 15.5 dpc, testis development in both mutants had partly recovered. Mitotic analysis in vivo and in vitro suggested that the testicular phenotypes in these mutants arise in part through reduced proliferation of the gonadal supporting cells. These data raise the possibility that human FGF9 mutations causative for dominant skeletal conditions can also lead to loss of FGF9 function in the developing testis, at least in mice. Our data suggest that, in humans, testis development is largely tolerant of deleterious FGF9 mutations which lead to skeletal defects, thus offering an explanation as to why XY DSDs are rare in patients with pathogenic FGF9 variants.

Description: Alzheimer’s disease (AD) is the most common form of dementia and the most prevalent neurodegenerative disease. Genome-wide association studies have linked PICALM to AD risk. PICALM has been implicated in Aβ42 production and turnover, but whether it plays a direct role in modulating Aβ42 toxicity remains unclear. We found that increased expression of the Drosophila PICALM orthologue lap could rescue Aβ42 toxicity in an adult-onset model of AD, without affecting Aβ42 level. Imbalances in the glutamatergic system, leading to excessive, toxic stimulation, have been associated with AD. We found that Aβ42 caused the accumulation of presynaptic vesicular glutamate transporter (VGlut) and increased spontaneous glutamate release. Increased lap expression reversed these phenotypes back to control levels, suggesting that lap may modulate glutamatergic transmission. We also found that lap modulated the localization of amphiphysin (Amph), the homologue of another AD risk factor BIN1, and that Amph itself modulated postsynaptic glutamate receptor (GluRII) localization. We propose a model where PICALM modulates glutamatergic transmission, together with BIN1, to ameliorate synaptic dysfunction and disease progression.

Description: Aggregates of Aβ peptide and the microtubule-associated protein tau are key molecular hallmarks of Alzheimer’s disease (AD). However, the interaction between these two pathologies and the mechanisms underlying disease progression have remained unclear. Numerous failed clinical trials suggest the necessity for greater mechanistic understanding in order to refine strategies for therapeutic discovery and development. To this end, we have generated a transgenic Caenorhabditis elegans model expressing both human Aβ1-42 peptide and human tau protein pan-neuronally. We observed exacerbated behavioral dysfunction and age-dependent neurodegenerative changes in the Aβ;tau transgenic animals. Further, these changes occurred in the Aβ;tau transgenic animals at greater levels than worms harboring either the Aβ1-42 or tau transgene alone and interestingly without changes to the levels of tau expression, phosphorylation or aggregation. Functional changes were partially rescued with the introduction of a genetic suppressor of tau pathology. Taken together, the data herein support a synergistic role for both Aβ and tau in driving neuronal dysfunction seen in AD. Additionally, we believe that the utilization of the genetically tractable C. elegans model will provide a key resource for dissecting mechanisms driving AD molecular pathology.

Description: The theoretical conditions for small-angle inelastic scattering where the incident electron can effectively be treated as a particle moving in a uniform potential is examined. The motivation for this work is the recent development of a multislice method that combines plasmon energy losses with elastic scattering using Monte Carlo methods. Since plasmon excitation is delocalized, it was assumed that the Bloch wave nature of the incident electron in the crystal does not affect the scattering cross-section. It is shown here that for a delocalized excitation the mixed dynamic form factor term of the scattering cross-section is zero and the scattered intensities follow a Poisson distribution. These features are characteristic of particle-like scattering and validate the use of Monte Carlo methods to model plasmon losses in multislice simulations.

Description: This review examines methods of magnetic flux density measurements from the narrow grain boundary (GB) regions, the thickness of which is of the order of nanometers, produced in Nd–Fe–B-based sintered magnets. Despite of the complex crystallographic microstructure and the significant stray magnetic field of the sintered magnet, recent progress in electron holography allowed for the determination of the intrinsic magnetic flux density due to the GB which is embedded in the polycrystalline thin-foil. The methods appear to be useful as well for intensive studies about interface magnetism in a variety of systems.

Description: Some species of lactic acid bacteria used for the production of natural cheese produce exopolysaccharides (EPS). Electron microscopy is useful for analyzing the microstructure of EPS produced by lactic acid bacteria. However, pretreatments used to observe the microstructure of EPS by electron microscopy, such as dehydration and resin embedding, can result in EPS flowing out easily from the cell. Therefore, in this study, the Tokuyasu method was conducted on cryosection to reduce EPS outflow. Two types of observation method, namely, using lectin and ruthenium red, were conducted in an attempt to observe EPS produced by Lactobacillus helveticus SBT2171. Observation using the lectin method confirmed that colloidal gold particles conjugated with a lectin recognizing β-galactoside were present in the capsule. Structures that appeared to be β-galactoside-containing slime polysaccharides that were released from the cell wall were also observed. Observation using ruthenium red showed that capsular polysaccharides (CPS) in the capsule were present as a net-like structure. Colloidal gold conjugation with an anti-β-lactoglobulin antibody, in addition to ruthenium red staining, allowed the identification of slime polysaccharides released from the cell wall in the milk protein network derived from the culture medium. Based on these results, the Tokuyasu method was considered to be a useful pretreatment method to clarify and observe the presence of EPS. In particular, both CPS in the capsule and slime exopolysaccharides released from the cell wall were visualized.

Description: Transmission electron microscopy (TEM) and scanning transmission electron microscopy (STEM) enable the visualization of three-dimensional (3D) microstructures ranging from atomic to micrometer scales using 3D reconstruction techniques based on computed tomography algorithms. This 3D microscopy method is called electron tomography (ET) and has been utilized in the fields of materials science and engineering for more than two decades. Although atomic resolution is one of the current topics in ET research, the development and deployment of intermediate-resolution (non-atomic-resolution) ET imaging methods have garnered considerable attention from researchers. This research trend is probably not irrelevant due to the fact that the spatial resolution and functionality of 3D imaging methods of scanning electron microscopy (SEM) and X-ray microscopy have come to overlap with those of ET. In other words, there may be multiple ways to carry out 3D visualization using different microscopy methods for nanometer-scale objects in materials. From the above standpoint, this review paper aims to (i) describe the current status and issues of intermediate-resolution ET with regard to enhancing the effectiveness of TEM/STEM imaging and (ii) discuss promising applications of state-of-the-art intermediate-resolution ET for materials research with a particular focus on diffraction contrast ET for crystalline microstructures (superlattice domains and dislocations) including a demonstration of in situ dislocation tomography.

Description: Image contrast between carbon nanotubes (CNTs) and polytetrafluoroethylene (PTFE) in a CNT/PTFE composite film, which is difficult to obtain by conventional backscattered electron (BSE) imaging, was optimized to better elucidate the distribution of CNT in the film. Ultra-low landing energy condition (0.3 keV in this study) was used to prevent specimen damage due to electron beam irradiation. Signal acceptance maps, which represent the distributions of energy and take-off angle, were calculated to evaluate the features of the signal detection system used in this study. SEM images of this composite film were taken under several sets of conditions and analyzed using these acceptance maps. CNT and PTFE in the composite film can be clearly distinguished with material and topographic contrasts using the BSE signal under optimized energy and take-off angle ranges, even at ultra-low landing energy conditions.

Description: We have developed a new certified reference material (CRM) for image sharpness evaluation and magnification calibration for scanning electron microscopy (SEM). Designed to be suitable for the image sharpness evaluation by the derivative method, the CRM has nanoscale tungsten dot-array structure fabricated on silicon substrate, which gives steep contrast transition from the dot to the substrate in SEM image. The pitch of the dot-array was SI-traceably measured as a specified value with relative expanded uncertainty ($k=2$) of ~1.3%, which can be utilized for the magnification calibration of SEM. Since specimens, as one of the image formation parameters, easily affect image sharpness value, our CRM, as a ‘pinned specimen’, will play an important role to achieve robust and stable image sharpness measurement system.

Description: Here, we describe a method for tracking intracellular processing of small interfering RNA (siRNA) containing complexes using automated microscopy controls and image acquisition to minimize user effort and time. This technique uses fluorescence colocalization to monitor dual-labeled fluorescent siRNAs delivered by silica nanoparticles in different intracellular locations, including the early/late endosomes, fast/slow recycling endosomes, lysosomes and the endoplasmic reticulum. Combining the temporal association of siRNAs with each intracellular location, we reconstructed the intracellular pathways used in siRNA processing, and demonstrate how these pathways vary based on the chemical composition of the delivery vehicle.

Description: Candida auris is a multiresistant pathogenic yeast commonly isolated from bloodstream infections in immunocompromised patients. In this work, we infected G. mellonella larvae with 105 CFU of a reference strains and two clinical isolates of C. albicans and C. auris and we compared the outcomes of infection between both species. Larvae were evaluated every 24 h for a total of 120 h following the G. mellonella Health Index Scoring System, and survival, activity, melanization and cocoon formation were monitored. Our results showed that clinical isolates were significantly more pathogenic than reference strains independently of the tested species, producing lower survival and activity scores and higher melanization scores, and being C. albicans strains more virulent than C. auris strains. We did not find differences in mortality between aggregative and non-aggregative C. auris strains, although non-aggregative strains produced significantly lower activity scores and higher melanization scores than aggregative ones. Survival assays using Galleria mellonella have been previously employed to examine and classify strains of this and other microbial species based on their virulence before scaling the experiments to a mammal model. Taken together, these results show how a more complete evaluation of the model can improve the study of C. auris isolates.

Description: Cryptococcosis is the second most common invasive fungal infection reported in renal transplant recipients. Tissue granulomatous inflammation is necessary to contain Cryptococcus infection. This study aims to analyze the granuloma patterns and in situ expression of regulatory T (Treg) immune response in tissue samples from 12 renal transplant recipients with cryptococcosis. Fungal isolates were molecularly identified as Cryptococcus neoformans species complex. A detailed characterization of granulomas in tissue samples from 12 kidney transplant recipients with cryptococcosis was described by checking six lung and six skin biopsies by conventional histology and for immunohistochemical detection of CD4 and Treg markers: forkhead box P3 (FoxP3), interleukin (IL)-10 and transforming-growth factor (TGF)-β. Granulomas were classified as compact, loose or mixed. Patients with mixed (n = 4) and compact (n = 3) granulomatous inflammation patterns were associated with a better prognosis and presented a higher number of CD4+FoxP3+T cells compared to the group of patients with loose granulomas. In counterpart, three out of five patients with loose granulomas died with cryptococcosis. We suggest that Treg may have a protective role in the tissue response to Cryptococcus infection given its association with compact and mixed granulomas in patients with better clinical outcomes.

Description: The RNA chaperone Hfq regulates diverse processes in numerous bacteria. In this study, we compared phenotypes (growth rate, adherence, response to different stress conditions, and virulence in Galleria mellonella) of wild-type (WT) and isogenic hfq mutants of three serovars (1, 8 and 15) of the porcine pathogen A. pleuropneumoniae. Similar growth in rich broth was seen for all strains except Ap1∆hfq, which showed slightly reduced growth throughout the 24 hour time course, and the complemented Ap8∆hfqC mutant had a prolonged lag phase. Differences were seen between the three serovar WT strains regarding adherence, stress response and virulence in G. mellonella, and deletion of hfq affected some, but not all of these phenotypes, depending on serovar. Complementation by expression of cloned hfq from an endogenous promoter only restored some WT phenotypes, indicating that complex regulatory networks may be involved, and that levels of Hfq may be as important as presence/absence of the protein regarding its contribution to gene regulation. Our results support that Hfq is a pleiotropic global regulator in A. pleuropneumoniae, but serovar-related differences exist. These results highlight the importance of testing multiple strains/serovars within a given species when determining contributions of global regulators, such as Hfq, to expression of complex phenotypes.

Description: Multi-drug resistant (MDR) Klebsiella pneumoniae represent a global threat to healthcare due to lack of effective treatments and high mortality rates. The aim of this research was to explore the potential of administering zidovudine (AZT) in combination with an existing antibiotic to treat resistant K. pneumoniae infections. Two MDR K. pneumoniae strains were employed, producing either the NDM-1 or KPC-3 carbapenemase. Efficacy of combinations of AZT with meropenem were compared with monotherapies against infections in Galleria mellonella larvae by measuring larval mortality and bacterial burden. The effect of the same combinations in vitro was determined via checkerboard and time-kill assays. In vitro, both K. pneumoniae strains were resistant to meropenem but were susceptible to AZT. In G. mellonella, treatment with either AZT or meropenem alone offered minimal therapeutic benefit against infections with either strain. In contrast, combination therapy of AZT with meropenem presented significantly enhanced efficacy compared to monotherapies. This was correlated with prevention of bacterial proliferation within the larvae but not elimination. Checkerboard assays showed that the interaction between AZT and meropenem was not synergistic but indifferent. In summary, combination therapy of AZT with meropenem represents a potential treatment for carbapenemase-producing MDR K. pneumoniae and merits further investigation.

Description: The immune response of insects displays many structural and functional similarities to the innate immune response of mammals. As a result of these conserved features, insects may be used for evaluating microbial virulence or for testing the in vivo efficacy and toxicity of antimicrobial compounds and results show strong similarities to those from mammals. Galleria mellonella larvae are widely used in this capacity and have the advantage of being easy to use, inexpensive to purchase and house, and being free from the ethical and legal restrictions that relate to the use of mammals in these tests. G. mellonella larvae may be used to assess the in vivo toxicity and efficacy of novel antimicrobial compounds. A wide range of antibacterial and antifungal therapies have been evaluated in G. mellonella larvae and results have informed subsequent experiments in mammals. While insect larvae are a convenient and reproducible model to use, care must be taken in their use to ensure accuracy of results. The objective of this review is to provide a comprehensive account of the use of G. mellonella larvae for assessing the in vivo toxicity and efficacy of a wide range of antibacterial and antifungal agents.

Description: Persistent infections with the bacterial group-I carcinogen Helicobacter pylori (H. pylori) have been associated with a broad range of gastric disorders, including gastritis, ulceration, gastric cancer or mucosa-associated lymphoid tissue (MALT) lymphoma. Pathogenesis of H. pylori requires a balance between immune tolerance and defense. Although H. pylori induces inflammatory responses, the immune system cannot eliminate the pathogen. The detailed molecular mechanisms of how H. pylori interferes with cells of the immune system, in particular infiltrated B cells, are not well investigated. Previously, it was shown that the bacterial effector and oncoprotein cytotoxin-associated gene A (CagA) is delivered into B cells followed by its tyrosine-phosphorylation. To investigate the functional consequences in B cells colonized by CagA-positive H. pylori, we analyzed the global transcriptome of H. pylori-infected Mec-1 cells by RNA sequencing. We found 889 differentially expressed genes (DEGs) and validated JUN, FOSL2, HSPA1B, SRC, CXCR3, TLR-4, TNF-α, CXCL8, CCL2, CCL4, MHC class I and MHC class II molecules by qPCR, western blot, flow cytometry and ELISA assays. The H. pylori-specific mRNA expression signature reveals a downregulation of inflammation- and migration-associated genes, whereas central signal transduction regulators of cell survival and death are upregulated.

Description: The larva of the greater wax moth Galleria mellonella is an increasingly popular model for assessing the virulence of bacterial pathogens, and the effectiveness of antimicrobial agents. In this review, we discuss details of the components of the G. mellonella larval immune system that underpin its use as an alternative infection model, and provide an updated overview of the state of the art of research with G. mellonella- infection models to study bacterial virulence, and in the evaluation of antimicrobial efficacy. Emphasis is given to virulence studies with relevant human and veterinary pathogens, especially Escherichia coli and bacteria of the ESKAPE group. In addition, we make practical recommendations for larval rearing and testing, and overcoming potential limitations of the use of the model, that facilitate intra- and inter-laboratory reproducibility.

Description: Although numerous studies have clarified the synergistic pathogenesis in mouse models of influenza A virus (IAV)-associated dual infections, fewer studies have investigated the influence of intranasal liquid administration on the disease. This study explored the effects of intranasal PBS administration in mouse models of mimic IAV dual infection and the infectious dose of IAV that caused equivalent pathogenesis in different dual infection models. Weights, survival rates, virus loads, lung indexes and lung pathology were compared. We demonstrated that intranasal PBS administration following H1N1 or H3N2 infection increased weight loss, mortality, virus replication and lung damage. No difference was observed if the order was reversed or PBS was given simultaneously with IAV. To induce equivalent virulence, a 20-fold difference in the infectious dose was needed when the H3N2–PBS superinfection and H3N2–PBS coinfection or PBS–H3N2 superinfection groups were compared. Our study demonstrated that the unfavourable effect of intranasal liquid administration should not be neglected and that both the strain and infectious dose of IAV should be considered to avoid an illusion of synergistic pathogenicity when establishing IAV-associated dual infection model. A 20-fold lower dose than that of coinfection may be a better choice for secondary infection following IAV.

Description: Background Pseudomonas aeruginosa possesses an array of virulence genes ensuring successful infection development. A two-partner secretion system Exolysin BA (ExlBA) is expressed in the PA7-like genetic outliers consisting of ExlA, a pore-forming toxin and ExlB transporter protein. Presence of exlBA in multidrug-resistant (MDR) strains has not been investigated, particularly in the strains isolated from wounded soldiers. Methods We screened whole genome sequences of 2439 MDR- P. aeruginosa strains for the presence of exlBA. We compiled all exlBA positive strains and compared them with a diversity set for demographics, antimicrobial profiles and phenotypic characteristics: surface motility, biofilm formation, pyocyanin production and hemolysis. We compared the virulence of strains with comparable phenotypic characteristics in Galleria mellonella. Results We identified 33 exlBA-positive strains (1.5%). These strains have increased antibiotic resistance, they are more motile, produce more robust biofilms and have comparable pyocianin production with the diversity set despite the phenotypic differences within the group. In in vivo infection models, these strains were less virulent than Type III Secretion System (T3SS) positive counterparts. Conclusions exlBA-positive strains are wide spread among the PA7-like outliers. While not as virulent as strains possessing T3SS, these strains exhibit phenotypic features associated with virulence and are still lethal in vivo.

Description: Intestinal mucosal barriers help the body resist many intestinal inflammatory diseases, such as inflammatory bowel disease (IBD). In this study, we identified a novel bacterium promoting the repair of intestinal mucosa and investigated the potential mechanisms underlying its activity. Culture supernatant of Bacillus subtilis RZ001 upregulated the expression of mucin 2 (MUC2) and tight junction (TJ) proteins in HT-29 cells in vitro. Oral administration of B. subtilis RZ001 may have significantly reduced symptoms such as the dextran sulfate sodium (DSS)-induced decrease in body weight, shortening of colon length and overproduction of proinflammatory factors. The number of goblet cells and levels of MUC2 and TJ proteins were significantly increased in adult mice fed with B. subtilis RZ001. B. subtilis RZ001 cells upregulated the levels of MUC2 in the intestinal organoids. Furthermore, culture supernatant of B. subtilis RZ001 could suppress the Notch signalling pathway and activate the expression of atonal homolog 1 (Atoh1). The transcription factor Atoh1 is required for intestinal secretory cell differentiation and activates transcription of MUC2 via binding to E-boxes on the MUC2 promoter. Taken together, B. subtilis strain RZ001 has the potential for treating IBD. The present study is helpful to elucidate the mechanisms of B. subtilis action.

Description: Motivation Recent technological advances produce a wealth of high-dimensional descriptions of biological processes, yet extracting meaningful insight and mechanistic understanding from these data remains challenging. For example, in developmental biology, the dynamics of differentiation can now be mapped quantitatively using single-cell RNA sequencing, yet it is difficult to infer molecular regulators of developmental transitions. Here, we show that discovering informative features in the data is crucial for statistical analysis as well as making experimental predictions. Results We identify features based on their ability to discriminate between clusters of the data points. We define a class of problems in which linear separability of clusters is hidden in a low-dimensional space. We propose an unsupervised method to identify the subset of features that define a low-dimensional subspace in which clustering can be conducted. This is achieved by averaging over discriminators trained on an ensemble of proposed cluster configurations. We then apply our method to single-cell RNA-seq data from mouse gastrulation, and identify 27 key transcription factors (out of 409 total), 18 of which are known to define cell states through their expression levels. In this inferred subspace, we find clear signatures of known cell types that eluded classification prior to discovery of the correct low-dimensional subspace. Availability and implementation https://github.com/smelton/SMD. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation Polymerase chain reaction (PCR) has been a revolutionary biomedical advancement. However, for PCR to be appropriately used, one must spend a significant amount of effort on PCR primer design. Carefully designed PCR primers not only increase sensitivity and specificity, but also decrease effort spent on experimental optimization. Computer software removes the human element by performing and automating the complex and rigorous calculations required in PCR primer design. Classification and review of the available software options and their capabilities should be a valuable resource for any PCR application. Results This paper focuses on currently available free PCR primer design software and their major functions (https://pcrprimerdesign.github.io/). The software are classified according to their PCR applications, such as Sanger sequencing, reverse transcription quantitative PCR, single nucleotide polymorphism detection, splicing variant detection, methylation detection, microsatellite detection, multiplex PCR and targeted next generation sequencing, and conserved/degenerate primers to clone orthologous genes from related species, new gene family members in the same species, or to detect a group of related pathogens. Each software is summarized to provide a technical review of their capabilities and utilities.

Description: Motivation Network diffusion and label propagation are fundamental tools in computational biology, with applications like gene-disease association, protein function prediction and module discovery. More recently, several publications have introduced a permutation analysis after the propagation process, due to concerns that network topology can bias diffusion scores. This opens the question of the statistical properties and the presence of bias of such diffusion processes in each of its applications. In this work, we characterised some common null models behind the permutation analysis and the statistical properties of the diffusion scores. We benchmarked seven diffusion scores on three case studies: synthetic signals on a yeast interactome, simulated differential gene expression on a protein-protein interaction network and prospective gene set prediction on another interaction network. For clarity, all the datasets were based on binary labels, but we also present theoretical results for quantitative labels. Results Diffusion scores starting from binary labels were affected by the label codification, and exhibited a problem-dependent topological bias that could be removed by the statistical normalisation. Parametric and non-parametric normalisation addressed both points by being codification-independent and by equalising the bias. We identified and quantified two sources of bias -mean value and variance- that yielded performance differences when normalising the scores. We provided closed formulae for both and showed how the null covariance is related to the spectral properties of the graph. Despite none of the proposed scores systematically outperformed the others, normalisation was preferred when the sought positive labels were not aligned with the bias. We conclude that the decision on bias removal should be problem and data-driven, i.e. based on a quantitative analysis of the bias and its relation to the positive entities. Availability The code is publicly available at https://github.com/b2slab/diffuBench Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation MicroRNAs (miRNAs) are a class of non-coding RNAs that play critical roles in various biological processes. Many studies have shown that miRNAs are closely related to the occurrence, development and diagnosis of human diseases. Traditional biological experiments are costly and time consuming. As a result, effective computational models have become increasingly popular for predicting associations between miRNAs and diseases, which could effectively boost human disease diagnosis and prevention. Results We propose a novel computational framework, called AEMDA, to identify associations between miRNAs and diseases. AEMDA applies a learning-based method to extract dense and high-dimensional representations of diseases and miRNAs from integrated disease semantic similarity, miRNA functional similarity and heterogeneous related interaction data. In addition, AEMDA adopts a deep autoencoder that does not need negative samples to retrieve the underlying associations between miRNAs and diseases. Furthermore, the reconstruction error is used as a measurement to predict disease-associated miRNAs. Our experimental results indicate that AEMDA can effectively predict disease-related miRNAs and outperforms state-of-the-art methods. Availability and implementation The source code and data are available at https://github.com/CunmeiJi/AEMDA. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation Different from traditional linear RNAs (containing 5′ and 3′ ends), circular RNAs (circRNAs) are a special type of RNAs that have a closed ring structure. Accumulating evidence has indicated that circRNAs can directly bind proteins and participate in a myriad of different biological processes. Results For identifying the interaction of circRNAs with 37 different types of circRNA-binding proteins (RBPs), we develop an ensemble neural network, termed PASSION, which is based on the concatenated artificial neural network (ANN) and hybrid deep neural network frameworks. Specifically, the input of the ANN is the optimal feature subset for each RBP, which has been selected from six types of feature encoding schemes through incremental feature selection and application of the XGBoost algorithm. In turn, the input of the hybrid deep neural network is a stacked codon-based scheme. Benchmarking experiments indicate that the ensemble neural network reaches the average best area under the curve (AUC) of 0.883 across the 37 circRNA datasets when compared with XGBoost, k-nearest neighbor, support vector machine, random forest, logistic regression and Naive Bayes. Moreover, each of the 37 RBP models is extensively tested by performing independent tests, with the varying sequence similarity thresholds of 0.8, 0.7, 0.6 and 0.5, respectively. The corresponding average AUC obtained are 0.883, 0.876, 0.868 and 0.883, respectively, highlighting the effectiveness and robustness of PASSION. Extensive benchmarking experiments demonstrate that PASSION achieves a competitive performance for identifying binding sites between circRNA and RBPs, when compared with several state-of-the-art methods. Availability and implementation A user-friendly web server of PASSION is publicly accessible at http://flagship.erc.monash.edu/PASSION/. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Summary The COVID-19 crisis has elicited a global response by the scientific community that has led to a burst of publications on the pathophysiology of the virus. However, without coordinated efforts to organize this knowledge, it can remain hidden away from individual research groups. By extracting and formalizing this knowledge in a structured and computable form, as in the form of a knowledge graph, researchers can readily reason and analyze this information on a much larger scale. Here, we present the COVID-19 Knowledge Graph, an expansive cause-and-effect network constructed from scientific literature on the new coronavirus that aims to provide a comprehensive view of its pathophysiology. To make this resource available to the research community and facilitate its exploration and analysis, we also implemented a web application and released the KG in multiple standard formats. Availability The COVID-19 Knowledge Graph is publicly available under CC-0 license at https://github.com/covid19kg and https://bikmi.covid19-knowledgespace.de. Supplementary information Supplementary data are available online.

Description: Motivation One of the most important problems in drug discovery research is to precisely predict a new indication for an existing drug, i.e. drug repositioning. Recent recommendation system-based methods have tackled this problem using matrix completion models. The models identify latent factors contributing to known drug-disease associations, and then infer novel drug-disease associations by the correlations between latent factors. However, these models have not fully considered the various drug data sources and the sparsity of the drug-disease association matrix. In addition, using the global structure of the drug-disease association data may introduce noise, and consequently limit the prediction power. Results In this work, we propose a novel drug repositioning approach by using Bayesian inductive matrix completion (DRIMC). First, we embed four drug data sources into a drug similarity matrix and two disease data sources in a disease similarity matrix. Then, for each drug or disease, its feature is described by similarity values between it and its nearest neighbors, and these features for drugs and diseases are mapped onto a shared latent space. We model the association probability for each drug-disease pair by inductive matrix completion, where the properties of drugs and diseases are represented by projections of drugs and diseases, respectively. As the known drug-disease associations have been manually verified, they are more trustworthy and important than the unknown pairs. We assign higher confidence levels to known association pairs compared with unknown pairs. We perform comprehensive experiments on three benchmark datasets, and DRIMC improves prediction accuracy compared with six stat-of-the-art approaches. Availability and implementation Source code and datasets are available at https://github.com/linwang1982/DRIMC. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation Neural methods to extract drug-drug interactions (DDIs) from literature require a large number of annotations. In this study, we propose a novel method to effectively utilize external drug database information as well as information from large-scale plain text for DDI extraction. Specifically, we focus on drug description and molecular structure information as the drug database information. Results We evaluated our approach on the DDIExtraction 2013 shared task data set. We obtained the following results. First, large-scale raw text information can greatly improve the performance of extracting DDIs when combined with the existing model and it shows the state-of-the-art performance. Second, each of drug description and molecular structure information is helpful to further improve the DDI performance for some specific DDI types. Finally, the simultaneous use of the drug description and molecular structure information can significantly improve the performance on all the DDI types. We showed that the plain text, the drug description information, and molecular structure information are complementary and their effective combination are essential for the improvement. Availability https://github.com/tticoin/DESC_MOL-DDIE

Description: Summary Bulk RNA sequencing studies have demonstrated that human leukocyte antigen (HLA) genes may be expressed in a cell type-specific and allele-specific fashion. Single-cell gene expression assays have the potential to further resolve these expression patterns, but currently available methods do not perform allele-specific quantification at the molecule level. Here, we present scHLAcount, a post-processing workflow for single-cell RNA-seq data that computes allele-specific molecule counts of the HLA genes based on a personalized reference constructed from the sample’s HLA genotypes. Availability and implementation scHLAcount is available under the MIT license at https://github.com/10XGenomics/scHLAcount. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation Protein domains are subunits that can fold and function independently. Correct domain boundary assignment is thus a critical step toward accurate protein structure and function analyses. There is, however, no efficient algorithm available for accurate domain prediction from sequence. The problem is particularly challenging for proteins with discontinuous domains, which consist of domain segments that are separated along the sequence. Results We developed a new algorithm, FUpred, which predicts protein domain boundaries utilizing contact maps created by deep residual neural networks coupled with coevolutionary precision matrices. The core idea of the algorithm is to retrieve domain boundary locations by maximizing the number of intra-domain contacts, while minimizing the number of inter-domain contacts from the contact maps. FUpred was tested on a large-scale dataset consisting of 2549 proteins and generated correct single- and multi-domain classifications with a Matthew’s correlation coefficient of 0.799, which was 19.1% (or 5.3%) higher than the best machine learning (or threading)-based method. For proteins with discontinuous domains, the domain boundary detection and normalized domain overlapping scores of FUpred were 0.788 and 0.521, respectively, which were 17.3% and 23.8% higher than the best control method. The results demonstrate a new avenue to accurately detect domain composition from sequence alone, especially for discontinuous, multi-domain proteins. Availability and implementation https://zhanglab.ccmb.med.umich.edu/FUpred. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation Protein structure and function are essentially determined by how the side-chain atoms interact with each other. Thus, accurate protein side-chain packing (PSCP) is a critical step toward protein structure prediction and protein design. Despite the importance of the problem, however, the accuracy and speed of current PSCP programs are still not satisfactory. Results We present FASPR for fast and accurate PSCP by using an optimized scoring function in combination with a deterministic searching algorithm. The performance of FASPR was compared with four state-of-the-art PSCP methods (CISRR, RASP, SCATD and SCWRL4) on both native and non-native protein backbones. For the assessment on native backbones, FASPR achieved a good performance by correctly predicting 69.1% of all the side-chain dihedral angles using a stringent tolerance criterion of 20°, compared favorably with SCWRL4, CISRR, RASP and SCATD which successfully predicted 68.8%, 68.6%, 67.8% and 61.7%, respectively. Additionally, FASPR achieved the highest speed for packing the 379 test protein structures in only 34.3 s, which was significantly faster than the control methods. For the assessment on non-native backbones, FASPR showed an equivalent or better performance on I-TASSER predicted backbones and the backbones perturbed from experimental structures. Detailed analyses showed that the major advantage of FASPR lies in the optimal combination of the dead-end elimination and tree decomposition with a well optimized scoring function, which makes FASPR of practical use for both protein structure modeling and protein design studies. Availability and implementation The web server, source code and datasets are freely available at https://zhanglab.ccmb.med.umich.edu/FASPR and https://github.com/tommyhuangthu/FASPR. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation Next-generation sequencing is rapidly improving diagnostic rates in rare Mendelian diseases, but even with whole genome or whole exome sequencing, the majority of cases remain unsolved. Increasingly, RNA sequencing is being used to solve many cases that evade diagnosis through sequencing alone. Specifically, the detection of aberrant splicing in many rare disease patients suggests that identifying RNA splicing outliers is particularly useful for determining causal Mendelian disease genes. However, there is as yet a paucity of statistical methodologies to detect splicing outliers. Results We developed LeafCutterMD, a new statistical framework that significantly improves the previously published LeafCutter in the context of detecting outlier splicing events. Through simulations and analysis of real patient data, we demonstrate that LeafCutterMD has better power than the state-of-the-art methodology while controlling false-positive rates. When applied to a cohort of disease-affected probands from the Mayo Clinic Center for Individualized Medicine, LeafCutterMD recovered all aberrantly spliced genes that had previously been identified by manual curation efforts. Availability and implementation The source code for this method is available under the opensource Apache 2.0 license in the latest release of the LeafCutter software package available online at http://davidaknowles.github.io/leafcutter. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation Programmed DNA elimination (PDE) plays a crucial role in the transitions between germline and somatic genomes in diverse organisms ranging from unicellular ciliates to multicellular nematodes. However, software specific for the detection of DNA splicing events is scarce. In this paper, we describe Accurate Deletion Finder (ADFinder), an efficient detector of PDEs using high-throughput sequencing data. ADFinder can predict PDEs with relatively low sequencing coverage, detect multiple alternative splicing forms in the same genomic location and calculate the frequency for each splicing event. This software will facilitate research of PDEs and all down-stream analyses. Results By analyzing genome-wide DNA splicing events in two micronuclear genomes of Oxytricha trifallax and Tetrahymena thermophila, we prove that ADFinder is effective in predicting large scale PDEs. Availability and implementation The source codes and manual of ADFinder are available in our GitHub website: https://github.com/weibozheng/ADFinder. Supplementary information Supplementary data are available at Bioinformatics online.

Description:   The development of new drugs is costly, time consuming, and often accompanied with safety issues. Drug repurposing can avoid the expensive and lengthy process of drug development by finding new uses for already approved drugs. In order to repurpose drugs effectively, it is useful to know which proteins are targeted by which drugs. Computational models that estimate the interaction strength of new drug–target pairs have the potential to expedite drug repurposing. Several models have been proposed for this task. However, these models represent the drugs as strings, which is not a natural way to represent molecules. We propose a new model called GraphDTA that represents drugs as graphs and uses graph neural networks to predict drug–target affinity. We show that graph neural networks not only predict drug–target affinity better than non-deep learning models, but also outperform competing deep learning methods. Our results confirm that deep learning models are appropriate for drug–target binding affinity prediction, and that representing drugs as graphs can lead to further improvements. Availability of data and materials The proposed models are implemented in Python. Related data, pre-trained models, and source code are publicly available at https://github.com/thinng/GraphDTA. All scripts and data needed to reproduce the post-hoc statistical analysis are available from https://doi.org/10.5281/zenodo.3603523.

Description: Summary Fully realizing the promise of personalized medicine will require rapid and accurate classification of pathogenic human variation. Multiplexed assays of variant effect (MAVEs) can experimentally test nearly all possible variants in selected gene targets. Planning a MAVE study involves identifying target genes with clinical impact, and identifying scalable functional assays for that target. Here, we describe MaveQuest, a web-based resource enabling systematic variant effect mapping studies by identifying potential functional assays, disease phenotypes and clinical relevance for nearly all human protein-coding genes. Availability and implementation MaveQuest service: https://mavequest.varianteffect.org/. MaveQuest source code: https://github.com/kvnkuang/mavequest-front-end/. Supplementary information Supplementary data are available at Bioinformatics online.

Description: Motivation Nearly 40% of the genes in sequenced genomes have no experimentally- or computationally-derived functional annotations. To fill this gap, we seek to develop methods for network-based gene function prediction that can integrate heterogeneous data for multiple species with experimentally-based functional annotations and systematically transfer them to newly-sequenced organisms on a genomewide scale. However, the large sizes of such networks pose a challenge for the scalability of current methods. Results We develop a label propagation algorithm called FastSinkSource. By formally bounding its rate of progress, we decrease the running time by a factor of 100 without sacrificing accuracy. We systematically evaluate many approaches to construct multi-species bacterial networks and apply FastSinkSource and other state-of-the-art methods to these networks. We find that the most accurate and efficient approach is to pre-compute annotation scores for species with experimental annotations, and then to transfer them to other organisms. In this manner, FastSinkSource runs in under three minutes for 200 bacterial species. Availability and Implementation An implementation of our framework and all data used in this research are available at https://github.com/Murali-group/multi-species-GOA-prediction. Contact murali@cs.vt.edu Supplementary Information Supplementary information is available at Bioinformatics online.

Description: Motivation Infection with strains of different subtypes and the subsequent crossover reading between the two strands of genomic RNAs by host cells’ reverse transcriptase are the main causes of the vast HIV-1 sequence diversity. Such inter-subtype genomic recombinants can become circulating recombinant forms (CRFs) after widespread transmissions in a population. Complete prediction of all the subtype sources of a CRF strain is a complicated machine learning problem. It is also difficult to understand whether a strain is an emerging new subtype and if so, how to accurately identify the new components of the genetic source. Results We introduce a multi-label learning algorithm for the complete prediction of multiple sources of a CRF sequence as well as the prediction of its chronological number. The prediction is strengthened by a voting of various multi-label learning methods to avoid biased decisions. In our steps, frequency and position features of the sequences are both extracted to capture signature patterns of pure subtypes and CRFs. The method was applied to 7185 HIV-1 sequences, comprising 5530 pure subtype sequences and 1655 CRF sequences. Results have demonstrated that the method can achieve very high accuracy (reaching 99%) in the prediction of the complete set of labels of HIV-1 recombinant forms. A few wrong predictions are actually incomplete predictions, very close to the complete set of genuine labels. Availability https://github.com/Runbin-tang/The-source-of-HIV-CRFs-prediction Contact yuzuguo@aliyun.com;jinyan.li@uts.edu.au Supplementary information Supplementary data are available at Bioinformatics online.