Assignment Of Gene Function Research Articles

Inferring Gene Regulatory Networks from Single-Cell Time-Course Data Based on Temporal Convolutional Networks

Background: Time-course single-cell RNA sequencing (scRNA-seq) data represent dynamic gene expression values that change over time, which can be used to infer causal relationships between genes and construct dynamic gene regulatory networks (GRNs). However, most of the existing methods are designed for bulk RNA sequencing (bulk RNA-seq) data and static scRNA-seq data, and only a few methods, such as CNNC and DeepDRIM can be directly applied to time-course scRNA-seq data. Objective: This work aims to infer causal relationships between genes and construct dynamic gene regulatory networks using time-course scRNA-seq data. Methods: We propose an analytical method for inferring GRNs from single-cell time-course data based on temporal convolutional networks (scTGRN), which provides a supervised learning approach to infer causal relationships among genes. scTGRN constructs a 4D tensor representing gene expression features for each gene pair, then inputs the constructed 4D tensor into the temporal convolutional network to train and infer the causal relationship between genes. Results: We validate the performance of scTGRN on five real datasets and four simulated datasets, and the experimental results show that scTGRN outperforms existing models in constructing GRNs. In addition, we test the performance of scTGRN on gene function assignment, and scTGRN outperforms other models. Conclusion: The analysis shows that scTGRN can not only accurately identify the causal relationship between genes, but also can be used to achieve gene function assignment.

Gene essentiality in the solventogenic Clostridium acetobutylicum DSM 792.

Clostridium acetobutylicum is a leading candidate to synthesize valuable compounds like three and four carbons alcohols. Its ability to convert carbohydrates into a mixture of acetone, butanol, and ethanol as well as other chemicals of interest upon genetic engineering makes it an advantageous organism for the valorization of lignocellulose-derived sugar mixtures. Since, genetic optimization depends on the fundamental insights supplied by accurate gene function assignment, gene essentiality analysis is of great interest as it can shed light on the function of many genes whose functions are still to be confirmed. The data obtained in this study will be of great value for the research community aiming to develop C. acetobutylicum as a platform organism for the production of chemicals of interest.

Determining glycosyltransferase functional order via lethality due to accumulated O-antigen intermediates, exemplified with Shigella flexneri O-antigen biosynthesis.

The O antigen (OAg) polysaccharide is one of the most diverse surface molecules of Gram-negative bacterial pathogens. The structural classification of OAg, based on serological typing and sequence analysis, is important in epidemiology and the surveillance of outbreaks of bacterial infections. Despite the diverse chemical structures of OAg repeating units (RUs), the genetic basis of RU assembly remains poorly understood and represents a major limitation in assigning gene functions in polysaccharide biosynthesis. Here, we describe a genetic approach to interrogate the functional order of glycosyltransferases (GTs). Using Shigella flexneri as a model, we established an initial glycosyltransferase (IT)-controlled system, which allows functional order allocation of the subsequent GT in a 2-fold manner as follows: (i) first, by reporting the growth defects caused by the sequestration of UndP through disruption of late GTs and (ii) second, by comparing the molecular sizes of stalled OAg intermediates when each putative GT is disrupted. Using this approach, we demonstrate that for RfbF and RfbG, the GT involved in the assembly of S. flexneri backbone OAg RU, RfbG, is responsible for both the committed step of OAg synthesis and the third transferase for the second L-Rha. We also show that RfbF functions as the last GT to complete the S. flexneri OAg RU backbone. We propose that this simple and effective genetic approach can be also extended to define the functional order of enzymatic synthesis of other diverse polysaccharides produced both by Gram-negative and Gram-positive bacteria.IMPORTANCEThe genetic basis of enzymatic assembly of structurally diverse O antigen (OAg) repeating units (RUs) in Gram-negative pathogens is poorly understood, representing a major limitation in our understanding of gene functions for the synthesis of bacterial polysaccharides. We present a simple genetic approach to confidently assign glycosyltransferase (GT) functions and the order in which they act during assembly of the OAg RU. We employed this approach to determine the functional order of GTs involved in Shigella flexneri OAg assembly. This approach can be generally applied in interrogating GT functions encoded by other bacterial polysaccharides to advance our understanding of diverse gene functions in the biosynthesis of polysaccharides, key knowledge in advancing biosynthetic polysaccharide production.

