Identification of High-Confidence Structural Variants in Domesticated Rainbow Trout Using Whole-Genome Sequencing

Genomic structural variants (SVs) are a major source of genetic and phenotypic variation but have not been investigated systematically in rainbow trout (Oncorhynchus mykiss), an important aquaculture species of cold freshwater. The objectives of this study were 1) to identify and validate high-confidence SVs in rainbow trout using whole-genome re-sequencing; and 2) to examine the contribution of transposable elements (TEs) to SVs in rainbow trout. A total of 96 rainbow trout, including 11 homozygous lines and 85 outbred fish from three breeding populations, were whole-genome sequenced with an average genome coverage of 17.2×. Putative SVs were identified using the program Smoove which integrates LUMPY and other associated tools into one package. After rigorous filtering, 13,863 high-confidence SVs were identified. Pacific Biosciences long-reads of Arlee, one of the homozygous lines used for SV detection, validated 98% (3,948 of 4,030) of the high-confidence SVs identified in the Arlee homozygous line. Based on principal component analysis, the 85 outbred fish clustered into three groups consistent with their populations of origin, further indicating that the high-confidence SVs identified in this study are robust. The repetitive DNA content of the high-confidence SV sequences was 86.5%, which is much higher than the 57.1% repetitive DNA content of the reference genome, and is also higher than the repetitive DNA content of Atlantic salmon SVs reported previously. TEs thus contribute substantially to SVs in rainbow trout as TEs make up the majority of repetitive sequences. Hundreds of the high-confidence SVs were annotated as exon-loss or gene-fusion variants, and may have phenotypic effects. The high-confidence SVs reported in this study provide a foundation for further rainbow trout SV studies.

Introduction

Structural variants (SVs) refer to sequence variants greater than 50 bp in size (Alkan et al., 2011; Ho et al., 2020). The most common SVs include deletions, duplications, inversions, insertions and translocations. SVs are a major source of genetic variation, and affect a larger proportion of the human genome than single nucleotide polymorphism (SNP) and any other genetic variants (Sudmant et al., 2015). SVs can influence the expression of genes (Chiang et al., 2017; Spielmann et al., 2018), and impact numerous phenotypes in both plant and mammals (Spielmann and Mundlos, 2013; Deane-Coe et al., 2018; Zhou et al., 2019; Alonge et al., 2020; Liu et al., 2020). Therefore, it is essential to study SVs in order to explore the full spectrum of genetic variation.

Structural variants are poorly characterized in most species due to both financial and technical constraints. Only a few SV studies were reported in aquaculture species. Based on whole-genome sequences of 492 Atlantic salmon (Salmo salar L.), 15,483 high-confidence SVs were identified (Bertolotti et al., 2020). Using a combination of reduced representation sequencing and sequence capture enrichment, SNPs and SVs were identified in American lobster (Homarus americanus) (Dorant et al., 2020). This study also reported that 46 SVs were significantly associated with the annual variance of sea surface temperature while SNPs failed to reveal population adaption to local environments.

Repetitive sequences, especially transposable elements (TEs), are a major component of plant and animal genomes (Lisch, 2013; Sotero-Caio et al., 2017), and TEs are an important source for SVs (Payer and Burns, 2019). Initially, specialized microarrays were used to detect genome-wide TE insertions in the human genome (Huang et al., 2010). With the recent advancement of genome sequencing technology, both short-reads and long-reads can be used to detect TE insertions (Zhou et al., 2020; Zook et al., 2020). Computer programs specialized for the detection of mobile element insertions (MEIs) from whole-genome sequences have been developed (Chu et al., 2020), and were used to detect MEIs in many species such as human (Sudmant et al., 2015), grape (Zhou et al., 2019), rice (Kou et al., 2020), and tomato (Alonge et al., 2020; Dominguez et al., 2020).

Rainbow trout (Oncorhynchus mykiss) is one of the most widely cultured cold freshwater fish, with production on every continent except Antarctica (Food and Agriculture Organization (FAO), 2021). SVs are a largely unexplored feature in rainbow trout. Classical cytogenetic studies revealed chromosomal number variation from 2N = 58 to 2N = 64 in rainbow trout (Thorgaard, 1983). Recently, comparison of the de novo genome assembly and SNP linkage maps revealed a 55 Mb double inversion associated with a migratory trait in rainbow trout on chromosome Omy05 and another 14 Mb inversion on chromosome Omy20 (Pearse et al., 2019). However, a systematic study of SVs in rainbow trout has not been reported to date. The objectives of this study were 1) to identify and validate SVs in rainbow trout using whole-genome re-sequencing; and 2) to examine the contribution of TEs to SVs in rainbow trout. The high-confidence SVs reported in this study are useful resources to understand and explore the genetic variation in rainbow trout, and provide a foundation for further rainbow trout SV studies.

Materials and Methods

Ethics Statement

The fish used in this study were sampled according to our standard operating procedures of care and use of research animals, approved by the Institutional Animal Care and Use Committee, National Center for Cool and Cold Water Aquaculture, Agricultural Research Service, United States Department of Agriculture, United States of America. All efforts were made to minimize suffering and to ensure fish welfare.

Whole-Genome Sequencing and Read Mapping

