Multiscale spatial mapping of cell populations across anatomical sites in healthy human skin and basal cell carcinoma.

Our understanding of how human skin cells differ according to anatomical site and tumour formation is limited. To address this, we have created a multiscale spatial atlas of healthy skin and basal cell carcinoma (BCC), incorporating in vivo optical coherence tomography, single-cell RNA sequencing, spatial global transcriptional profiling, and in situ sequencing. Computational spatial deconvolution and projection revealed the localisation of distinct cell populations to specific tissue contexts. Although cell populations were conserved between healthy anatomical sites and in BCC, mesenchymal cell populations including fibroblasts and pericytes retained signatures of developmental origin. Spatial profiling and in silico lineage tracing support a hair follicle origin for BCC and demonstrate that cancer-associated fibroblasts are an expansion of a POSTN+ subpopulation associated with hair follicles in healthy skin. RGS5+ pericytes are also expanded in BCC suggesting a role in vascular remodelling. We propose that the identity of mesenchymal cell populations is regulated by signals emanating from adjacent structures and that these signals are repurposed to promote the expansion of skin cancer stroma. The resource we have created is publicly available in an interactive format for the research community.

Bioinformatics of germline variant discovery for rare disease diagnostics: current approaches and remaining challenges.

Next-generation sequencing (NGS) has revolutionized the field of rare disease diagnostics. Whole exome and whole genome sequencing are now routinely used for diagnostic purposes; however, the overall diagnosis rate remains lower than expected. In this work, we review current approaches used for calling and interpretation of germline genetic variants in the human genome, and discuss the most important challenges that persist in the bioinformatic analysis of NGS data in medical genetics. We describe and attempt to quantitatively assess the remaining problems, such as the quality of the reference genome sequence, reproducible coverage biases, or variant calling accuracy in complex regions of the genome. We also discuss the prospects of switching to the complete human genome assembly or the human pan-genome and important caveats associated with such a switch. We touch on arguably the hardest problem of NGS data analysis for medical genomics, namely, the annotation of genetic variants and their subsequent interpretation. We highlight the most challenging aspects of annotation and prioritization of both coding and non-coding variants. Finally, we demonstrate the persistent prevalence of pathogenic variants in the coding genome, and outline research directions that may enhance the efficiency of NGS-based disease diagnostics.

Pathogenic diversification of the gut commensal Providencia alcalifaciens via acquisition of a second type III secretion system.

Providencia alcalifaciens is a Gram-negative bacterium found in various water and land environments and organisms, including insects and mammals. Some P. alcalifaciens strains encode gene homologs of virulence factors found in pathogenic Enterobacterales members, such as Salmonella enterica serovar Typhimurium and Shigella flexneri. Whether these genes are pathogenic determinants in P. alcalifaciens is not known. In this study, we investigated P. alcalifaciens-host interactions at the cellular level, focusing on the role of two type III secretion systems (T3SS) belonging to the Inv-Mxi/Spa family. T3SS1b is widespread in Providencia spp. and encoded on the chromosome. A large plasmid that is present in a subset of P. alcalifaciens strains, primarily isolated from diarrheal patients, encodes for T3SS1a. We show that P. alcalifaciens 205/92 is internalized into eukaryotic cells, lyses its internalization vacuole, and proliferates in the cytosol. This triggers caspase-4-dependent inflammasome responses in gut epithelial cells. The requirement for the T3SS1a in entry, vacuole lysis, and cytosolic proliferation is host cell type-specific, playing a more prominent role in intestinal epithelial cells than in macrophages or insect cells. In a bovine ligated intestinal loop model, P. alcalifaciens colonizes the intestinal mucosa and induces mild epithelial damage with negligible fluid accumulation in a T3SS1a- and T3SS1b-independent manner. However, T3SS1b was required for the rapid killing of Drosophila melanogaster. We propose that the acquisition of two T3SS has allowed P. alcalifaciens to diversify its host range, from a highly virulent pathogen of insects to an opportunistic gastrointestinal pathogen of animals.

