Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control.

Rapid advances in high-throughput DNA sequencing technologies have enabled large-scale whole genome sequencing (WGS) studies. Before performing association analysis between phenotypes and genotypes, preprocessing and quality control (QC) of the raw sequence data need to be performed. Because many biostatisticians have not been working with WGS data so far, we first sketch Illumina's short-read sequencing technology. Second, we explain the general preprocessing pipeline for WGS studies. Third, we provide an overview of important QC metrics, which are applied to WGS data: on the raw data, after mapping and alignment, after variant calling, and after multisample variant calling. Fourth, we illustrate the QC with the data from the GENEtic SequencIng Study Hamburg-Davos (GENESIS-HD), a study involving more than 9000 human whole genomes. All samples were sequenced on an Illumina NovaSeq 6000 with an average coverage of 35× using a PCR-free protocol. For QC, one genome in a bottle (GIAB) trio was sequenced in four replicates, and one GIAB sample was successfully sequenced 70 times in different runs. Fifth, we provide empirical data on the compression of raw data using the DRAGEN original read archive (ORA). The most important quality metrics in the application were genetic similarity, sample cross-contamination, deviations from the expected Het/Hom ratio, relatedness, and coverage. The compression ratio of the raw files using DRAGEN ORA was 5.6:1, and compression time was linear by genome coverage. In summary, the preprocessing, joint calling, and QC of large WGS studies are feasible within a reasonable time, and efficient QC procedures are readily available.

Variant effect predictors: a systematic review and practical guide.

Large-scale association analyses using whole-genome sequence data have become feasible, but understanding the functional impacts of these associations remains challenging. Although many tools are available to predict the functional impacts of genetic variants, it is unclear which tool should be used in practice. This work provides a practical guide to assist in selecting appropriate tools for variant annotation. We conducted a MEDLINE search up to November 10, 2023, and included tools that are applicable to a broad range of phenotypes, can be used locally, and have been recently updated. Tools were categorized based on the types of variants they accept and the functional impacts they predict. Sequence Ontology terms were used for standardization. We identified 118 databases and software packages, encompassing 36 variant types and 161 functional impacts. Combining only three tools, namely SnpEff, FAVOR, and SparkINFERNO, allows predicting 99 (61%) distinct functional impacts. Thirty-seven tools predict 89 functional impacts that are not supported by any other tool, while 75 tools predict pathogenicity and can be used within the ACMG/AMP guidelines in a clinical context. We launched a website allowing researchers to select tools based on desired variants and impacts. In summary, more than 100 tools are already available to predict approximately 160 functional impacts. About 60% of the functional impacts can be predicted by the combination of three tools. Unexpectedly, recent tools do not predict more impacts than older ones. Future research should allow predicting the functionality of so far unsupported variant types, such as gene fusions.URL: https://cardio-care.shinyapps.io/VEP_Finder/ .Registration: OSF Registries on November 10, 2023, https://osf.io/s2gct .

TransposonUltimate: software for transposon classification, annotation and detection.

Most genomes harbor a large number of transposons, and they play an important role in evolution and gene regulation. They are also of interest to clinicians as they are involved in several diseases, including cancer and neurodegeneration. Although several methods for transposon identification are available, they are often highly specialised towards specific tasks or classes of transposons, and they lack common standards such as a unified taxonomy scheme and output file format. We present TransposonUltimate, a powerful bundle of three modules for transposon classification, annotation, and detection of transposition events. TransposonUltimate comes as a Conda package under the GPL-3.0 licence, is well documented and it is easy to install through https://github.com/DerKevinRiehl/TransposonUltimate. We benchmark the classification module on the large TransposonDB covering 891,051 sequences to demonstrate that it outperforms the currently best existing solutions. The annotation and detection modules combine sixteen existing softwares, and we illustrate its use by annotating Caenorhabditis elegans, Rhizophagus irregularis and Oryza sativa subs. japonica genomes. Finally, we use the detection module to discover 29 554 transposition events in the genomes of 20 wild type strains of C. elegans. Databases, assemblies, annotations and further findings can be downloaded from (https://doi.org/10.5281/zenodo.5518085).

Cysteine synthases CYSL-1 and CYSL-2 mediate C. elegans heritable adaptation to P. vranovensis infection.

