Scaling read aligners to hundreds of threads on general-purpose processors

Challenges and Approaches to Genotyping Repetitive DNA

Individuals within a species can exhibit vast variation in copy number of repetitive DNA elements. This variation may contribute to complex traits such as lifespan and disease, yet it is only infrequently considered in genotype-phenotype associations. Although the possible importance of copy number variation is widely recognized, accurate copy number quantification remains challenging. Here, we assess the technical reproducibility of several major methods for copy number estimation as they apply to the large repetitive ribosomal DNA array (rDNA). rDNA encodes the ribosomal RNAs and exists as a tandem gene array in all eukaryotes. Repeat units of rDNA are kilobases in size, often with several hundred units comprising the array, making rDNA particularly intractable to common quantification techniques. We evaluate pulsed-field gel electrophoresis, droplet digital PCR, and Nextera-based whole genome sequencing as approaches to copy number estimation, comparing techniques across model organisms and spanning wide ranges of copy numbers. Nextera-based whole genome sequencing, though commonly used in recent literature, produced high error. We explore possible causes for this error and provide recommendations for best practices in rDNA copy number estimation. We present a resource of high-confidence rDNA copy number estimates for a set of S. cerevisiae and C. elegans strains for future use. We furthermore explore the possibility for FISH-based copy number estimation, an alternative that could potentially characterize copy number on a cellular level.

A native chromatin immunoprecipitation (ChIP) protocol for studying histone modifications in strawberry fruits

BACKGROUND

Covalent modifications of histones and histone variants have great influence on chromatin structure, which is involved in the transcriptional regulation of gene expression. Chromatin immunoprecipitation (ChIP) is a powerful tool for studying in vivo DNA-histone interactions. Strawberry is a model for Rosaceae and non-climacteric fruits, in which histone modifications have been implicated to affect fruit development and ripening. However, a validated ChIP method has not been reported in strawberry, probably due to its high levels of polysaccharides which affect the quality of prepared chromatin and the efficiency of immunoprecipitation.

RESULTS

We describe a native chromatin immunoprecipitation (N-ChIP) protocol suitable for strawberry by optimizing the parameters for nuclei isolation, chromatin extraction, DNA fragmentation and validation analysis using quantitative real-time PCR (qRT-PCR). The qRT-PCR results show that both the active mark H3K36me3 and the silent mark H3K9me2 are efficiently immunoprecipitated for the enriched regions. Compared to X-ChIP (cross-linked chromatin followed by immunoprecipitation), our optimized N-ChIP procedure has a higher signal-to-noise ratio and a lower background for both the active and the silent histone modifications. Furthermore, high-throughput sequencing following N-ChIP demonstrates that nearly 90% of the enriched H3K9/K14ac peaks are overlapped between biological replicates, indicating its remarkable consistency and reproducibility.

CONCLUSIONS

An N-ChIP method suitable for the fleshy fruit tissues of woodland strawberry Fragaria vesca is described in this study. The efficiency and reproducibility of our optimized N-ChIP protocol are validated by both qRT-PCR and high-throughput sequencing. We conclude that N-ChIP is a more suitable method for strawberry fruit tissues relative to X-ChIP, which could be combined with high-throughput sequencing to investigate the impact of histone modifications in strawberry and potentially in other fruits with high content of polysaccharides.

Small-Molecule Antiviral β-d-N4-Hydroxycytidine Inhibits a Proofreading-Intact Coronavirus with a High Genetic Barrier to Resistance