Intracellular niche-specific profiling reveals transcriptional adaptations required for the cytosolic lifestyle of Salmonella enterica.

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a zoonotic pathogen that causes diarrheal disease in humans and animals. During salmonellosis, S. Typhimurium colonizes epithelial cells lining the gastrointestinal tract. S. Typhimurium has an unusual lifestyle in epithelial cells that begins within an endocytic-derived Salmonella-containing vacuole (SCV), followed by escape into the cytosol, epithelial cell lysis and bacterial release. The cytosol is a more permissive environment than the SCV and supports rapid bacterial growth. The physicochemical conditions encountered by S. Typhimurium within the epithelial cytosol, and the bacterial genes required for cytosolic colonization, remain largely unknown. Here we have exploited the parallel colonization strategies of S. Typhimurium in epithelial cells to decipher the two niche-specific bacterial virulence programs. By combining a population-based RNA-seq approach with single-cell microscopic analysis, we identified bacterial genes with cytosol-induced or vacuole-induced expression signatures. Using these genes as environmental biosensors, we defined that Salmonella is exposed to oxidative stress and iron and manganese deprivation in the cytosol and zinc and magnesium deprivation in the SCV. Furthermore, iron availability was critical for optimal S. Typhimurium replication in the cytosol, as well as entC, fepB, soxS, mntH and sitA. Virulence genes that are typically associated with extracellular bacteria, namely Salmonella pathogenicity island 1 (SPI1) and SPI4, showed increased expression in the cytosol compared to vacuole. Our study reveals that the cytosolic and vacuolar S. Typhimurium virulence gene programs are unique to, and tailored for, residence within distinct intracellular compartments. This archetypical vacuole-adapted pathogen therefore requires extensive transcriptional reprogramming to successfully colonize the mammalian cytosol.

Systematic benchmark of state-of-the-art variant calling pipelines identifies major factors affecting accuracy of coding sequence variant discovery.

BACKGROUND: Accurate variant detection in the coding regions of the human genome is a key requirement for molecular diagnostics of Mendelian disorders. Efficiency of variant discovery from next-generation sequencing (NGS) data depends on multiple factors, including reproducible coverage biases of NGS methods and the performance of read alignment and variant calling software. Although variant caller benchmarks are published constantly, no previous publications have leveraged the full extent of available gold standard whole-genome (WGS) and whole-exome (WES) sequencing datasets. RESULTS: In this work, we systematically evaluated the performance of 4 popular short read aligners (Bowtie2, BWA, Isaac, and Novoalign) and 9 novel and well-established variant calling and filtering methods (Clair3, DeepVariant, Octopus, GATK, FreeBayes, and Strelka2) using a set of 14 "gold standard" WES and WGS datasets available from Genome In A Bottle (GIAB) consortium. Additionally, we have indirectly evaluated each pipeline's performance using a set of 6 non-GIAB samples of African and Russian ethnicity. In our benchmark, Bowtie2 performed significantly worse than other aligners, suggesting it should not be used for medical variant calling. When other aligners were considered, the accuracy of variant discovery mostly depended on the variant caller and not the read aligner. Among the tested variant callers, DeepVariant consistently showed the best performance and the highest robustness. Other actively developed tools, such as Clair3, Octopus, and Strelka2, also performed well, although their efficiency had greater dependence on the quality and type of the input data. We have also compared the consistency of variant calls in GIAB and non-GIAB samples. With few important caveats, best-performing tools have shown little evidence of overfitting. CONCLUSIONS: The results show surprisingly large differences in the performance of cutting-edge tools even in high confidence regions of the coding genome. This highlights the importance of regular benchmarking of quickly evolving tools and pipelines. We also discuss the need for a more diverse set of gold standard genomes that would include samples of African, Hispanic, or mixed ancestry. Additionally, there is also a need for better variant caller assessment in the repetitive regions of the coding genome.

Adding function to the genome of African Salmonella Typhimurium ST313 strain D23580.

