Search

Abstract Background Attine ants live in an intensely studied tripartite mutualism with the fungus Leucoagaricus gongylophorus, which provides food to the ants, and with antibiotic-producing actinomycete bacteria. One hypothesis suggests that bacteria from the genus Pseudonocardia are the sole, co-evolved mutualists of attine ants and are transmitted vertically by the queens. A recent study identified a Pseudonocardia-produced antifungal, named dentigerumycin, associated with the lower attine Apterostigma dentigerum consistent with the idea that co-evolved Pseudonocardia make novel antibiotics. An alternative possibility is that attine ants sample actinomycete bacteria from the soil, selecting and maintaining those species that make useful antibiotics. Consistent with this idea, a Streptomyces species associated with the higher attine Acromyrmex octospinosus was recently shown to produce the well-known antifungal candicidin. Candicidin production is widespread in environmental isolates of Streptomyces, so this could either be an environmental contaminant or evidence of recruitment of useful actinomycetes from the environment. It should be noted that the two possibilities for actinomycete acquisition are not necessarily mutually exclusive. Results In order to test these possibilities we isolated bacteria from a geographically distinct population of A. octospinosus and identified a candicidin-producing Streptomyces species, which suggests that they are common mutualists of attine ants, most probably recruited from the environment. We also identified a Pseudonocardia species in the same ant colony that produces an unusual polyene antifungal, providing evidence for co-evolution of Pseudonocardia with A. octospinosus. Conclusions Our results show that a combination of co-evolution and environmental sampling results in the diversity of actinomycete symbionts and antibiotics associated with attine ants.

Abstract Background The Brassica species, related to Arabidopsis thaliana, include an important group of crops and represent an excellent system for studying the evolutionary consequences of polyploidy. Previous studies have led to a proposed structure for an ancestral karyotype and models for the evolution of the B. rapa genome by triplication and segmental rearrangement, but these have not been validated at the sequence level. Results We developed computational tools to analyse the public collection of B. rapa BAC end sequence, in order to identify candidates for representing collinearity discontinuities between the genomes of B. rapa and A. thaliana. For each putative discontinuity, one of the BACs was sequenced and analysed for collinearity with the genome of A. thaliana. Additional BAC clones were identified and sequenced as part of ongoing efforts to sequence four chromosomes of B. rapa. Strikingly few of the 19 inter-chromosomal rearrangements corresponded to the set of collinearity discontinuities anticipated on the basis of previous studies. Our analyses revealed numerous instances of newly detected collinearity blocks. For B. rapa linkage group A8, we were able to develop a model for the derivation of the chromosome from the ancestral karyotype. We were also able to identify a rearrangement event in the ancestor of B. rapa that was not shared with the ancestor of A. thaliana, and is represented in triplicate in the B. rapa genome. In addition to inter-chromosomal rearrangements, we identified and analysed 32 BACs containing the end points of segmental inversion events. Conclusion Our results show that previous studies of segmental collinearity between the A. thaliana, Brassica and ancestral karyotype genomes, although very useful, represent over-simplifications of their true relationships. The presence of numerous cryptic collinear genome segments and the frequent occurrence of segmental inversions mean that inference of the positions of genes in B. rapa based on the locations of orthologues in A. thaliana can be misleading. Our results will be of relevance to a wide range of plants that have polyploid genomes, many of which are being considered according to a paradigm of comprising conserved synteny blocks with respect to sequenced, related genomes.