Parental exposure to pathogens can prime offspring immunity in diverse organisms. The mechanisms by which this heritable priming occurs are largely unknown. Here we report that the soil bacteria Pseudomonas vranovensis is a natural pathogen of the nematode Caenorhabditis elegans and that parental exposure of animals to P. vranovensis promotes offspring resistance to infection. Furthermore, we demonstrate a multigenerational enhancement of progeny survival when three consecutive generations of animals are exposed to P. vranovensis. By investigating the mechanisms by which animals heritably adapt to P. vranovensis infection, we found that parental infection by P. vranovensis results in increased expression of the cysteine synthases cysl-1 and cysl-2 and the regulator of hypoxia inducible factor rhy-1 in progeny, and that these three genes are required for adaptation to P. vranovensis. These observations establish a CYSL-1, CYSL-2, and RHY-1 dependent mechanism by which animals heritably adapt to infection.

Genome sequence of the root-knot nematode Meloidogyne luci.

Root-knot nematodes from the genus Meloidogyne are polyphagous plant endoparasites and agricultural pests of global importance. Here, we report the high-quality genome sequence of Meloidogyne luci population SI-Smartno V13. The resulting genome assembly of M. luci SI-Smartno V13 consists of 327 contigs, with an N50 contig length of 1,711,905 bp and a total assembly length of 209.16 Mb.Root-knot nematodes from the genus Meloidogyne are polyphagous plant endoparasites and agricultural pests of global importance. Here, we report the high-quality genome sequence of Meloidogyne luci population SI-Smartno V13. The resulting genome assembly of M. luci SI-Smartno V13 consists of 327 contigs, with an N50 contig length of 1,711,905 bp and a total assembly length of 209.16 Mb.

Optogenetic Analysis of Depolarization-Dependent Glucagonlike Peptide-1 Release.

Incretin hormones play an important role in the regulation of food intake and glucose homeostasis. Glucagonlike peptide-1 (GLP-1)-secreting cells have been demonstrated to be electrically excitable and to fire action potentials (APs) with increased frequency in response to nutrient exposure. However, nutrients can also be metabolized or activate G-protein-coupled receptors, thus potentially stimulating GLP-1 secretion independent of their effects on the plasma membrane potential. Here we used channelrhodopsins to manipulate the membrane potential of GLUTag cells, a well-established model of GLP-1-secreting enteroendocrine L cells. Using channelrhodopsins with fast or slow on/off kinetics (CheTA and SSFO, respectively), we found that trains of light pulses could trigger APs and calcium elevation in GLUTag cells stably expressing either CheTA or SSFO. Tetrodotoxin reduced light-triggered AP frequency but did not impair calcium responses, whereas further addition of the calcium-channel blockers nifedipine and ω-conotoxin GVIA abolished both APs and calcium transients. Light pulse trains did not trigger GLP-1 secretion from CheTA-expressing cells under basal conditions but were an effective stimulus when cyclic adenosine monophosphate (cAMP) concentrations were elevated by forskolin plus 3-isobutyl 1-methylxanthine. In SSFO-expressing cells, light-stimulated GLP-1 release was observed at resting and elevated cAMP concentrations and was blocked by nifedipine plus ω-conotoxin GVIA but not tetrodotoxin. We conclude that cAMP elevation or cumulative membrane depolarization triggered by SSFO enhances the efficiency of light-triggered action potential firing, voltage-gated calcium entry, and GLP-1 secretion.

Function of bacteriophage G7C esterase tailspike in host cell adsorption.

Bacteriophages recognize and bind to their hosts with the help of receptor-binding proteins (RBPs) that emanate from the phage particle in the form of fibers or tailspikes. RBPs show a great variability in their shapes, sizes, and location on the particle. Some RBPs are known to depolymerize surface polysaccharides of the host while others show no enzymatic activity. Here we report that both RBPs of podovirus G7C - tailspikes gp63.1 and gp66 - are essential for infection of its natural host bacterium E. coli 4s that populates the equine intestinal tract. We characterize the structure and function of gp63.1 and show that unlike any previously described RPB, gp63.1 deacetylates surface polysaccharides of E. coli 4s leaving the backbone of the polysaccharide intact. We demonstrate that gp63.1 and gp66 form a stable complex, in which the N-terminal part of gp66 serves as an attachment site for gp63.1 and anchors the gp63.1-gp66 complex to the G7C tail. The esterase domain of gp63.1 as well as domains mediating the gp63.1-gp66 interaction is widespread among all three families of tailed bacteriophages.