Three groups of rainbow trout (Supplementary Table S1), including 11 double haploid (DH) lines, 39 outbred fish from the National Center for Cool and Cold Water Aquaculture (NCCCWA) (Leeds et al., 2010, 2016; Wiens et al., 2018) and 46 outbred fish from the Troutlodge May spawning (TLUM) population, were used in this study. DNA was extracted from fin clips following the manufacturer’s recommended protocols for AutoGenprep 965 (Autogen, Holliston, MA, United States). Sequencing and read mapping of the 11 DH lines were described in our previous study (Gao et al., 2018). For both the NCCCWA and TLUM samples, whole-genome DNA sequencing libraries were prepared using the KAPA HyperPrep kit (KAPA Biosystems, Wilmington, MA, United States), and were sequenced in paired-end (2 × 150 bp) mode on an Illumina HiSeq X sequencer. The sequence reads were mapped to the rainbow trout reference genome GCF_002163495.1 (Pearse et al., 2019) using BWA-MEM algorithm (Li, 2013), and alignments were converted to BAM (Binary sequence Alignment/Map) format using SAMtools v1.11 (Li et al., 2009). PCR duplicates were marked and removed using Picard¹ v2.18.2. Only reads mapped to the 29 chromosomes were used for SV detection as described below, and scaffolds that cannot be assigned to the 29 chromosomes and mitochondria sequence were excluded from SV detection.

SV Detection and Genotyping

We chose to use LUMPY (Layer et al., 2014) to detect SVs in this study because it integrates multiple signals across multiple samples to detect SVs, and several associated tools, as mentioned below, are available to facilitate SV detection. To reduce false positive SV calls, reads overlapping with problematic regions, including extreme high read coverage regions, low complexity regions and gaps in the reference genome, were excluded from SV detection. To identify extreme high-read coverage regions, mosdepth v0.2.9 (Pedersen and Quinlan, 2018) was used to estimate read coverage using the quantize option, and regions with coverage equal or above 200× were converted into a BED (Brower Extensible Data) file. RepeatMasker v4.1.0 (Smit et al., 2020) was used to mask the reference genome GCF_002163495.1, and the low complexity regions were extracted from the output file of RepeatMasker. The low complexity regions separated by less than 1,000 bp were merged using Bedtools v2.27.1 (Quinlan and Hall, 2010), and intervals greater than 200 bp were then retained. The BED file for gaps in the reference genome was generated using a Python script (Bertolotti et al., 2020). The BED files for the three problematic regions were combined, and regions separated by less than 1,000 bp were merged using Bedtools. As recommended in the user instructions of LUMPY, program Smoove² v0.2.5, which integrates the best practices of LUMPY and other associated tools, was used to detect SVs. We followed the four SV calling steps at the population level as described in the Smoove user guide. Briefly, SVs were called initially in each sample, and then SVtools v0.5.1 (Larson et al., 2019) was used to merge the SV calls across all samples. The SVs were genotyped using SVtyper v0.7.1 (Chiang et al., 2015) and annotated with read coverage using Duphold v0.2.1 (Pedersen and Quinlan, 2019).

To identify high-confidence SVs, we used the following steps to filter the putative SVs detected above: 1) SVs with size less than 300 bp were removed using BCFtools v1.9; 2) SVs overlapping with the extreme high-read coverage (equal or greater than 200×) regions or the gaps of the reference genome were excluded using Bedtools; 3) SVs were filtered based on Duphold annotation. Deletions with DHFFC (fold change for the variant depth relative to flanking regions) > 0.7, duplications with DHFFC < 1.3 and inversions with DHFFC > 1.3 or DHFFC < 0.7 in at least one samples were excluded using BCFtools; 4) SVs with a heterozygous genotype in at least one of the 11 DH lines were also excluded; 5) Variant allele frequencies (AF) were calculated using SVtools subcommand afreq, and SVs with AF < 0.05 or AF > 0.95 were excluded; 6) Genotype call rate for each remaining SV was calculated using the computer program KING v2.2.5 (Manichaikul et al., 2010), and SVs with genotype call rates less than 0.8 were filtered out.

To examine the repetitive DNA content of the SV regions, we extracted the SV sequences from the reference genome GCF_002163495.1 using Bedtools. The rainbow trout repetitive sequence database reported by Pearse et al. (2019) was used to mask the SV sequences using RepeatMasker v4.1.0 with default parameters.

SV Validation Using Long-Reads

Deep-coverage (111.6×) Pacific Biosciences (PacBio) long-reads of the Arlee homozygous line (unpublished data) were used to validate the high-confidence SVs present in the Arlee line. The long-reads were mapped to the reference genome GCF_002163495.1 using minimap2 v2.17-r941 (Li, 2018). The Python scripts and methods reported by Abel et al. (2020) were then used for SV validation. Briefly, the Python scripts were used to generate a BEDPE (paired-end BED) file to represent the breakpoints of split reads (based on alignment with the reference genome) and reads with insertion or deletion (based on cigar strings) larger than 50 bp, respectively. To support a putative SV call, the long-reads should have breakpoints overlapping with the putative SVs based on Bedtools (“pairtopair -type both -sloop 100 -is” or “intersect -r -f 0.9”). Also, the long-reads and the putative SVs should have the same type of breakpoints such as deletions or duplications. Only putative SVs supported by at least two long-reads were considered as validated SVs.

