Genome-wide mapping and assembly of structural variant breakpoints in the mouse genome

Behavioral and transcriptional effects of repeated electroconvulsive seizures in the neonatal MK-801-treated rat model of schizophrenia

RATIONALE

Electroconvulsive therapy (ECT) is an effective treatment modality for schizophrenia. However, its antipsychotic-like mechanism remains unclear.

OBJECTIVES

To gain insight into the antipsychotic-like actions of ECT, this study investigated how repeated treatments of electroconvulsive seizure (ECS), an animal model for ECT, affect the behavioral and transcriptomic profile of a neurodevelopmental animal model of schizophrenia.

METHODS

Two injections of MK-801 or saline were administered to rats on postnatal day 7 (PN7), and either repeated ECS treatments (E10X) or sham shock was conducted daily from PN50 to PN59. Ultimately, the rats were divided into vehicle/sham (V/S), MK-801/sham (M/S), vehicle/ECS (V/E), and MK-801/ECS (M/E) groups. On PN59, prepulse inhibition and locomotor activity were tested. Prefrontal cortex transcriptomes were analyzed with mRNA sequencing and network and pathway analyses, and quantitative real-time polymerase chain reaction (qPCR) analyses were subsequently conducted.

RESULTS

Prepulse inhibition deficit was induced by MK-801 and normalized by E10X. In M/S vs. M/E model, Egr1, Mmp9, and S100a6 were identified as center genes, and interleukin-17 (IL-17), nuclear factor kappa B (NF-κB), and tumor necrosis factor (TNF) signaling pathways were identified as the three most relevant pathways. In the V/E vs. V/S model, mitophagy, NF-κB, and receptor for advanced glycation end products (RAGE) pathways were identified. qPCR analyses demonstrated that Igfbp6, Btf3, Cox6a2, and H2az1 were downregulated in M/S and upregulated in M/E.

CONCLUSIONS

E10X reverses the behavioral changes induced by MK-801 and produces transcriptional changes in inflammatory, insulin, and mitophagy pathways, which provide mechanistic insight into the antipsychotic-like mechanism of ECT.

Incorporating CNV analysis improves the yield of exome sequencing for rare monogenic disorders-an important consideration for resource-constrained settings

Exome sequencing (ES) is a recommended first-tier diagnostic test for many rare monogenic diseases. It allows for the detection of both single-nucleotide variants (SNVs) and copy number variants (CNVs) in coding exonic regions of the genome in a single test, and this dual analysis is a valuable approach, especially in limited resource settings. Single-nucleotide variants are well studied; however, the incorporation of copy number variant analysis tools into variant calling pipelines has not been implemented yet as a routine diagnostic test, and chromosomal microarray is still more widely used to detect copy number variants. Research shows that combined single and copy number variant analysis can lead to a diagnostic yield of up to 58%, increasing the yield with as much as 18% from the single-nucleotide variant only pipeline. Importantly, this is achieved with the consideration of computational costs only, without incurring any additional sequencing costs. This mini review provides an overview of copy number variant analysis from exome data and what the current recommendations are for this type of analysis. We also present an overview on rare monogenic disease research standard practices in resource-limited settings. We present evidence that integrating copy number variant detection tools into a standard exome sequencing analysis pipeline improves diagnostic yield and should be considered a significantly beneficial addition, with relatively low-cost implications. Routine implementation in underrepresented populations and limited resource settings will promote generation and sharing of CNV datasets and provide momentum to build core centers for this niche within genomic medicine.

Identification and characterisation of de novo germline structural variants in two commercial pig lines using trio-based whole genome sequencing

BACKGROUND

De novo mutations arising in the germline are a source of genetic variation and their discovery broadens our understanding of genetic disorders and evolutionary patterns. Although the number of de novo single nucleotide variants (dnSNVs) has been studied in a number of species, relatively little is known about the occurrence of de novo structural variants (dnSVs). In this study, we investigated 37 deeply sequenced pig trios from two commercial lines to identify dnSVs present in the offspring. The identified dnSVs were characterised by identifying their parent of origin, their functional annotations and characterizing sequence homology at the breakpoints.

RESULTS

We identified four swine germline dnSVs, all located in intronic regions of protein-coding genes. Our conservative, first estimate of the swine germline dnSV rate is 0.108 (95% CI 0.038-0.255) per generation (one dnSV per nine offspring), detected using short-read sequencing. Two detected dnSVs are clusters of mutations. Mutation cluster 1 contains a de novo duplication, a dnSNV and a de novo deletion. Mutation cluster 2 contains a de novo deletion and three de novo duplications, of which one is inverted. Mutation cluster 2 is 25 kb in size, whereas mutation cluster 1 (197 bp) and the other two individual dnSVs (64 and 573 bp) are smaller. Only mutation cluster 2 could be phased and is located on the paternal haplotype. Mutation cluster 2 originates from both micro-homology as well as non-homology mutation mechanisms, where mutation cluster 1 and the other two dnSVs are caused by mutation mechanisms lacking sequence homology. The 64 bp deletion and mutation cluster 1 were validated through PCR. Lastly, the 64 bp deletion and the 573 bp duplication were validated in sequenced offspring of probands with three generations of sequence data.

CONCLUSIONS

Our estimate of 0.108 dnSVs per generation in the swine germline is conservative, due to our small sample size and restricted possibilities of dnSV detection from short-read sequencing. The current study highlights the complexity of dnSVs and shows the potential of breeding programs for pigs and livestock species in general, to provide a suitable population structure for identification and characterisation of dnSVs.

Chromosome-level genome of the three-spot damselfish, Dascyllus trimaculatus

Damselfishes (Family: Pomacentridae) are a group of ecologically important, primarily coral reef fishes that include over 400 species. Damselfishes have been used as model organisms to study recruitment (anemonefishes), the effects of ocean acidification (spiny damselfish), population structure, and speciation (Dascyllus). The genus Dascyllus includes a group of small-bodied species, and a complex of relatively larger bodied species, the Dascyllus trimaculatus species complex that is comprised of several species including D. trimaculatus itself. The three-spot damselfish, D. trimaculatus, is a widespread and common coral reef fish species found across the tropical Indo-Pacific. Here, we present the first-genome assembly of this species. This assembly contains 910 Mb, 90% of the bases are in 24 chromosome-scale scaffolds, and the Benchmarking Universal Single-Copy Orthologs score of the assembly is 97.9%. Our findings confirm previous reports of a karyotype of 2n = 47 in D. trimaculatus in which one parent contributes 24 chromosomes and the other 23. We find evidence that this karyotype is the result of a heterozygous Robertsonian fusion. We also find that the D. trimaculatus chromosomes are each homologous with single chromosomes of the closely related clownfish species, Amphiprion percula. This assembly will be a valuable resource in the population genomics and conservation of Damselfishes, and continued studies of the karyotypic diversity in this clade.

In vivo interrogation of regulatory genomes reveals extensive quasi-insufficiency in cancer evolution

In contrast to mono- or biallelic loss of tumor-suppressor function, effects of discrete gene dysregulations, as caused by non-coding (epi)genome alterations, are poorly understood. Here, by perturbing the regulatory genome in mice, we uncover pervasive roles of subtle gene expression variation in cancer evolution. Genome-wide screens characterizing 1,450 tumors revealed that such quasi-insufficiency is extensive across entities and displays diverse context dependencies, such as distinct cell-of-origin associations in T-ALL subtypes. We compile catalogs of non-coding regions linked to quasi-insufficiency, show their enrichment with human cancer risk variants, and provide functional insights by engineering regulatory alterations in mice. As such, kilo-/megabase deletions in a Bcl11b-linked non-coding region triggered aggressive malignancies, with allele-specific tumor spectra reflecting gradual gene dysregulations through modular and cell-type-specific enhancer activities. Our study constitutes a first survey toward a systems-level understanding of quasi-insufficiency in cancer and gives multifaceted insights into tumor evolution and the tissue-specific effects of non-coding mutations.

