Parameters for accurate genome alignment

Improved assembly of the Pungitius pungitius reference genome

The nine-spined stickleback (Pungitius pungitius) has been increasingly used as a model system in studies of local adaptation and sex chromosome evolution but its current reference genome assembly is far from perfect, lacking distinct sex chromosomes. We generated an improved assembly of the nine-spined stickleback reference genome (98.3% BUSCO completeness) with the aid of linked-read mapping. While the new assembly (v8) was of similar size as the earlier version (v7), we were able to assign 4.4 times more contigs to the linkage groups and improve the contiguity of the genome. Moreover, the new assembly contains a ∼22.8 Mb Y-linked scaffold (LG22) consisting mainly of previously assigned X-contigs, putative Y-contigs, putative centromere contigs, and highly repetitive elements. The male individual showed an even mapping depth on LG12 (pseudo X chromosome) and LG22 (Y-linked scaffold) in the segregating sites, suggesting near-pure X and Y representation in the v8 assembly. A total of 26,803 genes were annotated, and about 33% of the assembly was found to consist of repetitive elements. The high proportion of repetitive elements in LG22 (53.10%) suggests it can be difficult to assemble the complete sequence of the species' Y chromosome. Nevertheless, the new assembly is a significant improvement over the previous version and should provide a valuable resource for genomic studies of stickleback fishes.

An FPGA-based hardware accelerator supporting sensitive sequence homology filtering with profile hidden Markov models

BACKGROUND

Sequence alignment lies at the heart of genome sequence annotation. While the BLAST suite of alignment tools has long held an important role in alignment-based sequence database search, greater sensitivity is achieved through the use of profile hidden Markov models (pHMMs). Here, we describe an FPGA hardware accelerator, called HAVAC, that targets a key bottleneck step (SSV) in the analysis pipeline of the popular pHMM alignment tool, HMMER.

RESULTS

The HAVAC kernel calculates the SSV matrix at 1739 GCUPS on a ∼  $3000 Xilinx Alveo U50 FPGA accelerator card, ∼  227× faster than the optimized SSV implementation in nhmmer. Accounting for PCI-e data transfer data processing, HAVAC is 65× faster than nhmmer's SSV with one thread and 35× faster than nhmmer with four threads, and uses ∼  31% the energy of a traditional high end Intel CPU.

CONCLUSIONS

HAVAC demonstrates the potential offered by FPGA hardware accelerators to produce dramatic speed gains in sequence annotation and related bioinformatics applications. Because these computations are performed on a co-processor, the host CPU remains free to simultaneously compute other aspects of the analysis pipeline.

Chromosome-level genome assembly of Plagiognathops microlepis based on PacBio HiFi and Hi-C sequencing

Plagiognathops microlepis is an economic freshwater fish in the subfamily Xenocyprinae of Cyprinidae. It is widely distributed in the freshwater ecosystem of China, with moderate economic value and broad development prospects. However, the lack of genomic resources has limited our understanding on the genetic basis, phylogenetic status and adaptive evolution strategies of this fish. Here, we assembled a chromosome-level reference genome of P. microlepis by integrating Pacbio HiFi long-reads, Illumina short-reads and Hi-C sequencing data. The size of this genome is 1004.34 Mb with a contig N50 of 38.80 Mb. Using Hi-C sequencing data, 99.59% of the assembled sequences were further anchored to 24 chromosomes. A total of 578.91 Mb repeat sequences and 28,337 protein-coding genes were predicted in the current genome, of which, 26,929 genes were functionally annotated. This genome provides valuable information for investigating the phylogeny and evolutionary history of cyprinid fishes, as well as the genetic basis of adaptive strategies and special traits in P. microlepis.

ULTRA-Effective Labeling of Repetitive Genomic Sequence

In the age of long read sequencing, genomics researchers now have access to accurate repetitive DNA sequence (including satellites) that, due to the limitations of short read sequencing, could previously be observed only as unmappable fragments. Tools that annotate repetitive sequence are now more important than ever, so that we can better understand newly uncovered repetitive sequences, and also so that we can mitigate errors in bioinformatic software caused by those repetitive sequences. To that end, we introduce the 1.0 release of our tool for identifying and annotating locally-repetitive sequence, ULTRA (ULTRA Locates Tandemly Repetitive Areas). ULTRA is fast enough to use as part of an efficient annotation pipeline, produces state-of-the-art reliable coverage of repetitive regions containing many mutations, and provides interpretable statistics and labels for repetitive regions. It released under an open license, and available for download at https://github.com/TravisWheelerLab/ULTRA.

Mechanobiological stimulation in organ-on-a-chip systems reduces hepatic drug metabolic capacity in favor of regenerative specialization

Hepatic physiology depends on the liver's complex structural composition which among others, provides high oxygen supply rates, locally differential oxygen tension, endothelial paracrine signaling, as well as residual hemodynamic shear stress to resident hepatocytes. While functional improvements were shown by implementing these factors into hepatic culture systems, direct cause-effect relationships are often not well characterized-obfuscating their individual contribution in more complex microphysiological systems. By comparing increasingly complex hepatic in vitro culture systems that gradually implement these parameters, we investigate the influence of the cellular microenvironment to overall hepatic functionality in pharmacological applications. Here, hepatocytes were modulated in terms of oxygen tension and supplementation, endothelial coculture, and exposure to fluid shear stress delineated from oxygen influx. Results from transcriptomic and metabolomic evaluation indicate that particularly oxygen supply rates are critical to enhance cellular functionality-with cellular drug metabolism remaining comparable to physiological conditions after prolonged static culture. Endothelial signaling was found to be a major contributor to differential phenotype formation known as metabolic zonation, indicated by WNT pathway activity. Lastly, oxygen-delineated shear stress was identified to direct cellular fate towards increased hepatic plasticity and regenerative phenotypes at the cost of drug metabolic functionality - in line with regenerative effects observed in vivo. With these results, we provide a systematic evaluation of critical parameters and their impact in hepatic systems. Given their adherence to physiological effects in vivo, this highlights the importance of their implementation in biomimetic devices, such as organ-on-a-chip systems. Considering recent advances in basic liver biology, direct translation of physiological structures into in vitro models is a promising strategy to expand the capabilities of pharmacological models.

Ten common issues with reference sequence databases and how to mitigate them

Metagenomic sequencing has revolutionized our understanding of microbiology. While metagenomic tools and approaches have been extensively evaluated and benchmarked, far less attention has been given to the reference sequence database used in metagenomic classification. Issues with reference sequence databases are pervasive. Database contamination is the most recognized issue in the literature; however, it remains relatively unmitigated in most analyses. Other common issues with reference sequence databases include taxonomic errors, inappropriate inclusion and exclusion criteria, and sequence content errors. This review covers ten common issues with reference sequence databases and the potential downstream consequences of these issues. Mitigation measures are discussed for each issue, including bioinformatic tools and database curation strategies. Together, these strategies present a path towards more accurate, reproducible and translatable metagenomic sequencing.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

Secondary Contact, Introgressive Hybridization, and Genome Stabilization in Sticklebacks

