Optimized minimal genome-wide human sgRNA library.

Genome-wide clustered regularly interspaced short palindromic repeats (CRISPR)-based knockout screening is revolting the genetic analysis of a cellular or molecular phenotype in question but is challenged by the large size of single-guide RNA (sgRNA) library. Here we designed a minimal genome-wide human sgRNA library, H-mLib, which is composed of 21,159 sgRNA pairs assembled based on a dedicated selection strategy from all potential SpCas9/sgRNAs in the human genome. These sgRNA pairs were cloned into a dual-gRNA vector each targeting one gene, resulting in a compact library size nearly identical to the number of human protein-coding genes. The performance of the H-mLib was benchmarked to other CRISPR libraries in a proliferation screening conducted in K562 cells. We also identified groups of core essential genes and cell-type specific essential genes by comparing the screening results from the K562 and Jurkat cells. Together, the H-mLib exemplified high specificity and sensitivity in identifying essential genes while containing minimal library complexity, emphasizing its advantages and applications in CRISPR screening with limited cell numbers.

Ultralow-Input Genome Library Preparation for Nanopore Sequencing with Droplet MDA.

Genome sequencing of small species, such as those of meinofauna, can be challenging due to the extremely low input of genomic DNA. While nanopore sequencing is a promising technology for genome assembly due to its limitless long reads, recommended input of 1 µg for the Ligation Sequencing Kit often precludes the use of this technology. Here, I detail an unbiased droplet-based multiple displacement amplification of picogram order of DNA to realize nanopore sequencing with ultralow input of genomic DNA. For this purpose, a microfluidic chip of 10X Genomics Chromium Controller is utilized. With this method, over 10 µg of unbiased amplicons around 10 kbp in length can be obtained from as low as 50 µg of input DNA, which is enough for the construction of multiple sequencing libraries, or for the size selection of longer DNA fragments.

Microsatellite DNA markers for the commercially important fish Hypophthalmus donascimientoi (Pisces: Siluriformes: Pimelodidae) in the Amazon basin: isolation, characterization and amplification in congeneric species.

BACKGROUND: The genus Hypophthalmus comprises six species (H. edentatus, H. marginatus, H. fimbriatus, and H. oremaculatus), and the recently described: H. donascimientoi and H. celiae. The popular name for Hypophthalmus spp. in Brazil is mapará, this name refers to the six species. This group of fish has commercial importance for the states of Amazonas and Pará and, for this reason, requires studies to identify fish stocks. One approach is to use molecular markers, which have been very useful in studies with identification and population analysis of fish. Microsatellite molecular markers (SSRs) are one of the most informative markers for this purpose. There is little populations study of Hypophthalmus using SSRs, and there are less than six loci for the species Hypophthalmus marginatus available in the literature. With the construction of a genomic library of H. donascimientoi, we aimed to isolate and characterize SSRS markers and evaluate the extent of interspecific amplification. METHODS AND RESULTS: A genomic library was constructed with regions enriched of microsatellite for Hypophthalmus donascimientoi. A total of 126 contigs with 42 SSRs were used to design flanking primers for 39 microsatellites. Fifteen loci were characterized in three locations of the Solimões/Amazonas Rivers. The number of alleles ranged from one to 17 with a total of 126 alleles. The mean observed heterozygosity (HO) and expected heterozygosity (HE) were 0.721-0.692, respectively (S.d. HO 0.061 and HE 0.060). Two loci showed significant deviation in the HWE. The PIC ranged from 0.375 to 0.908. Such identified, 12 highly informative loci, and two moderately informative loci. Among the fifteen loci characterized, seven were successfully amplified in four other species of the genus. CONCLUSIONS: The microsatellite showed promise for estimating the genetic variability of H. donascimientoi and can be used as an efficient tool in population analyses of this species and in congeneric species analyzed.

Development of novel microsatellite markers for population differentiation and detection of natural selection in wild populations of butter catfish, Ompok bimaculatus (Bloch, 1794).