The oxidative phosphorylation inhibitor IM156 suppresses B-cell activation by regulating mitochondrial membrane potential and contributes to the mitigation of systemic lupus erythematosus

Current treatment strategies for autoimmune diseases may not sufficiently control aberrant metabolism in B-cells. To address this concern, we investigated a biguanide derivative, IM156, as a potential regulator for B-cell metabolism in vitro and in vivo on overactive B-cells stimulated by the pro-inflammatory receptor TLR-9 agonist CpG oligodeoxynucleotide, a mimic of viral/bacterial DNA. Using RNA sequencing, we analyzed the B-cell transcriptome expression, identifying the major molecular pathways affected by IM156 in vivo. We also evaluated the anti-inflammatory effects of IM156 in lupus-prone NZB/W F1 mice. CD19⁺B-cells exhibited higher mitochondrial mass and mitochondrial membrane potential compared to T-cells and were more susceptible to IM156-mediated oxidative phosphorylation inhibition. In vivo, IM156 inhibited mitochondrial oxidative phosphorylation, cell cycle progression, plasmablast differentiation, and activation marker levels in CpG oligodeoxynucleotide-stimulated mouse spleen B-cells. Interestingly, IM156 treatment significantly increased overall survival, reduced glomerulonephritis and inhibited B-cell activation in the NZB/W F1 mice. Thus, our data indicated that IM156 suppressed the mitochondrial membrane potentials of activated B-cells in mice, contributing to the mitigation of lupus activity. Hence, IM156 may represent a therapeutic alternative for autoimmune disease mediated by B-cell hyperactivity.

Towards mouse genetic-specific RNA-sequencing read mapping

Genetic variations affect behavior and cause disease but understanding how these variants drive complex traits is still an open question. A common approach is to link the genetic variants to intermediate molecular phenotypes such as the transcriptome using RNA-sequencing (RNA-seq). Paradoxically, these variants between the samples are usually ignored at the beginning of RNA-seq analyses of many model organisms. This can skew the transcriptome estimates that are used later for downstream analyses, such as expression quantitative trait locus (eQTL) detection. Here, we assessed the impact of reference-based analysis on the transcriptome and eQTLs in a widely-used mouse genetic population: the BXD panel of recombinant inbred lines. We highlight existing reference bias in the transcriptome data analysis and propose practical solutions which combine available genetic variants, genotypes, and genome reference sequence. The use of custom BXD line references improved downstream analysis compared to classical genome reference. These insights would likely benefit genetic studies with a transcriptomic component and demonstrate that genome references need to be reassessed and improved.

To understand how genetic variations affect behavior and cause disease it is common to quantify expression of transcripts by sequencing. Transcripts are extracted, fragmented, and the sequence of the fragments read. An important step for their quantification is to virtually assign the different fragments to the transcript they originate from using a reference genome. Reference genomes are costly to build, so usually only one high-quality reference per animal model species is available. When comparing genetically different individuals, using a single reference may introduce a bias because it might be more similar to some individuals than to others. Paradoxically, the variations at the core of genetic studies are thus ignored at the start of the analysis. We built customized references with known genetic variants for each of the mouse lines we had and quantified the impact of the reference at different levels of the bioinformatic analysis. We found that using customized references reduced the bias compared to using a single reference. Our study uses publicly available data and tools, so others can easily implement this improvement in their analyses.

Comparative genomics of Nocardia seriolae reveals recent importation and subsequent widespread dissemination in mariculture farms in the South Central Coast region, Vietnam

Between 2010 and 2015, nocardiosis outbreaks caused by Nocardia seriolae affected many permit farms throughout Vietnam, causing mass fish mortalities. To understand the biology, origin and epidemiology of these outbreaks, 20 N. seriolae strains collected from farms in four provinces in the South Central Coast region of Vietnam, along with two Taiwanese strains, were analysed using genetics and genomics. PFGE identified a single cluster amongst all Vietnamese strains that was distinct from the Taiwanese strains. Like the PFGE findings, phylogenomic and SNP genotyping analyses revealed that all Vietnamese N. seriolae strains belonged to a single, unique clade. Strains fell into two subclades that differed by 103 SNPs, with almost no diversity within clades (0–5 SNPs). There was no association between geographical origin and subclade placement, suggesting frequent N. seriolae transmission between Vietnamese mariculture facilities during the outbreaks. The Vietnamese strains shared a common ancestor with strains from Japan and China, with the closest strain, UTF1 from Japan, differing by just 220 SNPs from the Vietnamese ancestral node. Draft Vietnamese genomes range from 7.55 to 7.96 Mbp in size, have an average G+C content of 68.2 % and encode 7 602–7958 predicted genes. Several putative virulence factors were identified, including genes associated with host cell adhesion, invasion, intracellular survival, antibiotic and toxic compound resistance, and haemolysin biosynthesis. Our findings provide important new insights into the epidemiology and pathogenicity of N. seriolae and will aid future vaccine development and disease management strategies, with the ultimate goal of nocardiosis-free aquaculture.

Pierre CL, Williams A, Sohni A, Wilkinson MF, Gymrek M, Sebat J, Palmer AA. SNPs, short tandem repeats, and structural variants are responsible for differential gene expression across C57BL/6 and C57BL/10 substrains

Mouse substrains are an invaluable model for understanding disease. We compared C57BL/6J, which is the most commonly used inbred mouse strain, with eight C57BL/6 and five C57BL/10 closely related inbred substrains. Whole-genome sequencing and RNA-sequencing analysis yielded 352,631 SNPs, 109,096 indels, 150,344 short tandem repeats (STRs), 3,425 structural variants (SVs), and 2,826 differentially expressed genes (DE genes) among these 14 strains; 312,981 SNPs (89%) distinguished the B6 and B10 lineages. These SNPs were clustered into 28 short segments that are likely due to introgressed haplotypes rather than new mutations. Outside of these introgressed regions, we identified 53 SVs, protein-truncating SNPs, and frameshifting indels that were associated with DE genes. Our results can be used for both forward and reverse genetic approaches and illustrate how introgression and mutational processes give rise to differences among these widely used inbred substrains.

Progress in Methods for Copy Number Variation Profiling

Copy number variations (CNVs) are the predominant class of structural genomic variations involved in the processes of evolutionary adaptation, genomic disorders, and disease progression. Compared with single-nucleotide variants, there have been challenges associated with the detection of CNVs owing to their diverse sizes. However, the field has seen significant progress in the past 20–30 years. This has been made possible due to the rapid development of molecular diagnostic methods which ensure a more detailed view of the genome structure, further complemented by recent advances in computational methods. Here, we review the major approaches that have been used to routinely detect CNVs, ranging from cytogenetics to the latest sequencing technologies, and then cover their specific features.

The coral symbiont Candidatus Aquarickettsia is variably abundant in threatened Caribbean acroporids and transmitted horizontally