Advances in genomic studies have revealed that hybridization in nature is pervasive and raised questions about the dynamics of different genetic and evolutionary factors following the initial hybridization event. While recent research has proposed that the genomic outcomes of hybridization might be predictable to some extent, many uncertainties remain. With comprehensive whole-genome sequence data, we investigated the genetic introgression between 2 divergent lineages of 9-spined sticklebacks (Pungitius pungitius) in the Baltic Sea. We found that the intensity and direction of selection on the introgressed variation has varied across different genomic elements: while functionally important regions displayed reduced rates of introgression, promoter regions showed enrichment. Despite the general trend of negative selection, we identified specific genomic regions that were enriched for introgressed variants, and within these regions, we detected footprints of selection, indicating adaptive introgression. Geographically, we found the selection against the functional changes to be strongest in the vicinity of the secondary contact zone and weaken as a function of distance from the initial contact. Altogether, the results suggest that the stabilization of introgressed variation in the genomes is a complex, multistage process involving both negative and positive selection. In spite of the predominance of negative selection against introgressed variants, we also found evidence for adaptive introgression variants likely associated with adaptation to Baltic Sea environmental conditions.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet

BACKGROUND

Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult.

RESULTS

Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPβ for methylated motifs and the preference of c-Myc for unmethylated E-box motifs.

CONCLUSIONS

Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

Comparative Genomic Analysis of Bacterial Data in BV-BRC: An Example Exploring Antimicrobial Resistance

As genomic and related data continue to expand, research biologists are often hampered by the computational hurdles required to analyze their data. The National Institute of Allergy and Infectious Diseases (NIAID) established the Bioinformatics Resource Centers (BRC) to assist researchers with their analysis of genome sequence and other omics-related data. Recently, the PAThosystems Resource Integration Center (PATRIC), the Influenza Research Database (IRD), and the Virus Pathogen Database and Analysis Resource (ViPR) BRCs merged to form the Bacterial and Viral Bioinformatics Resource Center (BV-BRC) at https://www.bv-brc.org/ . The combined BV-BRC leverages the functionality of the original resources for bacterial and viral research communities with a unified data model, enhanced web-based visualization and analysis tools, and bioinformatics services. Here we demonstrate how antimicrobial resistance data can be analyzed in the new resource.

An evolutionarily distinct ringed seal in the Ilulissat Icefjord

The Earth's polar regions are low rates of inter- and intraspecific diversification. An extreme mammalian example is the Arctic ringed seal (Pusa hispida hispida), which is assumed to be panmictic across its circumpolar Arctic range. Yet, local Inuit communities in Greenland and Canada recognize several regional variants; a finding supported by scientific studies of body size variation. It is however unclear whether this phenotypic variation reflects plasticity, morphs or distinct ecotypes. Here, we combine genomic, biologging and survey data, to document the existence of a unique ringed seal ecotype in the Ilulissat Icefjord (locally 'Kangia'), Greenland; a UNESCO World Heritage site, which is home to the most productive marine-terminating glacier in the Arctic. Genomic analyses reveal a divergence of Kangia ringed seals from other Arctic ringed seals about 240 kya, followed by secondary contact since the Last Glacial Maximum. Despite ongoing gene flow, multiple genomic regions appear under strong selection in Kangia ringed seals, including candidate genes associated with pelage coloration, growth and osmoregulation, potentially explaining the Kangia seal's phenotypic and behavioural uniqueness. The description of 'hidden' diversity and adaptations in yet another Arctic species merits a reassessment of the evolutionary processes that have shaped Arctic diversity and the traditional view of this region as an evolutionary freezer. Our study highlights the value of indigenous knowledge in guiding science and calls for efforts to identify distinct populations or ecotypes to understand how these might respond differently to environmental change.

How genomic insights into the evolutionary history of clouded leopards inform their conservation

Clouded leopards (Neofelis spp.), a morphologically and ecologically distinct lineage of big cats, are severely threatened by habitat loss and fragmentation, targeted hunting, and other human activities. The long-held poor understanding of their genetics and evolution has undermined the effectiveness of conservation actions. Here, we report a comprehensive investigation of the whole genomes, population genetics, and adaptive evolution of Neofelis. Our results indicate the genus Neofelis arose during the Pleistocene, coinciding with glacial-induced climate changes to the distributions of savannas and rainforests, and signatures of natural selection associated with genes functioning in tooth, pigmentation, and tail development, associated with clouded leopards' unique adaptations. Our study highlights high-altitude adaptation as the main factor driving nontaxonomic population differentiation in Neofelis nebulosa. Population declines and inbreeding have led to reduced genetic diversity and the accumulation of deleterious variation that likely affect reproduction of clouded leopards, highlighting the urgent need for effective conservation efforts.

Homoeologous evolution of the allotetraploid genome of Poa annua L

BACKGROUND

Poa annua (annual bluegrass) is an allotetraploid turfgrass, an agronomically significant weed, and one of the most widely dispersed plant species on earth. Here, we report the chromosome-scale genome assemblies of P. annua's diploid progenitors, P. infirma and P. supina, and use multi-omic analyses spanning all three species to better understand P. annua's evolutionary novelty.

RESULTS

We find that the diploids diverged from their common ancestor 5.5 - 6.3 million years ago and hybridized to form P. annua ≤ 50,000 years ago. The diploid genomes are similar in chromosome structure and most notably distinguished by the divergent evolutionary histories of their transposable elements, leading to a 1.7 × difference in genome size. In allotetraploid P. annua, we find biased movement of retrotransposons from the larger (A) subgenome to the smaller (B) subgenome. We show that P. annua's B subgenome is preferentially accumulating genes and that its genes are more highly expressed. Whole-genome resequencing of several additional P. annua accessions revealed large-scale chromosomal rearrangements characterized by extensive TE-downsizing and evidence to support the Genome Balance Hypothesis.

CONCLUSIONS

The divergent evolutions of the diploid progenitors played a central role in conferring onto P. annua its remarkable phenotypic plasticity. We find that plant genes (guided by selection and drift) and transposable elements (mostly guided by host immunity) each respond to polyploidy in unique ways and that P. annua uses whole-genome duplication to purge highly parasitized heterochromatic sequences. The findings and genomic resources presented here will enable the development of homoeolog-specific markers for accelerated weed science and turfgrass breeding.

Genomic Epidemiological Analysis of Antimicrobial-Resistant Bacteria with Nanopore Sequencing

Antimicrobial-resistant (AMR) bacterial infections caused by clinically important bacteria, including ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) and mycobacteria (Mycobacterium tuberculosis and nontuberculous mycobacteria), have become a global public health threat. Their epidemic and pandemic clones often accumulate useful accessory genes in their genomes, such as AMR genes (ARGs) and virulence factor genes (VFGs). This process is facilitated by horizontal gene transfer among microbial communities via mobile genetic elements (MGEs), such as plasmids and phages. Nanopore long-read sequencing allows easy and inexpensive analysis of complex bacterial genome structures, although some aspects of sequencing data calculation and genome analysis methods are not systematically understood. Here we describe the latest and most recommended experimental and bioinformatics methods available for the construction of complete bacterial genomes from nanopore sequencing data and the detection and classification of genotypes of bacterial chromosomes, ARGs, VFGs, plasmids, and other MGEs based on their genomic sequences for genomic epidemiological analysis of AMR bacteria.