Abstract Background The Brassica species include an important group of crops and provide opportunities for studying the evolutionary consequences of polyploidy. They are related to Arabidopsis thaliana, for which the first complete plant genome sequence was obtained and their genomes show extensive, although imperfect, conserved synteny with that of A. thaliana. A large number of EST sequences, derived from a range of different Brassica species, are available in the public database, but no public microarray resource has so far been developed for these species. Results We assembled unigenes using ~800,000 EST sequences, mainly from three species: B. napus, B. rapa and B. oleracea. The assembly was conducted with the aim of co-assembling ESTs of orthologous genes (including homoeologous pairs of genes in B. napus from each of the A and C genomes), but resolving assemblies of paralogous, or paleo-homoeologous, genes (i.e. the genes related by the ancestral genome triplication observed in diploid Brassica species). 90,864 unique sequence assemblies were developed. These were incorporated into the BAC sequence annotation for the Brassica rapa Genome Sequencing Project, enabling the identification of cognate genomic sequences for a proportion of them. A 60-mer oligo microarray comprising 94,558 probes was developed using the unigene sequences. Gene expression was analysed in reciprocal resynthesised B. napus lines and the B. oleracea and B. rapa lines used to produce them. The analysis showed that significant expression could consistently be detected in leaf tissue for 35,386 unigenes. Expression was detected across all four genotypes for 27,355 unigenes, genome-specific expression patterns were observed for 7,851 unigenes and 180 unigenes displayed other classes of expression pattern. Principal component analysis (PCA) clearly resolved the individual microarray datasets for B. rapa, B. oleracea and resynthesised B. napus. Quantitative differences in expression were observed between the resynthesised B. napus lines for 98 unigenes, most of which could be classified into non-additive expression patterns, including 17 that showed cytoplasm-specific patterns. We further characterized the unigenes for which A genome-specific expression was observed and cognate genomic sequences could be identified. Ten of these unigenes were found to be Brassica-specific sequences, including two that originate from complex loci comprising gene clusters. Conclusion We succeeded in developing a Brassica community microarray resource. Although expression can be measured for the majority of unigenes across species, there were numerous probes that reported in a genome-specific manner. We anticipate that some proportion of these will represent species-specific transcripts and the remainder will be the consequence of variation of sequences within the regions represented by the array probes. Our studies demonstrated that the datasets obtained from the arrays can be used for typical analyses, including PCA and the analysis of differential expression. We have also demonstrated that Brassica-specific transcripts identified in silico in the sequence assembly of public EST database accessions are indeed reported by the array. These would not be detectable using arrays designed using A. thaliana sequences.

In the past three decades, dietary and lifestyle changes worldwide have resulted in a global increase in the prevalence of obesity in both adults and children. Known to be highly influenced by genetic, environmental and lifestyle factors, obesity is characterized by a low-grade chronic inflammation that contributes to the development of other metabolic diseases such as diabetes and cardiovascular disease. Recently, the gut microbiome has been added as a cause/contributor to the development of obesity. As differences in the microbiome between obese and normoweight individuals have been observed, we set out to determine whether infants harbor an obesogenic microbiome early on and whether the pre-pregnancy status of the mother (obese or normoweight) is correlated to their infant's microbiome composition. Using shotgun sequencing, we analyzed stool samples throughout the first year of life from infants born to obese ( n  = 23 participants, m  = 104 samples) and normoweight ( n  = 23 participants, m  = 99 samples) mothers. We found that the infants' microbiome diversity at taxonomic and functional levels was significantly influenced by time (ANOVA p  < 0.001) but not by the mother's pre-pregnancy status. Overall, no deterministic succession of taxa or functions was observed. However, infants born to obese mothers were found to have a significantly higher Bacillota/Bacteroidota ratio ( p  = 0.02) at six months, were significantly depleted from six months old of the well-established obesity biomarkers Akkermansia municiphila and Faecalibacterium prausnitzii ( p  < 0.01), and were at one week old, significantly enriched in pathways such as the UDP-N-acetyl-D-glucosamine biosynthesis II ( p  = 0.02) involved in leptin production, suggesting perhaps that there may exist some underlying mechanisms that dictate the development of an obesogenic microbiota early on., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Vallès, Arshad, Abdalbaqi, Inman, Ahmad, Drou, Gunsalus, Ali, Tahlak and Abdulle.)

Ferroptosis is a unique iron-dependent form of non-apoptotic cell death characterized by devastating lipid peroxidation. Whilst growing evidence suggests that ferroptosis is a type of autophagy-dependent cell death, the underlying molecular mechanisms regulating ferroptosis are largely unknown. In this study, through an unbiased RNA-sequencing screening, we demonstrate the activation of a multi-faceted tumor-suppressor protein Par-4/PAWR during ferroptosis. Functional studies reveal that genetic depletion of Par-4 effectively blocks ferroptosis, whereas Par-4 overexpression sensitizes cells to undergo ferroptosis. More importantly, we have determined that Par-4-triggered ferroptosis is mechanistically driven by the autophagic machinery. Upregulation of Par-4 promotes activation of ferritinophagy (autophagic degradation of ferritin) via the nuclear receptor co-activator 4 (NCOA4), resulting in excessive release of free labile iron and, hence, enhanced lipid peroxidation and ferroptosis. Inhibition of Par-4 dramatically suppresses the NCOA4-mediated ferritinophagy signaling axis. Our results also establish that Par-4 activation positively correlates with reactive oxygen species (ROS) production, which is critical for ferritinophagy-mediated ferroptosis. Furthermore, Par-4 knockdown effectively blocked ferroptosis-mediated tumor suppression in the mouse xenograft models. Collectively, these findings reveal that Par-4 has a crucial role in ferroptosis, which could be further exploited for cancer therapy., (© 2024. The Author(s).)