Principal Component Analysis

The SV genotypes of the NCCCWA and TLUM fish were used for principal analysis using PLINK v1.9 (Chang et al., 2015) with the top 10 principal components. The 11 DH lines were excluded from the principal component analysis because they were derived from diverse populations.

SV Annotation

We used SnpEff v5.0 (Cingolani et al., 2012) to annotate the high-confidence SVs. The GTF file of the rainbow trout reference genome, GCF_002163495.1, was downloaded from NCBI, and was used to build a reference database. The high-confidence SVs were then annotated using the default parameters of SnpEff.

Results

Identification and Validation of SVs in Rainbow Trout

After running the four SV calling steps at the population level using Smoove, a total of 298,410 putative SVs including 253,325 deletions, 37,648 duplications and 7,434 inversions were identified among the 96 fish used in this study. After rigorous filtering as described in the methods, 13,863 high-confidence SVs including 13,138 deletions, 719 duplications and 6 inversions were retained (Table 1). The average spacing across all chromosomes was 126.6 kb per high-confidence SV (Table 1). The number of high-confidence SVs per sample ranged from 2,618 to 5,924 (Supplementary Table S1) with a median of 5,283 high-confidence SVs per sample. The SV sizes ranged from 300 to 27,641 bp, and 90% of the SVs are shorter than 3,000 bp. The SV size distributions for deletions and duplications are shown in Figures 1A and B, respectively.

TABLE 1

Table 1. Number of the high-confidence SVs identified in rainbow trout.*

FIGURE 1

Figure 1. Size distributions of the high-confidence SVs identified in rainbow trout. (A) Deletions; (B) Duplications.

Among the 13,863 high-confidence SVs identified, 4,030 SVs were homozygous for the non-reference alleles in the Arlee line. We used the PacBio long-reads of the Arlee line to validate the high-confidence SVs, and 3,948 SVs (98.0%) were validated (Table 2).

TABLE 2

Table 2. Number of validated SVs using long-reads of the Arlee homozygous line.

Substantial TE Contribution to SVs in Rainbow Trout

Several peaks were observed in the size distribution of the high-confidence deletions (Figure 1A). The tallest peak was in the 1,610–1,630 bp size range. Among the 1,072 sequences in this size range, 1,030 (96%) sequences had high sequence similarity (identity > 95%) with the full sequence of omyk_TCE_37 (sequence length 1,621 bp, Supplementary Note), the most abundant Tc1-mariner transposon family in rainbow trout (Pearse et al., 2019). This deletion size spike was not unexpected as the Tc1-mariner family omyk_TCE_37 has been active recently and makes up more than 2% of the rainbow trout genome (Pearse et al., 2019).

The next most prominent peak was in the 370–390 bp size range. There were 528 high-confidence deletions in this size range. BLASTN searches against the repeat database of rainbow trout (Pearse et al., 2019) revealed that 250 (47%) of those sequences aligned well (identity > 80%) with the long terminal repeat (LTR) of the Gypsy retrotransposon Omy_1322 (Supplementary Note). The third most prominent peak was in the 1,460–1,470 bp size range. Among the 181 deletion sequences in this size range, 148 (82%) of those sequences were homologous to an unclassified TE sequence (Supplementary Note) in rainbow trout.

Since most sequences in the three peak regions were homologous to TE sequences in rainbow trout, we decided to examine the proportion of repetitive sequences among all high confident SVs reported in this study. As shown in Table 3, 88.3% nucleotides of the deletions were masked by RepeatMasker. For the duplications and inversions, 73.0 and 74.0% nucleotides were masked, respectively. Overall, 86.5% nucleotides of the high-confidence SVs were masked by RepeatMasker, which is much higher than 57.1% repetitive DNA content in the rainbow trout genome (Pearse et al., 2019). Since most repetitive sequences were TEs, TEs contributed substantially to the high-confidence SVs reported in this study.

TABLE 3

Table 3. Repetitive DNA content of the high-confidence SV sequences.*

Principal Component Analysis

True SVs can capture the population structure of animals (Bertolotti et al., 2020) and plants (Zhou et al., 2019). The 85 outbred fish used in this study were selected from three distinct rainbow trout populations, NCCCWA even-year population, NCCCWA odd-year population and Troutlodge odd-year May spawning population. Principal component analysis was performed to test whether we can group the fish into three clusters consistent with their populations of origin. As shown in Figure 2, the first principal component clearly separated the Troutlodge population from the two NCCCWA populations, and the second principal component separated the two NCCCWA year-class populations. This result further confirm that the high-confidence SVs reported in this study are robust.

FIGURE 2

Figure 2. Scatter plot of the 85 outbred fish based on the top two principal components. TLUM fish are represented by yellow dots. Red squares and purple triangles represent the NCCCWA odd-year and even-year fish, respectively.

Annotation of the High-Confidence SVs in Rainbow Trout