Coronaviruses (CoVs) have emerged from animal reservoirs to cause severe and lethal disease in humans, but there are currently no FDA-approved antivirals to treat the infections. One class of antiviral compounds, nucleoside analogues, mimics naturally occurring nucleosides to inhibit viral replication. While these compounds have been successful therapeutics for several viral infections, mutagenic nucleoside analogues, such as ribavirin and 5-fluorouracil, have been ineffective at inhibiting CoVs. This has been attributed to the proofreading activity of the viral 3'-5' exoribonuclease (ExoN). β-d-N4-Hydroxycytidine (NHC) (EIDD-1931; Emory Institute for Drug Development) has recently been reported to inhibit multiple viruses. Here, we demonstrate that NHC inhibits both murine hepatitis virus (MHV) (50% effective concentration [EC50] = 0.17 μM) and Middle East respiratory syndrome CoV (MERS-CoV) (EC50 = 0.56 μM) with minimal cytotoxicity. NHC inhibited MHV lacking ExoN proofreading activity similarly to wild-type (WT) MHV, suggesting an ability to evade or overcome ExoN activity. NHC inhibited MHV only when added early during infection, decreased viral specific infectivity, and increased the number and proportion of G:A and C:U transition mutations present after a single infection. Low-level NHC resistance was difficult to achieve and was associated with multiple transition mutations across the genome in both MHV and MERS-CoV. These results point to a virus-mutagenic mechanism of NHC inhibition in CoVs and indicate a high genetic barrier to NHC resistance. Together, the data support further development of NHC for treatment of CoVs and suggest a novel mechanism of NHC interaction with the CoV replication complex that may shed light on critical aspects of replication.IMPORTANCE The emergence of coronaviruses (CoVs) into human populations from animal reservoirs has demonstrated their epidemic capability, pandemic potential, and ability to cause severe disease. However, no antivirals have been approved to treat these infections. Here, we demonstrate the potent antiviral activity of a broad-spectrum ribonucleoside analogue, β-d-N4-hydroxycytidine (NHC), against two divergent CoVs. Viral proofreading activity does not markedly impact sensitivity to NHC inhibition, suggesting a novel interaction between a nucleoside analogue inhibitor and the CoV replicase. Further, passage in the presence of NHC generates only low-level resistance, likely due to the accumulation of multiple potentially deleterious transition mutations. Together, these data support a mutagenic mechanism of inhibition by NHC and further support the development of NHC for treatment of CoV infections.

Global transcriptome analysis of different stages of preimplantation embryo development in river buffalo

Background

Water buffalo (Bubalus bubalis) are divided into river buffalo and swamp buffalo subspecies and are essential livestock for agriculture and the local economy. Studies on buffalo reproduction have primarily focused on optimal fertility and embryonic mortality. There is currently limited knowledge on buffalo embryonic development, especially during the preimplantation period. Assembly of the river buffalo genome offers a reference for omics studies and facilitates transcriptomic analysis of preimplantation embryo development (PED).

Methods

We revealed transcriptomic profile of four stages (2-cell, 8-cell, Morula and Blastocyst) of PED via RNA-seq (Illumina HiSeq4000). Each stage comprised three biological replicates. The data were analyzed according to the basic RNA-seq analysis process. Ingenuity analysis of cell lineage control, especially transcription factor (TF) regulatory networks, was also performed.

Results

A total of 21,519 expressed genes and 67,298 transcripts were predicted from approximately 81.94 Gb of raw data. Analysis of transcriptome-wide expression, gene coexpression networks, and differentially expressed genes (DEGs) allowed for the characterization of gene-specific expression levels and relationships for each stage. The expression patterns of TFs, such as POU5F1, TEAD4, CDX4 and GATAs, were elucidated across diverse time series; most TF expression levels were increased during the blastocyst stage, during which time cell differentiation is initiated. All of these TFs were involved in the composition of the regulatory networks that precisely specify cell fate. These findings offer a deeper understanding of PED at the transcriptional level in the river buffalo.

Long-term stability and Red Queen-like strain dynamics in marine viruses

Viruses that infect microorganisms dominate marine microbial communities numerically, with impacts ranging from host evolution to global biogeochemical cycles^1,2. However, virus community dynamics, necessary for conceptual and mechanistic model development, remains difficult to assess. Here, we describe the long-term stability of a viral community by analysing the metagenomes of near-surface 0.02-0.2 μm samples from the San Pedro Ocean Time-series³ that were sampled monthly over 5 years. Of 19,907 assembled viral contigs (>5 kb, mean 15 kb), 97% were found in each sample (by >98% ID metagenomic read recruitment) to have relative abundances that ranged over seven orders of magnitude, with limited temporal reordering of rank abundances along with little change in richness. Seasonal variations in viral community composition were superimposed on the overall stability; maximum community similarity occurred at 12-month intervals. Despite the stability of viral genotypic clusters that had 98% sequence identity, viral sequences showed transient variations in single-nucleotide polymorphisms (SNPs) and constant turnover of minor population variants, each rising and falling over a few months, reminiscent of Red Queen dynamics⁴. The rise and fall of variants within populations, interpreted through the perspective of known virus-host interactions⁵, is consistent with the hypothesis that fluctuating selection acts on a microdiverse cloud of strains, and this succession is associated with ever-shifting virus-host defences and counterdefences. This results in long-term virus-host coexistence that is facilitated by perpetually changing minor variants.

Improved metagenomic analysis with Kraken 2

Although Kraken's k-mer-based approach provides a fast taxonomic classification of metagenomic sequence data, its large memory requirements can be limiting for some applications. Kraken 2 improves upon Kraken 1 by reducing memory usage by 85%, allowing greater amounts of reference genomic data to be used, while maintaining high accuracy and increasing speed fivefold. Kraken 2 also introduces a translated search mode, providing increased sensitivity in viral metagenomics analysis.