CURE on yeast genes of unknown function increases students' bioinformatics proficiency and research confidence.

Course-based undergraduate research experiences (CUREs) can reduce barriers to research opportunities while increasing student knowledge and confidence. However, the number of widely adopted, easily transferable CUREs is relatively small. Here, we describe a CURE aimed at determining the function of poorly characterized Saccharomyces cerevisiae genes. More than 20 years after sequencing of the yeast genome, nearly 10% of open reading frames (ORFs) still have at least one uncharacterized Gene Ontology (GO) term. We refer to these genes as "ORFans" and formed a consortium aimed at assigning functions to them. Specifically, over 70 faculty members attended summer workshops to learn the bioinformatics workflow and basic laboratory techniques described herein. Ultimately, this CURE was adapted for implementation at 34 institutions, resulting in over 1,300 students conducting course-based research on ORFans. Pre-/post-tests confirmed that students gained both (i) an understanding of gene ontology and (ii) knowledge regarding the use of bioinformatics to assign gene function. After using these data to craft their own hypotheses, then testing their predictions by constructing and phenotyping deletion strains, students self-reported significant gains in several areas, including computer modeling and exposure to a project where no one knows the outcome. Interestingly, most net gains self-reported by ORFan Gene Project participants were greater than published findings for CUREs assessed with the same survey instrument. The surprisingly strong impact of this CURE may be due to the incoming lack of experience of ORFan Project participants and/or the independent thought required to develop testable hypotheses from complex data sets.

AnnotaPipeline: An integrated tool to annotate eukaryotic proteins using multi-omics data.

Assignment of gene function has been a crucial, laborious, and time-consuming step in genomics. Due to a variety of sequencing platforms that generates increasing amounts of data, manual annotation is no longer feasible. Thus, the need for an integrated, automated pipeline allowing the use of experimental data towards validation of in silico prediction of gene function is of utmost relevance. Here, we present a computational workflow named AnnotaPipeline that integrates distinct software and data types on a proteogenomic approach to annotate and validate predicted features in genomic sequences. Based on FASTA (i) nucleotide or (ii) protein sequences or (iii) structural annotation files (GFF3), users can input FASTQ RNA-seq data, MS/MS data from mzXML or similar formats, as the pipeline uses both transcriptomic and proteomic information to corroborate annotations and validate gene prediction, providing transcription and expression evidence for functional annotation. Reannotation of the available Arabidopsis thaliana, Caenorhabditis elegans, Candida albicans, Trypanosoma cruzi, and Trypanosoma rangeli genomes was performed using the AnnotaPipeline, resulting in a higher proportion of annotated proteins and a reduced proportion of hypothetical proteins when compared to the annotations publicly available for these organisms. AnnotaPipeline is a Unix-based pipeline developed using Python and is available at: https://github.com/bioinformatics-ufsc/AnnotaPipeline.

DEFECTIVE EMBRYO AND MERISTEMS1 (DEM1) Is Essential for Cell Proliferation and Cell Differentiation in Tomato.

Most flowering plant species contain at least two copies of the DEFECTIVE EMBRYO AND MERISTEMS (DEM) gene with the encoded DEM proteins lacking homology to proteins of known biochemical function. In tomato (Sl; Solanum lycopersicum), stable mutations in the SlDEM1 locus result in shoot and root meristem defects with the dem1 mutant failing to progress past the cotyledon stage of seedling development. Generation of a Somatic Mutagenesis of DEM1 (SMD) transformant line in tomato allowed for the characterization of SlDEM1 gene function past the seedling stage of vegetative development with SMD plants displaying a range of leaf development abnormalities. Further, the sectored or stable in planta expression of specific regions of the SlDEM1 coding sequence also resulted in the generation of tomato transformants that displayed a range of vegetative development defects, which when considered together with the dem1 mutant seedling and SMD transformant line phenotypic data, allowed for the assignment of SlDEM1 gene function to early embryo development, adaxial epidermis cell development, lateral leaf blade expansion, and mesophyll cell proliferation and differentiation.

Mapping information-rich genotype-phenotype landscapes with genome-scale Perturb-seq