Stepwise use of genomics and transcriptomics technologies increases diagnostic yield in Mendelian disorders

Purpose: Multi-omics offer worthwhile and increasingly accessible technologies to diagnostic laboratories seeking potential second-tier strategies to help patients with unresolved rare diseases, especially patients clinically diagnosed with a rare OMIM (Online Mendelian Inheritance in Man) disease. However, no consensus exists regarding the optimal diagnostic care pathway to adopt after negative results with standard approaches. Methods: In 15 unsolved individuals clinically diagnosed with recognizable OMIM diseases but with negative or inconclusive first-line genetic results, we explored the utility of a multi-step approach using several novel omics technologies to establish a molecular diagnosis. Inclusion criteria included a clinical autosomal recessive disease diagnosis and single heterozygous pathogenic variant in the gene of interest identified by first-line analysis (60%-9/15) or a clinical diagnosis of an X-linked recessive or autosomal dominant disease with no causative variant identified (40%-6/15). We performed a multi-step analysis involving short-read genome sequencing (srGS) and complementary approaches such as mRNA sequencing (mRNA-seq), long-read genome sequencing (lrG), or optical genome mapping (oGM) selected according to the outcome of the GS analysis. Results: SrGS alone or in combination with additional genomic and/or transcriptomic technologies allowed us to resolve 87% of individuals by identifying single nucleotide variants/indels missed by first-line targeted tests, identifying variants affecting transcription, or structural variants sometimes requiring lrGS or oGM for their characterization. Conclusion: Hypothesis-driven implementation of combined omics technologies is particularly effective in identifying molecular etiologies. In this study, we detail our experience of the implementation of genomics and transcriptomics technologies in a pilot cohort of previously investigated patients with a typical clinical diagnosis without molecular etiology.

Ancient homomorphy of molluscan sex chromosomes sustained by reversible sex-biased genes and sex determiner translocation

Contrary to classic theory prediction, sex-chromosome homomorphy is prevalent in the animal kingdom but it is unclear how ancient homomorphic sex chromosomes avoid chromosome-scale degeneration. Molluscs constitute the second largest, Precambrian-originated animal phylum and have ancient, uncharacterized homomorphic sex chromosomes. Here, we profile eight genomes of the bivalve mollusc family of Pectinidae in a phylogenetic context and show 350 million years sex-chromosome homomorphy, which is the oldest known sex-chromosome homomorphy in the animal kingdom, far exceeding the ages of well-known heteromorphic sex chromosomes such as 130-200 million years in mammals, birds and flies. The long-term undifferentiation of molluscan sex chromosomes is potentially sustained by the unexpected intertwined regulation of reversible sex-biased genes, together with the lack of sexual dimorphism and occasional sex chromosome turnover. The pleiotropic constraint of regulation of reversible sex-biased genes is widely present in ancient homomorphic sex chromosomes and might be resolved in heteromorphic sex chromosomes through gene duplication followed by subfunctionalization. The evolutionary dynamics of sex chromosomes suggest a mechanism for 'inheritance' turnover of sex-determining genes that is mediated by translocation of a sex-determining enhancer. On the basis of these findings, we propose an evolutionary model for the long-term preservation of homomorphic sex chromosomes.

In Vitro Studies Reveal that Pseudomonas, from Odontotermes obesus Colonies, can Function as a Defensive Mutualist as it Prevents the Weedy Fungus While Keeping the Crop Fungus Unaffected

Insects that farm monocultures of fungi are canonical examples of nutritional symbiosis as well as independent evolution of agriculture in non-human animals. But just like in human agriculture, these fungal crops face constant threat of invasion by weeds which, if unchecked, take over the crop fungus. In fungus-growing termites, the crop fungus (Termitomyces) faces such challenges from the weedy fungus Pseudoxylaria. The mechanism by which Pseudoxylaria is suppressed is not known. However, evidence suggests that some bacterial secondary symbionts can serve as defensive mutualists by preventing the growth of Pseudoxylaria. However, such secondary symbionts must possess the dual, yet contrasting, capabilities of suppressing the weedy fungus while keeping the growth of the crop fungus unaffected. This study describes the isolation, identification, and culture-dependent estimation of the roles of several such putative defensive mutualists from the colonies of the wide-spread fungus-growing termite from India, Odontotermes obesus. From the 38 bacterial cultures tested, a strain of Pseudomonas showed significantly greater suppression of the weedy fungus than the crop fungus. Moreover, a 16S rRNA pan-microbiome survey, using the Nanopore platform, revealed Pseudomonas to be a part of the core microbiota of O. obesus. A meta-analysis of microbiota composition across different species of Odontotermes also confirms the widespread prevalence of Pseudomonas within this termite. These lines of evidence indicate that Pseudomonas could be playing the role of defensive mutualist within Odontotermes.

Application of MinION Amplicon Sequencing to Buccal Swab Samples for Improving Resolution and Throughput of Rumen Microbiota Analysis

The Illumina MiSeq platform has been widely used as a standard method for studying the rumen microbiota. However, the low resolution of taxonomic identification is the only disadvantage of MiSeq amplicon sequencing, as it targets a part of the 16S rRNA gene. In the present study, we performed three experiments to establish a high-resolution and high-throughput rumen microbial profiling approach using a combination of MinION platform and buccal swab sample, which is a proxy for rumen contents. In experiment 1, rumen contents and buccal swab samples were collected simultaneously from cannulated cattle (n = 6) and used for microbiota analysis using three different analytical workflows: amplicon sequencing of the V3–V4 region of the 16S rRNA gene using MiSeq and amplicon sequencing of near full-length 16S rRNA gene using MinION or PacBio Sequel II. All reads derived from the MinION and PacBio platforms were classified at the species-level. In experiment 2, rumen fluid samples were collected from beef cattle (n = 28) and used for 16S rRNA gene amplicon sequencing using the MinION platform to evaluate this sequencing platform for rumen microbiota analysis. We confirmed that the MinION platform allowed species-level taxa assignment for the predominant bacterial groups, which were previously identified at the family- and genus-level using the MiSeq platform. In experiment 3, buccal swab samples were collected from beef cattle (n = 30) and used for 16S rRNA gene amplicon sequencing using the MinION platform to validate the applicability of a combination of the MinION platform and buccal swab samples for rumen microbiota analysis. The distribution of predominant bacterial taxa in the buccal swab samples was similar to that in the rumen samples observed in experiment 2. Based on these results, we concluded that the combination of the MinION platform and buccal swab samples may be potentially applied for rumen microbial analysis in large-scale studies.

QitanTech Nanopore Long-Read Sequencing Enables Rapid Resolution of Complete Genomes of Multi-Drug Resistant Pathogens

