Samuel W. Lukowski

The Serum Resistome of a Globally Disseminated Multidrug Resistant Uropathogenic Escherichia coli Clone

Escherichia coli ST131 is a globally disseminated, multidrug resistant clone responsible for a high proportion of urinary tract and bloodstream infections. The rapid emergence and successful spread of E. coli ST131 is strongly associated with antibiotic resistance; however, this phenotype alone is unlikely to explain its dominance amongst multidrug resistant uropathogens circulating worldwide in hospitals and the community. Thus, a greater understanding of the molecular mechanisms that underpin the fitness of E. coli ST131 is required. In this study, we employed hyper-saturated transposon mutagenesis in combination with multiplexed transposon directed insertion-site sequencing to define the essential genes required for in vitro growth and the serum resistome (i.e. genes required for resistance to human serum) of E. coli EC958, a representative of the predominant E. coli ST131 clonal lineage. We identified 315 essential genes in E. coli EC958, 231 (73%) of which were also essential in E. coli K-12. The serum resistome comprised 56 genes, the majority of which encode membrane proteins or factors involved in lipopolysaccharide (LPS) biosynthesis. Targeted mutagenesis confirmed a role in serum resistance for 46 (82%) of these genes. The murein lipoprotein Lpp, along with two lipid A-core biosynthesis enzymes WaaP and WaaG, were most strongly associated with serum resistance. While LPS was the main resistance mechanism defined for E. coli EC958 in serum, the enterobacterial common antigen and colanic acid also impacted on this phenotype. Our analysis also identified a novel function for two genes, hyxA and hyxR, as minor regulators of O-antigen chain length. This study offers novel insight into the genetic make-up of E. coli ST131, and provides a framework for future research on E. coli and other Gram-negative pathogens to define their essential gene repertoire and to dissect the molecular mechanisms that enable them to survive in the bloodstream and cause disease.

Ferrets exclusively synthesize Neu5Ac and express naturally humanized influenza A virus receptors

Mammals express the sialic acids N-acetylneuraminic acid (Neu5Ac) and N-glycolylneuraminic acid (Neu5Gc) on cell surfaces, where they act as receptors for pathogens, including influenza A virus (IAV). Neu5Gc is synthesized from Neu5Ac by the enzyme cytidine monophosphate-N-acetylneuraminic acid hydroxylase (CMAH). In humans, this enzyme is inactive and only Neu5Ac is produced. Ferrets are susceptible to human-adapted IAV strains and have been the dominant animal model for IAV studies. Here we show that ferrets, like humans, do not synthesize Neu5Gc. Genomic analysis reveals an ancient, nine-exon deletion in the ferret CMAH gene that is shared by the Pinnipedia and Musteloidia members of the Carnivora. Interactions between two human strains of IAV with the sialyllactose receptor (sialic acid—α2,6Gal) confirm that the type of terminal sialic acid contributes significantly to IAV receptor specificity. Our results indicate that exclusive expression of Neu5Ac contributes to the susceptibility of ferrets to human-adapted IAV strains.

Gene Age Predicts the Strength of Purifying Selection Acting on Gene Expression Variation in Humans

Gene expression levels can be subject to selection. We hypothesized that the age of gene origin is associated with expression constraints, given that it affects the level of gene integration into the functional cellular environment. By studying the genetic variation affecting gene expression levels (cis expression quantitative trait loci [cis-eQTLs]) and protein levels (cis protein QTLs [cis-pQTLs]), we determined that young, primate-specific genes are enriched in cis-eQTLs and cis-pQTLs. Compared to cis-eQTLs of old genes originating before the zebrafish divergence, cis-eQTLs of young genes have a higher effect size, are located closer to the transcription start site, are more significant, and tend to influence genes in multiple tissues and populations. These results suggest that the expression constraint of each gene increases throughout its lifespan. We also detected a positive correlation between expression constraints (approximated by cis-eQTL properties) and coding constraints (approximated by Ka/Ks) and observed that this correlation might be driven by gene age. To uncover factors associated with the increase in gene-age-related expression constraints, we demonstrated that gene connectivity, gene involvement in complex regulatory networks, gene haploinsufficiency, and the strength of posttranscriptional regulation increase with gene age. We also observed an increase in heritability of gene expression levels with age, implying a reduction of the environmental component. In summary, we show that gene age shapes key gene properties during evolution and is therefore an important component of genome function.