The computer program SnpEff (Cingolani et al., 2012) was used to annotate the high-confidence SVs reported in this study. Four impact groups, high, moderate, low and modifier, were used to annotate the putative SV effects. Although 95.53% of the SV effects fell into the modifier group, 3.43% of the SV effects were in the high impact group, and 0.15% SVs were annotated with moderate effect. The rest 0.89% of the SV effects were in the low impact group. Hundreds of SVs were annotated as exon loss variants (1.5%) or gene fusion variants (0.4%) although most SVs were annotated as intron variants (45.5%) or intergenic variants (25.3%) (Supplementary Table S2). Therefore, the high-confidence SVs reported in this study provide a useful genomic resource for further investigation of SV impacts on phenotypes in rainbow trout.

Discussion

SVs in Rainbow Trout

Structural variants are a major source of genetic and phenotypic variation but are largely unexplored in many species, including rainbow trout. In this study, 13,863 high-confidence SVs were identified in rainbow trout using whole-genome sequencing. Both the high SV validation rate and principal component analysis indicated that the high-confidence SVs are robust. Although most SVs are in introns and intergenic regions, hundreds were annotated as high impact variants. Thus, the high-confidence SVs reported in this study provide a useful genomic resource for future SV studies in rainbow trout. Moreover, the outbred fish used in this study were from three breeding populations. Therefore, the high-confidence SVs reported in this study are relevant for the identification of SVs associated with traits targeted by selective breeding. However, there are some limitations with the approaches used in this study, and should be improved in the further studies. Although LUMPY is a widely used SV callers, it has some limitations regarding to SV sensitivity and accuracy just like any other SV callers (Layer et al., 2014; Cameron et al., 2019). Also, rigorous filtering steps were used to remove false-positive SVs, some true SVs might be filtered out. Furthermore, DNA sequencing technologies and SV detection algorithms are evolving rapidly (Ho et al., 2020). Thus, the full SV spectrum in rainbow trout remains to be explored.

Each DH line is formed from the duplication of a single gametic cell and hence contains two nearly identical copies of maternal or paternal genome (Parsons and Thorgaard, 1985). It is not surprising that more high-confidence SVs were obtained in each outbred fish (ranged from 4,809 to 5,924) than in each DH line (ranged from 2,618 to 4,149) (Supplementary Table S1). The SVs in each DH line were derived from only one haploid genome, whereas the SVs in each outbred diploid fish came from both the maternal and paternal genomes. The number of SVs in each DH line is largely determined by its genetic distance from the Swanson DH line, which was used for the rainbow trout reference genome assembly. Over 4,000 SVs were observed in each of the Arlee, Hot creek, OSU and Golden lines, which is much more than the number of SVs in each of the other DH lines (2,618–3,783). This observation is consistent with the previously published phylogenic tree of the DH lines based on SNP genotypes (Palti et al., 2014). The four DH lines are more distantly related to the Swanson line than the other DH lines used in this study.

It is essential to remove false positive SVs to retain only high-confidence SVs (Cameron et al., 2019; Ho et al., 2020). Among the 4,030 high-confidence SVs present in the Arlee DH line, 98% SVs were validated using PacBio long-reads. The high SV validation rate is likely due to multiple reasons. Both read mapping and SV detection are much easier for the Arlee DH line, which contains only homozygous genotypes. Also, it is unlikely that true SVs could not be validated due to a lack of long-read coverage in the SV regions because deep genome coverage (111.6×) of the PacBio long-reads were used for SV validation. Additional factors also contributed to this high SV validation rate. First, extreme high read coverage regions and low complexity regions were excluded from SV discovery as described in the section “Materials and Methods,” and putative SVs spanning extreme high read coverage regions were also removed from the putative SVs. The high read coverage regions are known for a high false discovery rate (FDR), which was as high as 99.2% in Atlantic salmon (Bertolotti et al., 2020). Second, the program Smoove, which was used to detect and filter SVs, integrates LUMPY, a widely used SV caller, and Duphold, a program to annotate SVs with the depth of sequence coverage, into one package. As described in the section “Materials and Methods,” low-confidence SVs were filtered out using the sequence depth annotation provided by the program Duphold. Third, additional stringent filtering steps were used to exclude low-confidence SVs. Because small size SVs have higher FDR (Cameron et al., 2019), putative SVs smaller than 300 bp were removed from the dataset. Furthermore, SVs with heterozygous genotypes in any of the 11 DH lines were removed. At last, SVs with AF < 0.05 or AF > 0.95 or SVs with genotype call rates <0.8 were also filtered out. Thus, the high-confidence SVs reported in this study are robust, and are useful resources for understanding and exploring the genetic variation in rainbow trout.

TEs Contribute Substantially to SVs in Rainbow Trout