Advancement of novel sequencing technologies facilitates modern life science and medicine unprecedentedly. Exploring complete genome sequences of bacteria by long-read sequencing technology is significant for microbial genomics research. However, third-generation long-read sequencing technologies are available with limited choices, which generate technological barrier to scientific research. Recently, a novel QitanTech nanopore long-read sequencing technology has emerged in China, but the potential application and performance were unexplored. Herein, we comprehensively evaluated the feasibility of the emerging sequencing technology in assembling complete genomes of MDR pathogens. The results showed that 500 Mbp QitanTech nanopore sequencing data could be generated within 8 h in one flow cell with the standard library preparation method. The mean read length, longest read length, and mean read-level accuracy of QitanTech sequencing data were 6,041 bp, 57,037 bp, and 81.50% (LAST)/81.40% (Minimap2), respectively. Two routine assembly strategies including long-read assembly and hybrid assembly enable the achievement of complete bacterial genomes. The accuracy of assembled draft bacterial genomes with QitanTech long-read data could be improved up to 99.9% dramatically by polishing using accurate short-read data. Furthermore, the assembled bacterial genomes cover accurate structures of complex resistance plasmids harboring critical resistance genes such as tet(X), tmexCD-toprJ, and bla_VIM–2, even the complex fusion MDR plasmid generated from homologous recombination. In conclusion, QitanTech nanopore sequencing, as a nanopore long-read sequencing technology launched in China, could be a good option for investigation of complex bacterial genomes. More potential applications based on this novel platform warrant investigations.

Characterization of UDP-glycosyltransferase family members reveals how major flavonoid glycoside accumulates in the roots of Scutellaria baicalensis

Background

Flavonoid glycosides extracted from roots of Scutellaria baicalensis exhibit strong pharmaceutical antitumor, antioxidative, anti-inflammatory, and antiviral activities. UDP glycosyltransferase (UGT) family members are responsible for the transfer of a glycosyl moiety from UDP sugars to a wide range of acceptor flavonoids. Baicalin is the major flavonoid glycoside found in S. baicalensis roots, and its aglycone baicalein is synthesized from a specially evolved pathway that has been elucidated. However, it is necessary to carry out a genome-wide study of genes involved in 7-O-glucuronidation, the final biosynthesis step of baicalin, which might elucidate the relationship between the enzymes and the metabolic accumulation patterns in this medicinal plant.

Results

We reported the phylogenetic analysis, tissue-specific expression, biochemical characterization and evolutionary analysis of glucosyltransferases (SbUGTs) and glucuronosyltransferases (SbUGATs) genes based on the recently released genome of S. baicalensis. A total of 124 UGTs were identified, and over one third of them were highly expressed in roots. In vitro enzyme assays showed that 6 SbUGTs could use UDP-glucose as a sugar donor and convert baicalein to oroxin A (baicalein 7-O-glucoside), while 4 SbUGATs used only UDP-glucuronic acid as the sugar donor and catalyzed baicalein to baicalin. SbUGAT4 and SbUGT2 are the most highly expressed SbUGAT and SbUGT genes in root tissues, respectively. Kinetic measurements revealed that SbUGAT4 had a lower Km value and higher Vmax/Km ratio to baicalein than those of SbUGT2. Furthermore, tandem duplication events were detected in SbUGTs and SbUGATs.

Conclusions

This study demonstrated that glucosylation and glucuronidation are two major glycosylated decorations in the roots of S. baicalensis. Higher expression level and affinity to substrate of SbUGAT4, and expansion of this gene family contribute high accumulation of baicalin in the root of S. baicalensis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08391-1.

MycoSNP: A Portable Workflow for Performing Whole-Genome Sequencing Analysis of Candida auris

Candida auris is an urgent public health threat characterized by high drug-resistant rates and rapid spread in healthcare settings worldwide. As part of the C. auris response, molecular surveillance has helped public health officials track the global spread and investigate local outbreaks. Here, we describe whole-genome sequencing analysis methods used for routine C. auris molecular surveillance in the United States; methods include reference selection, reference preparation, quality assessment and control of sequencing reads, read alignment, and single-nucleotide polymorphism calling and filtration. We also describe the newly developed pipeline MycoSNP, a portable workflow for performing whole-genome sequencing analysis of fungal organisms including C. auris.

Concerted genomic and epigenomic changes accompany stabilization of Arabidopsis allopolyploids

During evolution successful allopolyploids must overcome 'genome shock' between hybridizing species but the underlying process remains elusive. Here, we report concerted genomic and epigenomic changes in resynthesized and natural Arabidopsis suecica (TTAA) allotetraploids derived from Arabidopsis thaliana (TT) and Arabidopsis arenosa (AA). A. suecica shows conserved gene synteny and content with more gene family gain and loss in the A and T subgenomes than respective progenitors, although A. arenosa-derived subgenome has more structural variation and transposon distributions than A. thaliana-derived subgenome. These balanced genomic variations are accompanied by pervasive convergent and concerted changes in DNA methylation and gene expression among allotetraploids. The A subgenome is hypomethylated rapidly from F1 to resynthesized allotetraploids and convergently to the T-subgenome level in natural A. suecica, despite many other methylated loci being inherited from F1 to all allotetraploids. These changes in DNA methylation, including small RNAs, in allotetraploids may affect gene expression and phenotypic variation, including flowering, silencing of self-incompatibility and upregulation of meiosis- and mitosis-related genes. In conclusion, concerted genomic and epigenomic changes may improve stability and adaptation during polyploid evolution.

Divergent Genomic Adaptations in the Microbiomes of Arctic Subzero Sea-Ice and Cryopeg Brines

Subzero hypersaline brines are liquid microbial habitats within otherwise frozen environments, where concentrated dissolved salts prevent freezing. Such extreme conditions presumably require unique microbial adaptations, and possibly altered ecologies, but specific strategies remain largely unknown. Here we examined prokaryotic taxonomic and functional diversity in two seawater-derived subzero hypersaline brines: first-year sea ice, subject to seasonally fluctuating conditions; and ancient cryopeg, under relatively stable conditions geophysically isolated in permafrost. Overall, both taxonomic composition and functional potential were starkly different. Taxonomically, sea-ice brine communities (∼10⁵ cells mL^–1) had greater richness, more diversity and were dominated by bacterial genera, including Polaribacter, Paraglaciecola, Colwellia, and Glaciecola, whereas the more densely inhabited cryopeg brines (∼10⁸ cells mL^–1) lacked these genera and instead were dominated by Marinobacter. Functionally, however, sea ice encoded fewer accessory traits and lower average genomic copy numbers for shared traits, though DNA replication and repair were elevated; in contrast, microbes in cryopeg brines had greater genetic versatility with elevated abundances of accessory traits involved in sensing, responding to environmental cues, transport, mobile elements (transposases and plasmids), toxin-antitoxin systems, and type VI secretion systems. Together these genomic features suggest adaptations and capabilities of sea-ice communities manifesting at the community level through seasonal ecological succession, whereas the denser cryopeg communities appear adapted to intense bacterial competition, leaving fewer genera to dominate with brine-specific adaptations and social interactions that sacrifice some members for the benefit of others. Such cryopeg genomic traits provide insight into how long-term environmental stability may enable life to survive extreme conditions.

SVLR: Genome Structural Variant Detection Using Long-Read Sequencing Data