The symbiont "Candidatus Aquarickettsia rohweri" infects a diversity of aquatic hosts. In the threatened Caribbean coral, Acropora cervicornis, Aquarickettsia proliferates in response to increased nutrient exposure, resulting in suppressed growth and increased disease susceptibility and mortality of coral. This study evaluated the extent, as well as the ecology and evolution of Aquarickettsia infecting threatened corals, Ac. cervicornis, and Ac. palmata and their hybrid ("Ac. prolifera"). Aquarickettsia was found in all acroporids, with coral host and geographic location impacting the infection magnitude. Phylogenomic and genome-wide single-nucleotide variant analysis of Aquarickettsia found phylogenetic clustering by geographic region, not by coral taxon. Analysis of Aquarickettsia fixation indices suggests multiple sequential infections of the same coral colony are unlikely. Furthermore, relative to other Rickettsiales species, Aquarickettsia is undergoing positive selection, with Florida populations experiencing greater positive selection relative to other Caribbean locations. This may be due in part to Aquarickettsia proliferating in response to greater nutrient stress in Florida, as indicated by greater in situ replication rates in these corals. Aquarickettsia was not found to significantly codiversify with either the coral animal or the coral's algal symbiont (Symbiodinium "fitti"). Quantitative PCR analysis showed that gametes, larvae, recruits, and juveniles from susceptible, captive-reared coral genets were not infected with Aquarickettsia. Thus, horizontal transmission of Aquarickettsia via coral mucocytes or an unidentified host is more likely. The prevalence of Aquarickettsia in Ac. cervicornis and its high abundance in the Florida coral population suggests that coral disease mitigation efforts focus on preventing early infection via horizontal transmission.

Diverse tumorigenic consequences of human papillomavirus integration in primary oropharyngeal cancers

Human papillomavirus (HPV) causes 5% of all cancers and frequently integrates into host chromosomes. The HPV oncoproteins E6 and E7 are necessary but insufficient for cancer formation, indicating that additional secondary genetic events are required. Here, we investigate potential oncogenic impacts of virus integration. Analysis of 105 HPV-positive oropharyngeal cancers by whole-genome sequencing detects virus integration in 77%, revealing five statistically significant sites of recurrent integration near genes that regulate epithelial stem cell maintenance (i.e., SOX2, TP63, FGFR, MYC) and immune evasion (i.e., CD274). Genomic copy number hyperamplification is enriched 16-fold near HPV integrants, and the extent of focal host genomic instability increases with their local density. The frequency of genes expressed at extreme outlier levels is increased 86-fold within ±150 kb of integrants. Across 95% of tumors with integration, host gene transcription is disrupted via intragenic integrants, chimeric transcription, outlier expression, gene breaking, and/or de novo expression of noncoding or imprinted genes. We conclude that virus integration can contribute to carcinogenesis in a large majority of HPV-positive oropharyngeal cancers by inducing extensive disruption of host genome structure and gene expression.

Massive QTL analysis identifies pleiotropic genetic determinants for stress resistance, aroma formation, and ethanol, glycerol and isobutanol production in Saccharomyces cerevisiae

BACKGROUND

The brewer's yeast Saccharomyces cerevisiae is exploited in several industrial processes, ranging from food and beverage fermentation to the production of biofuels, pharmaceuticals and complex chemicals. The large genetic and phenotypic diversity within this species offers a formidable natural resource to obtain superior strains, hybrids, and variants. However, most industrially relevant traits in S. cerevisiae strains are controlled by multiple genetic loci. Over the past years, several studies have identified some of these QTLs. However, because these studies only focus on a limited set of traits and often use different techniques and starting strains, a global view of industrially relevant QTLs is still missing.

RESULTS

Here, we combined the power of 1125 fully sequenced inbred segregants with high-throughput phenotyping methods to identify as many as 678 QTLs across 18 different traits relevant to industrial fermentation processes, including production of ethanol, glycerol, isobutanol, acetic acid, sulfur dioxide, flavor-active esters, as well as resistance to ethanol, acetic acid, sulfite and high osmolarity. We identified and confirmed several variants that are associated with multiple different traits, indicating that many QTLs are pleiotropic. Moreover, we show that both rare and common variants, as well as variants located in coding and non-coding regions all contribute to the phenotypic variation.

CONCLUSIONS

Our findings represent an important step in our understanding of the genetic underpinnings of industrially relevant yeast traits and open new routes to study complex genetics and genetic interactions as well as to engineer novel, superior industrial yeasts. Moreover, the major role of rare variants suggests that there is a plethora of different combinations of mutations that can be explored in genome editing.

Gv1, a Zinc Finger Gene Controlling Endogenous MLV Expression

The genomes of inbred mice harbor around 50 endogenous murine leukemia virus (MLV) loci, although the specific complement varies greatly between strains. The Gv1 locus is known to control the transcription of endogenous MLVs and to be the dominant determinant of cell-surface presentation of MLV envelope, the G_IX antigen. Here, we identify a single Krüppel-associated box zinc finger protein (ZFP) gene, Zfp998, as Gv1 and show it to be necessary and sufficient to determine the GIX+ phenotype. By long-read sequencing of bacterial artificial chromosome clones from 129 mice, the prototypic GIX+ strain, we reveal the source of sufficiency and deficiency as splice-acceptor variations and highlight the varying origins of the chromosomal region encompassing Gv1. Zfp998 becomes the second identified ZFP gene responsible for epigenetic suppression of endogenous MLVs in mice and further highlights the prominent role of this gene family in control of endogenous retroviruses.

De novo structural mutation rates and gamete-of-origin biases revealed through genome sequencing of 2,396 families

Each human genome includes de novo mutations that arose during gametogenesis. While these germline mutations represent a fundamental source of new genetic diversity, they can also create deleterious alleles that impact fitness. Whereas the rate and patterns of point mutations in the human germline are now well understood, far less is known about the frequency and features that impact de novo structural variants (dnSVs). We report a family-based study of germline mutations among 9,599 human genomes from 33 multigenerational CEPH-Utah families and 2,384 families from the Simons Foundation Autism Research Initiative. We find that de novo structural mutations detected by alignment-based, short-read WGS occur at an overall rate of at least 0.160 events per genome in unaffected individuals, and we observe a significantly higher rate (0.206 per genome) in ASD-affected individuals. In both probands and unaffected samples, nearly 73% of de novo structural mutations arose in paternal gametes, and we predict most de novo structural mutations to be caused by mutational mechanisms that do not require sequence homology. After multiple testing correction, we did not observe a statistically significant correlation between parental age and the rate of de novo structural variation in offspring. These results highlight that a spectrum of mutational mechanisms contribute to germline structural mutations and that these mechanisms most likely have markedly different rates and selective pressures than those leading to point mutations.

Novel Genetic Rearrangements in Hepatitis B Virus: Complex Structural Variations and Structural Variation Polymorphisms

Chronic hepatitis B virus (HBV) causes serious clinical problems, such as liver cirrhosis and hepatocellular carcinoma. Current antiviral treatments suppress HBV; however, the clinical cure rate remains low. Basic research on HBV is indispensable to eradicate and cure HBV. Genetic alterations are defined by nucleotide substitutions and canonical forms of structural variations (SVs), such as insertion, deletion and duplication. Additionally, genetic changes inconsistent with the canonical forms have been reported, and these have been termed complex SVs. Detailed analyses of HBV using bioinformatical applications have detected complex SVs in HBV genomes. Sequence gaps and low sequence similarity have been observed in the region containing complex SVs. Additionally, insertional motif sequences have been observed in HBV strains with complex SVs. Following the analyses of complex SVs in the HBV genome, the role of SVs in the genetic diversity of orthohepadnavirus has been investigated. SV polymorphisms have been detected in comparisons of several species of orthohepadnaviruses. As mentioned, complex SVs are composed of multiple SVs. On the contrary, SV polymorphisms are observed as insertions of different SVs. Up to a certain point, nucleotide substitutions cause genetic differences. However, at some point, the nucleotide sequences are split into several particular patterns. These SVs have been observed as polymorphic changes. Different species of orthohepadnaviruses possess SVs which are unique and specific to a certain host of the virus. Studies have shown that SVs play an important role in the HBV genome. Further studies are required to elucidate their virologic and clinical roles.