Repetitive sequences, especially TEs, have been shown to be a major component of genome architecture, diversity and evolution (Lisch, 2013; Sotero-Caio et al., 2017). The rainbow trout genome contains 57.1% repetitive DNA (Pearse et al., 2019). In this study the repetitive DNA content is much higher in the SV sequences as 85.6% of their nucleotides were masked by RepeatMasker. Since most of the repeat sequences are TEs, TEs contribute substantially to SVs in rainbow trout. Consistent with the results of this study, TE enrichment in SVs was also observed in cichlids species (Faber-Hammond et al., 2019; Penso-Dolfin et al., 2020). However, the extent of TE contribution to SVs might vary across species. Although deletions caused by a recently active DNA transposon were observed in Atlantic salmon, 65% of the high-confidence SVs identified in Atlantic salmon contained no repeat sequences (Bertolotti et al., 2020). Therefore, rainbow trout and Atlantic salmon likely have quite different TE content in their high-confidence SV regions despite their close genetic relationship (Macqueen and Johnston, 2014) and similar repetitive DNA content (Lien et al., 2016; Pearse et al., 2019). Among the phase 3 release of human SVs based on the 1000 Genome Project, about 24% SVs belonged to the MEI class (Sudmant et al., 2015). In the latest release of human SVs (Abel et al., 2020), MEIs accounted for 27% SVs. In pig, Zhao et al. (2016) reported that 54.4% deletions were caused by TEs. In tomato, Alonge et al. (2020) reported that 84% deletions and 76% insertions larger than 100 bp contained at least one TE.

The exact amount of TE contribution to SVs in rainbow trout needs to be further investigated. Due to the technical constraints such as short-reads and read coverage and the SV detection method used in this study, the SV type and size biases are unavoidable. However, the contribution of TEs to SVs in rainbow trout might be underestimated due to the following reasons. First, extreme high read coverage regions which are enriched with TEs were excluded from the SV discovery, and putative SVs spanning these regions were excluded during the SV filtering process. Furthermore, there are still thousands of gaps in the current reference genome assembly (GCF_002163495.1) due in large part to the difficulty to assemble repeat and TE rich regions using Illumina short-reads. Thus, SVs located in the assembly gaps cannot be identified, which likely leads to underestimated TE content of the SV sequences. Second, the program LUMPY used for SV identification cannot detect sequence insertions. So, non-reference MEIs cannot be identified in this study; Third, low frequency SVs were filtered out. It is known that the bulk of MEIs identified in human have low variant allele frequencies (Sudmant et al., 2015; Wildschutte et al., 2016). It was possible that many SVs containing TEs were filtered out due to low allele frequencies; Lastly, long-read DNA sequencing was shown to be more effective than short-reads for detecting SVs that contain TE sequences (Zhou et al., 2020). However, short-reads were used for SV discovery in this study as the short-read sequencing technology is still much more affordable for whole-genome sequencing of multiple samples. Future studies using long-read sequencing technology will likely allow for more accurate quantification and characterization of the impact of TEs on SVs in rainbow trout.

The contribution of TEs to SVs in rainbow trout is not limited to a specific TE class. Both retrotransposon and DNA transposons derived SVs were identified in this study. The most prominent peak with a deletion size range 1,610–1,630 bp (Figure 1A) was due to a DNA transposon. Near half of the deletion sequences in the size range 370–390 bp aligned well with the LTR of Gypsy retrotransposon Omy_1322. These SV sequences were not full-length Gypsy retrotransposons, and they were solo LTRs. Solo LTRs are derived from ectopic recombination between LTRs (Thomas et al., 2018), and about 90% LTR retrotransposons in human genome are represented by solo LTRs (Subramanian et al., 2011). Since we did not observe deletions matching the full length of the Gypsy retrotransposon Omy_1322, we speculate that the majority of this retrotransposon sequences in rainbow trout are solo LTRs. Furthermore, among all the high-confidence SVs reported in this study, 25.5% nucleotides were masked by retrotransposons, 41.2% nucleotides were masked by DNA transposons, and 18.9% nucleotides were masked by unclassified TEs.

Conclusion

Structural variants are increasingly being recognized for their importance in evolutionary processes and impact on variation of economically important traits in plants and animals. In this paper we report for the first time on the discovery and validation of high-confidence SVs in rainbow trout which provides a foundation for further studies on their impact on phenotypic variation in this species. Our report also includes the finding of substantial TE contribution to SVs in rainbow trout, which provides an opening for further research on the role of TEs in the evolution of salmonid species given the smaller TE content in the high-confidence SVs in Atlantic salmon reported previously.

Data Availability Statement

The raw sequences of the 85 outbred fish were submitted to NCBI SRA (BioProject PRJNA681179), and the raw sequences of the 11 double haploid lines are also available in NCBI SRA (BioProject PRJNA386519). The PacBio long-reads of Arlee line are available in NCBI SRA under BioProject PRJNA623027. The VCF file for the high-confidence SVs reported in this study was submitted to the European Variation Archive under project PRJEB41774.

Ethics Statement

The fish used in this study were sampled according to our standard operating procedures of care and use of research animals, approved by the Institutional Animal Care and Use Committee, National Center for Cool and Cold Water Aquaculture, Agricultural Research Service, United States Department of Agriculture, United States of America. All efforts were made to minimize suffering and to ensure fish welfare.

Author Contributions

SL and YP conceived and planned the experiment. SL performed the data analysis and wrote the manuscript. GG mapped the short-reads to the rainbow trout reference genome. RL and GG contributed to SV identification and validation. GT provided the sequences of the double haploid lines. TL, KM, and GW were involved in selection of the outbred fish used for sequencing. All authors read and approved the final manuscript.

Funding

This study was supported by Agricultural Research Service CRIS project 8082-31000-013-00D.

Conflict of Interest