Salmonella Typhimurium sequence type (ST) 313 causes invasive nontyphoidal Salmonella (iNTS) disease in sub-Saharan Africa, targeting susceptible HIV+, malarial, or malnourished individuals. An in-depth genomic comparison between the ST313 isolate D23580 and the well-characterized ST19 isolate 4/74 that causes gastroenteritis across the globe revealed extensive synteny. To understand how the 856 nucleotide variations generated phenotypic differences, we devised a large-scale experimental approach that involved the global gene expression analysis of strains D23580 and 4/74 grown in 16 infection-relevant growth conditions. Comparison of transcriptional patterns identified virulence and metabolic genes that were differentially expressed between D23580 versus 4/74, many of which were validated by proteomics. We also uncovered the S. Typhimurium D23580 and 4/74 genes that showed expression differences during infection of murine macrophages. Our comparative transcriptomic data are presented in a new enhanced version of the Salmonella expression compendium, SalComD23580: http://bioinf.gen.tcd.ie/cgi-bin/salcom_v2.pl. We discovered that the ablation of melibiose utilization was caused by three independent SNP mutations in D23580 that are shared across ST313 lineage 2, suggesting that the ability to catabolize this carbon source has been negatively selected during ST313 evolution. The data revealed a novel, to our knowledge, plasmid maintenance system involving a plasmid-encoded CysS cysteinyl-tRNA synthetase, highlighting the power of large-scale comparative multicondition analyses to pinpoint key phenotypic differences between bacterial pathovariants.

Catching hidden variation: systematic correction of reference minor allele annotation in clinical variant calling.

PurposeWe comprehensively assessed the influence of reference minor alleles (RMAs), one of the inherent problems of the human reference genome sequence.MethodsThe variant call format (VCF) files provided by the 1000 Genomes and Exome Aggregation Consortium (ExAC) consortia were used to identify RMA sites. All coding RMA sites were checked for concordance with UniProt and the presence of same codon variants. RMA-corrected predictions of functional effect were obtained with SIFT, PolyPhen-2, and PROVEAN standalone tools and compared with dbNSFP v2.9 for consistency.ResultsWe systematically characterized the problem of RMAs and identified several possible ways in which RMA could interfere with accurate variant discovery and annotation. We have discovered a systematic bias in the automated variant effect prediction at the RMA loci, as well as widespread switching of functional consequences for variants located in the same codon as the RMA. As a convenient way to address the problem of RMAs we have developed a simple bioinformatic tool that identifies variation at RMA sites and provides correct annotations for all such substitutions. The tool is available free of charge at http://rmahunter.bioinf.me.ConclusionCorrection of RMA annotation enhances the accuracy of next-generation sequencing-based methods in clinical practice.

Long-Range Signaling in MutS and MSH Homologs via Switching of Dynamic Communication Pathways.

Allostery is conformation regulation by propagating a signal from one site to another distal site. This study focuses on the long-range communication in DNA mismatch repair proteins MutS and its homologs where intramolecular signaling has to travel over 70 Å to couple lesion detection to ATPase activity and eventual downstream repair. Using dynamic network analysis based on extensive molecular dynamics simulations, multiple preserved communication pathways were identified that would allow such long-range signaling. The pathways appear to depend on the nucleotides bound to the ATPase domain as well as the type of DNA substrate consistent with previously proposed functional cycles of mismatch recognition and repair initiation by MutS and homologs. A mechanism is proposed where pathways are switched without major conformational rearrangements allowing for efficient long-range signaling and allostery.

Targeted chromatin profiling reveals novel enhancers in Ig H and Ig L chain Loci.

The assembly and expression of mouse Ag receptor genes are controlled by a collection of cis-acting regulatory elements, including transcriptional promoters and enhancers. Although many powerful enhancers have been identified for Ig (Ig) and TCR (Tcr) loci, it remained unclear whether additional regulatory elements remain undiscovered. In this study, we use chromatin profiling of pro-B cells to define 38 epigenetic states in mouse Ag receptor loci, each of which reflects a distinct regulatory potential. One of these chromatin states corresponds to known transcriptional enhancers and identifies a new set of candidate elements in all three Ig loci. Four of the candidates were subjected to functional assays, and all four exhibit enhancer activity in B but not in T lineage cells. The new regulatory elements identified by focused chromatin profiling most likely have important functions in the creation, refinement, and expression of Ig repertoires.