Detecting Causal Variants in Mendelian Disorders Using Whole-Genome Sequencing

Increasingly affordable sequencing technologies are revolutionizing the field of genomic medicine. It is now feasible to interrogate all major classes of variation in an individual across the entire genome for less than $1000 USD. While the generation of patient sequence information using these technologies has become routine, the analysis and interpretation of this data remains the greatest obstacle to widespread clinical implementation. This chapter summarizes the steps to identify, annotate, and prioritize variant information required for clinical report generation. We discuss methods to detect each variant class and describe strategies to increase the likelihood of detecting causal variant(s) in Mendelian disease. Lastly, we describe a sample workflow for synthesizing large amount of genetic information into concise clinical reports.

Mating system evolution and genetic structure of diploid sexual populations of Cyrtomium falcatum in Japan

Evolution of mating systems has become one of the most important research areas in evolutionary biology. Cyrtomium falcatum is a homosporous fern species native to eastern Asia. Two subspecies belonging to a sexual diploid race of C. falcatum are recognized: subsp. littorale and subsp. australe. Subspecies littorale shows intermediate selfing rates, while subsp. australe is an obligate outcrosser. We aimed to evaluate the process of mating system evolution and divergence for the two subspecies using restriction site associated DNA sequencing (RAD-seq). The results showed that subsp. littorale had lower genetic diversity and stronger genetic drift than subsp. australe. Fluctuations in the effective population size over time were evaluated by extended Bayesian skyline plot and Stairway plot analyses, both of which revealed a severe population bottleneck about 20,000 years ago in subsp. littorale. This bottleneck and the subsequent range expansion after the LGM appear to have played an important role in the divergence of the two subspecies and the evolution of selfing in subsp. littorale. These results shed new light on the relationship between mating system evolution and past demographic change in fern species.

Sequencing facility and DNA source associated patterns of virus-mappable reads in whole-genome sequencing data

Numerous viral sequences have been reported in the whole-genome sequencing (WGS) data of human blood. However, it is not clear to what degree the virus-mappable reads represent true viral sequences rather than random-mapping or noise originating from sample preparation, sequencing processes, or other sources. Identification of patterns of virus-mappable reads may generate novel indicators for evaluating the origins of these viral sequences. We characterized paired-end unmapped reads and reads aligned to viral references in human WGS datasets, then compared patterns of the virus-mappable reads among DNA sources and sequencing facilities which produced these datasets. We then examined potential origins of the source- and facility-associated viral reads. The proportions of clean unmapped reads among the seven sequencing facilities were significantly different (P < 2 × 10-16). We identified 260,339 reads that were mappable to a total of 99 viral references in 2535 samples. The majority (86.7%) of these virus-mappable reads (corresponding to 47 viral references), which can be classified into four groups based on their distinct patterns, were strongly associated with sequencing facility or DNA source (adjusted P value <0.01). Possible origins of these reads include artificial sequences in library preparation, recombinant vectors in cell culture, and phages co-contaminated with their host bacteria. The sequencing facility-associated virus-mappable reads and patterns were repeatedly observed in other datasets produced in the same facilities. We have constructed an analytic framework and profiled the unmapped reads mappable to viral references. The results provide a new understanding of sequencing facility- and DNA source-associated batch effects in deep sequencing data and may facilitate improved bioinformatics filtering of reads.

Performance of copy number variants detection based on whole-genome sequencing by DNBSEQ platforms

BACKGROUND

DNBSEQ™ platforms are new massively parallel sequencing (MPS) platforms that use DNA nanoball technology. Use of data generated from DNBSEQ™ platforms to detect single nucleotide variants (SNVs) and small insertions and deletions (indels) has proven to be quite effective, while the feasibility of copy number variants (CNVs) detection is unclear.

RESULTS

Here, we first benchmarked different CNV detection tools based on Illumina whole-genome sequencing (WGS) data of NA12878 and then assessed these tools in CNV detection based on DNBSEQ™ sequencing data from the same sample. When the same tool was used, the CNVs detected based on DNBSEQ™ and Illumina data were similar in quantity, length and distribution, while great differences existed within results from different tools and even based on data from a single platform. We further estimated the CNV detection power based on available CNV benchmarks of NA12878 and found similar precision and sensitivity between the DNBSEQ™ and Illumina platforms. We also found higher precision of CNVs shorter than 1 kbp based on DNBSEQ™ platforms than those based on Illumina platforms by using Pindel, DELLY and LUMPY. We carefully compared these two available benchmarks and found a large proportion of specific CNVs between them. Thus, we constructed a more complete CNV benchmark of NA12878 containing 3512 CNV regions.

CONCLUSIONS

We assessed and benchmarked CNV detections based on WGS with DNBSEQ™ platforms and provide guidelines for future studies.

Reconstruction of clone- and haplotype-specific cancer genome karyotypes from bulk tumor samples

Many cancer genomes are extensively rearranged with aberrant chromosomal karyotypes. Deriving these karyotypes from high-throughput DNA sequencing of bulk tumor samples is complicated because most tumors are a heterogeneous mixture of normal cells and subpopulations of cancer cells, or clones, that harbor distinct somatic mutations. We introduce a new algorithm, Reconstructing Cancer Karyotypes (RCK), to reconstruct haplotype-specific karyotypes of one or more rearranged cancer genomes from DNA sequencing data from a bulk tumor sample. RCK leverages evolutionary constraints on the somatic mutational process in cancer to reduce ambiguity in the deconvolution of admixed sequencing data into multiple haplotype-specific cancer karyotypes. RCK models mixtures containing an arbitrary number of derived genomes and allows the incorporation of information both from short-read and long-read DNA sequencing technologies. We compare RCK to existing approaches on 17 primary and metastatic prostate cancer samples. We find that RCK infers cancer karyotypes that better explain the DNA sequencing data and conform to a reasonable evolutionary model. RCK's reconstructions of clone- and haplotype-specific karyotypes will aid further studies of the role of intra-tumor heterogeneity in cancer development and response to treatment. RCK is freely available as open source software.

Mapping and characterization of structural variation in 17,795 human genomes

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

ERVcaller: identifying polymorphic endogenous retrovirus and other transposable element insertions using whole-genome sequencing data

MOTIVATION

Approximately 8% of the human genome is derived from endogenous retroviruses (ERVs). In recent years, an increasing number of human diseases have been found to be associated with ERVs. However, it remains challenging to accurately detect the full spectrum of polymorphic (unfixed) ERVs using whole-genome sequencing (WGS) data.

RESULTS