KM is employed by Troutlodge Inc.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to thank Roseanna Long and Kristy Shewbridge for preparing the DNA samples used in this study. SL is also grateful to Dr. Ning Jiang at the Michigan State University for helpful discussion about transposable elements. We would also like to thank Dr. Ben Koop and David Minkley at the University of Victoria for sharing the repetitive sequences of rainbow trout. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture (USDA). USDA is an equal opportunity provider and employer.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.639355/full#supplementary-material

Supplementary Table 1 | Details of the samples used in this study including sample names, populations, genetic lines, read coverage and number of the high-confidence SVs per sample.

Supplementary Table 2 | SV effects by type. The high-confidence SVs were annotated using SnpEff.

Supplementary Note | The sequences of three transposable elements highlighted in this study.

Footnotes

1. ^ http://broadinstitute.github.io/picard/
2. ^ https://github.com/brentp/smoove

References

Abel, H. J., Larson, D. E., Regier, A. A., Chiang, C., Das, I., Kanchi, K. L., et al. (2020). Mapping and characterization of structural variation in 17,795 human genomes. Nature 583, 83–89. doi: 10.1038/s41586-020-2371-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Alkan, C., Coe, B. P., and Eichler, E. E. (2011). Genome structural variation discovery and genotyping. Nat. Rev. Genet. 12, 363–376. doi: 10.1038/nrg2958

PubMed Abstract | CrossRef Full Text | Google Scholar

Alonge, M., Wang, X. G., Benoit, M., Soyk, S., Pereira, L., Zhang, L., et al. (2020). Major impacts of widespread structural variation on gene expression and crop improvement in tomato. Cell 182, 145–161. doi: 10.1016/j.cell.2020.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertolotti, A. C., Layer, R. M., Gundappa, M. K., Gallagher, M. D., Pehlivanoglu, E., Nome, T., et al. (2020). The structural variation landscape in 492 Atlantic salmon genomes. Nat. Commun. 11. doi: 10.1038/s41467-020-18972-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cameron, D. L., Di Stefano, L., and Papenfuss, A. T. (2019). Comprehensive evaluation and characterisation of short read general-purpose structural variant calling software. Nat. Commun. 10:3240.

Google Scholar

Chang, C. C., Chow, C. C., Tellier, L. C. A. M., Vattikuti, S., Purcell, S. M., and Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4:7. doi: 10.1186/S13742-015-0047-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiang, C., Layer, R. M., Faust, G. G., Lindberg, M. R., Rose, D. B., Garrison, E. P., et al. (2015). SpeedSeq: ultra-fast personal genome analysis and interpretation. Nat. Methods 12, 966–968. doi: 10.1038/Nmeth.3505

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiang, C., Scott, A. J., Davis, J. R., Tsang, E. K., Li, X., Kim, Y., et al. (2017). The impact of structural variation on human gene expression. Nat. Genet. 49, 692–699. doi: 10.1038/ng.3834

PubMed Abstract | CrossRef Full Text | Google Scholar

Chu, C., Zhao, B., Park, P. J., and Lee, E. A. (2020). Identification and genotyping of transposable element insertions from genome sequencing data. Curr. Protoc. Hum. Genet. 107:e102. doi: 10.1002/cphg.102

PubMed Abstract | CrossRef Full Text | Google Scholar

Cingolani, P., Platts, A., Wang, L. L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w(1118); iso-2; iso-3. Fly 6, 80–92. doi: 10.4161/fly.19695

PubMed Abstract | CrossRef Full Text | Google Scholar

Deane-Coe, P. E., Chu, E. T., Slavney, A., Boyko, A. R., and Sams, A. J. (2018). Direct-to-consumer DNA testing of 6,000 dogs reveals 98.6-kb duplication associated with blue eyes and heterochromia in Siberian Huskies. PLoS Genet. 14:e1007648. doi: 10.1371/journal.pgen.1007648

PubMed Abstract | CrossRef Full Text | Google Scholar

Dominguez, M., Dugas, E., Benchouaia, M., Leduque, B., Jimenez-Gomez, J. M., Colot, V., et al. (2020). The impact of transposable elements on tomato diversity. Nat. Commun. 11:4058. doi: 10.1038/S41467-020-17874-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorant, Y., Cayuela, H., Wellband, K., Laporte, M., Rougemont, Q., Mérot, C., et al. (2020). Copy number variants outperform SNPs to reveal genotype-temperature association in a marine species. Mol. Ecol. 29, 4765–4782. doi: 10.1111/mec.15565

PubMed Abstract | CrossRef Full Text | Google Scholar

Faber-Hammond, J. J., Bezault, E., Lunt, D. H., Joyce, D. A., and Renn, S. C. P. (2019). The genomic substrate for adaptive radiation: copy number variation across 12 tribes of African cichlid species. Genome Biol. Evol. 11, 2856–2874. doi: 10.1093/gbe/evz185

PubMed Abstract | CrossRef Full Text | Google Scholar

Food and Agriculture Organization (FAO) (2021). Fisheries and Aquaculture [Online]. Available online at: http://www.fao.org/fishery/culturedspecies/Oncorhynchus_mykiss/en (accessed January 20, 2021).

Google Scholar