Differential mismatch recognition specificities of eukaryotic MutS homologs, MutSα and MutSß.

In eukaryotes, the recognition of the DNA postreplication errors and initiation of the mismatch repair is carried out by two MutS homologs: MutSα and MutSß. MutSα recognizes base mismatches and 1 to 2 unpaired nucleotides whereas MutSß recognizes longer insertion-deletion loops (IDLs) with 1 to 15 unpaired nucleotides as well as certain mismatches. Results from molecular dynamics simulations of native MutSß:IDL-containing DNA and MutSα:mismatch DNA complexes as well as complexes with swapped DNA substrates provide mechanistic insight into how the differential substrate specificities are achieved by MutSα and MutSß, respectively. Our simulations results suggest more extensive interactions between MutSß and IDL-DNA and between MutSα and mismatch-containing DNA that suggest corresponding differences in stability. Furthermore, our simulations suggest more expanded mechanistic details involving a different degree of bending when DNA is bound to either MutSα or MutSß and a more likely opening of the clamp domains when noncognate substrates are bound. The simulation results also provide detailed information on key residues in MutSß and MutSα that are likely involved in recognizing IDL-DNA and mismatch-containing DNA, respectively.

Modulation of Salmonella virulence by a novel SPI-2 injectisome effector that interacts with the dystrophin-associated protein complex.

The injectisome encoded by Salmonella pathogenicity island 2 (SPI-2) had been thought to translocate 28 effectors. Here, we used a proteomic approach to characterize the secretome of a clinical strain of invasive non-typhoidal Salmonella enterica serovar Enteritidis that had been mutated to cause hyper-secretion of the SPI-2 injectisome effectors. Along with many known effectors, we discovered the novel SseM protein. sseM is widely distributed among the five subspecies of Salmonella enterica, is found in many clinically relevant serovars, and is co-transcribed with pipB2, a SPI-2 effector gene. The translocation of SseM required a functional SPI-2 injectisome. Following expression in human cells, SseM interacted with five components of the dystrophin-associated protein complex (DAPC), namely, ß-2-syntrophin, utrophin/dystrophin, α-catulin, α-dystrobrevin, and ß-dystrobrevin. The interaction between SseM and ß-2-syntrophin and α-dystrobrevin was verified in Salmonella Typhimurium-infected cells and relied on the postsynaptic density-95/discs large/zonula occludens-1 (PDZ) domain of ß-2-syntrophin and a sequence corresponding to a PDZ-binding motif (PBM) in SseM. A ΔsseM mutant strain had a small competitive advantage over the wild-type strain in the S. Typhimurium/mouse model of systemic disease. This phenotype was complemented by a plasmid expressing wild-type SseM from S. Typhimurium or S. Enteritidis and was dependent on the PBM of SseM. Therefore, a PBM within a Salmonella effector mediates interactions with the DAPC and modulates the systemic growth of bacteria in mice. Furthermore, the ΔsseM mutant strain displayed enhanced replication in bone marrow-derived macrophages, demonstrating that SseM restrains intracellular bacterial growth to modulate Salmonella virulence. IMPORTANCE: In Salmonella enterica, the injectisome machinery encoded by Salmonella pathogenicity island 2 (SPI-2) is conserved among the five subspecies and delivers proteins (effectors) into host cells, which are required for Salmonella virulence. The identification and functional characterization of SPI-2 injectisome effectors advance our understanding of the interplay between Salmonella and its host(s). Using an optimized method for preparing secreted proteins and a clinical isolate of the invasive non-typhoidal Salmonella enterica serovar Enteritidis strain D24359, we identified 22 known SPI-2 injectisome effectors and one new effector-SseM. SseM modulates bacterial growth during murine infection and has a sequence corresponding to a postsynaptic density-95/discs large/zonula occludens-1 (PDZ)-binding motif that is essential for interaction with the PDZ-containing host protein ß-2-syntrophin and other components of the dystrophin-associated protein complex (DAPC). To our knowledge, SseM is unique among Salmonella effectors in containing a functional PDZ-binding motif and is the first bacterial protein to target the DAPC.