Genome structural variants (SVs) have great impacts on human phenotype and diversity, and have been linked to numerous diseases. Long-read sequencing technologies arise to make it possible to find SVs of as long as 10,000 nucleotides. Thus, long read-based SV detection has been drawing attention of many recent research projects, and many tools have been developed for long reads to detect SVs recently. In this article, we present a new method, called SVLR, to detect SVs based on long-read sequencing data. Comparing with existing methods, SVLR can detect three new kinds of SVs: block replacements, block interchanges, and translocations. Although these new SVs are structurally more complicated, SVLR achieves accuracies that are comparable with those of the classic SVs. Moreover, for the classic SVs that can be detected by state-of-the-art methods (e.g., SVIM and Sniffles), our experiments demonstrate recall improvements of up to 38% without harming the precisions (i.e., >78%). We also point out three directions to further improve SV detection in the future. Source codes: https://github.com/GWYSDU/SVLR.

Insights into amphicarpy from the compact genome of the legume Amphicarpaea edgeworthii

Amphicarpy (seed heteromorphy) is a unique and fascinating reproductive strategy wherein a single plant produces both aerial and subterranean fruits. This strategy is believed to be an adaptation to life under stressful or uncertain environments. Here, we sequenced and de novo assembled a chromosome-level genome assembly of the legume Amphicarpaea edgeworthii Benth. The 299-Mb A. edgeworthii genome encodes 27 899 protein-coding genes and is the most compact sequenced legume genome reported until date. Its reduced genome size may be attributed to the reduced long-terminal repeat retrotransposon content, which stems from the unequal homologous recombination. Gene families related to immunity and stress resistance have been contracted in A. edgeworthii, which is consistent with the notion that the amphicarpic reproductive strategy may be a complementary mechanism for its weak environmental-adaptation ability. We demonstrated the 'ABCE' model for the differentiation of chasmogamous and cleistogamous flowers. In addition, the characteristics of aerial and subterranean seeds in hard-seededness were explored. Thus, we suggest that the A. edgeworthii genome, which is the first of an amphicarpic plant, offers significant insights into its unusual reproductive strategy that is a key resource towards comprehending the evolution of angiosperms.

The genome of Tripterygium wilfordii and characterization of the celastrol biosynthesis pathway

Tripterygium wilfordii is a vine from the Celastraceae family that is used in traditional Chinese medicine (TCM). The active ingredient, celastrol, is a friedelane-type pentacyclic triterpenoid with putative roles as an antitumor, immunosuppressive, and anti-obesity agent. Here, we report a reference genome assembly of T. wilfordii with high-quality annotation using a hybrid sequencing strategy. The total genome size obtained is 340.12 Mb, with a contig N50 value of 3.09 Mb. We successfully anchored 91.02% of sequences into 23 pseudochromosomes using high-throughput chromosome conformation capture (Hi–C) technology. The super-scaffold N50 value was 13.03 Mb. We also annotated 31,593 structural genes, with a repeat percentage of 44.31%. These data demonstrate that T. wilfordii diverged from Malpighiales species approximately 102.4 million years ago. By integrating genome, transcriptome and metabolite analyses, as well as in vivo and in vitro enzyme assays of two cytochrome P450 (CYP450) genes, TwCYP712K1 and TwCYP712K2, it is possible to investigate the second biosynthesis step of celastrol and demonstrate that this was derived from a common ancestor. These data provide insights and resources for further investigation of pathways related to celastrol, and valuable information to aid the conservation of resources, as well as understand the evolution of Celastrales.

Million-year-old DNA sheds light on the genomic history of mammoths

Temporal genomic data hold great potential for studying evolutionary processes, including speciation. However, sampling across speciation events would in many cases require genomic time series that stretch well into the Early Pleistocene (>1 million years). Although theoretical models suggest that DNA should survive on this timescale¹, the oldest genomic data recovered so far is from a 560-780 ka old horse specimen². Here we report the recovery of genome-wide data from three Early and Middle Pleistocene mammoth specimens, two of which are more than one million years old. We find that two distinct mammoth lineages were present in eastern Siberia during the Early Pleistocene. One of these gave rise to the woolly mammoth, whereas the other represents a previously unrecognised lineage that was ancestral to the first mammoths to colonise North America. Our analyses reveal that the North American Columbian mammoth traces its ancestry to a Middle Pleistocene hybridisation between these two lineages, with roughly equal admixture proportions. Finally, we show that the majority of protein-coding changes associated with cold adaptation in woolly mammoths were present already a million years ago. These findings highlight the potential of deep time palaeogenomics to expand our understanding of speciation and long-term adaptive evolution.

Comparative analysis of nuclear, chloroplast, and mitochondrial genomes of watermelon and melon provides evidence of gene transfer

During plant evolution, there is genetic communication between organelle and nuclear genomes. A comparative analysis was performed on the organelle and nuclear genomes of the watermelon and melon. In the watermelon, chloroplast-derived sequences accounted for 7.6% of the total length of the mitochondrial genome. In the melon, chloroplast-derived sequences accounted for approximately 2.73% of the total mitochondrial genome. In watermelon and melon, the chloroplast-derived small-fragment sequences are either a subset of large-fragment sequences or appeared multiple times in the mitochondrial genome, indicating that these fragments may have undergone multiple independent migration integrations or emerged in the mitochondrial genome after migration, replication, and reorganization. There was no evidence of migration from the mitochondria to chloroplast genome. A sequence with a total length of about 73 kb (47%) in the watermelon chloroplast genome was homologous to a sequence of about 313 kb in the nuclear genome. About 33% of sequences in the watermelon mitochondrial genome was homologous with a 260 kb sequence in the nuclear genome. A sequence with a total length of about 38 kb (25%) in the melon chloroplast genome was homologous with 461 sequences in the nuclear genome, with a total length of about 301 kb. A 3.4 Mb sequence in the nuclear genome was homologous with a melon mitochondrial sequence. These results indicate that, during the evolution of watermelon and melon, a large amount of genetic material was exchanged between the nuclear genome and the two organelle genomes in the cytoplasm.

Analysis of a phase-variable restriction modification system of the human gut symbiont Bacteroides fragilis

The genomes of gut Bacteroidales contain numerous invertible regions, many of which contain promoters that dictate phase-variable synthesis of surface molecules such as polysaccharides, fimbriae, and outer surface proteins. Here, we characterize a different type of phase-variable system of Bacteroides fragilis, a Type I restriction modification system (R-M). We show that reversible DNA inversions within this R-M locus leads to the generation of eight specificity proteins with distinct recognition sites. In vitro grown bacteria have a different proportion of specificity gene combinations at the expression locus than bacteria isolated from the mammalian gut. By creating mutants, each able to produce only one specificity protein from this region, we identified the R-M recognition sites of four of these S-proteins using SMRT sequencing. Transcriptome analysis revealed that the locked specificity mutants, whether grown in vitro or isolated from the mammalian gut, have distinct transcriptional profiles, likely creating different phenotypes, one of which was confirmed. Genomic analyses of diverse strains of Bacteroidetes from both host-associated and environmental sources reveal the ubiquity of phase-variable R-M systems in this phylum.