Gao, G. T., Nome, T., Pearse, D. E., Moen, T., Naish, K. A., Thorgaard, G. H., et al. (2018). A new single nucleotide polymorphism database for rainbow trout generated through whole genome resequencing. Front. Genet. 9:147. doi: 10.3389/Fgene.2018.00147

PubMed Abstract | CrossRef Full Text | Google Scholar

Ho, S. V. S., Urban, A. E., and Mills, R. E. (2020). Structural variation in the sequencing era. Nat. Rev. Genet. 21, 171–189. doi: 10.1038/s41576-019-0180-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C. R. L., Schneider, A. M., Lu, Y. Q., Niranjan, T., Shen, P. L., Robinson, M. A., et al. (2010). Mobile interspersed repeats are major structural variants in the human genome. Cell 141, 1171–1182. doi: 10.1016/j.cell.2010.05.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Kou, Y., Liao, Y., Toivainen, T., Lv, Y., Tian, X., Emerson, J. J., et al. (2020). Evolutionary genomics of structural variation in Asian rice (Oryza sativa) domestication. Mol. Biol. Evol. 37, 3507–3524. doi: 10.1093/molbev/msaa185

CrossRef Full Text | Google Scholar

Larson, D. E., Abel, H. J., Chiang, C., Badve, A., Das, I., Eldred, J. M., et al. (2019). svtools: population-scale analysis of structural variation. Bioinformatics 35, 4782–4787. doi: 10.1093/bioinformatics/btz492

PubMed Abstract | CrossRef Full Text | Google Scholar

Layer, R. M., Chiang, C., Quinlan, A. R., and Hall, I. M. (2014). LUMPY: a probabilistic framework for structural variant discovery. Genome Biol. 15:R84. doi: 10.1186/Gb-2014-15-6-R84

PubMed Abstract | CrossRef Full Text | Google Scholar

Leeds, T. D., Silverstein, J. T., Weber, G. M., Vallejo, R. L., Palti, Y., Rexroad, C. E., et al. (2010). Response to selection for bacterial cold water disease resistance in rainbow trout. J. Anim. Sci. 88, 1936–1946. doi: 10.2527/jas.2009-2538

PubMed Abstract | CrossRef Full Text | Google Scholar

Leeds, T. D., Vallejo, R. L., Weber, G. M., Gonzalez-Pena, D., and Silverstein, J. T. (2016). Response to five generations of selection for growth performance traits in rainbow trout (Oncorhynchus mykiss). Aquaculture 465, 341–351. doi: 10.1016/j.aquaculture.2016.08.036

CrossRef Full Text | Google Scholar

Li, H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:1303.3997v2

Google Scholar

Li, H. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094–3100. doi: 10.1093/bioinformatics/bty191

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352

PubMed Abstract | CrossRef Full Text | Google Scholar

Lien, S., Koop, B. F., Sandve, S. R., Miller, J. R., Kent, M. P., Nome, T., et al. (2016). The Atlantic salmon genome provides insights into rediploidization. Nature 533, 200–205. doi: 10.1038/nature17164

PubMed Abstract | CrossRef Full Text | Google Scholar

Lisch, D. (2013). How important are transposons for plant evolution? Nat. Rev. Genet. 14, 49–61. doi: 10.1038/nrg3374

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y. C., Du, H. L., Li, P. C., Shen, Y. T., Peng, H., Liu, S. L., et al. (2020). Pan-genome of wild and cultivated soybeans. Cell 182, 162–176.e13. doi: 10.1016/j.cell.2020.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Macqueen, D. J., and Johnston, I. A. (2014). A well-constrained estimate for the timing of the salmonid whole genome duplication reveals major decoupling from species diversification. Proc. R. Soc. B Biol. Sci. 281:20132881. doi: 10.1098/Rspb.2013.2881

PubMed Abstract | CrossRef Full Text | Google Scholar

Manichaikul, A., Mychaleckyj, J. C., Rich, S. S., Daly, K., Sale, M., and Chen, W. M. (2010). Robust relationship inference in genome-wide association studies. Bioinformatics 26, 2867–2873. doi: 10.1093/bioinformatics/btq559

PubMed Abstract | CrossRef Full Text | Google Scholar

Palti, Y., Gao, G. T., Miller, M. R., Vallejo, R. L., Wheeler, P. A., Quillet, E., et al. (2014). A resource of single-nucleotide polymorphisms for rainbow trout generated by restriction-site associated DNA sequencing of doubled haploids. Mol. Ecol. Resour. 14, 588–596. doi: 10.1111/1755-0998.12204

PubMed Abstract | CrossRef Full Text | Google Scholar

Parsons, J. E., and Thorgaard, G. H. (1985). Production of androgenetic diploid rainbow trout. J. Hered. 76, 177–181. doi: 10.1093/oxfordjournals.jhered.a110060

PubMed Abstract | CrossRef Full Text | Google Scholar