BACKGROUND: Butter catfish (Ompok bimaculatus) is a preferred species in South East Asia, with huge aquaculture potential. However, there is limited information about genetic stock composition due to insufficient markers. The goal of this study was to develop de novo microsatellite markers. METHODS AND RESULTS: For sequencing, genomic SMRT bell libraries (1.5 Kbp size) were prepared for O. bimaculatus. A total of 114 SSR containing sequences were used for primer designing. Polymorphic loci were validated by genotyping 83 individuals from four distant riverine populations, viz., Brahmaputra, Bichiya, Gomti and Kaveri. A total of 30 microsatellite loci were polymorphic, of which five were found to be associated with functional genes and eight (four positive and four negative) loci were found to be under selection pressure. A total of 115 alleles were detected in all loci and PIC ranged from 0.539 to 0.927 and pair-wise FST values from 0.1267 to 0.26002 (p < 0.001), with an overall FST value of 0.17047, indicating the presence of population sub-structure. Cross-species transferability of 29 loci (96.67%) was successful in congener species, Ompok pabda. CONCLUSION: The novel SSR markers developed in this study would facilitate stock characterization of natural populations, to be used in future selection breeding programs and planning conservation strategies in these species. Identified non-neutral markers will give insights into the effect of local adaptation on genetic differentiation in the natural population of this species.

Systematic discovery of recombinases for efficient integration of large DNA sequences into the human genome.

Large serine recombinases (LSRs) are DNA integrases that facilitate the site-specific integration of mobile genetic elements into bacterial genomes. Only a few LSRs, such as Bxb1 and PhiC31, have been characterized to date, with limited efficiency as tools for DNA integration in human cells. In this study, we developed a computational approach to identify thousands of LSRs and their DNA attachment sites, expanding known LSR diversity by >100-fold and enabling the prediction of their insertion site specificities. We tested their recombination activity in human cells, classifying them as landing pad, genome-targeting or multi-targeting LSRs. Overall, we achieved up to seven-fold higher recombination than Bxb1 and genome integration efficiencies of 40-75% with cargo sizes over 7 kb. We also demonstrate virus-free, direct integration of plasmid or amplicon libraries for improved functional genomics applications. This systematic discovery of recombinases directly from microbial sequencing data provides a resource of over 60 LSRs experimentally characterized in human cells for large-payload genome insertion without exposed DNA double-stranded breaks.

A pipeline for assembling low copy nuclear markers from plant genome skimming data for phylogenetic use.

