Introduction
The Decapoda represent a highly speciose and diverse group of crustaceans that provide important ecosystem services across mainly marine and freshwater environments, with many of the larger species being of significant commercial importance for fisheries and aquaculture industries. The value of generating and using genomic resources for key crustacean species and groups is widely appreciated for a variety of purposes including phylogenetic and population studies, selective breeding programs and broodstock management (Tan et al., 2016; Grandjean et al., 2017; Wolfe et al., 2019; Zenger et al., 2019; Zhang et al., 2019). However, comprehensive genomic studies are particularly limited for decapod crustaceans (Zenger et al., 2019; Zhang et al., 2019), compared to other aquatic invertebrate groups. Most genomic-related studies on decapod crustaceans have focused on mitogenome recovery and transcriptomics studies, with few genome assemblies attempted and most of these on commercially-important species (Song et al., 2016; Yuan et al., 2018; Zhang et al., 2019). Of the decapod species for which genome assemblies have been published and are publicly available, the quality is highly variable and BUSCO completeness poor with the exception of Procambarus virginalis (87.0%) and Litopeneaus vannamei (84.6%) (Table 1). It is now apparent that the assembly and annotation of decapod genomes are highly challenging due to the need to deal with large repetitive genomes, often with high heterozygosity, and should not to be undertaken by the faint-hearted or on a shoe-string budget (Zhang et al., 2019).
Table 1
| Order | Sub/Infraorder | Species | Assembly size (bp) | Number of scaffolds | Scaffold N50 length (bp) | Longest scaffold (bp) | Sequencing technology | Source | References |
|---|---|---|---|---|---|---|---|---|---|
| Class Branchiopoda | |||||||||
| Diplostraca | Anomopoda | Daphnia pulex | 197,206,209 | 5,186 | 642,089 | 4,193,030 | SS | Ensembl (V1.0) | Colbourne et al. (2011) |
| Diplostraca | N/A | Eulimnadia texana | 120,535,642 | 108 | 18,070,303 | 42,684,797 | PE, PB | NCBI (NKDA01) | Baldwin-Brown et al. (2017) |
| Class Malacostraca | |||||||||
| Amphipoda | Talitrida | Parhyale hawaiensis | 2,752,560,740 | 278,189 | 20,228,728 | 75,825,039 | PE, MP | NCBI (LQNS02) | Kao et al. (2016) |
| Decapoda | Astacidea | Cherax quadricarinatus | 3,236,648,033 | 508,682 | 33,235 | 970,867 | PE, ONT | NCBI (VSFE01) | This study |
| Decapoda | Astacidea | Procambarus virginalis | 3,338,655,684 | 1,980,964 | 37,475 | 717,999 | PE, LJD | NCBI (MRZY01) | Gutekunst et al. (2018) |
| Decapoda | Brachyura | Eriocheir sinensis | 1,118,179,535 | 1,768,652 | 111,755 | 2,002,076 | PE | GigaScience repo | Song et al. (2016) |
| Decapoda | Caridea | Exopalaemon carinicauda | 6,699,723,695 | 9,470,451 | 962 | 135,963 | PE | NCBI (QUOF01) | Li et al. (2019) |
| Decapoda | Caridea | Neocaridina denticulata | 1,284,468,468 | 3,346,358 | 400 | 124,746 | PE | requested from author | Kenny et al. (2014) |
| Decapoda | Dendrobranchiata | Litopenaeus vannamei | 1,663,565,311 | 4,682 | 605,555 | 3,458,385 | PE, PB, BAC | NCBI (QCYY01) | Zhang et al. (2019) |
| Decapoda | Dendrobranchiata | Marsupenaeus japonicus | 1,660,270,162 | 2,434,740 | 912 | 29,048 | PE | NCBI (NIUR01) | Yuan et al. (2018) |
| Decapoda | Dendrobranchiata | Penaeus monodon | 1,447,415,504 | 2,525,346 | 769 | 26,545 | PE | NCBI (NIUS01) | Yuan et al. (2018) |
Genomes used in this comparative study.
Sequencing Technology abbreviations.
PE, Paired-end short reads (Illumina).
MP, Mate-pair short reads (Illumina).
LJD, Long jumping distance (Illumina).
PB, Pacific Biosciences long reads.
ONT, Oxford Nanopore Technology long reads.
SS, Paired-end Sanger sequencing (plasmid, fosmid libraries).
The lack of genomic resources coupled with limited understanding of the molecular basis of gene expression and phenotypic variation will continue to limit advances in aquaculture-based productivity of decapods. Understanding the molecular basis of phenotypic variation and gene function is therefore important for selective breeding programs for traits such as increased growth and disease resistance (Zenger et al., 2019). Similarly, whole genome assemblies support genome-wide association studies (GWAS) to identify trait-specific loci and for genomic-based selective breeding.
Freshwater crayfishes make up an important group of Decapod crustaceans (Crandall and De Grave, 2017) naturally occurring on all continents, with the exception of Antarctica and Africa (Toon et al., 2010; Bracken-Grissom et al., 2014). They reach their greatest species richness in Australia and North America. A number of species are subject to aquaculture and recreational activities, indigenous fishers for food security and bait fisheries for anglers, and others are of conservation concern due to a range of activities threatening vulnerable freshwater environments (Piper, 2000; Richman et al., 2015; CABI, 2019). The genomic architecture and karyotype evolution of freshwater crayfish is also of interest as they have some of the largest chromosome numbers (2n = 200) recorded for invertebrate species (Tan et al., 2004). In Australia, freshwater crayfish species from the genus Cherax Erichson, known as smooth yabbies, include the largest commercial crayfish species in the world (Austin, 1987, 1996; Austin and Knott, 1996). Of these, the economically-important and best known is the red claw crayfish, Cherax quadricarinatus von Martens, distributed widely across northern Australia (Austin, 1996; FAO, 2019). This species is capable of growing to greater than 200 grams making it the second largest commercial species of Cherax (Austin, 1987) behind the marron (C. cainii Austin) (Austin and Ryan, 2002). The popularity of the red claw as an aquaculture and ornamental species and its adaptability has resulted in widespread translocation, resulting in an extensive global distribution and it is acknowledged as a major invasive species of inland aquatic ecosystems in the tropics (Ahyong and Yeo, 2007; Larson and Olden, 2013; Saoud et al., 2013). Due to its large size and ease with which it can be maintained in captivity, C. quadricarinatus is also increasingly being used as a model organism to address fundamental questions relevant to molecular biology, physiology, functional genomics and cell biology (Fernández et al., 2012; Pamuru et al., 2012; Ventura et al., 2019).
The genetics of C. quadricarinatus has been well studied using PCR-based approaches and recently, next generation sequencing has been used for mitogenomic and transcriptomic studies (Baker et al., 2008; Gan et al., 2016; Tan et al., 2016). In this paper, we present the first genome assembly for a southern hemisphere crayfish using a hybrid assembly approach utilizing long Nanopore reads and short Illumina reads, which has proven to be an efficient and effective approach for the assembly and annotation of species with large and repetitive genomes (Austin et al., 2017; Tan et al., 2018; Gan et al., 2019; Sánchez-Herrero et al., 2019), but has only been minimally explored to support decapod assemblies (Van Quyen et al., 2020).
We then benchmark our assembly against seven other published decapod genome assemblies and present a preliminary phylogenomic analysis for decapod crustaceans. Lastly, we use these data to investigate the occurrence and diversification of nuclear mitochondrial DNA (NUMTs)1 in decapod genomes, which is directly relevant to debates on the veracity of the most common approach to the DNA barcoding of life for molecular-based taxonomic identification and biodiversity assessment based on the mitochondrial cytochrome oxidase 1 (cox1) gene region (Hebert et al., 2003; Cristescu, 2014; deWaard et al., 2018).
Materials and Methods
Data Generation
Tissue samples from two male C. quadricarinatus individuals were collected from Northern Territory in Australia, DWN1 from Rapid Creek and M2R2 from Charles Darwin University Aquaculture Centre, both located in Darwin following institutionally endorsed ethical, biosecurity and safety guidelines. Genomic DNA was extracted using E.Z.N.A.® Tissue DNA Kit (Omega Bio-tek) from tail muscle tissue. One microgram of DWN1 was sent to Macrogen, Korea for PCR-free library preparation and sequencing on the Illumina HiSeq platform. Additional Illumina sequencing on the Illumina NovaSeq 6000 system was also performed at the Deakin Genomics Centre using a PCR-based library constructed with NEBNext Ultra Illumina library preparation kit. Using tissue samples from DWN1 and M2R2 individuals, for each Nanopore library, ~1 μg of gDNA as measured by Qubit Fluorometer (Invitrogen) was processed using the LSK108 library preparation kit followed by sequencing on the R9.4 MinION flowcell. Nanopore data were base called with Albacore versions compatible with kits and flow cells used (discontinued, was available on https://community.nanoporetech.com) and adapter-trimmed with the Porechop tool (Wick, 2017).
De novo Assembly and Scaffolding of the Crayfish Genome
Raw reads were pre-processed prior to assembly. The fastp tool (–poly_g_min_len 1) (Chen et al., 2018) was used to trim polyG sequences from the 3' end of reads generated on the Illumina NovaSeq platform. All short reads were trimmed with Trimmomatic (Bolger et al., 2014) to eliminate adapters (ILLUMINACLIP:2:30:10) and low quality sequences (AVGQUAL:20), and subsequently assembled with the Platanus assembler (Kajitani et al., 2014) using default parameters. The repeat content and large size of this genome necessitated scaffolding of the assembly with several data types; these included the use of short and long reads in order to span gaps of various sizes, in addition to using data from a previous C. quadricarinatus transcriptome project (Tan et al., 2016) to stitch together scaffolds containing coding genes. A first level of scaffolding was performed with Platanus, which scaffolds the assembled contigs and further closes gaps with insert size and sequence information from short paired-end libraries. Subsequently, as second level of scaffolding was accomplished with long Nanopore reads using LINKS (Warren et al., 2015), which was run for 20 iterations with multiple distance values (d) and step of sliding window (t), in addition to set values for k-mer (-k 19), minimum number of links (-l 5) and maximum link ratio between two best pairs (-a 0.3). Detailed parameters used in each iteration is available as Data Sheet 1. Resulting scaffolds were polished with pilon v1.22 (–fix all) (Walker et al., 2014) for five iterations. A third scaffolding step was performed using C. quadricarinatus RNA-seq reads previously generated (ERP004477) (Tan et al., 2016). HISAT2 was used to align RNA-seq reads to the assembly (Kim et al., 2019), information from which was used by BESST (Sahlin et al., 2014) for further scaffolding. All assemblies were evaluated for “completeness” based on the detection of single copy conserved genes (arthropoda_odb9) with BUSCO (Simão et al., 2015).
Genome Size Estimation
Jellyfish (Marçais and Kingsford, 2011) was used to count the occurrence of 19-, 21- and 25-mers in pre-processed short reads, generating k-mer frequency histograms that were uploaded to the GenomeScope webserver (max kmer coverage disabled) (Vurture et al., 2017) for estimations of haploid genome size, heterozygosity, and repetitive content.
Gene Prediction and Functional Annotation
A first iteration of the MAKER v2.31.10 annotation pipeline (Cantarel et al., 2008) was run to produce initial gene models based on transcript and protein hints (est2genome = 1; protein2genome = 1) from C. quadricarinatus transcriptome sequences (Tan et al., 2016), protein sequences from UniProt-SwissProt (Consortium, 2018) and an additional set of proteins from four other published decapod genomes: Eriocheir sinensis (Song et al., 2016), Litopenaeus vannamei (Zhang et al., 2019), Marsupenaeus japonicas, and Penaeus monodon (Yuan et al., 2018). These gene models were used to train two ab initio gene predictors, SNAP (Korf, 2004) and AUGUSTUS (Stanke et al., 2006), which were provided in a second MAKER iteration. Gene models resulting from this iteration were again used to retrain SNAP and AUGUSTUS before a third and final iteration of MAKER. Genes with Annotation Edit Distance values (AED) ≤ 0.5 were retained (Eilbeck et al., 2009). A small AED value points to a lesser degree of variance between the predicted gene and evidences/hints used in prediction. For functional annotation, DIAMOND (Buchfink et al., 2014) was used to find homology between the predicted genes and known proteins in the UniProtKB (Swiss-Prot and TrEMBL) database (Consortium, 2018). The predicted protein sequences were also scanned for motifs, signatures, and protein domains using InterProScan (Jones et al., 2014).
Comparative Analyses of Published Decapod Genomes
Publicly-available genome assemblies for seven published decapod species (Eriocheir sinensis, Exopalaemon carinicauda, Litopenaeus vannamei, Marsupenaeus japonicus, Neocaridina denticulata, Penaeus monodon, Procambarus virginalis) (Kenny et al., 2014; Song et al., 2016; Gutekunst et al., 2018; Yuan et al., 2018; Li et al., 2019; Zhang et al., 2019) and three other non-decapod crustaceans (Daphnia pulex, Eulimnadia texana, Parhyale hawaiensis) (Colbourne et al., 2011; Kao et al., 2016; Baldwin-Brown et al., 2017) were downloaded (sources in Table 1) and compared in terms of genome and assembly sizes, repeat content and completeness. “Completeness” of a genome was assessed with BUSCO (Simão et al., 2015) based on the arthropoda_odb9 databases. To obtain repeat profiles, repeat families were identified de novo using RepeatModeler (Smit and Hubley, 2019), which uses RepeatScout (Price et al., 2005) and RECON (Bao and Eddy, 2002) to scan for repeats. The output from this step contains a set of consensus sequences from each identified repeat family that was used by RepeatMasker (Smit et al., 2019) to mask repeats within the assembly.
Preliminary Detection of NUMTs
The eight assembled genomes for C. quadricarinatus, Er. sinensis, Ex. carinicauda, L. vannamei, M. japonicus, N. denticulata, Pe. monodon, and Pr. virginalis were also scanned for the presence of putative NUMTs. The blastn tool (Altschul et al., 1990) was used to align the mitochondrial cox1 gene of each species to its corresponding genome (exception: Procambarus clarkii cox1 was used since the full cox1 sequence for Pr. virginalis is unavailable on NCBI). Each output was filtered to eliminate alignments based on the following criteria:
Alignment length <100 nucleotides.
Alignment length that spans 95% of a genomic scaffold.
Alignment at the edge of a genomic scaffold.
An illustration of how these filters were applied is provided in Data Sheet 2. Sequences of these putative NUMTs were extracted based on alignment start and stop coordinates and were aligned with MAFFT (Katoh and Standley, 2013). IQ-TREE (Nguyen et al., 2014) was used to perform model testing and construct a maximum-likelihood (ML) phylogenetic tree (-m TESTMERGE).
Phylogenetic Analyses
A ML tree was also constructed based on orthologous protein sequences predicted by BUSCO. This tree consists of six decapods (C. quadricarinatus, Er. sinensis, L. vannamei, M. japonicus, Pe. monodon, Pr. virginalis) with D. pulex, Eu. texana and Pa. hawaiensis as non-decapod outgroup species, rooted with Drosophila melanogaster. Since this analysis used sequences predicted by BUSCO, N. denticulata and Ex. carinicauda with low BUSCO completeness were excluded from this analysis. Briefly, protein sequences within each orthologous group (based on metazoa_odb9 and arthropoda_odb9) were aligned with MAFFT (Katoh and Standley, 2013) and trimmed with Gblocks (Castresana, 2000). Only multiple sequence alignments containing sequences from all 10 species were retained and concatenated into a supermatrix with FASconCat-G (Kück and Longo, 2014). IQ-TREE (Nguyen et al., 2014) was used for initial model testing and construction of the ML tree (-m TESTMERGE), with support values indicated by the SH-like approximate likelihood ratio test (SH-aLRT) (Guindon et al., 2010) and ultrafast bootstrap (UFBoot) (Minh et al., 2013) (-alrt 1000 -bb 1000).
Results and Discussion
A Giant Genome for a Large Freshwater Crayfish
This study produced 191 Gbp of HiSeq data (2 × 100 bp) and 774 Gbp of NovaSeq data (2 × 150 bp) in addition to 36 Gbp of Nanopore reads (average: 3,419 bp) to assist with scaffolding. Kmer-based methods estimated a 5 Gbp size for the Cherax quadricarinatus genome (Data Sheet 3), a value that is within the 3.82 to 6.06 Gbp size range that has been reported for species in the Infraorder Astacidea (Gregory, 2020) (Figure 1A). Based on this estimate, data generated in this study yields sequencing depths of 193× and 7× of short and long reads, respectively. Sequence data is available in the Sequence Read Archive (SRA) on NCBI (BioProject: PRJNA559771).
Figure 1
Characteristics of the C. quadricarinatus Genome
Hybrid assembly of this genome resulted in a 3.24 Gbp final assembly contained in 508,682 scaffolds, with a N50 length of 33,235 bp (Table 1). Details of each assembly at each scaffolding level are available in Data Sheet 4. The BUSCO tool reports the presence of 81.3% and 12.8% of complete and fragmented arthropod BUSCOs, respectively. These values are comparable to the completeness evaluation of other recently sequenced decapod genomes that are relatively much smaller (L. vannamei, Pr. virginalis, Er. sinensis) (Table 1, Figure 1B). GenomeScope profiles display a small shoulder on k-mer distributions, indicating some level of heterozygosity within the genome and estimated an average heterozygosity level of 0.44%, which translates into approximately 1 mutation in 230 bp (Data Sheet 3). This is on par with the heterozygosity level of 0.53% reported for the marbled crayfish, Pr. virginalis (Gutekunst et al., 2018), and is higher than in prawns with both M. japonicus and Pe. monodon, both with 0.19% and 0.21% heterozygosity, respectively (Yuan et al., 2018). In addition, 33.73% of this assembly was masked for repeats, with a large proportion of interspersed repeats being long interspersed nuclear elements (LINEs), representing a repeat profile that is similar to that of Pr. virginalis, another crayfish species (Figure 1C). While the report of only 5.9% missing BUSCO genes in this C. quadricarinatus assembly suggests that majority of the gene space is present, the assembly reported in this study makes up only 64.8% of the expected 5 Gbp genome with 85.2% of short reads successfully mapped back to the assembly, indicating that there remains regions of this large genome that are missing from the assembly. This is not unusual for decapod genomes as they are rarely assembled to their full size, as seen from data points in Figure 1A, which mostly fall well below the “line of identity” between assembly vs. estimated genome size. The exception being Pr. virginalis with an assembled size close to its 3.5 Gbp genome size. This study tackles the ambitious assembly of a gigantic genome riddled with large repetitive structures and heterozygous regions that pose serious challenges to computational and bioinformatic resources, but there are clear benefits to hybrid assembly as an efficient method to add to the limited pool of genomic resources available for the order Decapoda. This assembly is available on NCBI (BioProject: PRJNA559771, WGS: VSFE00000000). It is noteworthy that this study has incorporated the largest dataset of long Nanopore reads to date for a decapod genome assembly (36 Gbp, ~7×), the other being the recent updated assembly of the black tiger prawn genome that incorporated the use of 2.5 Gbp of long Nanopore reads (~1.25×) (Van Quyen et al., 2020). The only other study to use long reads in decapod whole genome research is by Zhang et al. (2019) who utilized PacBio reads to support the assembly of the genome for the Pacific white shrimp L. vannamei.
Predicted Genes and Functional
An average of 95.4% of RNA-seq reads from five C. quadricarinatus tissue types in Tan et al. (2016) were successfully aligned to the assembly. The annotation process predicted a total of 19,494 protein-coding genes, with an average gene length of 9,768 bp containing an average of 4.4 exons per gene. This number of predicted genes is within the range of that typically reported for other decapod crustaceans (Er. sinensis: 14436, M. japonicus: 16716, Pe. monodon: 18100, Pr. virginalis: 21772, L. vannamei: 25596). Of the annotated protein-coding genes, 88% are functionally annotated based on homology to existing UniProt protein sequences or through identification of protein domains and signatures. Protein and transcript sequences, BLAST homology alignments and InterProScan results are available as Data Sheet 5.
Decapod Evolution Based on Nuclear Genes
Alignment and trimming of 97 orthologous genes resulted in a supermatrix of 16,679 amino acid characters. Maximum-likelihood analysis based on these protein sequences generated the tree shown in Figure 1D. Rooted with the fruit fly, the phylogeny recovers branchiopods and amphipods as outgroup species to the decapod species. Shrimp and prawn species (L. vannamei, M. japonicus, Pe. monodon) are clustered in a clade as the Dendrobranchiata, which in turn forms a sister relationship with a second clade consisting of the other three decapod species. In this latter clade, crayfish species are placed as sister taxa (C. quadricarinatus and Pr. virginalis), consistent with their taxonomic status as species within Astacidea, but with surprisingly weak nodal support (SH-aLRT 55.1%, UFBoot 66.0%). These relationships are consistent with findings reported in mitogenome- and nuclear-based studies (Bracken et al., 2009; Shen et al., 2013; Tan et al., 2019; Wolfe et al., 2019). Alignment and phylogenetic tree produced in this analysis are available as Data Sheet 6.
Integration of the Mitochondrial cox1 Gene into Decapod Genomes
The concept of DNA barcoding, by which a short universal DNA sequence is used to discriminate among species, is being very widely used to support taxonomic identification, detection of cryptic species and for establishing a reference database to support ecological, biodiversity and conservation-related studies (Ratnasingham and Hebert, 2007). A relatively recent and potentially powerful application of barcoding is the use of the rapidly increasing COI database as a reference for supporting environmental DNA metabarcoding studies (Deiner et al., 2017; Günther et al., 2018). However, DNA barcoding has its critics (Moritz and Cicero, 2004; Collins and Cruickshank, 2013), with a persistent concern being the impact of mitochondrial DNA copies in the nuclear genome (Bensasson et al., 2001; Hazkani-Covo et al., 2010) on the veracity of the DNA barcoding methodology, potentially making DNA barcoding unreliable for certain taxonomic groups (Sorenson and Quinn, 1998), including crustaceans (Song et al., 2008).
NUMTs have now been reported in a diversity of animal groups including crustaceans, with cytochrome oxidase 1 (cox1) NUMTs found in planktonic copepods (Bucklin et al., 1999) and snapping shrimp (Williams and Knowlton, 2001). Nguyen et al. (2002) and Munasinghe et al. (2003) reported the first NUMTs for crayfish, from the genus Cherax, and Song et al. (2008) for species of Orconectes. Our study extends the findings of Song et al. (2008) in that we identify cox1 pseudogenes in each of the decapod genome assemblies we have examined, but find the number to be widely variable among taxa. This analysis identified the most cox1 pseudogenes in C. quadricarinatus (34 insertion sites), followed by Pe. monodon (11), Pr. virginalis (8), M. japonicus (6), Ex. carinicauda (5), N. denticulata (3), Er. sinensis (2), and L. vannamei (1) (Figure 1E). NUMTs identified from these decapod species form monophyletic groups by species in the phylogeny, suggesting that the cox1 mitochondrial gene was integrated into each nuclear genome subsequent to the divergence of species in this analysis, with either multiple independent transfer events or further evolution of NUMTs within each species through duplication events. The NUMT phylogenetic tree is available in Data Sheet 6.
In a broad-based study of NUMTs in 85 eukaryotic genomes, Hazkani-Covo et al. (2010) found a correlation with genome size. While finding the highest number of NUMTs in the red claw crayfish, with the second largest decapod genome assembled, is consistent with this observation, the evidence for a similar overall trend is not so clear across all the decapod species. In this context, it should be noted that we only have a relatively small sampling of decapod genomes and the varying assembly quality and volume of sequence data among studies makes testing this hypothesis fraught at this stage. There are also several caveats associated with this analysis. The identification process is dependent on the search methods and implemented filters, which if too conservative, can limit detection and exclude true NUMTs with highly divergent and variable sequences, especially if the time of insertion from the mitogenome occurred millions of years ago (Tsuji et al., 2012). It is also unknown whether other published decapod studies have been post-processed to exclude scaffolds containing mitochondrial genes during the submission process to databases. Nevertheless, results from this analysis provide an insight into the prevalence of NUMTs in decapod genomes, encouraging further exploration into the evolution of NUMTs in future studies.
Conclusion
Following on from mitogenomic and transcriptomic studies, we present the first draft genome for the red claw crayfish (Cherax quadricarinatus) based on relatively large volumes of short and long genomic reads from Illumina and Oxford Nanopore (ONT) platforms. While the assembly is relatively fragmented, it is much better than many of those currently available for decapod species, and the quality of the annotation is equivalent to other more recently-sequenced decapod genomes. Due to the very large size and repetitive structures of the red claw genome, the assembly was highly challenging. However, we demonstrated the value of long Nanopore reads and a hybrid assembly approach for improving an assembly based on short reads alone (Data Sheet 4), which gives encouragement for tackling other crayfish and crustacean taxa with large and complex genomes. This draft genome will be an important and valuable resource to support ongoing comparative genomic, phylogenomics and molecular-based breeding studies for aquaculture, conservation and biodiversity-related studies and can be approved upon over time, with the generation of additional long read data. However, a major challenge still remains in relation to the computational resources needed to assemble large repetitive genomes from predominately short reads, even when aided with long reads (Lewin et al., 2018). Computationally, assembly, scaffolding and polishing processes to achieve the final draft genome took almost 85 processor-weeks, and annotation required another 100 processor-weeks. A worthwhile next step would be to investigate the efficacy of a long-read led crayfish genome assembly, now feasible as result of declining costs and improved accuracy of long read sequencing, and which should be less expensive and lead to greatly reduced processor time for assembly tasks.
Statements
Data availability statement
The datasets generated for this study can be found in the NCBI BioProject PRJNA559771, BioSamples SAMN12558911 and SAMN12885547 (WGS accession: VSFE00000000; SRA accessions: SRR10214688, SRR10416208, SRR10416210, SRR10445782, SRR10445783, SRR10445784, SRR10484712).
Author contributions
CA, HG, and LC conceived the project. CA collected the samples. HG and YL performed the sequencing. MT analyzed the data. MT and CA wrote the manuscript. FG and LC modified the manuscript. All authors read and approved the final version of the manuscript.
Acknowledgments
This research utilized computational resources and services provided by the National Computational Infrastructure (NCI), which was supported by the Australian Government. We thank the Charles Darwin University Aquaculture Centre and the Museum and Art Gallery of the Northern Territory for assistance in obtaining samples. We would also like to thank Robert Ruge of Deakin University for assistance and use of the SIT HPC Cluster.
Conflict of interest
The 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00201/full#supplementary-material
Footnotes
1.^The natural integration of DNA from mitochondria into the nuclear genome resulting in nuclear copies of mitochondrial pseudogenes (NUMTs).
References
1
AhyongS. T.YeoD. C. J. (2007). Feral populations of the Australian Red-Claw crayfish (Cherax quadricarinatus von Martens) in water supply catchments of Singapore. Biol. Invasions9, 943–946. 10.1007/s1s0530-007-90944-0
2
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215, 403–410. 10.1016/S0S022-2836(05)803600-2
3
AustinC. (1987). Diversity of Australian freshwater crayfish increases potential for aquaculture. Aust. Fish.46, 30–31.
4
AustinC. (1996). Systematics of the freshwater crayfish genus Cherax erichson (Decapoda: Parastacidae) in Northern and Eastern Australia: electrophoretic and morphological variation. Aust. J. Zool.44, 259–296. 10.1071/ZO9O960259
5
AustinC.KnottB. (1996). Systematics of the freshwater crayfish genus Cherax erichson (Decapoda: Parastacidae) in South-Western Australia: electrophoretic, morphological and habitat variation. Aust. J. Zool.44, 223–258. 10.1071/ZO9O960223
6
AustinC. M.RyanS. G. (2002). Allozyme evidence for a new species of freshwater crayfish of the genus Cherax Erichson (Decapoda : Parastacidae) from the south-west of Western Australia. Invertebr. Syst.16, 357–367. 10.1071/IT0T1010
7
AustinC. M.TanM. H.HarrissonK. A.LeeY. P.CroftL. J.SunnucksP.et al. (2017). De novo genome assembly and annotation of Australia's largest freshwater fish, the Murray cod (Maccullochella peelii), from illumina and nanopore sequencing read. Gigascience6, 1–6. 10.1093/gigascience/gix0x63
8
BakerN.de BruynM.MatherP. B. (2008). Patterns of molecular diversity in wild stocks of the redclaw crayfish (Cherax quadricarinatus) from northern Australia and Papua New Guinea: impacts of Plio-Pleistocene landscape evolution. Freshw. Biol.53, 1592–1605. 10.1111/j.13655-2427.2008.01996.x
9
Baldwin-BrownJ. G.WeekWeeksS. C.LongA. D. (2017). A new standard for crustacean genomes: the highly contiguous, annotated genome assembly of the clam shrimp Eulimnadia texana reveals HOX gene order and identifies the sex chromosome. Genome Biol. Evol.10, 143–156. 10.1093/gbe/evx2x80
10
BaoZ.EddyS. R. (2002). Automated de novo identification of repeat sequence families in sequenced genomes. Genome Res.12, 1269, 1269–1276. 10.1101/gr.88502
11
BensassonD.ZhangD.-X.HartlD. L.HewittG. M. (2001). Mitochondrial pseudogenes: evolution's misplaced witnesses. Trends Ecol. Evol.16, 314–321. 10.1016/S0S169-5347(01)021511-6
12
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics30, 2114–2120. 10.1093/bioinformatics/btu1u70
13
BrackenH. D.ToonA.FelderD. L.MartinJ. W.FinleyM.RasmussenJ.et al. (2009). The decapod tree of life: compiling the data and moving toward a consensus of decapod evolution. Arthropod Syst. Phylogeny67, 99–116.
14
Bracken-GrissomH. D.AhyongS. T.WilkinsonR. D.FeldmannR. M.SchweitzerC. E.BreinholtJ. W.et al. (2014). The emergence of lobsters: phylogenetic relationships, morphological evolution and divergence time comparisons of an ancient group (Decapoda: Achelata, Astacidea, Glypheidea, Polychelida). Syst. Biol.63, 457–479. 10.1093/sysbio/syu0u08
15
BuchfinkB.XieC.HusonD. H. (2014). Fast and sensitive protein alignment using DIAMOND. Nat. Methods12, 59–60. 10.1038/nmeth.3176
16
BucklinA.GuarnieriM.HillR.BentleyA.KaartvedtS. (1999). “Taxonomic and systematic assessment of planktonic copepods using mitochondrial COI sequence variation and competitive, species-specific PCR,” in Molecular Ecology of Aquatic Communities, eds J. P. Zehr and M. A. Voytek (Dordrecht: Springer), 239–254.
17
CABI (2019). Datasheet on Cherax quadricarinatus.
18
CantarelB. L.KorfI.RobbS. M. C.ParraG.RossE.MooreB.et al. (2008). MAKER: An easy-to-use annotation pipeline designed for emerging model organism genomes. Genome Res.18, 188–196. 10.1101/gr.6743907
19
CastresanaJ. (2000). Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol.17, 540–552. 10.1093/oxfordjournals.molbev.a0a26334
20
ChenS.ZhouY.ChenY.GuJ. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics34, i8i84–i8i90. 10.1093/bioinformatics/bty5y60
21
ColbourneJ. K.PfrenderM. E.GilbertD.ThomasW. K.TuckerA.OakleyT. H.et al. (2011). The ecoresponsive genome of Daphnia pulex. Science331, 555–561. 10.1126/science.1197761
22
CollinsR. A.CruickshankR. H. (2013). The seven deadly sins of DNA barcoding. Mol. Ecol. Resour.13, 969–975. 10.1111/17555-0998.12046
23
ConsortiumT. U. (2018). UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res.47, D5D06–D5D15. 10.1093/nar/gky1y049
24
CrandallK. A.De GraveS. (2017). An updated classification of the freshwater crayfishes (Decapoda: Astacidea) of the world, with a complete species list. J. Crustacean Biol.37, 615–653. 10.1093/jcbiol/rux0x70
25
CristescuM. E. (2014). From barcoding single individuals to metabarcoding biological communities: towards an integrative approach to the study of global biodiversity. Trends Ecol. Evol.29, 566–571. 10.1016/j.tree.2014.08.001
26
DeinerK.BikH. M.MächlerE.SeymourM.Lacoursière-RousselA.AltermattF.et al. (2017). Environmental DNA metabarcoding: Transforming how we survey animal and plant communities. Mol. Ecol.26, 5872–5895. 10.1111/mec.14350
27
deWaardJ. R.Levesque-BeaudinV.deWaardS. L.IvanovaN. V.McKeownJ. T. A.MiskieR.et al. (2018). Expedited assessment of terrestrial arthropod diversity by coupling Malaise traps with DNA barcoding. Genome62, 85–95. 10.1139/gen-20188-0093
28
EilbeckK.MooreB.HoltC.YandellM. (2009). Quantitative measures for the management and comparison of annotated genomes. BMC Bioinformatics10:67. 10.1186/14711-21055-10-67
29
FAO (2019). Cherax quadricarinatus (von Martens, 1868).
30
FernándezM. S.BustosC.LuquetG.SaezD.Neira-CarrilloA.CorneillatM.et al. (2012). Proteoglycan occurrence in gastrolith of the crayfish Cherax quadricarinatus (Malacostraca: Decapoda). J. Crustacean Biol.32, 802–815. 10.1163/193724012X22X6X49804
31
GanH. M.TanM. H.AustinC. M. (2016). The complete mitogenome of the red claw crayfish Cherax quadricarinatus (Von Martens, 1868) (Crustacea: Decapoda: Parastacidae). Mitochondrial DNA Part A27, 385–386. 10.3109/19401736.2014.895997
32
GanH. M.TanM. H.AustinC. M.ShermanC. D. H.WongY. T.StrugnellJ.et al. (2019). Best foot forward: Nanopore long reads, hybrid meta-Assembly, and haplotig purging optimizes the first genome assembly for the Southern Hemisphere blacklip abalone (Haliotis rubra). Front. Genet. 10:889. 10.3389/fgene.2019.00889
33
GrandjeanF.TanM. H.GanH. M.LeeY. P.KawaiT.DistefanoR. J.et al. (2017). Rapid recovery of nuclear and mitochondrial genes by genome skimming from Northern Hemisphere freshwater crayfish. Zool. Scripta46, 718–728. 10.1111/zsc.12247
34
GregoryT. R. (2020). “Animal Genome Size Database”. Available online at: http://www.genomesize.com
35
GuindonS.DufayardJ.-F.LefortV.AnisimovaM.HordijkW.GascuelO. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol.59, 307–321. 10.1093/sysbio/syq0q10
36
GüntherB.KnebelsbergerT.NeumannH.LaakmannS.Martínez ArbizuP. (2018). Metabarcoding of marine environmental DNA based on mitochondrial and nuclear genes. Sci. Rep.8:1482. 10.1038/s4s1598-018-32917-x
37
GutekunstJ.AndriantsoaR.FalckenhaynC.HannaK.SteinW.RasamyJ.et al. (2018). Clonal genome evolution and rapid invasive spread of the marbled crayfish. Nat. Ecol. Evol.2, 567–573. 10.1038/s4s1559-018-04677-9
38
Hazkani-CovoE.ZellerR. M.MartinW. (2010). Molecular poltergeists: mitochondrial DNA copies (numts) in sequenced nuclear genomes. PLoS Genet.6:e1e000834. 10.1371/journal.pgen.1000834
39
HebertP. D. N.CywinskaA.BallS. L.deWaardJ. R. (2003). Biological identifications through DNA barcodes. Proc. Biol. Sci.270, 313–321. 10.1098/rspb.2002.2218
40
JonesP.BinnsD.ChangH.-Y.FraserM.LiW.McAnullaC.et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics30, 1236–1240. 10.1093/bioinformatics/btu0u31
41
KajitaniR.ToshimotoK.NoguchiH.ToyodaA.OguraY.OkunoM.et al. (2014). Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res.24, 1384–1395. 10.1101/gr.170720.113
42
KaoD.LaiA. G.StamatakiE.RosicS.KonstantinidesN.JarvisE.et al. (2016). The genome of the crustacean Parhyale hawaiensis, a model for animal development, regeneration, immunity and lignocellulose digestion. Elife5:e2e0062. 10.7554/eLife.20062.057
43
KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol.30, 772–780. 10.1093/molbev/mst0t10
44
KennyN.SinY.ShenX.ZheQ.WangW.ChanT.et al. (2014). Genomic sequence and experimental tractability of a new decapod shrimp model, Neocaridina denticulata. Mar. Drugs12, 1419–1437. 10.3390/md1d2031419
45
KimD.PaggiJ. M.ParkC.BennettC.SalzbergS. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol, 37, 907–915.
46
KorfI. (2004). Gene finding in novel genomes. BMC Bioinformatics5:59. 10.1186/14711-21055-5-59
47
KückP.LongoG. C. (2014). FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studies. Front. Zool.11:81. 10.1186/s1s2983-014-0081-x
48
LarsonE. R.OldenJ. D. (2013). Crayfish occupancy and abundance in lakes of the Pacific Northwest, USA. Freshw. Sci.32, 94–107. 10.1899/12-051.1
49
LewinH. A.RobinsonG. E.KressW. J.BakerW. J.CoddingtonJ.CrandallK. A.et al. (2018). Earth BioGenome Project: sequencing life for the future of life. PNAS115, 4325–4333. 10.1073/pnas.1720115115
50
LiJ.LvJ.LiuP.ChenP.WangJ.LiJ. (2019). Genome survey and high-resolution backcross genetic linkage map construction of the ridgetail white prawn Exopalaemon carinicauda applications to QTL mapping of growth traits. BMC Genomics20:598. 10.1186/s1s2864-019-5981-x
51
MarçaisG.KingsfordC. (2011). A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics27, 764–770. 10.1093/bioinformatics/btr0r11
52
MinhB. Q.NguyenM. A. T.von HaeselerA. (2013). Ultrafast approximation for phylogenetic bootstrap. Mol. Biol. Evol.30, 1188–1195. 10.1093/molbev/mst0t24
53
MoritzC.CiceroC. (2004). DNA Barcoding: Promise and Pitfalls. PLoS Biol.2:e3e54. 10.1371/journal.pbio.0020354
54
MunasingheD. H. N.MurphyN. P.AustinC. M. (2003). Utility of mitochondrial DNA sequences from four gene regions for systematic studies of Australian freshwater crayfish of the genus Cherax (Decapoda: Parastacidae). J. Crustacean Biol.23, 402–417. 10.1163/200219755-99990350
55
NguyenL.-T.SchmidtH. A.von HaeselerA.MinhB. Q. (2014). IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol.32, 268–274. 10.1093/molbev/msu3u00
56
NguyenT. T. T.MurphyN. P.AustinC. M. (2002). Amplification of multiple copies of mitochondrial cytochrome b gene fragments in the Australian freshwater crayfish, Cherax destructor Clark (Parastacidae: Decapoda). Anim. Genet.33, 304–308. 10.1046/j.13655-2052.2002.00867.x
57
PamuruR. R.RosenO.ManorR.ChungJ. S.ZmoraN.GlazerL.et al. (2012). Stimulation of molt by RNA interference of the molt-inhibiting hormone in the crayfish Cherax quadricarinatus. Gen. Comp. Endocrinol.178, 227–236. 10.1016/j.ygcen.2012.05.007
58
PiperL. (2000). Potential for Expansion of the Freshwater Crayfish Industry in Australia.Canberra: RIRDC, 19.
59
PriceA. L.JonesN. C.PevznerP. A. (2005). De novo identification of repeat families in large genomes. Bioinformatics21(Suppl. 1), i3i51–i3i58. 10.1093/bioinformatics/bti1i018
60
RatnasinghamS.HebertP. D. N. (2007). bold: the barcode of life data system (http://www.barcodinglife.org). Mol. Ecol. Notes 7, 355–364. 10.1111/j.14711-8286.2007.01678.x
61
RichmanN. I.BöhmM.AdamsS. B.AlvarezF.BergeyE. A.BunnJ. J. S.et al. (2015). Multiple drivers of decline in the global status of freshwater crayfish (Decapoda: Astacidea). Philos. T. R. Soc. B370:20140060. 10.1098/rstb.2014.0060
62
SahlinK.VezziF.NystedtB.LundebergJ.ArvestadL. (2014). BESST - efficient scaffolding of large fragmented assemblies. BMC Bioinformatics15:281. 10.1186/14711-21055-15-281
63
Sánchez-HerreroJ. F.Frías-LópezC.EscuerP.Hinojosa-AlvarezS.ArnedoM. A.Sánchez-GraciaA.et al. (2019). The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): a valuable resource for functional and evolutionary genomic studies in chelicerates. Gigascience8:giz0z99. 10.1093/gigascience/giz0z99
64
SaoudI. P.GhanawiJ.ThompsonK. R.WebsterC. D. (2013). A review of the culture and diseases of redclaw crayfish Cherax quadricarinatus (Von Martens 1868). J. World Aquac. Soc.44, 1–29. 10.1111/jwas.12011
65
ShenH.BrabandA.ScholtzG. (2013). Mitogenomic analysis of decapod crustacean phylogeny corroborates traditional views on their relationships. Mol. Phylogenet. Evol.66, 776–789. 10.1016/j.ympev.2012.11.002
66
SimãoF. A.WaterhouseR. M.IoannidisP.KriventsevaE. V.ZdobnovE. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics31, 3210–3212. 10.1093/bioinformatics/btv3v51
67
SmitA.HubleyR. (2019). RepeatModeler-1.0.11. Institute for Systems Biology.
68
SmitA.HubleyR.GreenP. (2019). RepeatMasker Open-4.0. 2013–2015. RepeatMasker Home Page.
69
SongH.BuhayJ. E.WhitingM. F.CrandallK. A. (2008). Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. PNAS105, 13486–13491. 10.1073/pnas.0803076105
70
SongL.BianC.LuoY.WangL.YouX.LiJ.et al. (2016). Draft genome of the Chinese mitten crab, Eriocheir sinensis. Gigascience5:5. 10.1186/s1s3742-016-0112-y
71
SorensonM. D.QuinnT. W. (1998). Numts: a challenge for avian systematics and population biology. Auk115, 214–221. 10.2307/4089130
72
StankeM.KellerO.GunduzI.HayesA.WaackS.MorgensternB. (2006). AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res.34(Suppl. 2), W4W35–W4W39. 10.1093/nar/gkl2l00
73
TanM. H.AustinC. M.HammerM. P.LeeY. P.CroftL. J.GanH. M. (2018). Finding Nemo: hybrid assembly with Oxford Nanopore and Illumina reads greatly improves the clownfish (Amphiprion ocellaris) genome assembly. Gigascience7, 1–6. 10.1093/gigascience/gix1x37
74
TanM. H.GanH. M.GanH. Y.LeeY. P.CroftL. J.SchultzM. B.et al. (2016). First comprehensive multi-tissue transcriptome of Cherax quadricarinatus (Decapoda: Parastacidae) reveals unexpected diversity of endogenous cellulase. Org. Divers. Evol.16, 185–200. 10.1007/s1s3127-015-02377-3
75
TanM. H.GanH. M.LeeY. P.Bracken-GrissomH.ChanT.-Y.MillerA. D.et al. (2019). Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition. Sci. Rep.9:1075. 10.1038/s4s1598-019-471455-0
76
TanX.QinJ. G.ChenB.ChenL.LiX. (2004). Karyological analyses on redclaw crayfish Cherax quadricarinatus (Decapoda: Parastacidae). Aquaculture234, 65–76. 10.1016/j.aquaculture.2003.12.020
77
ToonA.Pérez-LosadaM.SchweitzerC. E.FeldmannR. M.CarlsonM.CrandallK. A. (2010). Gondwanan radiation of the Southern Hemisphere crayfishes (Decapoda: Parastacidae): evidence from fossils and molecules. J. Biogeogr.37, 2275–2290. 10.1111/j.13655-2699.2010.02374.x
78
TsujiJ.FrithM. C.TomiiK.HortonP. (2012). Mammalian NUMT insertion is non-random. Nucleic Acids Res.40, 9073–9088. 10.1093/nar/gks4s24
79
Van QuyenD.GanH. M.LeeY. P.NguyenD. D.NguyenT. H.TranX. T.et al. (2020). Improved genomic resources for the black tiger prawn (Penaeus monodon). Mar. Genomics100751. 10.1016/j.margen.2020.100751
80
VenturaT.StewartM. J.ChandlerJ. C.RotgansB.ElizurA.HewittA. W. (2019). Molecular aspects of eye development and regeneration in the Australian redclaw crayfish, Cherax quadricarinatus. Aquac. Res.4, 27–36. 10.1016/j.aaf.2018.04.001
81
VurtureG. W.SedlazeckF. J.NattestadM.UnderwoodC. J.FangH.GurtowskiJ.et al. (2017). GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics33, 2202–2204. 10.1093/bioinformatics/btx1x53
82
WalkerB. J.AbeelT.SheaT.PriestM.AbouellielA.SakthikumarS.et al. (2014). Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE9:e1e12963. 10.1371/journal.pone.0112963
83
WarrenR. L.YangC.VandervalkB. P.BehsazB.LagmanA.JonesS. J. M.et al. (2015). LINKS: scalable, alignment-free scaffolding of draft genomes with long reads. Gigascience4:35. 10.1186/s1s3742-015-00766-3
84
WickR. (2017). “Porechop”. Github. Available online at: https://github.com/rrwick/Porechop
85
WilliamsS. T.KnowltonN. (2001). Mitochondrial pseudogenes are pervasive and often insidious in the snapping shrimp genus Alpheus. Mol. Biol. Evol.18, 1484–1493. 10.1093/oxfordjournals.molbev.a0a03934
86
WolfeJ. M.BreinholtJ. W.CrandallK. A.LemmonA. R.LemmonE. M.TimmL. E.et al. (2019). A phylogenomic framework, evolutionary timeline and genomic resources for comparative studies of decapod crustaceans. Proc. Biol. Sci.286:2019.0079. 10.1098/rspb.2019.0079
87
YuanJ.ZhangX.LiuC.YuY.WeiJ.LiF.et al. (2018). Genomic resources and comparative analyses of two economical penaeid shrimp species, Marsupenaeus japonicus and Penaeus monodon. Mar. Genomics39, 22–25. 10.1016/j.margen.2017.12.006
88
ZengerK. R.KhatkarM. S.JonesD. B.KhalilisamaniN.JerryD. R.RaadsmaH. W. (2019). Genomic selection in aquaculture: application, limitations and opportunities with special reference to marine shrimp and pearl oysters. Front. Genet.9:693. 10.3389/fgene.2018.00693
89
ZhangX.YuanJ.SunY.LiS.GaoY.YuY.et al. (2019). Penaeid shrimp genome provides insights into benthic adaptation and frequent molting. Nat. Commun.10:356. 10.1038/s4s1467-018-081977-4
Summary
Keywords
freshwater crayfish, Parastacidae, genome, hybrid assembly, aquaculture, Illumina, Oxford Nanopore
Citation
Tan MH, Gan HM, Lee YP, Grandjean F, Croft LJ and Austin CM (2020) A Giant Genome for a Giant Crayfish (Cherax quadricarinatus) With Insights Into cox1 Pseudogenes in Decapod Genomes. Front. Genet. 11:201. doi: 10.3389/fgene.2020.00201
Received
03 December 2019
Accepted
20 February 2020
Published
06 March 2020
Volume
11 - 2020
Edited by
Peng Xu, Xiamen University, China
Reviewed by
Xiaozhu Wang, Auburn University, United States; Xiaojun Zhang, Institute of Oceanology (CAS), China; Lisui Bao, University of Chicago, United States
Updates
Copyright
© 2020 Tan, Gan, Lee, Grandjean, Croft and Austin.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mun Hua Tan mun.tan@deakin.edu.au
This article was submitted to Livestock Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.