Payer, L. M., and Burns, K. H. (2019). Transposable elements in human genetic disease. Nat. Rev. Genet. 20, 760–772. doi: 10.1038/s41576-019-0165-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Pearse, D. E., Barson, N. J., Nome, T., Gao, G. T., Campbell, M. A., Abadia-Cardoso, A., et al. (2019). Sex-dependent dominance maintains migration supergene in rainbow trout. Nat. Ecol. Evol. 3, 1731–1742. doi: 10.1038/s41559-019-1044-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedersen, B. S., and Quinlan, A. R. (2018). Mosdepth: quick coverage calculation for genomes and exomes. Bioinformatics 34, 867–868. doi: 10.1093/bioinformatics/btx699

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedersen, B. S., and Quinlan, A. R. (2019). Duphold: scalable, depth-based annotation and curation of high-confidence structural variant calls. Gigascience 8:giZ040. doi: 10.1093/gigascience/giz040

PubMed Abstract | CrossRef Full Text | Google Scholar

Penso-Dolfin, L., Man, A., Mehta, T., Haerty, W., and Di Palma, F. (2020). Analysis of structural variants in four African cichlids highlights an association with developmental and immune related genes. BMC Evol. Biol. 20:69. doi: 10.1186/S12862-020-01629-0

CrossRef Full Text | Google Scholar

Quinlan, A. R., and Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi: 10.1093/bioinformatics/btq033

PubMed Abstract | CrossRef Full Text | Google Scholar

Smit, A. F. A., Hubley, R., and Green, P. (2020). Repeat Maskers [Online]. Available online at: http://www.repeatmasker.org/ (accessed Octobor 6, 2020).

Google Scholar

Sotero-Caio, C. G., Platt, R. N., Suh, A., and Ray, D. A. (2017). Evolution and diversity of transposable elements in vertebrate genomes. Genome Biol. Evol. 9, 161–177. doi: 10.1093/gbe/evw264

PubMed Abstract | CrossRef Full Text | Google Scholar

Spielmann, M., Lupianez, D. G., and Mundlos, S. (2018). Structural variation in the 3D genome. Nat. Rev. Genet. 19, 453–467. doi: 10.1038/s41576-018-0007-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Spielmann, M., and Mundlos, S. (2013). Structural variations, the regulatory landscape of the genome and their alteration in human disease. Bioessays 35, 533–543. doi: 10.1002/bies.201200178

PubMed Abstract | CrossRef Full Text | Google Scholar

Subramanian, R. P., Wildschutte, J. H., Russo, C., and Coffin, J. M. (2011). Identification, characterization, and comparative genomic distribution of the HERV-K (HML-2) group of human endogenous retroviruses. Retrovirology 8:90. doi: 10.1186/1742-4690-8-90

PubMed Abstract | CrossRef Full Text | Google Scholar

Sudmant, P. H., Rausch, T., Gardner, E. J., Handsaker, R. E., Abyzov, A., Huddleston, J., et al. (2015). An integrated map of structural variation in 2,504 human genomes. Nature 526, 75–81. doi: 10.1038/nature15394

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, J., Perron, H., and Feschotte, C. (2018). Variation in proviral content among human genomes mediated by LTR recombination. Mob. DNA 9:36. doi: 10.1186/s13100-018-0142-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Thorgaard, G. H. (1983). Chromosomal differences among rainbow trout populations. Copeia 1983, 650–662. doi: 10.2307/1444329

CrossRef Full Text | Google Scholar

Wiens, G. D., Palti, Y., and Leeds, T. D. (2018). Three generations of selective breeding improved rainbow trout (Oncorhynchus mykiss) disease resistance against natural challenge with Flavobacterium psychrophilum during early life-stage rearing. Aquaculture 497, 414–421. doi: 10.1016/j.aquaculture.2018.07.064

CrossRef Full Text | Google Scholar

Wildschutte, J. H., Williams, Z. H., Montesion, M., Subramanian, R. P., Kidd, J. M., and Coffin, J. M. (2016). Discovery of unfixed endogenous retrovirus insertions in diverse human populations. Proc. Natl. Acad. Sci. U.S.A. 113, E2326–E2334. doi: 10.1073/pnas.1602336113

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, P. J., Li, J. H., Kang, H. M., Wang, H. F., Fan, Z. Y., Yin, Z. J., et al. (2016). Structural variant detection by large-scale sequencing reveals new evolutionary evidence on breed divergence between Chinese and European pigs. Sci. Rep. 6:18501. doi: 10.1038/Sre18501

CrossRef Full Text | Google Scholar

Zhou, W. C., Emery, S. B., Flasch, D. A., Wang, Y. F., Kwan, K. Y., Kidd, J. M., et al. (2020). Identification and characterization of occult human-specific LINE-1 insertions using long-read sequencing technology. Nucleic Acids Res. 48, 1146–1163. doi: 10.1093/nar/gkz1173

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y. F., Minio, A., Massonnet, M., Solares, E., Lv, Y. D., Beridze, T., et al. (2019). The population genetics of structural variants in grapevine domestication. Nat. Plants 5, 965–979. doi: 10.1038/s41477-019-0507-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zook, J. M., Hansen, N. F., Olson, N. D., Chapman, L., Mullikin, J. C., Xiao, C. N., et al. (2020). A robust benchmark for detection of germline large deletions and insertions. Nat. Biotechnol. 38, 1347–1355. doi: 10.1038/s41587-020-0538-8

PubMed Abstract | CrossRef Full Text | Google Scholar