SUMMARYA central goal of genetics is to define the relationships between genotypes and phenotypes. High-content phenotypic screens such as Perturb-seq (CRISPR-based screens with single-cell RNA-sequencing readouts) enable massively parallel functional genomic mapping but, to date, have been used at limited scales. Here, we perform genome-scale Perturb-seq targeting all expressed genes with CRISPR interference (CRISPRi) across >2.5 million human cells. We use transcriptional phenotypes to predict the function of poorly characterized genes, uncovering new regulators of ribosome biogenesis (including CCDC86, ZNF236, and SPATA5L1), transcription (C7orf26), and mitochondrial respiration (TMEM242). In addition to assigning gene function, single-cell transcriptional phenotypes allow for in-depth dissection of complex cellular phenomena—from RNA processing to differentiation. We leverage this ability to systematically identify genetic drivers and consequences of aneuploidy and to discover an unanticipated layer of stress-specific regulation of the mitochondrial genome. Our information-rich genotype-phenotype map reveals a multidimensional portrait of gene and cellular function.

PPA-GCN: A Efficient GCN Framework for Prokaryotic Pathways Assignment.

With the rapid development of sequencing technology, completed genomes of microbes have explosively emerged. For a newly sequenced prokaryotic genome, gene functional annotation and metabolism pathway assignment are important foundations for all subsequent research work. However, the assignment rate for gene metabolism pathways is lower than 48% on the whole. It is even lower for newly sequenced prokaryotic genomes, which has become a bottleneck for subsequent research. Thus, the development of a high-precision metabolic pathway assignment framework is urgently needed. Here, we developed PPA-GCN, a prokaryotic pathways assignment framework based on graph convolutional network, to assist functional pathway assignments using KEGG information and genomic characteristics. In the framework, genomic gene synteny information was used to construct a network, and ideas of self-supervised learning were inspired to enhance the framework’s learning ability. Our framework is applicable to the genera of microbe with sufficient whole genome sequences. To evaluate the assignment rate, genomes from three different genera (Flavobacterium (65 genomes) and Pseudomonas (100 genomes), Staphylococcus (500 genomes)) were used. The initial functional pathway assignment rate of the three test genera were 27.7% (Flavobacterium), 49.5% (Pseudomonas) and 30.1% (Staphylococcus). PPA-GCN achieved excellence performance of 84.8% (Flavobacterium), 77.0% (Pseudomonas) and 71.0% (Staphylococcus) for assignment rate. At the same time, PPA-GCN was proved to have strong fault tolerance. The framework provides novel insights into assignment for metabolism pathways and is likely to inform future deep learning applications for interpreting functional annotations and extends to all prokaryotic genera with sufficient genomes.

Required Gene Set for Autotrophic Growth of Clostridium autoethanogenum.

ABSTRACTThe majority of the genes present in bacterial genomes remain poorly characterized, with up to one-third of those that are protein encoding having no definitive function. Transposon insertion sequencing represents a high-throughput technique that can help rectify this deficiency. The technology, however, can only be realistically applied to those species in which high rates of DNA transfer can be achieved. Here, we have developed a number of approaches that overcome this barrier in the autotrophic species Clostridium autoethanogenum by using a mariner-based transposon system. The inherent instability of such systems in the Escherichia coli conjugation donor due to transposition events was counteracted through the incorporation of a conditionally lethal codA marker on the plasmid backbone. Relatively low frequencies of transformation of the plasmid into C. autoethanogenum were circumvented through the use of a plasmid that is conditional for replication coupled with the routine implementation of an Illumina library preparation protocol that eliminates plasmid-based reads. A transposon library was then used to determine the essential genes needed for growth using carbon monoxide as the sole carbon and energy source.IMPORTANCE Although microbial genome sequences are relatively easily determined, assigning gene function remains a bottleneck. Consequently, relatively few genes are well characterized, leaving the function of many as either hypothetical or entirely unknown. High-throughput transposon sequencing can help remedy this deficiency, but is generally only applicable to microbes with efficient DNA transfer procedures. These exclude many microorganisms of importance to humankind either as agents of disease or as industrial process organisms. Here, we developed approaches to facilitate transposon insertion sequencing in the acetogen Clostridium autoethanogenum, a chassis being exploited to convert single-carbon waste gases CO and CO2 into chemicals and fuels at an industrial scale. This allowed the determination of gene essentiality under heterotrophic and autotrophic growth, providing insights into the utilization of CO as a sole carbon and energy source. The strategies implemented are translatable and will allow others to apply transposon insertion sequencing to other microbes where DNA transfer has until now represented a barrier to progress.