Oxeiptosis is a recently identified reactive oxygen species (ROS)-sensitive, caspase independent, non-inflammatory regulated cell death pathway. The activation of Kelch-like ECH-associated protein 1-Phosphoglycerate mutase 5-Apoptosis inducing factor mitochondria associated 1 (KEAP1-PGAM5-AIFM1) pathway is the key signaling event in the execution of oxeiptosis. In the present study, we demonstrate that sanguinarine (SNG), a quaternary benzophenanthridine alkaloid, induces oxeiptosis in human colorectal cancer (CRC) cells via ROS, specifically hydrogen peroxide (H 2 O 2 )-dependent activation of KEAP1-PGAM5-AIFM1 signaling axis. Whilst, knockdown of KEAP1, PGAM5, and AIFM1 largely abolishes SNG-induced oxeiptosis, hence reinforcing the importance of the role of this pathway in the SNG-mediated cytotoxicity. Moreover, extracellular addition of H 2 O 2 sensitizes SNG-induced oxeiptosis in CRC cells, while removal of intracellular ROS by ROS scavengers, not only alleviated the overproduction of ROS caused by SNG, but also reversed the biochemical events associated with oxeiptosis. Finally, in vivo study demonstrates that SNG effectively reduces the tumor growth in HT-29 xenograft mouse model through features associated with oxeiptosis. This study highlights oxeiptosis as a novel tumor suppressive mechanism and further investigation of the role of oxeiptosis in cancer treatment is warranted., (© 2023. The Author(s).)

The studies on the prevalence and genetic diversity of begomoviruses in Saudi Arabia are minimal. In this study, field-grown symptomatic tomato and muskmelon plants were collected, and initially, begomovirus infection was confirmed by the core coat protein sequences. Four tomato and two muskmelon plants with viral infections were further evaluated for Illumina MiSeq sequencing, and twelve sequences (2.7-2.8 kb) equivalent to the full-length DNA-A or DNA-B components of begomoviruses were obtained along with eight sequences (~1.3-1.4 kb) equivalent to the begomovirus-associated DNA-satellite components. Four begomovirus sequences obtained from tomato plants were variants of tomato yellow leaf curl virus (TYLCV) with nt sequence identities of 95.3-100%. Additionally, two tomato plants showed a mixed infection of TYLCV and cotton leaf curl Gezira virus (CLCuGeV), okra yellow crinkle Cameroon alphasatellite (OYCrCMA), and okra leaf curl Oman betasatellite (OLCuOMB). Meanwhile, from muskmelon plants, two sequences were closely related (99-99.6%) to the tomato leaf curl Palampur virus (ToLCPalV) DNA-A, whereas two other sequences showed 97.9-100% sequence identities to DNA-B of ToLCPalV, respectively. Complete genome sequences of CLCuGeV and associated DNA-satellites were also obtained from these muskmelon plants. The nt sequence identities of the CLCuGeV, OYCrCMA, and OLCuOMB isolates obtained were 98.3-100%, 99.5-100%, and 95.6-99.7% with their respective available variants. The recombination was only detected in TYLCV and OLCuOMB isolates. To our knowledge, this is the first identification of a mixed infection of bipartite and monopartite begomoviruses associated with DNA-satellites from tomato and muskmelon in Saudi Arabia. The begomovirus variants reported in this study were clustered with Iranian isolates of respective begomovirus components in the phylogenetic dendrogram. Thus, the Iranian agroecological route can be a possible introduction of these begomoviruses and/or their associated DNA-satellites into Saudi Arabia.

Exposure to microgravity affects astronauts' health in adverse ways. However, less is known about the extent to which fibroblast differentiation during the wound healing process is affected by the lack of gravity. One of the key steps of this process is the differentiation of fibroblasts into myofibroblasts, which contribute functionally through extracellular matrix production and remodeling. In this work, we utilized collagen-based three-dimensional (3D) matrices to mimic interstitial tissue and studied fibroblast differentiation under simulated microgravity (sµG). Our results demonstrated that alpha-smooth muscle actin (αSMA) expression and translocation of Smad2/3 into the cell nucleus were reduced upon exposure to sµG compared to the 1 g control, which suggests the impairment of fibroblast differentiation under sµG. Moreover, matrix remodeling and production were decreased under sµG, which is in line with the impaired fibroblast differentiation. We further investigated changes on a transcriptomic level using RNA sequencing. The results demonstrated that sµG has less effect on fibroblast transcriptomes, while sµG triggers changes in the transcriptome of myofibroblasts. Several genes and biological pathways found through transcriptome analysis have previously been reported to impair fibroblast differentiation. Overall, our data indicated that fibroblast differentiation, as well as matrix production and remodeling, are impaired in 3D culture under sµG conditions.