Karyotypic Evolution of Sauropsid Vertebrates Illuminated by Optical and Physical Mapping of the Painted Turtle and Slider Turtle Genomes

Recent sequencing and software enhancements have advanced our understanding of the evolution of genomic structure and function, especially addressing novel evolutionary biology questions. Yet fragmentary turtle genome assemblies remain a challenge to fully decipher the genetic architecture of adaptive evolution. Here, we use optical mapping to improve the contiguity of the painted turtle (Chrysemys picta) genome assembly and use de novo fluorescent in situ hybridization (FISH) of bacterial artificial chromosome (BAC) clones, BAC-FISH, to physically map the genomes of the painted and slider turtles (Trachemys scripta elegans). Optical mapping increased C. picta's N50 by ~242% compared to the previous assembly. Physical mapping permitted anchoring ~45% of the genome assembly, spanning 5544 genes (including 20 genes related to the sex determination network of turtles and vertebrates). BAC-FISH data revealed assembly errors in C. picta and T. s. elegans assemblies, highlighting the importance of molecular cytogenetic data to complement bioinformatic approaches. We also compared C. picta's anchored scaffolds to the genomes of other chelonians, chicken, lizards, and snake. Results revealed a mostly one-to-one correspondence between chromosomes of painted and slider turtles, and high homology among large syntenic blocks shared with other turtles and sauropsids. Yet, numerous chromosomal rearrangements were also evident across chelonians, between turtles and squamates, and between avian and non-avian reptiles.

Meta-omics Provides Insights into the Impact of Hydrocarbon Contamination on Microbial Mat Functioning

Photosynthetic microbial mats are stable, self-supported communities. Due to their coastal localization, these mats are frequently exposed to hydrocarbon contamination and are able to grow on it. To decipher how this contamination disturbs the functioning of microbial mats, we compared two mats: a contaminated mat exposed to chronic petroleum contamination and a reference mat. The taxonomic and metabolic structures of the mats in spring and fall were determined using metagenomic and metatranscriptomic approaches. Extremely high contamination disturbed the seasonal variations of the mat. ABC transporters, two-component systems, and type IV secretion system-related genes were overabundant in the contaminated mats. Xenobiotic degradation metabolism was minor in the metagenomes of both mats, and only the expression of genes involved in polycyclic aromatic hydrocarbon degradation was higher in the contaminated mat. Interestingly, the expression rates of genes involved in hydrocarbon activation decreased during the 1-year study period, concomitant with the decrease in easily degradable hydrocarbons, suggesting a transient effect of hydrocarbon contamination. Alteromonadales and Oceanospirillales hydrocarbonoclastic bacteria appeared to be key in hydrocarbon remediation in the contaminated mat. Overall, the contaminated microbial mat was able to cope with hydrocarbon contamination and displayed an adaptive functioning that modified seasonal behaviour.

Regional sequence expansion or collapse in heterozygous genome assemblies

High levels of heterozygosity present a unique genome assembly challenge and can adversely impact downstream analyses, yet is common in sequencing datasets obtained from non-model organisms. Here we show that by re-assembling a heterozygous dataset with variant parameters and different assembly algorithms, we are able to generate assemblies whose protein annotations are statistically enriched for specific gene ontology categories. While total assembly length was not significantly affected by assembly methodologies tested, the assemblies generated varied widely in fragmentation level and we show local assembly collapse or expansion underlying the enrichment or depletion of specific protein functional groups. We show that these statistically significant deviations in gene ontology groups can occur in seemingly high-quality assemblies, and result from difficult-to-detect local sequence expansion or contractions. Given the unpredictable interplay between assembly algorithm, parameter, and biological sequence data heterozygosity, we highlight the need for better measures of assembly quality than N50 value, including methods for assessing local expansion and collapse.

In the genomic era, genomes must be reconstructed from fragments using computational methods, or assemblers. How do we know that a new genome assembly is correct? This is important because errors in assembly can lead to downstream problems in gene predictions and these inaccurate results can contaminate databases, affecting later comparative studies. A particular challenge occurs when a diploid organism inherits two highly divergent genome copies from its parents. While it is widely appreciated that this type of data is difficult for assemblers to handle properly, here we show that the process is prone to more errors than previously appreciated. Specifically, we document examples of regional expansion and collapse, affecting downstream gene prediction accuracy, but without changing the overall genome assembly size or other metrics of accuracy. Our results suggest that assembly evaluation methods should be altered to identify whether regional expansions and collapses are present in the genome assembly.

Hayai-Annotation Plants: an ultra-fast and comprehensive functional gene annotation system in plants

Summary

Hayai-Annotation Plants is a browser-based interface for an ultra-fast and accurate functional gene annotation system for plant species using R. The pipeline combines the sequence-similarity searches, using USEARCH against UniProtKB (taxonomy Embryophyta), with a functional annotation step. Hayai-Annotation Plants provides five layers of annotation: i) protein name; ii) gene ontology terms consisting of its three main domains (Biological Process, Molecular Function and Cellular Component); iii) enzyme commission number; iv) protein existence level; and v) evidence type. It implements a new algorithm that gives priority to protein existence level to propagate GO and EC information and annotated Arabidopsis thaliana representative peptide sequences (Araport11) within 5 min at the PC level.

Availability and implementation

The software is implemented in R and runs on Macintosh and Linux systems. It is freely available at https://github.com/kdri-genomics/Hayai-Annotation-Plants under the GPLv3 license.

Supplementary information

Supplementary data are available at Bioinformatics online.

Accurate Detection of HPV Integration Sites in Cervical Cancer Samples Using the Nanopore MinION Sequencer Without Error Correction

During the carcinogenesis of cervical cancer, the DNA of human papillomavirus (HPV) is frequently integrated into the human genome, which might be a biomarker for the early diagnosis of cervical cancer. Although the detection sensitivity of virus infection status increased significantly through the Illumina sequencing platform, there were still disadvantages remain for further improvement, including the detection accuracy and the complex integrated genome structure identification, etc. Nanopore sequencing has been proven to be a fast yet accurate technique of detecting pathogens in clinical samples with significant longer sequencing length. However, the identification of virus integration sites, especially HPV integration sites was seldom carried out by using nanopore platform. In this study, we evaluated the feasibility of identifying HPV integration sites by nanopore sequencer. Specifically, we re-sequenced the integration sites of a previously published sample by both nanopore and Illumina sequencing. After analyzing the results, three points of conclusions were drawn: first, 13 out of 19 previously published integration sites were found from all three datasets (i.e., nanopore, Illumina, and the published data), indicating a high overlap rate and comparability among the three platforms; second, our pipeline of nanopore and Illumina data identified 66 unique integration sites compared with previous published paper with 13 of them being verified by Sanger sequencing, indicating the higher integration sites detection sensitivity of our results compared with published data; third, we established a pipeline which could be used in HPV integration site detection by nanopore sequencing data without doing error correction analysis. In summary, a new nanopore data analysis method was tested and proved to be reliable in integration sites detection compared with methods of existing Illumina data analysis pipeline with less sequencing data required. It provides a solid evidence and tool to support the potential application of nanopore in virus status identification.

Simultaneous detection and comprehensive analysis of HPV and microbiome status of a cervical liquid-based cytology sample using Nanopore MinION sequencing