Improved metagenomic analysis with Kraken 2

Although Kraken’s k-mer-based approach provides a fast taxonomic classification of metagenomic sequence data, its large memory requirements can be limiting for some applications. Kraken 2 improves upon Kraken 1 by reducing memory usage by 85%, allowing greater amounts of reference genomic data to be used, while maintaining high accuracy and increasing speed fivefold. Kraken 2 also introduces a translated search mode, providing increased sensitivity in viral metagenomics analysis.

The Indispensable Role of Histone Methyltransferase PoDot1 in Extracellular Glycoside Hydrolase Biosynthesis of Penicillium oxalicum

Histone methylation is associated with transcription regulation, but its role for glycoside hydrolase (GH) biosynthesis is still poorly understood. We identified the histone H3 lysine 79 (H3K79)-specific methyltransferase PoDot1 in Penicillium oxalicum. PoDot1 affects conidiation by regulating the transcription of key regulators (BrlA, FlbC, and StuA) of asexual development and is required in normal hyphae septum and branch formation by regulating the transcription of five septin-encoding genes, namely, aspA, aspB, aspC, aspD, and aspE. Tandem affinity purification/mass spectrometry showed that PoDot1 has no direct interaction with transcription machinery, but it affects the expressions of extracellular GH genes extensively. The expression of genes (amy15A, amy13A, cel7A/cbh1, cel61A, chi18A, cel3A/bgl1, xyn10A, cel7B/eg1, cel5B/eg2, and cel6A/cbh2) that encode the top 10 GHs was remarkably downregulated by Podot1 deletion (ΔPodot1). Consistent with the decrease in gene transcription level, the activities of amylases and cellulases were significantly decreased in ΔPodot1 mutants in agar (solid) and fermentation (liquid) media. The repression of GH gene expressions caused by PoDot1 deletion was not mediated by key transcription factors, such as AmyR, ClrB, CreA, and XlnR, but was accompanied by defects in global demethylated H3K79 (H3K79me2) and trimethylated H3K79 (H3K79me3). The impairment of H3K79me2 on specific GH gene loci was observed due to PoDot1 deletion. The results implies that defects of H3K79 methylation is the key reason of the downregulated transcription level of GH-encoding genes and reveals the indispensable role of PoDot1 in extracellular GH biosynthesis.

Alternative Sigma Factor B in Bovine Mastitis-Causing Staphylococcus aureus: Characterization of Its Role in Biofilm Formation, Resistance to Hydrogen Peroxide Stress, Regulon Members

This study examines treatments of the bacterial pathogen Staphylococcus aureus, namely, in the context of its being a major cause of subclinical bovine mastitis. Such infections caused by S. aureus among dairy cows are difficult to detect and can easily become chronic, leading to reduced productivity and large losses for dairy manufacturers. In this study, the role of alternative sigma factor B (σ^B), which has been shown to be a global regulator for S. aureus infections, was explored in a mastitis-causing S. aureus strain, RF122. For comparison with the wild-type strain, a sigB null (ΔsigB) mutant was constructed and analyzed for its phenotypes and transcriptome. Our study found that σ^B is essential for biofilm formation as the ΔsigB mutant strain produced significantly less biofilm than did the wild-type strain at 48 h. σ^B is involved in response to H₂O₂ stress. However, σ^B plays a minor or no role in resistance to antiseptics (e.g., povidone-iodine and chlorhexidine), resistance to tested antibiotics, hemolysin activity, and invasion ability. RNA sequencing identified 225 σ^B-dependent genes, of which 171 are positively regulated and 54 are negatively regulated. The identified genes are involved in stress response, pathogenesis, and metabolic mechanisms. Quantitative TaqMan RT-PCR was performed to verify the RNA sequencing results; i.e., σ^B is a positive regulator for asp23, sarA, katA, yabJ, sodA, SAB2006c, and nrdD expressions. In the RF122 strain, σ^B plays a role in biofilm formation, general stress response (e.g., H₂O₂), and regulation of virulence factors and virulence-associated genes.

Has gene expression neofunctionalization in the fire ant antennae contributed to queen discrimination behavior?