Unicellular eukaryotic phytoplankton, such as diatoms, rely on microbial communities for survival despite lacking specialized compartments to house microbiomes (e.g., animal gut). Microbial communities have been widely shown to benefit from diatom excretions that accumulate within the microenvironment surrounding phytoplankton cells, known as the phycosphere. However, mechanisms that enable diatoms and other unicellular eukaryotes to nurture specific microbiomes by fostering beneficial bacteria and repelling harmful ones are mostly unknown. We hypothesized that diatom exudates may tune microbial communities and employed an integrated multiomics approach using the ubiquitous diatom Asterionellopsis glacialis to reveal how it modulates its naturally associated bacteria. We show that A. glacialis reprograms its transcriptional and metabolic profiles in response to bacteria to secrete a suite of central metabolites and two unusual secondary metabolites, rosmarinic acid and azelaic acid. While central metabolites are utilized by potential bacterial symbionts and opportunists alike, rosmarinic acid promotes attachment of beneficial bacteria to the diatom and simultaneously suppresses the attachment of opportunists. Similarly, azelaic acid enhances growth of beneficial bacteria while simultaneously inhibiting growth of opportunistic ones. We further show that the bacterial response to azelaic acid is numerically rare but globally distributed in the world's oceans and taxonomically restricted to a handful of bacterial genera. Our results demonstrate the innate ability of an important unicellular eukaryotic group to modulate select bacteria in their microbial consortia, similar to higher eukaryotes, using unique secondary metabolites that regulate bacterial growth and behavior inversely across different bacterial populations., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Familial breast cancer is estimated to account for 15-20% of all cases of breast cancer. Surveillance for familial breast cancer is well-established world-wide. However, this service does not exist in Jordan, due to the scarcity of information with regard to the genetic profiling of these patients, and therefore lack of recommendations for policy-makers. As such, patients with very strong family history of breast or ovarian cancers are not screened routinely; leading to preventable delay in diagnosis. Whole coding sequencing for BCRA1/BCRA2 using next-generation sequencing (NGS)/Ion PGM System was performed. Sanger sequencing were then used to confirm the pathogenic variants detected by NGS. In this study, 192 breast cancer patients (and 8 ovarian cancer cases) were included. The prevalence of recurrent pathogenic mutations was 14.5%, while the prevalence of newly detected mutations was 3.5%. Two novel pathogenic mutations were identified in BRCA2 genes. The common mutations in the Ashkenazi population used for screening may not apply in the Jordanian population, as previously reported mutations were not prevalent, and other new mutations were identified. These data will aid to establish a specific screening test for BRCA 1/BRCA2 in the Jordanian population.

The mechanisms behind the ability of Plasmodium falciparum to evade host immune system are poorly understood and are a major roadblock in achieving malaria elimination. Here, we use integrative genomic profiling and a longitudinal pediatric cohort in Burkina Faso to demonstrate the role of post-transcriptional regulation in host immune response in malaria. We report a strong signature of miRNA expression differentiation associated with P. falciparum infection (127 out of 320 miRNAs, B-H FDR 5%) and parasitemia (72 miRNAs, B-H FDR 5%). Integrative miRNA-mRNA analysis implicates several infection-responsive miRNAs (e.g., miR-16-5p, miR-15a-5p and miR-181c-5p) promoting lymphocyte cell death. miRNA cis-eQTL analysis using whole-genome sequencing data identified 1,376 genetic variants associated with the expression of 34 miRNAs (B-H FDR 5%). We report a protective effect of rs114136945 minor allele on parasitemia mediated through miR-598-3p expression. These results highlight the impact of post-transcriptional regulation, immune cell death processes and host genetic regulatory control in malaria.