The cytochrome P450 Cyp6t3 is not required for ecdysone biosynthesis in Drosophila melanogaster.

The steroid hormone 20-hydroxyecdysone (20E) is essential for proper development and the timing of intermediary stage transitions in insects. As a result, there is intense interest in identifying and defining the roles of the enzymes and signaling pathways that regulate 20E production in the prothoracic gland (PG), the major endocrine organ of juvenile insect phases. Transcriptomics is one powerful tool that has been used to identify novel genes that are up- or down-regulated in the PG which may contribute to 20E regulation. Additional functional characterization of putative regulatory candidate genes typically involves qRT-PCR and/or RNAi mediated knockdown of the candidate mRNA in the PG to assess whether the gene’s expression shows temporal regulation in the PG and whether its expression is essential for proper 20E production and the correct timing of developmental transitions. While these methods have proved fruitful for identifying novel regulators of 20E production, characterizing the null phenotype of putative regulatory genes is the gold standard for assigning gene function since RNAi is known to generate various types of “off target” effects. Here we describe the genetic null mutant phenotype of the Drosophila melanogaster Cyp6t3 gene . Cyp6t3 was originally identified as a differentially regulated gene in a PG microarray screen and assigned a place in the “Black Box” step of the E biosynthetic pathway based on RNAi mediated knockdown phenotypes and rescue experiments involving feeding of various intermediate compounds of the E biosynthetic pathway. In contrast, we find that Crispr generated null mutations in Cyp6t3 are viable and have normal developmental timing. Therefore, we conclude that Cyp6t3 is not required for E production under typical lab growth conditions and therefore is not an obligate enzymatic component of the Black Box.

Large scale genomic rearrangements in selected Arabidopsis thaliana T-DNA lines are caused by T-DNA insertion mutagenesis

BackgroundExperimental proof of gene function assignments in plants is based on mutant analyses. T-DNA insertion lines provided an invaluable resource of mutants and enabled systematic reverse genetics-based investigation of the functions of Arabidopsis thaliana genes during the last decades.ResultsWe sequenced the genomes of 14 A. thaliana GABI-Kat T-DNA insertion lines, which eluded flanking sequence tag-based attempts to characterize their insertion loci, with Oxford Nanopore Technologies (ONT) long reads. Complex T-DNA insertions were resolved and 11 previously unknown T-DNA loci identified, resulting in about 2 T-DNA insertions per line and suggesting that this number was previously underestimated. T-DNA mutagenesis caused fusions of chromosomes along with compensating translocations to keep the gene set complete throughout meiosis. Also, an inverted duplication of 800 kbp was detected. About 10 % of GABI-Kat lines might be affected by chromosomal rearrangements, some of which do not involve T-DNA. Local assembly of selected reads was shown to be a computationally effective method to resolve the structure of T-DNA insertion loci. We developed an automated workflow to support investigation of long read data from T-DNA insertion lines. All steps from DNA extraction to assembly of T-DNA loci can be completed within days.ConclusionsLong read sequencing was demonstrated to be an effective way to resolve complex T-DNA insertions and chromosome fusions. Many T-DNA insertions comprise not just a single T-DNA, but complex arrays of multiple T-DNAs. It is becoming obvious that T-DNA insertion alleles must be characterized by exact identification of both T-DNA::genome junctions to generate clear genotype-to-phenotype relations.

An Overview of Antennal Esterases in Lepidoptera.

Lepidoptera are used as a model for the study of insect olfactory proteins. Among them, odorant degrading enzymes (ODEs), that degrade odorant molecules to maintain the sensitivity of antennae, have received less attention. In particular, antennal esterases (AEs; responsible for ester degradation) are crucial for intraspecific communication in Lepidoptera. Currently, transcriptomic and genomic studies have provided AEs in several species. However, efforts in gene annotation, classification, and functional assignment are still lacking. Therefore, we propose to combine evidence at evolutionary, structural, and functional level to update ODEs as well as key information into an easier classification, particularly of AEs. Finally, the kinetic parameters for putative inhibition of ODEs are discussed in terms of its role in future integrated pest management (IPM) strategies.