Human papillomavirus (HPV) is a major pathogen that causes cervical cancer and many other related diseases. HPV infection related cervical microbiome could be an induce factor of cervical cancer. However, it is uncommon to find a single test on the market that can simultaneously provide information on both HPV and the microbiome. Herein, a novel method was developed in this study to simultaneously detect HPV infection and microbiota composition promptly and accurately. It provides a new and simple way to detect vaginal pathogen situation and also provide valuable information for clinical diagnose. This approach combined multiplex PCR, which targeted both HPV16 E6E7 and full-length 16S rRNA, and Nanopore sequencing to generate enough information to understand the vagina condition of patients. One HPV positive liquid-based cytology (LBC) sample was sequenced and analyzed. After comparing with Illumina sequencing, the results from Nanopore showed a similar microbiome composition. An instant sequencing evaluation showed that 15 min sequencing is enough to identify the top 10 most abundant bacteria. Moreover, two HPV integration sites were identified and verified by Sanger sequencing. This approach has many potential applications in pathogen detection and can potentially aid in providing a more rapid clinical diagnosis.

The Impact of cDNA Normalization on Long-Read Sequencing of a Complex Transcriptome

Normalization of cDNA is widely used to improve the coverage of rare transcripts in analysis of transcriptomes employing next-generation sequencing. Recently, long-read technology has been emerging as a powerful tool for sequencing and construction of transcriptomes, especially for complex genomes containing highly similar transcripts and transcript-spliced isoforms. Here, we analyzed the transcriptome of sugarcane, a highly polyploidy plant genome, by PacBio isoform sequencing (Iso-Seq) of two different cDNA library preparations, with and without a normalization step. The results demonstrated that, while the two libraries included many of the same transcripts, many longer transcripts were removed, and many new generally shorter transcripts were detected by normalization. For the same input cDNA and data yield, the normalized library recovered more total transcript isoforms and number of predicted gene families and orthologous groups, resulting in a higher representation for the sugarcane transcriptome, compared to the non-normalized library. The non-normalized library, on the other hand, included a wider transcript length range with more longer transcripts above ∼1.25 kb and more transcript isoforms per gene family and gene ontology terms per transcript. A large proportion of the unique transcripts comprising ∼52% of the normalized library were expressed at a lower level than the unique transcripts from the non-normalized library, across three tissue types tested including leaf, stalk, and root. About 83% of the total 5,348 predicted long noncoding transcripts was derived from the normalized library, of which ∼80% was derived from the lowly expressed fraction. Functional annotation of the unique transcripts suggested that each library enriched different functional transcript fractions. This demonstrated the complementation of the two approaches in obtaining a complete transcriptome of a complex genome at the sequencing depth used in this study.

Tracing the Evolution of the Plastome and Mitogenome in the Chloropicophyceae Uncovered Convergent tRNA Gene Losses and a Variant Plastid Genetic Code

The tiny green algae belonging to the Chloropicophyceae play a key role in marine phytoplankton communities; this newly erected class of prasinophytes comprises two genera (Chloropicon and Chloroparvula) containing each several species. We sequenced the plastomes and mitogenomes of eight Chloropicon and five Chloroparvula species to better delineate the phylogenetic affinities of these taxa and to infer the suite of changes that their organelle genomes sustained during evolution. The relationships resolved in organelle-based phylogenomic trees were essentially congruent with previously reported rRNA trees, and similar evolutionary trends but distinct dynamics were identified for the plastome and mitogenome. Although the plastome sustained considerable changes in gene content and order at the time the two genera split, subsequently it remained stable and maintained a very small size. The mitogenome, however, was remodeled more gradually and showed more fluctuation in size, mainly as a result of expansions/contractions of intergenic regions. Remarkably, the plastome and mitogenome lost a common set of three tRNA genes, with the trnI(cau) and trnL(uaa) losses being accompanied with important variations in codon usage. Unexpectedly, despite the disappearance of trnI(cau) from the plastome in the Chloroparvula lineage, AUA codons (the codons recognized by this gene product) were detected in certain plastid genes. By comparing the sequences of plastid protein-coding genes from chloropicophycean and phylogenetically diverse chlorophyte algae with those of the corresponding predicted proteins, we discovered that the AUA codon was reassigned from isoleucine to methionine in Chloroparvula. This noncanonical genetic code has not previously been uncovered in plastids.

Toward on-site food authentication using nanopore sequencing

•MinION DNA metabarcoding is a promising tool for species identification in food.•MinION and Illumina MiSeq sequencing platforms perform equally accurate.•Species identification with MinION sequencing requires dedicated bioinformatics.

Whole-Genome Alignment

Whole-genome alignment (WGA) is the prediction of evolutionary relationships at the nucleotide level between two or more genomes. It combines aspects of both colinear sequence alignment and gene orthology prediction and is typically more challenging to address than either of these tasks due to the size and complexity of whole genomes. Despite the difficulty of this problem, numerous methods have been developed for its solution because WGAs are valuable for genome-wide analyses such as phylogenetic inference, genome annotation, and function prediction. In this chapter, we discuss the meaning and significance of WGA and present an overview of the methods that address it. We also examine the problem of evaluating whole-genome aligners and offer a set of methodological challenges that need to be tackled in order to make most effective use of our rapidly growing databases of whole genomes.

Identification and characterization of novel conserved RNA structures in Drosophila

BACKGROUND

Comparative genomics approaches have facilitated the discovery of many novel non-coding and structured RNAs (ncRNAs). The increasing availability of related genomes now makes it possible to systematically search for compensatory base changes - and thus for conserved secondary structures - even in genomic regions that are poorly alignable in the primary sequence. The wealth of available transcriptome data can add valuable insight into expression and possible function for new ncRNA candidates. Earlier work identifying ncRNAs in Drosophila melanogaster made use of sequence-based alignments and employed a sliding window approach, inevitably biasing identification toward RNAs encoded in the more conserved parts of the genome.

RESULTS

To search for conserved RNA structures (CRSs) that may not be highly conserved in sequence and to assess the expression of CRSs, we conducted a genome-wide structural alignment screen of 27 insect genomes including D. melanogaster and integrated this with an extensive set of tiling array data. The structural alignment screen revealed ∼30,000 novel candidate CRSs at an estimated false discovery rate of less than 10%. With more than one quarter of all individual CRS motifs showing sequence identities below 60%, the predicted CRSs largely complement the findings of sliding window approaches applied previously. While a sixth of the CRSs were ubiquitously expressed, we found that most were expressed in specific developmental stages or cell lines. Notably, most statistically significant enrichment of CRSs were observed in pupae, mainly in exons of untranslated regions, promotors, enhancers, and long ncRNAs. Interestingly, cell lines were found to express a different set of CRSs than were found in vivo. Only a small fraction of intergenic CRSs were co-expressed with the adjacent protein coding genes, which suggests that most intergenic CRSs are independent genetic units.

CONCLUSIONS

This study provides a more comprehensive view of the ncRNA transcriptome in fly as well as evidence for differential expression of CRSs during development and in cell lines.