Salmonella enterica serovar Typhimurium ST313 sublineage 2.2 has emerged in Malawi with a characteristic gene expression signature and a fitness advantage.

Invasive non-typhoidal Salmonella (iNTS) disease is a serious bloodstream infection that targets immune-compromised individuals, and causes significant mortality in sub-Saharan Africa. Salmonella enterica serovar Typhimurium ST313 causes the majority of iNTS in Malawi. We performed an intensive comparative genomic analysis of 608 S. Typhimurium ST313 isolates dating between 1996 and 2018 from Blantyre, Malawi. We discovered that following the arrival of the well-characterized S. Typhimurium ST313 lineage 2 in 1999, two multidrug-resistant variants emerged in Malawi in 2006 and 2008, designated sublineages 2.2 and 2.3, respectively. The majority of S. Typhimurium isolates from human bloodstream infections in Malawi now belong to sublineages 2.2 or 2.3. To understand the emergence of the prevalent ST313 sublineage 2.2, we studied two representative strains, D23580 (lineage 2) and D37712 (sublineage 2.2). The chromosome of ST313 lineage 2 and sublineage 2.2 only differed by 29 SNPs/small indels and a 3 kb deletion of a Gifsy-2 prophage region including the sseI pseudogene. Lineage 2 and sublineage 2.2 had distinctive plasmid profiles. The transcriptome was investigated in 15 infection-relevant in vitro conditions and within macrophages. During growth in physiological conditions that do not usually trigger S. Typhimurium SPI2 gene expression, the SPI2 genes of D37712 were transcriptionally active. We identified down-regulation of flagellar genes in D37712 compared with D23580. Following phenotypic confirmation of transcriptomic differences, we discovered that sublineage 2.2 had increased fitness compared with lineage 2 during mixed growth in minimal media. We speculate that this competitive advantage is contributing to the emergence of sublineage 2.2 in Malawi.

Acute response to pathogens in the early human placenta at single-cell resolution.

The placenta is a selective maternal-fetal barrier that provides nourishment and protection from infections. However, certain pathogens can attach to and even cross the placenta, causing pregnancy complications with potential lifelong impacts on the child's health. Here, we profiled at the single-cell level the placental responses to three pathogens associated with intrauterine complications-Plasmodium falciparum, Listeria monocytogenes, and Toxoplasma gondii. We found that upon exposure to the pathogens, all placental lineages trigger inflammatory responses that may compromise placental function. Additionally, we characterized the responses of fetal macrophages known as Hofbauer cells (HBCs) to each pathogen and propose that they are the probable niche for T. gondii. Finally, we revealed how P. falciparum adapts to the placental microenvironment by modulating protein export into the host erythrocyte and nutrient uptake pathways. Altogether, we have defined the cellular networks and signaling pathways mediating acute placental inflammatory responses that could contribute to pregnancy complications.

Pathogenic diversification of the gut commensal Providencia alcalifaciens via acquisition of a second type III secretion system.

Providencia alcalifaciens is a Gram-negative bacterium found in a wide variety of water and land environments and organisms. It has been isolated as part of the gut microbiome of animals and insects, as well as from stool samples of patients with diarrhea. Specific P. alcalifaciens strains encode gene homologs of virulence factors found in other pathogenic members of the same Enterobacterales order, such as Salmonella enterica serovar Typhimurium and Shigella flexneri. Whether these genes are also pathogenic determinants in P. alcalifaciens is not known. Here we have used P. alcalifaciens 205/92, a clinical isolate, with in vitro and in vivo infection models to investigate P. alcalifaciens -host interactions at the cellular level. Our particular focus was the role of two type III secretion systems (T3SS) belonging to the Inv-Mxi/Spa family. T3SS 1b is widespread in Providencia spp. and encoded on the chromosome. T3SS 1a is encoded on a large plasmid that is present in a subset of P. alcalifaciens strains, which are primarily isolates from diarrheal patients. Using a combination of electron and fluorescence microscopy and gentamicin protection assays we show that P. alcalifaciens 205/92 is internalized into eukaryotic cells, rapidly lyses its internalization vacuole and proliferates in the cytosol. This triggers caspase-4 dependent inflammasome responses in gut epithelial cells. The requirement for the T3SS 1a in entry, vacuole lysis and cytosolic proliferation is host-cell type specific, playing a more prominent role in human intestinal epithelial cells as compared to macrophages. In a bovine ligated intestinal loop model, P. alcalifaciens colonizes the intestinal mucosa, inducing mild epithelial damage with negligible fluid accumulation. No overt role for T3SS 1a or T3SS 1b was seen in the calf infection model. However, T3SS 1b was required for the rapid killing of Drosophila melanogaster . We propose that the acquisition of two T3SS by horizontal gene transfer has allowed P. alcalifaciens to diversify its host range, from a highly virulent pathogen of insects to an opportunistic gastrointestinal pathogen of animals.