Identification of key genes involved in post-traumatic stress disorder: Evidence from bioinformatics analysis.

BACKGROUNDPost-traumatic stress disorder (PTSD) is a serious stress-related disorder.AIMTo identify the key genes and pathways to uncover the potential mechanisms of PTSD using bioinformatics methods.METHODSGene expression profiles were obtained from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) were identified by using GEO2R. Gene functional annotation and pathway enrichment were then conducted. The gene-pathway network was constructed with Cytoscape software. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis was applied for validation, and text mining by Coremine Medical was used to confirm the connections among genes and pathways.RESULTSWe identified 973 DEGs including 358 upregulated genes and 615 downregulated genes in PTSD. A group of centrality hub genes and significantly enriched pathways (MAPK, Ras, and ErbB signaling pathways) were identified by using gene functional assignment and enrichment analyses. Six genes (KRAS, EGFR, NFKB1, FGF12, PRKCA, and RAF1) were selected to validate using qRT-PCR. The results of text mining further confirmed the correlation among hub genes and the enriched pathways. It indicated that these altered genes displayed functional roles in PTSD via these pathways, which might serve as key signatures in the pathogenesis of PTSD.CONCLUSIONThe current study identified a panel of candidate genes and important pathways, which might help us deepen our understanding of the underlying mechanism of PTSD at the molecular level. However, further studies are warranted to discover the critical regulatory mechanism of these genes via relevant pathways in PTSD.

Comparative Homology Modeling and Physicochemical Characterization of Cyprinus carpio Hsp70 Protein

Background: Expression of the heat shock proteins (Hsp) is responsible for the protection from adverse climatic changes particularly heat stress in Common carp (Cyprinus carpio). Although, with advancement of molecular techniques, Hsp70 protein has been isolated but this protein needs to be characterized by both physicochemically and structurally for the functional annotation of fish genome. So this current study was undertaken with aim of generation of various protein models and also for thorough physiochemical characterization of this protein.Methods: In this study, Hsp70 protein of common carp was characterized by both physiochemical and structurally through insilco platform and as the crystal structure of this protein is not available, protein models were created though homology modelling upon Modeller version 9.21, Phyre2 and Swiss-model and then the generated predicted models were evaluated through Rampage, Errat, Verify 3D, ProQ and ProSA analysis.Result: Our investigation showed that this protein is very stable, hydrophilic with no disulphide bonds and the protein models which were generated from this study, are of good quality with z value of - 9.58, -9.48 and -10.93 and quality factor of 82.56, 59.43 and 95.27 respectively. So this study was concluded that the generated Hsp70 protein models would provide an avenue for the other researchers for development of high-throughput gene function assignment in fish.

A Trypanosoma brucei ORFeome-Based Gain-of-Function Library Identifies Genes That Promote Survival during Melarsoprol Treatment.

Trypanosoma brucei is an early branching protozoan parasite that causes human and animal African trypanosomiasis. Forward genetics approaches are powerful tools for uncovering novel aspects of trypanosomatid biology, pathogenesis, and therapeutic approaches against trypanosomiasis. Here, we have generated a T. brucei cloned ORFeome consisting of >90% of the targeted 7,245 genes and used it to make an inducible gain-of-function parasite library broadly applicable to large-scale forward genetic screens. We conducted a proof-of-principle genetic screen to identify genes whose expression promotes survival in melarsoprol, a critical drug of last resort. The 57 genes identified as overrepresented in melarsoprol survivor populations included the gene encoding the rate-limiting enzyme for the biosynthesis of an established drug target (trypanothione), validating the tool. In addition, novel genes associated with gene expression, flagellum localization, and mitochondrion localization were identified, and a subset of those genes increased melarsoprol resistance upon overexpression in culture. These findings offer new insights into trypanosomatid basic biology, implications for drug targets, and direct or indirect drug resistance mechanisms. This study generated a T. brucei ORFeome and gain-of-function parasite library, demonstrated the library's usefulness in forward genetic screening, and identified novel aspects of melarsoprol resistance that will be the subject of future investigations. These powerful genetic tools can be used to broadly advance trypanosomatid research.IMPORTANCE Trypanosomatid parasites threaten the health of more than 1 billion people worldwide. Because their genomes are highly diverged from those of well-established eukaryotes, conservation is not always useful in assigning gene functions. However, it is precisely among the trypanosomatid-specific genes that ideal therapeutic targets might be found. Forward genetics approaches are an effective way to identify novel gene functions. We used an ORFeome approach to clone a large percentage of Trypanosoma brucei genes and generate a gain-of-function parasite library. This library was used in a genetic screen to identify genes that promote resistance to the clinically significant yet highly toxic drug melarsoprol. Hits arising from the screen demonstrated the library's usefulness in identifying known pathways and uncovered novel aspects of resistance mediated by proteins localized to the flagellum and mitochondrion. The powerful new genetic tools generated herein are expected to promote advances in trypanosomatid biology and therapeutic development in the years to come.