Disrupted posttranscriptional regulation of the cystic fibrosis transmembrane conductance regulator (CFTR) by a 5′UTR mutation is associated with a cftr‐related disease

Cystic fibrosis (CF) is characterized as a single‐gene disorder with a simple, autosomal recessive mode of inheritance. However, translation of cystic fibrosis transmembrane conductance regulator (CFTR) genotype into CF phenotype is influenced by nucleotide sequence variations at multiple genetic loci, and individuals heterozygous for CFTR mutations are predisposed to a range of CFTR‐related conditions, such as disseminated bronchiectasis. CF disease severity and CFTR‐related conditions are more akin to complex, multifactorial traits, which are increasingly being associated with mutations that perturb gene expression. We have identified a patient with disseminated bronchiectasis, who is heterozygous for a single‐nucleotide substitution in the CFTR 5′ untranslated region (UTR) (c.‐34C>T). The c.‐34C>T mutation creates an upstream AUG codon and upstream open reading frame that overlaps, and is out of frame with, the CFTR protein coding sequence. Using luciferase reporter constructs, we have shown that the c.‐34C>T mutation decreases gene expression by 85–99%, by reducing translation efficiency and mRNA stability. This is the first CFTR regulatory mutation shown to act at a posttranscriptional level that reduces the synthesis of normal CFTR (Class V), and reaffirms the importance of regulatory mutations as a genetic basis of multifactorial phenotypes. ©2011 Wiley‐Liss, Inc.

FGB mutations leading to congenital quantitative fibrinogen deficiencies: an update and report of four novel mutations

INTRODUCTION Causative mutations leading to congenital quantitative fibrinogen are frequently clustered in FGA encoding the fibrinogen Aα-chain. Mutations of FGB encoding the Bβ-chain are less common and of interest since the Bβ-chain is considered the rate-limiting factor in the hepatic production of the fibrinogen hexamer. METHOD Four novel FGB mutations were identified in two afibrinogenemic (one new-born and one 30 years old male) and hypofibrinogenemic (a 49 years old female) patient, with heterogeneous thrombotic and bleeding phenotype. The human fibrinogen beta chain precursor protein sequence (P02675) was obtained from the UniProt database. The resulting models were analysed in SwissPdbViewer 4.1 and POV-Ray 3.7. RESULTS The FGB c.895T>C p.Y299H (numbering from the initiator Met) and the FGB c.1415G>T p.G472V were predicted to be deleterious by SIFT analysis. The first replaces an uncharged aromatic amino acid side chain by a positively charged side chain modifying the balance in the distribution of hydrophobic and hydrophilic of the 10 Å neighbourhood residues. The second replaces one non-charged aliphatic side chain by another without any changes for the 10 Å surrounding region. The FGB c.352C>T p.Q118X leads to a severe premature termination codon and the FGB intron 4: IVS4-1G>C (c719-1G>C) leads to skipping of exon 5 or usage of a cryptic acceptor site located upstream or downstream of the normal site. CONCLUSIONS The continuous characterization of novel molecular defects responsible for fibrinogen deficiency combined with modelling of mutant proteins will continue to provide a better comprehension of the complexity of fibrinogen synthesis and physiology.

Hypofibrinogenemia and liver disease: a new case of Aguadilla fibrinogen and review of the literature.

Fibrinogen storage disease (FSD) is characterized by hypofibrinogenemia and hepatic inclusions due to impaired release of mutant fibrinogen which accumulates and aggregates in the hepatocellular endoplasmic reticulum. Liver disease is variable.

The genetic regulation of transcription in human endometrial tissue