We designed a new tool, ERVcaller, to detect and genotype transposable element (TE) insertions, including ERVs, in the human genome. We evaluated ERVcaller using both simulated and real benchmark WGS datasets. Compared to existing tools, ERVcaller consistently obtained both the highest sensitivity and precision for detecting simulated ERV and other TE insertions derived from real polymorphic TE sequences. For the WGS data from the 1000 Genomes Project, ERVcaller detected the largest number of TE insertions per sample based on consensus TE loci. By analyzing the experimentally verified TE insertions, ERVcaller had 94.0% TE detection sensitivity and 96.6% genotyping accuracy. Polymerase chain reaction and Sanger sequencing in a small sample set verified 86.7% of examined insertion statuses and 100% of examined genotypes. In conclusion, ERVcaller is capable of detecting and genotyping TE insertions using WGS data with both high sensitivity and precision. This tool can be applied broadly to other species.

AVAILABILITY AND IMPLEMENTATION

http://www.uvm.edu/genomics/software/ERVcaller.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONY: A Bayesian procedure for detecting copy number variations from sequencing read depths

Copy number variations (CNVs) are genomic structural mutations consisting of abnormal numbers of fragment copies. Next-generation sequencing of read-depth signals mirrors these variants. Some tools used to predict CNVs by depth have been published, but most of these tools can be applied to only a specific data type due to modeling limitations. We develop a tool for copy number variation detection by a Bayesian procedure, i.e., CONY, that adopts a Bayesian hierarchical model and an efficient reversible-jump Markov chain Monte Carlo inference algorithm for whole genome sequencing of read-depth data. CONY can be applied not only to individual samples for estimating the absolute number of copies but also to case-control pairs for detecting patient-specific variations. We evaluate the performance of CONY and compare CONY with competing approaches through simulations and by using experimental data from the 1000 Genomes Project. CONY outperforms the other methods in terms of accuracy in both single-sample and paired-samples analyses. In addition, CONY performs well regardless of whether the data coverage is high or low. CONY is useful for detecting both absolute and relative CNVs from read-depth data sequences. The package is available at https://github.com/weiyuchung/CONY.

Mapping and characterization of structural variation in 17,795 human genomes

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline¹ to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

The abiding relevance of mouse models of rare mutations to psychiatric neuroscience and therapeutics

Studies using powerful family-based designs aided by large scale case-control studies, have been instrumental in cracking the genetic complexity of the disease, identifying rare and highly penetrant risk mutations and providing a handle on experimentally tractable model systems. Mouse models of rare mutations, paired with analysis of homologous cognitive and sensory processing deficits and state-of-the-art neuroscience methods to manipulate and record neuronal activity have started providing unprecedented insights into pathogenic mechanisms and building the foundation of a new biological framework for understanding mental illness. A number of important principles are emerging, namely that degradation of the computational mechanisms underlying the ordered activity and plasticity of both local and long-range neuronal assemblies, the building blocks necessary for stable cognition and perception, might be the inevitable consequence and the common point of convergence of the vastly heterogeneous genetic liability, manifesting as defective internally- or stimulus-driven neuronal activation patterns and triggering the constellation of schizophrenia symptoms. Animal models of rare mutations have the unique potential to help us move from "which" (gene) to "how", "where" and "when" computational regimes of neural ensembles are affected. Linking these variables should improve our understanding of how symptoms emerge and how diagnostic boundaries are established at a circuit level. Eventually, a better understanding of pathophysiological trajectories at the level of neural circuitry in mice, aided by basic human experimental biology, should guide the development of new therapeutics targeting either altered circuitry itself or the underlying biological pathways.

Detecting inherited and novel structural variants in low-coverage parent-child sequencing data

Structural variants (SVs) are a class of genomic variation shared by members of the same species. Though relatively rare, they represent an increasingly important class of variation, as SVs have been associated with diseases and susceptibility to some types of cancer. Common approaches to SV detection require the sequencing and mapping of fragments from a test genome to a high-quality reference genome. Candidate SVs correspond to fragments with discordant mapped configurations. However, because errors in the sequencing and mapping will also create discordant arrangements, many of these predictions will be spurious. When sequencing coverage is low, distinguishing true SVs from errors is even more challenging. In recent work, we have developed SV detection methods that exploit genome information of closely related individuals - parents and children. Our previous approaches were based on the assumption that any SV present in a child's genome must have come from one of their parents. However, using this strict restriction may have resulted in failing to predict any rare but novel variants present only in the child. In this work, we generalize our previous approaches to allow the child to carry novel variants. We consider a constrained optimization approach where variants in the child are of two types either inherited - and therefore must be present in a parent - or novel. For simplicity, we consider only a single parent and single child each of which have a haploid genome. However, even in this restricted case, our approach has the power to improve variant prediction. We present results on both simulated candidate variant regions, parent-child trios from the 1000 Genomes Project, and a subset of the 17 Platinum Genomes.

Ectopic expression of the Stabilin2 gene triggered by an intracisternal A particle (IAP) element in DBA/2J strain of mice

Stabilin2 (Stab2) encodes a large transmembrane protein which is predominantly expressed in the liver sinusoidal endothelial cells (LSECs) and functions as a scavenger receptor for various macromolecules including hyaluronans (HA). In DBA/2J mice, plasma HA concentration is ten times higher than in 129S6 or C57BL/6J mice, and this phenotype is genetically linked to the Stab2 locus. Stab2 mRNA in the LSECs was significantly lower in DBA/2J than in 129S6, leading to reduced STAB2 proteins in the DBA/2J LSECs. We found a retrovirus-derived transposable element, intracisternal A particle (IAP), in the promoter region of Stab2^DBA which likely interferes with normal expression in the LSECs. In contrast, in other tissues of DBA/2J mice, the IAP drives high ectopic Stab2^DBA transcription starting within the 5′ long terminal repeat of IAP in a reverse orientation and continuing through the downstream Stab2^DBA. Ectopic transcription requires the Stab2-IAP element but is dominantly suppressed by the presence of loci on 59.7–73.0 Mb of chromosome (Chr) 13 from C57BL/6J, while the same region in 129S6 requires additional loci for complete suppression. Chr13:59.9–73 Mb contains a large number of genes encoding Krüppel-associated box-domain zinc-finger proteins that target transposable elements-derived sequences and repress their expression. Despite the high amount of ectopic Stab2^DBA transcript in tissues other than liver, STAB2 protein was undetectable and unlikely to contribute to the plasma HA levels of DBA/2J mice. Nevertheless, the IAP insertion and its effects on the transcription of the downstream Stab2^DBA exemplify that stochastic evolutional events could significantly influence susceptibility to complex but common diseases.

Analysis of unmapped regions associated with long deletions in Korean whole genome sequences based on short read data

While studies aimed at detecting and analyzing indels or single nucleotide polymorphisms within human genomic sequences have been actively conducted, studies on detecting long insertions/deletions are not easy to orchestrate. For the last 10 years, the availability of long read data of human genomes from PacBio or Nanopore platforms has increased, which makes it easier to detect long insertions/deletions. However, because long read data have a critical disadvantage due to their relatively high cost, many next generation sequencing data are produced mainly by short read sequencing machines. Here, we constructed programs to detect so-called unmapped regions (UMRs, where no reads are mapped on the reference genome), scanned 40 Korean genomes to select UMR long deletion candidates, and compared the candidates with the long deletion break points within the genomes available from the 1000 Genomes Project (1KGP). An average of about 36,000 UMRs were found in the 40 Korean genomes tested, 284 UMRs were common across the 40 genomes, and a total of 37,943 UMRs were found. Compared with the 74,045 break points provided by the 1KGP, 30,698 UMRs overlapped. As the number of compared samples increased from 1 to 40, the number of UMRs that overlapped with the break points also increased. This eventually reached a peak of 80.9% of the total UMRs found in this study. As the total number of overlapped UMRs could probably grow to encompass 74,045 break points with the inclusion of more Korean genomes, this approach could be practically useful for studies on long deletions utilizing short read data.