Background: Genome skimming is a popular method in plant phylogenomics that do not include a biased enrichment step, relying on random shallow sequencing of total genomic DNA. From these data the plastome is usually readily assembled and constitutes the bulk of phylogenetic information generated in these studies. Despite a few attempts to use genome skims to recover low copy nuclear loci for direct phylogenetic use, such endeavor remains neglected. Causes might include the trade-off between libraries with few reads and species with large genomes (i.e., missing data caused by low coverage), but also might relate to the lack of pipelines for data assembling. Methods: A pipeline and its companion R package designed to automate the recovery of low copy nuclear markers from genome skimming libraries are presented. Additionally, a series of analyses aiming to evaluate the impact of key assembling parameters, reference selection and missing data are presented. Results: A substantial amount of putative low copy nuclear loci was assembled and proved useful to base phylogenetic inference across the libraries tested (4 to 11 times more data than previously assembled plastomes from the same libraries). Discussion: Critical aspects of assembling low copy nuclear markers from genome skims include the minimum coverage and depth of a sequence to be used. More stringent values of these parameters reduces the amount of assembled data and increases the relative amount of missing data, which can compromise phylogenetic inference, in turn relaxing the same parameters might increase sequence error. These issues are discussed in the text, and parameter tuning through multiple comparisons tracking their effects on support and congruence is highly recommended when using this pipeline. The skimmingLoci pipeline (https://github.com/mreginato/skimmingLoci) might stimulate the use of genome skims to recover nuclear loci for direct phylogenetic use, increasing the power of genome skimming data to resolve phylogenetic relationships, while reducing the amount of sequenced DNA that is commonly wasted.

Generating Custom Pooled CRISPR Libraries for Genetic Dissection of Biological Pathways.

Genetic screens are a classic approach to dissecting biological pathways including membrane trafficking. Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 have enabled the utility of this approach in diploid models, including cultured mammalian cells. Here, we present detailed protocols for generating custom CRISPR libraries. These methods are useful for generating genome-wide libraries for new model organisms that lack an existing genome-wide library, and for generating smaller focused libraries.

Validation of a random Vibrio parahaemolyticus genomic library by selection of quinolone resistance in a heterologous host.

Vibrio parahaemolyticus is a shellfish-borne pathogen that is a highly prevalent causative agent of inflammatory gastroenteritis in humans. Genomic libraries have proven useful for the identification of novel gene functions in many bacterial species. In this study we prepared a library containing 40 kb fragments of randomly sheared V. parahaemolyticus genomic DNA and introduced this into Escherichia coli HB101 using a commercially available low copy cosmid system. In order to estimate coverage and suitability of the library and potentially identify novel antimicrobial resistance determinants, we screened for the acquisition of resistance to the fluoroquinolone norfloxacin - a phenotype exhibited by V. parahaemolyticus but not the heterologous E. coli host. Upon selection on solid medium containing norfloxacin, 0.52% of the library population was resistant, consistent with the selection of a single resistance locus. End-sequencing identified six distinct insert fragments. All clones displayed fourfold increased norfloxacin MIC compared with E. coli HB101 carrying an empty vector. The common locus contained within resistant clones included qnr, a previously described quinolone resistance gene. These results indicate that the library was unbiased, of sufficient coverage and that heterologous expression was possible. While we hope that this library proves useful for identifying the genetic determinants of complex phenotypes such as those related to virulence, not all norfloxacin resistance genes were detected in our screen. As such, we discuss the benefits and limitations of this approach for identifying the genetic basis of uncharacterized bacterial phenotypes.

Protocol variations in run-on transcription dataset preparation produce detectable signatures in sequencing libraries.

BACKGROUND: A variety of protocols exist for producing whole genome run-on transcription datasets. However, little is known about how differences between these protocols affect the signal within the resulting libraries. RESULTS: Using run-on transcription datasets generated from the same biological system, we show that a variety of GRO- and PRO-seq preparation methods leave identifiable signatures within each library. Specifically we show that the library preparation method results in differences in quality control metrics, as well as differences in the signal distribution at the 5 ' end of transcribed regions. These shifts lead to disparities in eRNA identification, but do not impact analyses aimed at inferring the key regulators involved in changes to transcription. CONCLUSIONS: Run-on sequencing protocol variations result in technical signatures that can be used to identify both the enrichment and library preparation method of a particular data set. These technical signatures are batch effects that limit detailed comparisons of pausing ratios and eRNAs identified across protocols. However, these batch effects have only limited impact on our ability to infer which regulators underlie the observed transcriptional changes.

Effect of short-term hindlimb immobilization on skeletal muscle atrophy and the transcriptome in a low compared with high responder to endurance training model.

Skeletal muscle atrophy is a physiological response to disuse, aging, and disease. We compared changes in muscle mass and the transcriptome profile after short-term immobilization in a divergent model of high and low responders to endurance training to identify biological processes associated with the early atrophy response. Female rats selectively bred for high response to endurance training (HRT) and low response to endurance training (LRT; n = 6/group; generation 19) underwent 3 day hindlimb cast immobilization to compare atrophy of plantaris and soleus muscles with line-matched controls (n = 6/group). RNA sequencing was utilized to identify Gene Ontology Biological Processes with differential gene set enrichment. Aerobic training performed prior to the intervention showed HRT improved running distance (+60.6 ± 29.6%), while LRT were unchanged (-0.3 ± 13.3%). Soleus atrophy was greater in LRT vs. HRT (-9.0 ±8.8 vs. 6.2 ±8.2%; P<0.05) and there was a similar trend in plantaris (-16.4 ±5.6% vs. -8.5 ±7.4%; P = 0.064). A total of 140 and 118 biological processes were differentially enriched in plantaris and soleus muscles, respectively. Soleus muscle exhibited divergent LRT and HRT responses in processes including autophagy and immune response. In plantaris, processes associated with protein ubiquitination, as well as the atrogenes (Trim63 and Fbxo32), were more positively enriched in LRT. Overall, LRT demonstrate exacerbated atrophy compared to HRT, associated with differential gene enrichments of biological processes. This indicates that genetic factors that result in divergent adaptations to endurance exercise, may also regulate biological processes associated with short-term muscle unloading.

Mapping of Nonhomologous End Joining-Mediated Integration Facilitates Genome-Scale Trackable Mutagenesis in Yarrowia lipolytica.

Genome-scale mutagenesis, phenotypic screening, and tracking the causal mutations is a powerful approach for genetic analysis. However, classic mutagenesis approaches require extensive effort to identify causal mutations. It is desirable to demonstrate a powerful approach for rapid trackable mutagenesis. Here, we mapped the distribution of nonhomologous end joining (NHEJ)-mediated integration for the first time and demonstrated that it can be used for constructing the genome-scale trackable mutagenesis library in Yarrowia lipolytica. The sequencing of 9.15 × 105 insertions showed that NHEJ-mediated integration inserted DNA randomly across the chromosomes, and the transcriptional regulatory regions exhibited integration preference. The insertions were located in both nucleosome-occupancy regions and nucleosome-free regions. Using NHEJ-mediated integration to construct the genome-scale mutagenesis library, the new targets that improved ß-carotene biosynthesis and acetic acid tolerance were identified rapidly. This mutagenesis approach is readily applicable to other organisms with strong NHEJ preference and will contribute to cell factory construction.

Using Genome Scale Mutant Libraries to Identify Essential Genes.

Identification of essential genes is key to understanding the required processes and gene products of organisms under one or more conditions. Transposon sequencing (Tn-seq) has been used to predict essential genes or ones that conditionally impact fitness in a wide variety of organisms. Here, we describe the generation of genome-scale mutant libraries and the analysis of Tn-seq data to identify essential genes from cultures grown in a single condition as well as those that are conditionally important by analyzing the behavior of these mutant libraries in different growth environments. While we illustrate the approach using data derived from Tn-seq analysis of the α-proteobacteria Rhodobacter sphaeroides and Zymomonas mobilis, the protocols and systems we describe should be generally applicable to a variety of organisms.

Development of highly polymorphic microsatellite markers for Campanula glomerata L. (Campanulaceae).

BACKGROUND: Campanula glomerata L. (Campanulaceae) is a self-incompatible, insect-pollinated herb occurring in calcareous grasslands, and is declining and (critically) endangered in many parts of its European distribution range. It often exists as small and isolated populations. A recovery plan of C. glomerata has been implemented in southern Belgium, involving plant translocations. METHODS AND RESULTS: We developed microsatellite markers using an enriched genomic library and characterized 16 loci in 111 individuals from eight populations. These 16 loci were highly polymorphic, with 11 to 31 alleles per locus for a total of 329 alleles, and expected heterozygosity (He) ranging from 0.470 to 0.938. CONCLUSIONS: These highly polymorphic loci constitute a promising tool for detailed genetic analyses: assigning individuals to distinct multilocus genotypes will allow quantifying pollen dispersal, clonal propagation and sexual recruitment and identifying admixed seed progeny and their pollen donors. Evaluating the genetic status of existing populations and a genetic monitoring of the translocated populations will contribute to optimize success in restoring viable and evolutionary resilient populations.

Genomic library preparation and hybridization capture of formalin-fixed tissues and allozyme supernatant for population genomics and considerations for combining capture- and RADseq-based single nucleotide polymorphism data sets.

Until recently many historical museum specimens were largely inaccessible to genomic inquiry, but high-throughput sequencing (HTS) approaches have allowed researchers to successfully sequence genomic DNA from dried and fluid-preserved museum specimens. In addition to preserved specimens, many museums contain large series of allozyme supernatant samples, but the amenability of these samples to HTS has not yet been assessed. Here, we compared the performance of a target-capture approach using alternative sources of genomic DNA from 10 specimens of spring salamanders (Plethodontidae: Gyrinophilus porphyriticus) collected between 1985 and 1990: allozyme supernatants, allozyme homogenate pellets and formalin-fixed tissues. We designed capture probes based on double-digest restriction-site associated sequencing (RADseq) derived loci from frozen blood samples available for seven of the specimens and assessed the success and consistency of capture and RADseq approaches. This study design enabled direct comparisons of data quality and potential biases among the different data sets for phylogenomic and population genomic analyses. We found that in phylogenetic analyses, all enrichment types for a given specimen clustered together. In principal component space all capture-based samples clustered together, but RADseq samples did not cluster with corresponding capture-based samples. Single nucleotide polymorphism calls were on average 18.3% different between enrichment types for a given individual, but these discrepancies were primarily due to differences in heterozygous/homozygous single nucleotide polymorphism calls. We demonstrate that both allozyme supernatant and formalin-fixed samples can be successfully used for population genomic analyses and we discuss ways to identify and reduce biases associated with combining capture and RADseq data.

Construction of Genomic Library and Screening of Edwardsiella tarda Immunogenic Proteins for Their Protective Efficacy Against Edwardsiellosis.

Edwardsiella tarda is a severe aquaculture pathogen that can infect many hosts including humans, animals, and fish. Timely diagnosis and treatment are crucial for the control of edwardsiellosis in the aqua industry. By using rabbit polyclonal antibody, an expression gene library of virulent Edwardsiella tarda strain ED-BDU 1 isolated in south India was constructed and screened. The identified immune expressive proteins were characterized, and the corresponding coding sequences were cloned, expressed, and the purified recombinant proteins were used as antigens. The identified immunoreactive proteins namely HflC, HflK, and YhcI were studied for their immune protective potential in vivo by challenge experiments. The protective efficacy of HflC, HflK, and YhcI showed that the clearance of Edwardsiella from the host with ~ 60% survivability. Further, the immunoreactive proteins induce a strong immune response upon infection and elicit the significant production of IL-10, IFN-γ, Th1, and Th2 mediated mRNA expression and were therefore effective in vaccine production for edwardsiellosis.

Susceptibility of Field-Collected Nyssorhynchus darlingi to Plasmodium spp. in Western Amazonian Brazil.

Mosquito susceptibility to Plasmodium spp. infection is of paramount importance for malaria occurrence and sustainable transmission. Therefore, understanding the genetic features underlying the mechanisms of susceptibility traits is pivotal to assessing malaria transmission dynamics in endemic areas. The aim of this study was to investigate the susceptibility of Nyssorhynchus darlingi-the dominant malaria vector in Brazil-to Plasmodium spp. using a reduced representation genome-sequencing protocol. The investigation was performed using a genome-wide association study (GWAS) to identify mosquito genes that are predicted to modulate the susceptibility of natural populations of the mosquito to Plasmodium infection. After applying the sequence alignment protocol, we generated the variant panel and filtered variants; leading to the detection of 202,837 SNPs in all specimens analyzed. The resulting panel was used to perform GWAS by comparing the pool of SNP variants present in Ny. darlingi infected with Plasmodium spp. with the pool obtained in field-collected mosquitoes with no evidence of infection by the parasite (all mosquitoes were tested separately using RT-PCR). The GWAS results for infection status showed two statistically significant variants adjacent to important genes that can be associated with susceptibility to Plasmodium infection: Cytochrome P450 (cyp450) and chitinase. This study provides relevant knowledge on malaria transmission dynamics by using a genomic approach to identify mosquito genes associated with susceptibility to Plasmodium infection in Ny. darlingi in western Amazonian Brazil.

Protocol for design, construction, and selection of genome phage (gPhage) display libraries.

This protocol describes the genomic phage (gPhage) display platform, a large-scale antigen and epitope mapping technique. We constructed a gPhage display peptide library of a eukaryotic organism, Trypanosoma cruzi (causative agent of Chagas disease), to map the antibody response landscape against the parasite. Here, we used an organism with a relatively large but intronless genome, although future applications could include other prevalent or (re)emerging infectious organisms carrying genomes with a limited number of introns. For complete details on the use and execution of this protocol, please refer to Teixeira et al. (2021).

Construction of a 5 immune-related lncRNA-based prognostic model of NSCLC via bioinformatics.

ABSTRACT: Participate in tumorigenic, oncogenic, and tumor suppressive pathways through gene expression regulation. We aimed to build an immune-related long noncoding RNA (lncRNA) prognostic model to enhance nonsmall cell lung cancer (NSCLC) prognostic prediction.The original data were collected from the cancer genome atlas database. Perl and R software were used for statistical analysis. The effects of lncRNAs expression on prognosis were analyzed by Gene Expression Profiling Interactive Analysis. Silico functional analysis were performed by DAVID Bioinformatics Resources.The median risk score as a dividing value separated patients into high- and low-risk groups. These 2 groups had different 5-year survival rates, median survival times, and immune statuses. The 5-lncRNA signature was validated as an independent prognostic factor with high accuracy (area under the receiver operating characteristic = 0.722). Silico functional analysis connected the lncRNAs with immune-related biological processes and pathways in carcinogenesis.The novel immune-related lncRNA prognostic model had significant clinical implication for enhancing lung adenocarcinoma outcome prediction and guiding the choice of treatment.

Reconstructing and counting genomic fragments through tagmentation-based haploid phasing.

Single-cell sequencing provides a new level of granularity in studying the heterogeneous nature of cancer cells. For some cancers, this heterogeneity is the result of copy number changes of genes within the cellular genomes. The ability to accurately determine such copy number changes is critical in tracing and understanding tumorigenesis. Current single-cell genome sequencing methodologies infer copy numbers based on statistical approaches followed by rounding decimal numbers to integer values. Such methodologies are sample dependent, have varying calling sensitivities which heavily depend on the sample's ploidy and are sensitive to noise in sequencing data. In this paper we have demonstrated the concept of integer-counting by using a novel bioinformatic algorithm built on our library construction chemistry in order to detect the discrete nature of the genome.