Queen discrimination behavior in the fire ant Solenopsis invicta maintains its two types of societies: colonies with one (monogyne) or many (polygyne) queens, yet the underlying genetic mechanism is poorly understood. This behavior is controlled by two supergene alleles, SB and Sb, with ~600 genes. Polygyne workers, having either the SB/SB or SB/Sb genotype, accept additional SB/Sb queens into their colonies but kill SB/SB queens. In contrast, monogyne workers, all SB/SB, reject all additional queens regardless of genotype. Because the SB and Sb alleles have suppressed recombination, determining which genes within the supergene mediate this differential worker behavior is difficult. We hypothesized that the alternate worker genotypes sense queens differently because of the evolution of differential expression of key genes in their main sensory organ, the antennae. To identify such genes, we sequenced RNA from four replicates of pooled antennae from three classes of workers: monogyne SB/SB, polygyne SB/SB, and polygyne SB/Sb. We identified 81 differentially expressed protein-coding genes with 13 encoding potential chemical metabolism or perception proteins. We focused on the two odorant perception genes: an odorant receptor SiOR463 and an odorant-binding protein SiOBP12. We found that SiOR463 has been lost in the Sb genome. In contrast, SiOBP12 has an Sb-specific duplication, SiOBP12b', which is expressed in the SB/Sb worker antennae, while both paralogs are expressed in the body. Comparisons with another fire ant species revealed that SiOBP12b' antennal expression is specific to S. invicta and suggests that queen discrimination may have evolved, in part, through expression neofunctionalization.

Carbon source regulates polysaccharide capsule biosynthesis in Streptococcus pneumoniae

The exopolysaccharide capsule of Streptococcus pneumoniae is an important virulence factor, but the mechanisms that regulate capsule thickness are not fully understood. Here, we investigated the effects of various exogenously supplied carbohydrates on capsule production and gene expression in several pneumococcal serotypes. Microscopy analyses indicated a near absence of the capsular polysaccharide (CPS) when S. pneumoniae was grown on fructose. Moreover, serotype 7F pneumococci produced much less CPS than strains of other serotypes (6B, 6C, 9V, 15, and 23F) when grown on glucose or sucrose. RNA-sequencing revealed carbon source-dependent regulation of distinct genes of WT strains and capsule-switch mutants of serotypes 6B and 7F, but could not explain the mechanism of capsule thickness regulation. In contrast, ³¹P NMR of whole-cell extract from capsule-knockout strains (Δcps) clearly revealed the accumulation or absence of capsule precursor metabolites when cells were grown on glucose or fructose, respectively. This finding suggests that fructose uptake mainly results in intracellular fructose 1-phosphate, which is not converted to CPS precursors. In addition, serotype 7F strains accumulated more precursors than did 6B strains, indicating less efficient conversion of precursor metabolites into the CPS in 7F, in line with its thinner capsule. Finally, isotopologue sucrose labeling and NMR analyses revealed that the uptake of the labeled fructose subunit into the capsule is <10% that of glucose. Our findings on the effects of carbon sources on CPS production in different S. pneumoniae serotypes may contribute to a better understanding of pneumococcal diseases and could inform future therapeutic approaches.

A genome-wide screen identifies IRF2 as a key regulator of caspase-4 in human cells

Caspase-4, the cytosolic LPS sensor, and gasdermin D, its downstream effector, constitute the non-canonical inflammasome, which drives inflammatory responses during Gram-negative bacterial infections. It remains unclear whether other proteins regulate cytosolic LPS sensing, particularly in human cells. Here, we conduct a genome-wide CRISPR/Cas9 screen in a human monocyte cell line to identify genes controlling cytosolic LPS-mediated pyroptosis. We find that the transcription factor, IRF2, is required for pyroptosis following cytosolic LPS delivery and functions by directly regulating caspase-4 levels in human monocytes and iPSC-derived monocytes. CASP4, GSDMD, and IRF2 are the only genes identified with high significance in this screen highlighting the simplicity of the non-canonical inflammasome. Upon IFN-γ priming, IRF1 induction compensates IRF2 deficiency, leading to robust caspase-4 expression. Deficiency in IRF2 results in dampened inflammasome responses upon infection with Gram-negative bacteria. This study emphasizes the central role of IRF family members as specific regulators of the non-canonical inflammasome.

Sustained Wolbachia-mediated blocking of dengue virus isolates following serial passage in Aedes aegypti cell culture