A decade of advances in transposon-insertion sequencing.

It has been 10 years since the introduction of modern transposon-insertion sequencing (TIS) methods, which combine genome-wide transposon mutagenesis with high-throughput sequencing to estimate the fitness contribution or essentiality of each genetic component in a bacterial genome. Four TIS variations were published in 2009: transposon sequencing (Tn-Seq), transposon-directed insertion site sequencing (TraDIS), insertion sequencing (INSeq) and high-throughput insertion tracking by deep sequencing (HITS). TIS has since become an important tool for molecular microbiologists, being one of the few genome-wide techniques that directly links phenotype to genotype and ultimately can assign gene function. In this Review, we discuss the recent applications of TIS to answer overarching biological questions. We explore emerging and multidisciplinary methods that build on TIS, with an eye towards future applications.

Heterologous reconstitution of the biosynthesis pathway for 4-demethyl-premithramycinone, the aglycon of antitumor polyketide mithramycin

BackgroundMithramycin is an anti-tumor compound of the aureolic acid family produced by Streptomyces argillaceus. Its biosynthesis gene cluster has been cloned and characterized, and several new analogs with improved pharmacological properties have been generated through combinatorial biosynthesis. To further study these compounds as potential new anticancer drugs requires their production yields to be improved significantly. The biosynthesis of mithramycin proceeds through the formation of the key intermediate 4-demethyl-premithramycinone. Extensive studies have characterized the biosynthesis pathway from this intermediate to mithramycin. However, the biosynthesis pathway for 4-demethyl-premithramycinone remains unclear.ResultsExpression of cosmid cosAR7, containing a set of mithramycin biosynthesis genes, in Streptomyces albus resulted in the production of 4-demethyl-premithramycinone, delimiting genes required for its biosynthesis. Inactivation of mtmL, encoding an ATP-dependent acyl-CoA ligase, led to the accumulation of the tricyclic intermediate 2-hydroxy-nogalonic acid, proving its essential role in the formation of the fourth ring of 4-demethyl-premithramycinone. Expression of different sets of mithramycin biosynthesis genes as cassettes in S. albus and analysis of the resulting metabolites, allowed the reconstitution of the biosynthesis pathway for 4-demethyl-premithramycinone, assigning gene functions and establishing the order of biosynthetic steps.ConclusionsWe established the biosynthesis pathway for 4-demethyl-premithramycinone, and identified the minimal set of genes required for its assembly. We propose that the biosynthesis starts with the formation of a linear decaketide by the minimal polyketide synthase MtmPKS. Then, the cyclase/aromatase MtmQ catalyzes the cyclization of the first ring (C7–C12), followed by formation of the second and third rings (C5–C14; C3–C16) catalyzed by the cyclase MtmY. Formation of the fourth ring (C1–C18) requires MtmL and MtmX. Finally, further oxygenation and reduction is catalyzed by MtmOII and MtmTI/MtmTII respectively, to generate the final stable tetracyclic intermediate 4-demethyl-premithramycinone. Understanding the biosynthesis of this compound affords enhanced possibilities to generate new mithramycin analogs and improve their production titers for bioactivity investigation.

Gametophytic Abortion in Heterozygotes but Not in Homozygotes: Implied Chromosome Rearrangement during T-DNA Insertion at the ASF1 Locus in Arabidopsis.