Background: As high-throughput sequencing applications continue to evolve, the rapid growth in quantity and variety of sequence-based data calls for the development of new software libraries and tools for data analysis and visualization. Often, effective use of these tools requires computational skills beyond those of many researchers. To ease this computational barrier, we have created a dynamic web-based platform, NASQAR (Nucleic Acid SeQuence Analysis Resource)., Results: NASQAR offers a collection of custom and publicly available open-source web applications that make extensive use of a variety of R packages to provide interactive data analysis and visualization. The platform is publicly accessible at http://nasqar.abudhabi.nyu.edu/ . Open-source code is on GitHub at https://github.com/nasqar/NASQAR , and the system is also available as a Docker image at https://hub.docker.com/r/aymanm/nasqarall . NASQAR is a collaboration between the core bioinformatics teams of the NYU Abu Dhabi and NYU New York Centers for Genomics and Systems Biology., Conclusions: NASQAR empowers non-programming experts with a versatile and intuitive toolbox to easily and efficiently explore, analyze, and visualize their Transcriptomics data interactively. Popular tools for a variety of applications are currently available, including Transcriptome Data Preprocessing, RNA-seq Analysis (including Single-cell RNA-seq), Metagenomics, and Gene Enrichment.

During neuronal development, β-actin serves an important role in growth cone mediated axon guidance. Consistent with this notion, in vivo ablation of the β-actin gene leads to abnormalities in the nervous system. However, whether β-actin is involved in the regulation of neuronal gene programs is not known. In this study, we directly reprogramed β-actin+/+ WT, β-actin+/- HET and β-actin-/- KO mouse embryonic fibroblast (MEFs) into chemically induced neurons (CiNeurons). Using RNA-seq analysis, we profiled the transcriptome changes among the CiNeurons. We discovered that induction of neuronal gene programs was impaired in KO CiNeurons in comparison to WT ones, whereas HET CiNeurons showed an intermediate levels of induction. ChIP-seq analysis of heterochromatin markers demonstrated that the impaired expression of neuronal gene programs correlated with the elevated H3K9 and H3K27 methylation levels at gene loci in β-actin deficient MEFs, which is linked to the loss of chromatin association of the BAF complex ATPase subunit Brg1. Together, our study shows that heterochromatin alteration in β-actin null MEFs impedes the induction of neuronal gene programs during direct reprograming. These findings are in line with the notion that H3K9Me3-based heterochromatin forms a major epigenetic barrier during cell fate change., Competing Interests: The authors have declared that no competing interests exist.

In eukaryotic cells, actin regulates both cytoplasmic and nuclear functions. However, whether actin-based structures are present in the mitochondria and are involved in mitochondrial functions has not been investigated. Here, using wild-type ?-actin +/+ and knockout (KO) ?-actin ?/? mouse embryonic fibroblasts we show evidence for the defect in maintaining mitochondrial membrane potential (MMP) in ?-actin-null cells. MMP defects were associated with impaired mitochondrial DNA (mtDNA) transcription and nuclear oxidative phosphorylation (OXPHOS) gene expression. Using super-resolution microscopy we provided direct evidence on the presence of ?-actin-containing structures inside mitochondria. Large aggregates of TFAM-stained nucleoids were observed in bulb-shaped mitochondria in KO cells, suggesting defects in mitochondrial nucleoid segregation without ?-actin. The observation that mitochondria-targeted ?-actin rescued mtDNA transcription and MMP suggests an indispensable functional role of a mitochondrial ?-actin pool necessary for mitochondrial quality control., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

During differentiation and development, cell fate and identity are established by waves of genetic reprogramming. Although the mechanisms are largely unknown, during these events, dynamic chromatin reorganization is likely to ensure that multiple genes involved in the same cellular functions are coregulated, depending on the nuclear environment. In this study, using high-content screening of embryonic fibroblasts from a β-actin knockout (KO) mouse, we found major chromatin rearrangements and changes in histone modifications, such as methylated histone (H)3-lysine-(K)9. Genome-wide H3K9 trimethylation-(Me)3 landscape changes correlate with gene up- and down-regulation in β-actin KO cells. Mechanistically, we found loss of chromatin association by the Brahma-related gene ( Brg)/Brahma-associated factor (BAF) chromatin remodeling complex subunit Brg1 in the absence of β-actin. This actin-dependent chromatin reorganization was concomitant with the up-regulation of sets of genes involved in angiogenesis, cytoskeletal organization, and myofibroblast features in β-actin KO cells. Some of these genes and phenotypes were gained in a β-actin dose-dependent manner. Moreover, reintroducing a nuclear localization signal-containing β-actin in the knockout cells affected nuclear features and gene expression. Our results suggest that, by affecting the genome-wide organization of heterochromatin through the chromatin-binding activity of the BAF complex, β-actin plays an essential role in the determination of gene expression programs and cellular identity.-Xie, X., Almuzzaini, B., Drou, N., Kremb, S., Yousif, A., Östlund Farrants, A.-K., Gunsalus, K., Percipalle, P. β-Actin-dependent global chromatin organization and gene expression programs control cellular identity.