Wolbachia is an intracellular endosymbiont of insects that inhibits the replication of a range of pathogens in its arthropod hosts. The release of Wolbachia into wild populations of mosquitoes is an innovative biocontrol effort to suppress the transmission of arthropod-borne viruses (arboviruses) to humans, most notably dengue virus. The success of the Wolbachia-based approach hinges upon the stable persistence of the ‘pathogen blocking’ effect, whose mechanistic basis is poorly understood. Evidence suggests that Wolbachia may affect viral replication via a combination of competition for host resources and activation of host immunity. The evolution of resistance against Wolbachia and pathogen blocking in the mosquito or the virus could reduce the public health impact of the symbiont releases. Here, we investigate if dengue 3 virus (DENV-3) is capable of accumulating adaptive mutations that improve its replicative capacity during serial passage in Wolbachia wMel-infected cells. During the passaging regime, viral isolates in Wolbachia-infected cells exhibited greater variation in viral loads compared to controls. The viral loads of these isolates declined rapidly during passaging due to the blocking effects of Wolbachia carriage, with several being lost all together and the remainder recovering to low but stable levels. We attempted to sequence the genomes of the surviving passaged isolates but, given their low abundance, were unable to obtain sufficient depth of coverage for evolutionary analysis. In contrast, viral loads in Wolbachia-free control cells were consistently high during passaging. The surviving isolates passaged in the presence of Wolbachia exhibited a reduced ability to replicate even in Wolbachia-free cells. These experiments demonstrate the challenge for dengue in evolving resistance to Wolbachia-mediated blocking.

Variants in KIAA0825 underlie autosomal recessive postaxial polydactyly

Postaxial polydactyly (PAP) is a common limb malformation that often leads to cosmetic and functional complications. Molecular evaluation of polydactyly can serve as a tool to elucidate genetic and signaling pathways that regulate limb development, specifically, the anterior-posterior specification of the limb. To date, only five genes have been identified for nonsyndromic PAP: FAM92A, GLI1, GLI3, IQCE and ZNF141. In this study, two Pakistani multiplex consanguineous families with autosomal recessive nonsyndromic PAP were clinically and molecularly evaluated. From both pedigrees, a DNA sample from an affected member underwent exome sequencing. In each family, we identified a segregating frameshift (c.591dupA [p.(Q198Tfs*21)]) and nonsense variant (c.2173A > T [p.(K725*)]) in KIAA0825 (also known as C5orf36). Although KIAA0825 encodes a protein of unknown function, it has been demonstrated that its murine ortholog is expressed during limb development. Our data contribute to the establishment of a catalog of genes important in limb patterning, which can aid in diagnosis and obtaining a better understanding of the biology of polydactyly.

Complete Genome Sequences of Human Astrovirus Prototype Strains (Types 1 to 8)

We report the complete genome sequences of the eight human astrovirus Oxford prototype strains. These sequences share 94.9% to 99.9% nucleotide identity with open reading frame 2 (ORF2) genes of astrovirus genomes previously deposited in GenBank and include the first complete genome of human astrovirus type 7.

We report the complete genome sequences of the eight human astrovirus Oxford prototype strains. These sequences share 94.9% to 99.9% nucleotide identity with open reading frame 2 (ORF2) genes of astrovirus genomes previously deposited in GenBank and include the first complete genome of human astrovirus type 7.

Complete Genome Sequence of the Lumpy Skin Disease Virus Recovered from the First Outbreak in the Northern Caucasus Region of Russia in 2015

We report here the complete genome sequence of a lumpy skin disease virus (LSDV) isolate obtained in the Northern Caucasus region of Russia in 2015. The LSDV/Russia/Dagestan/2015 genome sequence grouped with field LSDV isolates found in Serbia and Greece, suggesting the monophyletic origin of LSDV isolates that recently affected countries in the Northern Hemisphere.

Comparative effect of allopolyploidy on transposable element composition and gene expression between Gossypium hirsutum and its two diploid progenitors

An allopolyploidization event formed allotetraploid Gossypium species from an A-genome diploid species and a D-genome diploid species. To explore the responses of transposable elements (TEs) to allopolyploidy, we assembled parallel TE datasets from G. hirsutum, G. arboreum and G. raimondii and analyzed the TE types and the effects of TEs on orthologous gene expression in the three Gossypium genomes. Gypsy was the most abundant TE type and most TEs were located ∼500 bp from genes in all three genomes. In G. hirsutum, 35.6% of genes harbored TE insertions, whereas insertions were more frequent in G. arboreum and G. raimondii. G. hirsutum had the highest proportion of uniquely matching 24-nt small interfering RNAs (siRNAs) that targeted TEs. TEs, particularly those targeted by 24-nt siRNAs, were associated with reduced gene expression, but the effect of TEs on orthologous gene expression varied substantially among species. Orthologous gene expression levels in G. hirsutum were intermediate between those of G. arboreum and G. raimondii, which did not experience TE expansion or reduction resulting from allopolyploidization. This study underscores the diversity of TEs co-opted by host genes and provides insights into the roles of TEs in regulating gene expression in Gossypium.