Single-cell and spatially resolved interactomics of tooth-associated keratinocytes in periodontitis.

Periodontitis affects billions of people worldwide. To address relationships of periodontal niche cell types and microbes in periodontitis, we generated an integrated single-cell RNA sequencing (scRNAseq) atlas of human periodontium (34-sample, 105918-cell), including sulcular and junctional keratinocytes (SK/JKs). SK/JKs displayed altered differentiation states and were enriched for effector cytokines in periodontitis. Single-cell metagenomics revealed 37 bacterial species with cell-specific tropism. Fluorescence in situ hybridization detected intracellular 16 S and mRNA signals of multiple species and correlated with SK/JK proinflammatory phenotypes in situ. Cell-cell communication analysis predicted keratinocyte-specific innate and adaptive immune interactions. Highly multiplexed immunofluorescence (33-antibody) revealed peri-epithelial immune foci, with innate cells often spatially constrained around JKs. Spatial phenotyping revealed immunosuppressed JK-microniches and SK-localized tertiary lymphoid structures in periodontitis. Here, we demonstrate impacts on and predicted interactomics of SK and JK cells in health and periodontitis, which requires further investigation to support precision periodontal interventions in states of chronic inflammation.

Energetic and structural details of the trigger-loop closing transition in RNA polymerase II.

An evolutionarily conserved element in RNA polymerase II, the trigger loop (TL), has been suggested to play an important role in the elongation rate, fidelity of selection of the matched nucleoside triphosphate (NTP), catalysis of transcription elongation, and translocation in both eukaryotes and prokaryotes. In response to NTP binding, the TL undergoes large conformational changes to switch between distinct open and closed states to tighten the active site and avail catalysis. A computational strategy for characterizing the conformational transition pathway is presented to bridge the open and closed states of the TL. Information from a large number of independent all-atom molecular dynamics trajectories from Hamiltonian replica exchange and targeted molecular dynamics simulations is gathered together to assemble a connectivity map of the conformational transition. The results show that with a cognate NTP, TL closing should be a spontaneous process. One major intermediate state is identified along the conformational transition pathway, and the key structural features are characterized. The complete pathway from the open TL to the closed TL provides a clear picture of the TL closing.

Synthesis of [m.n]cyclophanes: regiochemistry transfer from vinyl halides to cyclophanes via Fischer carbene complexes.

The double benzannulation of bis-carbene complexes of chromium with α,ω-diynes generates [m.n]cyclophanes in which all three rings are generated in a single reaction. This triple annulation process is very flexible allowing for the construction of symmetrical [n.n]cyclophanes and unsymmetrical [m.n]cyclophanes as well as isomers in which the two benzene rings are both meta bridged or both para bridged, and isomers that contain both meta and para bridges. The connectivity patterns of the bridges in the cyclophanes can be controlled by regioselectivity transfer from the bis-vinyl carbene complexes in which the substitution pattern of the vinyl groups in the carbene complexes dictate the connectivity pattern in the [m.n]cyclophanes. This synthesis of [n.n]cyclophanes is quite flexible with regard to ring size and can be used with tether lengths ranging from n=2 to n=16 and thus to ring sizes with up to 40 member rings. The only limitation to regioselectivity transfer from the carbene complexes to the [m.n]cyclophanes was found in the synthesis of para,para-cyclophanes with four carbon tethers for which the loss of fidelity occurred with the unexpected formation of meta,para-cyclophanes.