To investigate the phenomic and genomic traits that allow green algae to survive in deserts, we characterized a ubiquitous species, Chloroidium sp. UTEX 3007 , which we isolated from multiple locations in the United Arab Emirates (UAE). Metabolomic analyses of Chloroidium sp. UTEX 3007 indicated that the alga accumulates a broad range of carbon sources, including several desiccation tolerance-promoting sugars and unusually large stores of palmitate. Growth assays revealed capacities to grow in salinities from zero to 60 g/L and to grow heterotrophically on >40 distinct carbon sources. Assembly and annotation of genomic reads yielded a 52.5 Mbp genome with 8153 functionally annotated genes. Comparison with other sequenced green algae revealed unique protein families involved in osmotic stress tolerance and saccharide metabolism that support phenomic studies. Our results reveal the robust and flexible biology utilized by a green alga to successfully inhabit a desert coastline.

Motivation: The de novo assembly of genomes from whole- genome shotgun sequence data is a computationally intensive, multi-stage task and it is not known a priori which methods and parameter settings will produce optimal results. In current de novo assembly projects, a popular strategy involves trying many approaches, using different tools and settings, and then comparing and contrasting the results in order to select a final assembly for publication., Results: Herein, we present RAMPART, a configurable workflow management system for de novo genome assembly, which helps the user identify combinations of third-party tools and settings that provide good results for their particular genome and sequenced reads. RAMPART is designed to exploit High performance computing environments, such as clusters and shared memory systems, where available., Availability and Implementation: RAMPART is available under the GPLv3 license at: https://github.com/TGAC/RAMPART., (© The Author 2015. Published by Oxford University Press.)

Natural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10(-3) substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.

Lactobacillus salivarius is part of the vertebrate indigenous microbiota of the gastrointestinal tract, oral cavity, and milk. The properties associated with some L. salivarius strains have led to their use as probiotics. Here we describe the draft genome of the pig isolate L. salivarius cp400, providing insights into host-niche specialization.

Modern sequencing platforms generate enormous quantities of data in ever-decreasing amounts of time. Additionally, techniques such as multiplex sequencing allow one run to contain hundreds of different samples. With such data comes a significant challenge to understand its quality and to understand how the quality and yield are changing across instruments and over time. As well as the desire to understand historical data, sequencing centres often have a duty to provide clear summaries of individual run performance to collaborators or customers. We present StatsDB, an open-source software package for storage and analysis of next generation sequencing run metrics. The system has been designed for incorporation into a primary analysis pipeline, either at the programmatic level or via integration into existing user interfaces. Statistics are stored in an SQL database and APIs provide the ability to store and access the data while abstracting the underlying database design. This abstraction allows simpler, wider querying across multiple fields than is possible by the manual steps and calculation required to dissect individual reports, e.g. "provide metrics about nucleotide bias in libraries using adaptor barcode X, across all runs on sequencer A, within the last month". The software is supplied with modules for storage of statistics from FastQC, a commonly used tool for analysis of sequence reads, but the open nature of the database schema means it can be easily adapted to other tools. Currently at The Genome Analysis Centre (TGAC), reports are accessed through our LIMS system or through a standalone GUI tool, but the API and supplied examples make it easy to develop custom reports and to interface with other packages.

For 10,000 years pigs and humans have shared a close and complex relationship. From domestication to modern breeding practices, humans have shaped the genomes of domestic pigs. Here we present the assembly and analysis of the genome sequence of a female domestic Duroc pig (Sus scrofa) and a comparison with the genomes of wild and domestic pigs from Europe and Asia. Wild pigs emerged in South East Asia and subsequently spread across Eurasia. Our results reveal a deep phylogenetic split between European and Asian wild boars ∼1 million years ago, and a selective sweep analysis indicates selection on genes involved in RNA processing and regulation. Genes associated with immune response and olfaction exhibit fast evolution. Pigs have the largest repertoire of functional olfactory receptor genes, reflecting the importance of smell in this scavenging animal. The pig genome sequence provides an important resource for further improvements of this important livestock species, and our identification of many putative disease-causing variants extends the potential of the pig as a biomedical model.