Mutational landscape of T-cell lymphoma in mice lacking the DNA mismatch repair gene Mlh1: no synergism with ionizing radiation

Biallelic germline mutations in the DNA mismatch repair gene MLH1 lead to constitutional mismatch repair-deficiency syndrome and an increased risk for childhood hematopoietic malignancies, including lymphoma and leukemia. To examine how Mlh1 dysfunction promotes lymphoma as well as the influence of ionizing radiation (IR) exposure, we used an Mlh1-/- mouse model and whole-exome sequencing to assess genomic alterations in 23 T-cell lymphomas, including 8 spontaneous and 15 IR-associated lymphomas. Exposure to IR accelerated T-cell lymphoma induction in the Mlh1-/- mice, and whole-exome sequencing revealed that IR exposure neither increased the number of mutations nor altered the mutation spectrum of the lymphomas. Frequent mutations were evident in genes encoding transcription factors (e.g. Ikzf1, Trp53, Bcl11b), epigenetic regulators (e.g. Suv420h1, Ep300, Kmt2d), transporters (e.g. Rangap1, Kcnj16), extracellular matrix (e.g. Megf6, Lrig1), cell motility (e.g. Argef19, Dnah17), protein kinase cascade (e.g. Ptpro, Marcks) and in genes involved in NOTCH (e.g. Notch1), and PI3K/AKT (e.g. Pten, Akt2) signaling pathways in both spontaneous and IR-associated lymphomas. Frameshift mutations in mononucleotide repeat sequences within the genes Trp53, Ep300, Kmt2d, Notch1, Pten and Marcks were newly identified in the lymphomas. The lymphomas also exhibited a few chromosomal abnormalities. The results establish a landscape of genomic alterations in spontaneous and IR-associated lymphomas that occur in the context of mismatch repair dysfunction and suggest potential targets for cancer treatment.

The genomic pool of standing structural variation outnumbers single nucleotide polymorphism by threefold in the marine teleost Chrysophrys auratus