T-DNA insertional mutations in Arabidopsis genes have conferred huge benefits to the research community, greatly facilitating gene function analyses. However, the insertion process can cause chromosomal rearrangements. Here, we show an example of a likely rearrangement following T-DNA insertion in the Anti-Silencing Function 1B (ASF1B) gene locus on Arabidopsis chromosome 5, so that the phenotype was not relevant to the gene of interest, ASF1B. ASF1 is a histone H3/H4 chaperone involved in chromatin remodeling in the sporophyte and during reproduction. Plants that were homozygous for mutant alleles asf1a or asf1b were developmentally normal. However, following self-fertilization of double heterozygotes (ASF1A/asf1a ASF1B/asf1b, hereafter AaBb), defects were visible in both male and female gametes. Half of the AaBb and aaBb ovules displayed arrested embryo sacs with functional megaspore identity. Similarly, half of the AaBb and aaBb pollen grains showed centromere defects, resulting in pollen abortion at the bi-cellular stage of the male gametophyte. However, inheritance of the mutant allele in a given gamete did not solely determine the abortion phenotype. Introducing functional ASF1B failed to rescue the AaBb- and aaBb- mediated abortion, suggesting that heterozygosity in the ASF1B gene causes gametophytic defects, rather than the loss of ASF1. The presence of reproductive defects in heterozygous mutants but not in homozygotes, and the characteristic all-or-nothing pollen viability within tetrads, were both indicative of commonly-observed T-DNA-mediated translocation activity for this allele. Our observations reinforce the importance of complementation tests in assigning gene function using reverse genetics.

Identification of genes required for the fitness of Streptococcus equi subsp. equi in whole equine blood and hydrogen peroxide.

The availability of next-generation sequencing techniques provides an unprecedented opportunity for the assignment of gene function. Streptococcus equi subspecies equi is the causative agent of strangles in horses, one of the most prevalent and important diseases of equids worldwide. However, the live attenuated vaccines that are utilized to control this disease cause adverse reactions in some animals. Here, we employ transposon-directed insertion-site sequencing (TraDIS) to identify genes that are required for the fitness of S. equi in whole equine blood or in the presence of H2O2 to model selective pressures exerted by the equine immune response during infection. We report the fitness values of 1503 and 1471 genes, representing 94.5 and 92.5 % of non-essential genes in S. equi , following incubation in whole blood and in the presence of H2O2, respectively. Of these genes, 36 and 15 were identified as being important to the fitness of S. equi in whole blood or H2O2, respectively, with 14 genes being important in both conditions. Allelic replacement mutants were generated to validate the fitness results. Our data identify genes that are important for S. equi to resist aspects of the immune response in vitro, which can be exploited for the development of safer live attenuated vaccines to prevent strangles.

An assessment of genome annotation coverage across the bacterial tree of life.

Although gene-finding in bacterial genomes is relatively straightforward, the automated assignment of gene function is still challenging, resulting in a vast quantity of hypothetical sequences of unknown function. But how prevalent are hypothetical sequences across bacteria, what proportion of genes in different bacterial genomes remain unannotated, and what factors affect annotation completeness? To address these questions, we surveyed over 27 000 bacterial genomes from the Genome Taxonomy Database, and measured genome annotation completeness as a function of annotation method, taxonomy, genome size, 'research bias' and publication date. Our analysis revealed that 52 and 79 % of the average bacterial proteome could be functionally annotated based on protein and domain-based homology searches, respectively. Annotation coverage using protein homology search varied significantly from as low as 14 % in some species to as high as 98 % in others. We found that taxonomy is a major factor influencing annotation completeness, with distinct trends observed across the microbial tree (e.g. the lowest level of completeness was found in the Patescibacteria lineage). Most lineages showed a significant association between genome size and annotation incompleteness, likely reflecting a greater degree of uncharacterized sequences in 'accessory' proteomes than in 'core' proteomes. Finally, research bias, as measured by publication volume, was also an important factor influencing genome annotation completeness, with early model organisms showing high completeness levels relative to other genomes in their own taxonomic lineages. Our work highlights the disparity in annotation coverage across the bacterial tree of life and emphasizes a need for more experimental characterization of accessory proteomes as well as understudied lineages.