A non-coding region near Follistatin controls head colour polymorphism in the Gouldian finch

Discrete colour morphs coexisting within a single population are common in nature. In a broad range of organisms, sympatric colour morphs often display major differences in other traits, including morphology, physiology or behaviour. Despite the repeated occurrence of this phenomenon, our understanding of the genetics that underlie multi-trait differences and the factors that promote the long-term maintenance of phenotypic variability within a freely interbreeding population are incomplete. Here, we investigated the genetic basis of red and black head colour in the Gouldian finch (Erythrura gouldiae), a classic polymorphic system in which naturally occurring colour morphs also display differences in aggressivity and reproductive success. We show that the candidate locus is a small (approx. 70 kb) non-coding region mapping to the Z chromosome near the Follistatin (FST) gene. Unlike recent findings in other systems where phenotypic morphs are explained by large inversions containing hundreds of genes (so-called supergenes), we did not identify any structural rearrangements between the two haplotypes using linked-read sequencing technology. Nucleotide divergence between the red and black alleles was high when compared to the remainder of the Z chromosome, consistent with their maintenance as balanced polymorphisms over several million years. Our results illustrate how pleiotropic phenotypes can arise from simple genetic variation, probably regulatory in nature.

ULTRA: A Model Based Tool to Detect Tandem Repeats

In biological sequences, tandem repeats consist of tens to hundreds of residues of a repeated pattern, such as atgatgatgatgatg ('atg' repeated), often the result of replication slippage. Over time, these repeats decay so that the original sharp pattern of repetition is somewhat obscured, but even degenerate repeats pose a problem for sequence annotation: when two sequences both contain shared patterns of similar repetition, the result can be a false signal of sequence homology. We describe an implementation of a new hidden Markov model for detecting tandem repeats that shows substantially improved sensitivity to labeling decayed repetitive regions, presents low and reliable false annotation rates across a wide range of sequence composition, and produces scores that follow a stable distribution. On typical genomic sequence, the time and memory requirements of the resulting tool (ULTRA) are competitive with the most heavily used tool for repeat masking (TRF). ULTRA is released under an open source license and lays the groundwork for inclusion of the model in sequence alignment tools and annotation pipelines.

Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviours

Strains of red fox (Vulpes vulpes) with markedly different behavioural phenotypes have been developed in the famous long-term selective breeding programme known as the Russian farm-fox experiment. Here we sequenced and assembled the red fox genome and re-sequenced a subset of foxes from the tame, aggressive and conventional farm-bred populations to identify genomic regions associated with the response to selection for behaviour. Analysis of the re-sequenced genomes identified 103 regions with either significantly decreased heterozygosity in one of the three populations or increased divergence between the populations. A strong positional candidate gene for tame behaviour was highlighted: SorCS1, which encodes the main trafficking protein for AMPA glutamate receptors and neurexins and suggests a role for synaptic plasticity in fox domestication. Other regions identified as likely to have been under selection in foxes include genes implicated in human neurological disorders, mouse behaviour and dog domestication. The fox represents a powerful model for the genetic analysis of affiliative and aggressive behaviours that can benefit genetic studies of behaviour in dogs and other mammals, including humans.

Nanopore sequencing and full genome de novo assembly of human cytomegalovirus TB40/E reveals clonal diversity and structural variations

BACKGROUND

Human cytomegalovirus (HCMV) has a double-stranded DNA genome of approximately 235 Kbp that is structurally complex including extended GC-rich repeated regions. Genomic recombination events are frequent in HCMV cultures but have also been observed in vivo. Thus, the assembly of HCMV whole genomes from technologies producing shorter than 500 bp sequences is technically challenging. Here we improved the reconstruction of HCMV full genomes by means of a hybrid, de novo genome-assembly bioinformatics pipeline upon data generated from the recently released MinION MkI B sequencer from Oxford Nanopore Technologies.

RESULTS

The MinION run of the HCMV (strain TB40/E) library resulted in ~ 47,000 reads from a single R9 flowcell and in ~ 100× average read depth across the virus genome. We developed a novel, self-correcting bioinformatics algorithm to assemble the pooled HCMV genomes in three stages. In the first stage of the bioinformatics algorithm, long contigs (N50 = 21,892) of lower accuracy were reconstructed. In the second stage, short contigs (N50 = 5686) of higher accuracy were assembled, while in the final stage the high quality contigs served as template for the correction of the longer contigs resulting in a high-accuracy, full genome assembly (N50 = 41,056). We were able to reconstruct a single representative haplotype without employing any scaffolding steps. The majority (98.8%) of the genomic features from the reference strain were accurately annotated on this full genome construct. Our method also allowed the detection of multiple alternative sub-genomic fragments and non-canonical structures suggesting rearrangement events between the unique (UL /US) and the repeated (T/IRL/S) genomic regions.

CONCLUSIONS

Third generation high-throughput sequencing technologies can accurately reconstruct full-length HCMV genomes including their low-complexity and highly repetitive regions. Full-length HCMV genomes could prove crucial in understanding the genetic determinants and viral evolution underpinning drug resistance, virulence and pathogenesis.

Multiplexed Targeted Sequencing for Oxford Nanopore MinION: A Detailed Library Preparation Procedure

MinION is a small form factor sequencer recently retailed by Oxford Nanopore technologies. This lighter-sized USB3.0-interfaced device uses innovative nanotechnology to generate extra-long reads from libraries prepared using only standard molecular biology lab equipment. The flexibility and the portability of the platform makes it ideal for point-of-interest and real-time surveillance applications. However, MinION's limited capacity is not enough for the study of specific targets within larger genomes. Apart from just PCR-amplifying regions of interest, the capture of long reads spanning the edges of known-unknown genomic regions is of great importance for structural studies, such as the identification of mobile elements' integrations sites, bridging over low complexity repetitive regions etc.In this study, using MinION-kit-included and commercially available reagents, we have developed an easy and versatile wet-lab procedure for the targeted enrichment of MinION libraries, capturing DNA fragments of interest before the ligation of the sensitive MinION sequencing-adapters. This method allows for simultaneous target-enrichment and barcode-multiplexing of up to 12 libraries, which can be loaded in the same sequencing run.

Picky comprehensively detects high-resolution structural variants in nanopore long reads

Acquired genomic structural variants (SVs) are major hallmarks of cancer genomes, but they are challenging to reconstruct from short-read sequencing data. Here we exploited the long reads of the nanopore platform using our customized pipeline, Picky ( https://github.com/TheJacksonLaboratory/Picky ), to reveal SVs of diverse architecture in a breast cancer model. We identified the full spectrum of SVs with superior specificity and sensitivity relative to short-read analyses, and uncovered repetitive DNA as the major source of variation. Examination of genome-wide breakpoints at nucleotide resolution uncovered micro-insertions as the common structural features associated with SVs. Breakpoint density across the genome is associated with the propensity for interchromosomal connectivity and was found to be enriched in promoters and transcribed regions of the genome. Furthermore, we observed an over-representation of reciprocal translocations from chromosomal double-crossovers through phased SVs. We demonstrate that Picky analysis is an effective tool for comprehensive detection of SVs in cancer genomes from long-read data.