Recent studies have highlighted an important role of structural variation (SV) in ecological and evolutionary processes, but few have studied nonmodel species in the wild. As part of our long-term research programme on the nonmodel teleost fish Australasian snapper (Chrysophrys auratus), we aim to build one of the first catalogues of genomic variants (SNPs and indels, and deletions, duplications and inversions) in fishes and evaluate overlap of genomic variants with regions under putative selection (Tajima's D and π), and coding sequences (genes). For this, we analysed six males and six females from three locations in New Zealand and generated a high-resolution genomic variation catalogue. We characterized 20,385 SVs and found they intersected with almost a third of all annotated genes. Together with small indels, SVs account for three times more variation in the genome in terms of bases affected compared to SNPs. We found that a sizeable portion of detected SVs was in the upper and lower genomic regions of Tajima's D and π, indicating that some of these have an effect on the phenotype. Together, these results shed light on the often neglected genomic variation that is produced by SVs and highlights the need to go beyond the mere measure of SNPs when investigating evolutionary processes, such as species diversification and adaptation.

Fluctuating fortunes: genomes and habitat reconstructions reveal global climate-mediated changes in bats' genetic diversity

Over the last approximately 2.6 Myr, Earth's climate has been dominated by cyclical ice ages that have profoundly affected species' population sizes, but the impact of impending anthropogenic climate change on species' extinction potential remains a worrying problem. We investigated 11 bat species from different taxonomic, ecological and geographical backgrounds using combined information from palaeoclimatic habitat reconstructions and genomes to analyse biotic impacts of historic climate change. We discover tightly correlated fluctuations between species' historic distribution and effective population size, identify frugivores as particularly susceptible to global warming, pinpoint large insectivores as having overall low effective population size and flag the onset of the Holocene (approx. 10-12 000 years ago) as the period with the generally lowest effective population sizes across the last approximately 1 Myr. Our study shows that combining genomic and palaeoclimatological approaches reveals effects of climatic shifts on genetic diversity and may help predict impacts of future climate change.

Novel Genetic Rearrangements Termed "Structural Variation Polymorphisms" Contribute to the Genetic Diversity of Orthohepadnaviruses

The genetic diversity of orthohepadnaviruses is not yet fully understood. This study was conducted to investigate the role of structural variations (SVs) in their diversity. Genetic sequences of orthohepadnaviruses were retrieved from databases. The positions of sequence gaps were investigated, since they were found to be related to SVs, and they were further used to search for SVs. Then, a combination of pair-wise and multiple alignment analyses was performed to analyze the genomic structure. Unique patterns of SVs were observed; genetic sequences at certain genomic positions could be separated into multiple patterns, such as no SV, SV pattern 1, SV pattern 2, and SV pattern 3, which were observed as polymorphic changes. We provisionally referred to these genetic changes as SV polymorphisms. Our data showed that higher frequency of sequence gaps and lower genetic identity were observed in the pre-S1-S2 region of various types of HBVs. Detailed examination of the genetic structure in the pre-S region by a combination of pair-wise and multiple alignment analyses showed that the genetic diversity of orthohepadnaviruses in the pre-S1 region could have been also induced by SV polymorphisms. Our data showed that novel genetic rearrangements provisionally termed SV polymorphisms were observed in various orthohepadnaviruses.

Variant calling and quality control of large-scale human genome sequencing data

Next-generation sequencing has allowed genetic studies to collect genome sequencing data from a large number of individuals. However, raw sequencing data are not usually interpretable due to fragmentation of the genome and technical biases; therefore, analysis of these data requires many computational approaches. First, for each sequenced individual, sequencing data are aligned and further processed to account for technical biases. Then, variant calling is performed to obtain information on the positions of genetic variants and their corresponding genotypes. Quality control (QC) is applied to identify individuals and genetic variants with sequencing errors. These procedures are necessary to generate accurate variant calls from sequencing data, and many computational approaches have been developed for these tasks. This review will focus on current widely used approaches for variant calling and QC.

Comprehensive evaluation and characterisation of short read general-purpose structural variant calling software

In recent years, many software packages for identifying structural variants (SVs) using whole-genome sequencing data have been released. When published, a new method is commonly compared with those already available, but this tends to be selective and incomplete. The lack of comprehensive benchmarking of methods presents challenges for users in selecting methods and for developers in understanding algorithm behaviours and limitations. Here we report the comprehensive evaluation of 10 SV callers, selected following a rigorous process and spanning the breadth of detection approaches, using high-quality reference cell lines, as well as simulations. Due to the nature of available truth sets, our focus is on general-purpose rather than somatic callers. We characterise the impact on performance of event size and type, sequencing characteristics, and genomic context, and analyse the efficacy of ensemble calling and calibration of variant quality scores. Finally, we provide recommendations for both users and methods developers.

A number of computational methods have been developed for calling structural variants (SVs) using short read sequencing data. Here, the authors perform a comprehensive benchmarking analysis comparing 10 general-purpose callers and provide recommendations for both users and methods developers.

AluMine: alignment-free method for the discovery of polymorphic Alu element insertions

Background

Recently, alignment-free sequence analysis methods have gained popularity in the field of personal genomics. These methods are based on counting frequencies of short k-mer sequences, thus allowing faster and more robust analysis compared to traditional alignment-based methods.

Results

We have created a fast alignment-free method, AluMine, to analyze polymorphic insertions of Alu elements in the human genome. We tested the method on 2,241 individuals from the Estonian Genome Project and identified 28,962 potential polymorphic Alu element insertions. Each tested individual had on average 1,574 Alu element insertions that were different from those in the reference genome. In addition, we propose an alignment-free genotyping method that uses the frequency of insertion/deletion-specific 32-mer pairs to call the genotype directly from raw sequencing reads. Using this method, the concordance between the predicted and experimentally observed genotypes was 98.7%. The running time of the discovery pipeline is approximately 2 h per individual. The genotyping of potential polymorphic insertions takes between 0.4 and 4 h per individual, depending on the hardware configuration.

Conclusions

AluMine provides tools that allow discovery of novel Alu element insertions and/or genotyping of known Alu element insertions from personal genomes within few hours.

Predicting Local Inversions Using Rectangle Clustering and Representative Rectangle Prediction

As a specific type of structural variation, inversions are enjoying particular traction as a result of their established role in evolution. Using third-generation sequencing technology to predict inversions is growing in interest, but many such methods focus on improving sensitivity, giving rise to either too many false positives or very long running times. In this paper, we propose a new framework for inversion detection based on a combination of two novel theoretical models: rectangle clustering and representative rectangle prediction. This combination can automatically filter out false positive inversion predictions while retaining correct ones, leading to a method that has both high sensitivity and high positive prediction values (PPV). Further, this new framework can run very fast on available data. Our software can be freely obtained at https://github.com/UTbioinf/RigInv.

The Genome of C57BL/6J "Eve", the Mother of the Laboratory Mouse Genome Reference Strain

Isogenic laboratory mouse strains enhance reproducibility because individual animals are genetically identical. For the most widely used isogenic strain, C57BL/6, there exists a wealth of genetic, phenotypic, and genomic data, including a high-quality reference genome (GRCm38.p6). Now 20 years after the first release of the mouse reference genome, C57BL/6J mice are at least 26 inbreeding generations removed from GRCm38 and the strain is now maintained with periodic reintroduction of cryorecovered mice derived from a single breeder pair, aptly named Adam and Eve. To provide an update to the mouse reference genome that more accurately represents the genome of today's C57BL/6J mice, we took advantage of long read, short read, and optical mapping technologies to generate a de novo assembly of the C57BL/6J Eve genome (B6Eve). Using these data, we have addressed recurring variants observed in previous mouse genomic studies. We have also identified structural variations, closed gaps in the mouse reference assembly, and revealed previously unannotated coding sequences. This B6Eve assembly explains discrepant observations that have been associated with GRCm38-based analyses, and will inform a reference genome that is more representative of the C57BL/6J mice that are in use today.

Structural variants in 3000 rice genomes

Investigation of large structural variants (SVs) is a challenging yet important task in understanding trait differences in highly repetitive genomes. Combining different bioinformatic approaches for SV detection, we analyzed whole-genome sequencing data from 3000 rice genomes and identified 63 million individual SV calls that grouped into 1.5 million allelic variants. We found enrichment of long SVs in promoters and an excess of shorter variants in 5′ UTRs. Across the rice genomes, we identified regions of high SV frequency enriched in stress response genes. We demonstrated how SVs may help in finding causative variants in genome-wide association analysis. These new insights into rice genome biology are valuable for understanding the effects SVs have on gene function, with the prospect of identifying novel agronomically important alleles that can be utilized to improve cultivated rice.

Mild inborn errors of metabolism in commonly used inbred mouse strains

Inbred mouse strains are a cornerstone of translational research but paradoxically many strains carry mild inborn errors of metabolism. For example, α-aminoadipic acidemia and branched-chain ketoacid dehydrogenase deficiency are known in C57BL/6J mice. Using RNA sequencing, we now reveal the causal variants in Dhtkd1 and Bckdhb, and the molecular mechanism underlying these metabolic defects. C57BL/6J mice have decreased Dhtkd1 mRNA expression due to a solitary long terminal repeat (LTR) in intron 4 of Dhtkd1. This LTR harbors an alternate splice donor site leading to a partial splicing defect and as a consequence decreased total and functional Dhtkd1 mRNA, decreased DHTKD1 protein and α-aminoadipic acidemia. Similarly, C57BL/6J mice have decreased Bckdhb mRNA expression due to an LTR retrotransposon in intron 1 of Bckdhb. This transposable element encodes an alternative exon 1 causing aberrant splicing, decreased total and functional Bckdhb mRNA and decreased BCKDHB protein. Using a targeted metabolomics screen, we also reveal elevated plasma C5-carnitine in 129 substrains. This biochemical phenotype resembles isovaleric acidemia and is caused by an exonic splice mutation in Ivd leading to partial skipping of exon 10 and IVD protein deficiency. In summary, this study identifies three causal variants underlying mild inborn errors of metabolism in commonly used inbred mouse strains.

PiggyBac transposon tools for recessive screening identify B-cell lymphoma drivers in mice

B-cell lymphoma (BCL) is the most common hematologic malignancy. While sequencing studies gave insights into BCL genetics, identification of non-mutated cancer genes remains challenging. Here, we describe PiggyBac transposon tools and mouse models for recessive screening and show their application to study clonal B-cell lymphomagenesis. In a genome-wide screen, we discover BCL genes related to diverse molecular processes, including signaling, transcriptional regulation, chromatin regulation, or RNA metabolism. Cross-species analyses show the efficiency of the screen to pinpoint human cancer drivers altered by non-genetic mechanisms, including clinically relevant genes dysregulated epigenetically, transcriptionally, or post-transcriptionally in human BCL. We also describe a CRISPR/Cas9-based in vivo platform for BCL functional genomics, and validate discovered genes, such as Rfx7, a transcription factor, and Phip, a chromatin regulator, which suppress lymphomagenesis in mice. Our study gives comprehensive insights into the molecular landscapes of BCL and underlines the power of genome-scale screening to inform biology.

A virome-wide clonal integration analysis platform for discovering cancer viral etiology

Oncoviral infection is responsible for 12%–15% of cancer in humans. Convergent evidence from epidemiology, pathology, and oncology suggests that new viral etiologies for cancers remain to be discovered. Oncoviral profiles can be obtained from cancer genome sequencing data; however, widespread viral sequence contamination and noncausal viruses complicate the process of identifying genuine oncoviruses. Here, we propose a novel strategy to address these challenges by performing virome-wide screening of early-stage clonal viral integrations. To implement this strategy, we developed VIcaller, a novel platform for identifying viral integrations that are derived from any characterized viruses and shared by a large proportion of tumor cells using whole-genome sequencing (WGS) data. The sensitivity and precision were confirmed with simulated and benchmark cancer data sets. By applying this platform to cancer WGS data sets with proven or speculated viral etiology, we newly identified or confirmed clonal integrations of hepatitis B virus (HBV), human papillomavirus (HPV), Epstein-Barr virus (EBV), and BK Virus (BKV), suggesting the involvement of these viruses in early stages of tumorigenesis in affected tumors, such as HBV in TERT and KMT2B (also known as MLL4) gene loci in liver cancer, HPV and BKV in bladder cancer, and EBV in non-Hodgkin's lymphoma. We also showed the capacity of VIcaller to identify integrations from some uncharacterized viruses. This is the first study to systematically investigate the strategy and method of virome-wide screening of clonal integrations to identify oncoviruses. Searching clonal viral integrations with our platform has the capacity to identify virus-caused cancers and discover cancer viral etiologies.

iMGEins: detecting novel mobile genetic elements inserted in individual genomes

Background

Recent advances in sequencing technology have allowed us to investigate personal genomes to find structural variations, which have been studied extensively to identify their association with the physiology of diseases such as cancer. In particular, mobile genetic elements (MGEs) are one of the major constituents of the human genomes, and cause genome instability by insertion, mutation, and rearrangement.

Result

We have developed a new program, iMGEins, to identify such novel MGEs by using sequencing reads of individual genomes, and to explore the breakpoints with the supporting reads and MGEs detected. iMGEins is the first MGE detection program that integrates three algorithmic components: discordant read-pair mapping, split-read mapping, and insertion sequence assembly. Our evaluation results showed its outstanding performance in detecting novel MGEs from simulated genomes, as well as real personal genomes. In detail, the average recall and precision rates of iMGEins are 96.67 and 100%, respectively, which are the highest among the programs compared. In the testing with real human genomes of the NA12878 sample, iMGEins shows the highest accuracy in detecting MGEs within 20 bp proximity of the breakpoints annotated.

Conclusion

In order to study the dynamics of MGEs in individual genomes, iMGEins was developed to accurately detect breakpoints and report inserted MGEs. Compared with other programs, iMGEins has valuable features of identifying novel MGEs and assembling the MGEs inserted.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5290-9) contains supplementary material, which is available to authorized users.