We report the annotation and analysis of the draft genome sequence of Brassica rapa accession Chiifu-401-42, a Chinese cabbage. We modeled 41,174 protein coding genes in the B. rapa genome, which has undergone genome triplication. We used Arabidopsis thaliana as an outgroup for investigating the consequences of genome triplication, such as structural and functional evolution. The extent of gene loss (fractionation) among triplicated genome segments varies, with one of the three copies consistently retaining a disproportionately large fraction of the genes expected to have been present in its ancestor. Variation in the number of members of gene families present in the genome may contribute to the remarkable morphological plasticity of Brassica species. The B. rapa genome sequence provides an important resource for studying the evolution of polyploid genomes and underpins the genetic improvement of Brassica oil and vegetable crops.

Streptomyces spp. are common symbionts of the leaf-cutting ant species Acromyrmex octospinosus, which feeds on basidiomycete fungus leaf matter and harvests the lipid- and carbohydrate-rich gongylidia as a food source. A. octospinosus and other ant genera use antifungal compounds produced by Streptomyces spp. and other actinomycetes in order to help defend their fungal gardens from parasitic fungi. Herein, we report the draft genome sequence of Streptomyces strain S4, an antifungal-producing symbiont of A. octospinosus.

H04402 065 is one of a very small group of strains of proteolytic Clostridium botulinum that form type A5 neurotoxin. Here, we report the complete 3.9-Mb genome sequence and annotation of strain H04402 065, which was isolated from a botulism patient in the United Kingdom in 2004.

Here we describe a virus discovery protocol for a range of different virus genera, that can be applied to biopsy-sized tissue samples. Our viral enrichment procedure, validated using canine and human liver samples, significantly improves viral read copy number and increases the length of viral contigs that can be generated by de novo assembly. This in turn enables the Illumina next generation sequencing (NGS) platform to be used as an effective tool for viral discovery from tissue samples., (© 2011 Daly et al.)

Background: The species Brassica rapa includes important vegetable and oil crops. It also serves as an excellent model system to study polyploidy-related genome evolution because of its paleohexaploid ancestry and its close evolutionary relationships with Arabidopsis thaliana and other Brassica species with larger genomes. Therefore, its genome sequence will be used to accelerate both basic research on genome evolution and applied research across the cultivated Brassica species., Results: We have determined and analyzed the sequence of B. rapa chromosome A3. We obtained 31.9 Mb of sequences, organized into nine contigs, which incorporated 348 overlapping BAC clones. Annotation revealed 7,058 protein-coding genes, with an average gene density of 4.6 kb per gene. Analysis of chromosome collinearity with the A. thaliana genome identified conserved synteny blocks encompassing the whole of the B. rapa chromosome A3 and sections of four A. thaliana chromosomes. The frequency of tandem duplication of genes differed between the conserved genome segments in B. rapa and A. thaliana, indicating differential rates of occurrence/retention of such duplicate copies of genes. Analysis of 'ancestral karyotype' genome building blocks enabled the development of a hypothetical model for the derivation of the B. rapa chromosome A3., Conclusions: We report the near-complete chromosome sequence from a dicotyledonous crop species. This provides an example of the complexity of genome evolution following polyploidy. The high degree of contiguity afforded by the clone-by-clone approach provides a benchmark for the performance of whole genome shotgun approaches presently being applied in B. rapa and other species with complex genomes.

Background: The Brassica species, related to Arabidopsis thaliana, include an important group of crops and represent an excellent system for studying the evolutionary consequences of polyploidy. Previous studies have led to a proposed structure for an ancestral karyotype and models for the evolution of the B. rapa genome by triplication and segmental rearrangement, but these have not been validated at the sequence level., Results: We developed computational tools to analyse the public collection of B. rapa BAC end sequence, in order to identify candidates for representing collinearity discontinuities between the genomes of B. rapa and A. thaliana. For each putative discontinuity, one of the BACs was sequenced and analysed for collinearity with the genome of A. thaliana. Additional BAC clones were identified and sequenced as part of ongoing efforts to sequence four chromosomes of B. rapa. Strikingly few of the 19 inter-chromosomal rearrangements corresponded to the set of collinearity discontinuities anticipated on the basis of previous studies. Our analyses revealed numerous instances of newly detected collinearity blocks. For B. rapa linkage group A8, we were able to develop a model for the derivation of the chromosome from the ancestral karyotype. We were also able to identify a rearrangement event in the ancestor of B. rapa that was not shared with the ancestor of A. thaliana, and is represented in triplicate in the B. rapa genome. In addition to inter-chromosomal rearrangements, we identified and analysed 32 BACs containing the end points of segmental inversion events., Conclusion: Our results show that previous studies of segmental collinearity between the A. thaliana, Brassica and ancestral karyotype genomes, although very useful, represent over-simplifications of their true relationships. The presence of numerous cryptic collinear genome segments and the frequent occurrence of segmental inversions mean that inference of the positions of genes in B. rapa based on the locations of orthologues in A. thaliana can be misleading. Our results will be of relevance to a wide range of plants that have polyploid genomes, many of which are being considered according to a paradigm of comprising conserved synteny blocks with respect to sequenced, related genomes.