A spatial human thymus cell atlas mapped to a continuous tissue axis.

T cells develop from circulating precursors, which enter the thymus and migrate throughout specialised sub-compartments to support maturation and selection. This process starts already in early fetal development and is highly active until the involution of the thymus in adolescence. To map the micro-anatomical underpinnings of this process in pre- vs. post-natal states, we undertook a spatially resolved analysis and established a new quantitative morphological framework for the thymus, the Cortico-Medullary Axis. Using this axis in conjunction with the curation of a multimodal single-cell, spatial transcriptomics and high-resolution multiplex imaging atlas, we show that canonical thymocyte trajectories and thymic epithelial cells are highly organised and fully established by post-conception week 12, pinpoint TEC progenitor states, find that TEC subsets and peripheral tissue genes are associated with Hassall's Corpuscles and uncover divergence in the pace and drivers of medullary entry between CD4 vs. CD8 T cell lineages. These findings are complemented with a holistic toolkit for spatial analysis and annotation, providing a basis for a detailed understanding of T lymphocyte development.

Transcriptomic profiling of Escherichia coli K-12 in response to a compendium of stressors.

Environmental perturbations impact multiple cellular traits, including gene expression. Bacteria respond to these stressful situations through complex gene interaction networks, thereby inducing stress tolerance and survival of cells. In this paper, we study the response mechanisms of E. coli when exposed to different environmental stressors via differential expression and co-expression analysis. Gene co-expression networks were generated and analyzed via Weighted Gene Co-expression Network Analysis (WGCNA). Based on the gene co-expression networks, genes with similar expression profiles were clustered into modules. The modules were analysed for identification of hub genes, enrichment of biological processes and transcription factors. In addition, we also studied the link between transcription factors and their differentially regulated targets to understand the regulatory mechanisms involved. These networks validate known gene interactions and provide new insights into genes mediating transcriptional regulation in specific stress environments, thus allowing for in silico hypothesis generation.

Data on RNA-seq analysis of the oviducts of five closely related species genus Littorina (Mollusca, Caenogastropoda): L. saxatilis, L. arcana, L. compressa, L. obtusata, L. fabalis.

In the evolution of invertebrates, the transition from egg-layers to brooders occurred many times. However, the molecular mechanisms underlying this transition are still not well understood. Recently diverged species genus Littorina (Mollusca, Gastropoda, Caenogastropoda, Littorinimorpha): Littorina saxatilis, L. arcana, L. compressa, L. obtusata and L. fabalis might be a fruitful model for elucidation of these mechanisms. All five species sympatrically inhabit an intertidal zone. Only L. saxatilis is ovoviviparous while the other four species form clutches. Although in L. saxatilis jelly gland of the pallial oviduct function as a brood pouch, it is not deeply modified at the morphological level in comparison to egg-laying relatives. Comparative analysis of transcriptomic profiles of the pallial oviducts of these closely related species might help to uncover the molecular mechanisms of the egg-laying to brooding transition. Unraveling of the mechanisms underlying this transition in L. saxatilis is important not only in aspects of reproduction biology and strategy, but also in a broader view as an example of relatively fast evolutionary transformations. We generated an RNA-seq dataset (224 104 446 clean reads) for oviducts of five species genus Littorina. Libraries of all five species were sequenced using Illumina HiSeq 2500; additional reads for L. arcana were obtained using Illumina NovaSeq 6000. Transcriptomic profiles were analyzed in pooled samples (of three individuals) with two biological replicates for each species (each biological replicate was prepared and sequenced as a separate library). The transcriptome was assembled de novo and annotated with five assembles corresponding to each species. The raw data were uploaded to the SRA database, the BioProject IDs are PRJNA662103 ("obtusata" group) and PRJNA707549 ("saxatilis" group).