Gene exchange between two divergent species of the fungal human pathogen, Coccidioides

The fungal genus Coccidioides is composed of two species, Coccidioides immitis and Coccidioides posadasii. These two species are the causal agents of coccidioidomycosis, a pulmonary disease also known as valley fever. The two species are thought to have shared genetic material due to gene exchange in spite of their long divergence. To quantify the magnitude of shared ancestry between them, we analyzed the genomes of a population sample from each species. Next, we inferred what is the expected size of shared haplotypes that might be inherited from the last common ancestor of the two species and find a cutoff to find what haplotypes have conclusively been exchanged between species. Finally, we precisely identified the breakpoints of the haplotypes that have crossed the species boundary and measure the allele frequency of each introgression in this sample. We find that introgressions are not uniformly distributed across the genome. Most, but not all, of the introgressions segregate at low frequency. Our results show that divergent species can share alleles, that species boundaries can be porous, and highlight the need for a systematic exploration of gene exchange in fungal species.

Characterization of hepatitis B virus with complex structural variations

BACKGROUND

Hepatitis B virus (HBV) infection is one of the most serious public health issues. Recent HBV genetic research has revealed novel genetic rearrangements termed complex structural variations (SVs), which are composed of combinations of SVs such as insertions, deletions, and duplications. An extensive search was made for complex SVs of HBV and their characteristics were analyzed.

RESULTS

Fifty-five HBV strains with complex SVs were identified by analyzing genetic sequences of HBV with bioinformatical tools. Along with 15 HBV strains with complex SVs in a previous report, a total of 70 HBV strains harboring complex SVs were analyzed. Complex SVs in the HBV genome were located frequently between nt 1500 and 2000. Insertions were observed in 65/70 (92.9%) of HBV strains with complex SVs. As insertional motif sequences, hepatocyte nuclear factor 1 binding site, a sequence complementary to part of box α in enhancer II, and insertions of unknown origins were observed. The complex SVs were classified into six groups, and combination of insertion and deletion was observed more frequently than other patterns.

CONCLUSION

Through an extensive search of HBV sequences, new strains with complex SVs were identified in this study. Characteristics of HBV with complex SVs were clarified by the analysis of 70 HBV strains harboring complex SVs. Further investigation is required to elucidate its role in pathogenesis of HBV-related liver disease.

Structural Variant Prediction in Extended Pedigrees Through Sparse Negative Binomial Genome Signal Recovery

Structural variants (SVs) are rearrangements, such as deletions, insertions, duplications, inversions, and translocations, in an individual's genome relative to a reference. SV detection is often marred by high false positive rates due to errors in sequencing and mapping. In previous work, we proposed a maximum likelihood approach to SV prediction that incorporated low-coverage sequencing data and coverage distribution. In particular, we developed a negative binomial framework to reflect a more realistic representation DNA fragment distributions sampled from an individual's genome. In this paper, we leverage relationships between an off spring and both parents, in addition to the negative binomial framework, to improve SV identification accuracy. We present numerical results on both simulated genomes as well as two sequenced parent-child trios from the 1000 Genomes Project.

Detection of large sequence insertions by a hybrid approach that combine de novo assembly and resequencing of medium-coverage genome sequences

Large sequence insertion (LSI) is one of the structural variations (SVs) that may cause phenotypic differences in plants. To identify the LSIs using medium-coverage sequencing data of four wild soybean (Glycine soja) genotypes, we designed a hybrid approach combining de novo assembly and read mapping. Total reads and reads with both ends unmapped were independently assembled into "ordinary contigs" and "orphan contigs", respectively, and subjected to pairwise alignment and stringent filtering. This approach predicted 24 LSIs averaging 2682 bp in size, with no overlap with SVs detected by Pindel, BreakDancer, or ScanIndel, and they were validated by PCR. Compared with the soybean (Glycine max) reference genome, 20 LSIs were located outside genic regions. One of the four LSIs within a genic region, LSI05, is located in the coding DNA sequence region of a protein kinase superfamily gene (Glyma.08G123500). It caused delayed translation initiation and loss of 24 amino acids in the wild soybean genotype CW12. LSI05 was more frequently observed in 29 G. soja accessions than in 34 G. max accessions. Identified LSIs would be genomic resources harboring novel gene contents for studying SVs and improving crops. Moreover, our cost-efficient approach may be applicable to other plant species.

Genome-Wide Analysis of Genetic Variations and the Detection of Rich Variants of NBS-LRR Encoding Genes in Common Wild Rice Lines

Common wild rice (Oryza rufipogon Griff.) is invaluable genetic resource for rice resistance breeding. Whole-genome re-sequencing was conducted to systematically analyze the variations in two new inbred lines (Huaye 3 and Huaye 4) developed from a common wild rice. A total of 4,841,127 SNPs, 1,170,479 InDels, 24,080 structural variations (SVs), and 298 copy number variations (CNVs) were identified in three materials. Approximately 16.24 and 5.64% of the total SNPs and InDels of Huaye 3 and Huaye 4 were located in genic regions, respectively. Together, 12,486 and 15,925 large-effect SNPs, and 12,417 and 14,513 large-effect InDels, which affect the integrity of the encoded protein, were identified in Huaye 3 and Huaye 4, respectively. The distribution map of 194 and 245 NBS-LRR encoding homologs was constructed across 12 rice chromosomes. Further, GO enrichment analysis of the homologs with identical genotype variations in Huaye 3 and Huaye 4 revealed 67, 82, and 58 homologs involved in cell death, response to stress, and both terms, respectively. Comparative analysis displayed that 550 out of 652 SNPs and 129 out of 147 InDels were present in a widely used blast-susceptible rice variety (LTH). Protein-protein interaction analysis revealed a strong interaction between NBS-LRR candidates and several known R genes. One homolog of disease resistance protein (RPM1) was involved in the plant-pathogen interaction pathway. Artificial inoculation of disease/insect displayed resistance phenotypes against rice blast and brown planthopper in two lines. The results will provide allele-specific markers for rice molecular breeding.