Homoeologous regions of Brassica genomes were analyzed at the sequence level. These represent segments of the Brassica A genome as found in Brassica rapa and Brassica napus and the corresponding segments of the Brassica C genome as found in Brassica oleracea and B. napus. Analysis of synonymous base substitution rates within modeled genes revealed a relatively broad range of times (0.12 to 1.37 million years ago) since the divergence of orthologous genome segments as represented in B. napus and the diploid species. Similar, and consistent, ranges were also identified for single nucleotide polymorphism and insertion-deletion variation. Genes conserved across the Brassica genomes and the homoeologous segments of the genome of Arabidopsis thaliana showed almost perfect collinearity. Numerous examples of apparent transduplication of gene fragments, as previously reported in B. oleracea, were observed in B. rapa and B. napus, indicating that this phenomenon is widespread in Brassica species. In the majority of the regions studied, the C genome segments were expanded in size relative to their A genome counterparts. The considerable variation that we observed, even between the different versions of the same Brassica genome, for gene fragments and annotated putative genes suggest that the concept of the pan-genome might be particularly appropriate when considering Brassica genomes.

Background: The Brassica species include an important group of crops and provide opportunities for studying the evolutionary consequences of polyploidy. They are related to Arabidopsis thaliana, for which the first complete plant genome sequence was obtained and their genomes show extensive, although imperfect, conserved synteny with that of A. thaliana. A large number of EST sequences, derived from a range of different Brassica species, are available in the public database, but no public microarray resource has so far been developed for these species., Results: We assembled unigenes using approximately 800,000 EST sequences, mainly from three species: B. napus, B. rapa and B. oleracea. The assembly was conducted with the aim of co-assembling ESTs of orthologous genes (including homoeologous pairs of genes in B. napus from each of the A and C genomes), but resolving assemblies of paralogous, or paleo-homoeologous, genes (i.e. the genes related by the ancestral genome triplication observed in diploid Brassica species). 90,864 unique sequence assemblies were developed. These were incorporated into the BAC sequence annotation for the Brassica rapa Genome Sequencing Project, enabling the identification of cognate genomic sequences for a proportion of them. A 60-mer oligo microarray comprising 94,558 probes was developed using the unigene sequences. Gene expression was analysed in reciprocal resynthesised B. napus lines and the B. oleracea and B. rapa lines used to produce them. The analysis showed that significant expression could consistently be detected in leaf tissue for 35,386 unigenes. Expression was detected across all four genotypes for 27,355 unigenes, genome-specific expression patterns were observed for 7,851 unigenes and 180 unigenes displayed other classes of expression pattern. Principal component analysis (PCA) clearly resolved the individual microarray datasets for B. rapa, B. oleracea and resynthesised B. napus. Quantitative differences in expression were observed between the resynthesised B. napus lines for 98 unigenes, most of which could be classified into non-additive expression patterns, including 17 that showed cytoplasm-specific patterns. We further characterized the unigenes for which A genome-specific expression was observed and cognate genomic sequences could be identified. Ten of these unigenes were found to be Brassica-specific sequences, including two that originate from complex loci comprising gene clusters., Conclusion: We succeeded in developing a Brassica community microarray resource. Although expression can be measured for the majority of unigenes across species, there were numerous probes that reported in a genome-specific manner. We anticipate that some proportion of these will represent species-specific transcripts and the remainder will be the consequence of variation of sequences within the regions represented by the array probes. Our studies demonstrated that the datasets obtained from the arrays can be used for typical analyses, including PCA and the analysis of differential expression. We have also demonstrated that Brassica-specific transcripts identified in silico in the sequence assembly of public EST database accessions are indeed reported by the array. These would not be detectable using arrays designed using A. thaliana sequences.