Study question Do genetic effects regulate gene expression in human endometrium? Summary answer This study demonstrated strong genetic effects on endometrial gene expression and some evidence for genetic regulation of gene expression in a menstrual cycle stage-specific manner. What is known already Genetic effects on expression levels for many genes are tissue specific. Endometrial gene expression varies across menstrual cycle stages and between individuals, but there are limited data on genetic control of expression in endometrium. Study design, size, duration We analysed genome-wide genotype and gene expression data to map cis expression quantitative trait loci (eQTL) in endometrium. Participants/materials, setting, methods We recruited 123 women of European ancestry. DNA samples from blood were genotyped on Illumina HumanCoreExome chips. Total RNA was extracted from endometrial tissues. Whole-transcriptome profiles were characterized using Illumina Human HT-12 v4.0 Expression Beadchips. We performed eQTL mapping with ~8 000 000 genotyped and imputed single nucleotide polymorphisms (SNPs) and 12 329 genes. Main results and the role of chance We identified a total of 18 595 cis SNP-probe associations at a study-wide level of significance (P < 1 × 10-7), which correspond to independent eQTLs for 198 unique genes. The eQTLs with the largest effect in endometrial tissue were rs4902335 for CHURC1 (P = 1.05 × 10-32) and rs147253019 for ZP3 (P = 8.22 × 10-30). We further performed a context-specific eQTL analysis to investigate if genetic effects on gene expression regulation act in a menstrual cycle-specific manner. Interestingly, five cis-eQTLs were identified with a significant stage-by-genotype interaction. The strongest stage interaction was the eQTL for C10ORF33 (PYROXD2) with SNP rs2296438 (P = 2.0 × 10-4), where we observe a 2-fold difference in the average expression levels of heterozygous samples depending on the stage of the menstrual cycle. Large scale data The summary eQTL results are publicly available to browse or download. Limitations, reasons for caution A limitation of the present study was the relatively modest sample size. It was not powered to identify trans-eQTLs and larger sample sizes will also be needed to provide better power to detect cis-eQTLs and cycle stage-specific effects, given the substantial changes in expression across the menstrual cycle for many genes. Wider implications of the findings Identification of endometrial eQTLs provides a platform for better understanding genetic effects on endometriosis risk and other endometrial-related pathologies. Study funding/competing interest(s) Funding for this work was provided by NHMRC Project Grants GNT1026033, GNT1049472, GNT1046880, GNT1050208, GNT1105321 and APP1083405. There are no competing interests.

ascend: R package for analysis of single cell RNA-seq data

Summary ascend is an R package comprised of fast, streamlined analysis functions optimized to address the statistical challenges of single cell RNA-seq. The package incorporates novel and established methods to provide a flexible framework to perform filtering, quality control, normalization, dimension reduction, clustering, differential expression and a wide-range of plotting. ascend is designed to work with scRNA-seq data generated by any high-throughput platform, and includes functions to convert data objects between software packages. Availability The R package and associated vignettes are freely available at https://github.com/IMB-Computational-Genomics-Lab/ascend. Contact joseph.powell@uq.edu.au Supplementary information An example dataset is available at ArrayExpress, accession number E-MTAB-6108

CRIg-expressing peritoneal macrophages are associated with disease severity in patients with cirrhosis and ascites

Infections are an important cause of morbidity and mortality in patients with decompensated cirrhosis and ascites. Hypothesizing that innate immune dysfunction contributes to susceptibility to infection, we assessed ascitic fluid macrophage phenotype and function. The expression of complement receptor of the immunoglobulin superfamily (CRIg) and CCR2 defined two phenotypically and functionally distinct peritoneal macrophage subpopulations. The proportion of CRIghi macrophages differed between patients and in the same patient over time, and a high proportion of CRIghi macrophages was associated with reduced disease severity (model for end-stage liver disease) score. As compared with CRIglo macrophages, CRIghi macrophages were highly phagocytic and displayed enhanced antimicrobial effector activity. Transcriptional profiling by RNA sequencing and comparison with human macrophage and murine peritoneal macrophage expression signatures highlighted similarities among CRIghi cells, human macrophages, and mouse F4/80hi resident peritoneal macrophages and among CRIglo macrophages, human monocytes, and mouse F4/80lo monocyte-derived peritoneal macrophages. These data suggest that CRIghi and CRIglo macrophages may represent a tissue-resident population and a monocyte-derived population, respectively. In conclusion, ascites fluid macrophage subset distribution and phagocytic capacity is highly variable among patients with chronic liver disease. Regulating the numbers and/or functions of these macrophage populations could provide therapeutic opportunities in cirrhotic patients.

Single cell RNA sequencing of stem cell-derived retinal ganglion cells

We used single cell sequencing technology to characterize the transcriptomes of 1,174 human embryonic stem cell-derived retinal ganglion cells (RGCs) at the single cell level. The human embryonic stem cell line BRN3B-mCherry (A81-H7), was differentiated to RGCs using a guided differentiation approach. Cells were harvested at day 36 and prepared for single cell RNA sequencing. Our data indicates the presence of three distinct subpopulations of cells, with various degrees of maturity. One cluster of 288 cells showed increased expression of genes involved in axon guidance together with semaphorin interactions, cell-extracellular matrix interactions and ECM proteoglycans, suggestive of a more mature RGC phenotype.