Characterization and Functional Analysis of Fads Reveals Δ5 Desaturation Activity during Long-Chain Polyunsaturated Fatty Acid Biosynthesis in Dwarf Surf Clam Mulinia lateralis.

Fatty acid desaturases (Fads), as key enzymes in the biosynthesis of long-chain polyunsaturated fatty acids (LC-PUFAs), catalyze the desaturation between defined carbons of fatty acyl chains and control the degree of unsaturation of fatty acids. In the present study, two Fads genes, designated MulFadsA and MulFadsB, were identified from the genome of the dwarf surf clam Mulinia lateralis (Mollusca, Mactridae), and their spatiotemporal expression was examined. MulFadsA and MulFadsB contained the corresponding conserved functional domains and clustered closely with their respective orthologs from other mollusks. Both genes were expressed in the developmental stages and all tested adult tissues of M. lateralis, with MulFadsA exhibiting significantly higher expression levels in adult tissues than MulFadsB. Subsequently, the effects of dietary microalgae on Fads expressions in the dwarf surf clam were investigated by feeding clams with two types of unialgal diets varying in fatty acid content, i.e., Chlorella pyrenoidosa (Cp) and Platymonas helgolandica (Ph). The results show that the expressions of MulFads were significantly upregulated among adult tissues in the Cp group compared with those in the Ph group. In addition, we observed the desaturation activity of MulFadsA via heterologous expression in yeasts, revealing Δ5 desaturation activity toward PUFA substrates. Taken together, these results provide a novel perspective on M. lateralis LC-PUFA biosynthesis, expanding our understanding of fatty acid synthesis in marine mollusks.

Sequencing-Based Transcriptome-Wide Targeted Genotyping for Evolutionary and Ecological Studies.

Transcriptome-wide targeted genotyping is highly attractive for evolutionary and ecological studies but, until recently, accomplishing this goal presented a major technical barrier for the study of non-model organisms. Our group has recently developed a high-throughput targeted genotyping approach (called HD-Marker) based on the high specificity and accuracy of oligo extension-ligation assays that facilitates the design of assays tailored to meet specific genotyping needs. HD-Marker allows for targeted genotyping of over 10 000 genes in a single tube, with strikingly high capture rate (98%-99%) and genotyping accuracy (97%-99%). With the remarkable advantages of cost-effectiveness and flexibility, we envision that HD-Marker has broad application potential in evolutionary and ecological studies.

Syngas to light olefins conversion with high olefin/paraffin ratio using ZnCrOx/AlPO-18 bifunctional catalysts.

Direct synthesis of light olefins from syngas (STO) using a bifunctional catalyst composed of oxide and zeolite has attracted extensive attention in both academia and industry. It is highly desirable to develop robust catalysts that could enhance the CO conversion while simultaneously maintain high selectivity to C2-C4 olefins. Herein, we report a bifunctional catalyst consisting of ZnCr binary oxide (ZnCrOx) and low-Si AlPO-18 zeolite, showing both satisfying selectivity to C2-C4 olefins of 45.0% (86.7%, CO2 free) and high olefin/paraffin ratio of 29.9 at the CO conversion of 25.2% under mild reaction conditions (4.0 MPa, 390 °C). By optimizing the reaction conditions, the CO conversion could be markedly increased to 49.3% with a slight drop in selectivity. CD3CN/CO-FTIR characterizations and theoretical calculations demonstrate that low-Si AlPO-18 zeolite has lower acid strength, and is therefore less reactive toward the hydride transfer in the STO reaction, leading to a higher olefin/paraffin ratio.

HD-Marker: a highly multiplexed and flexible approach for targeted genotyping of more than 10,000 genes in a single-tube assay.

Targeted genotyping of transcriptome-scale genetic markers is highly attractive for genetic, ecological, and evolutionary studies, but achieving this goal in a cost-effective manner remains a major challenge, especially for laboratories working on nonmodel organisms. Here, we develop a high-throughput, sequencing-based GoldenGate approach (called HD-Marker), which addresses the array-related issues of original GoldenGate methodology and allows for highly multiplexed and flexible targeted genotyping of more than 12,000 loci in a single-tube assay (in contrast to fewer than 3100 in the original GoldenGate assay). We perform extensive analyses to demonstrate the power and performance of HD-Marker on various multiplex levels (296, 795, 1293, and 12,472 genic SNPs) across two sequencing platforms in two nonmodel species (the scallops Chlamys farreri and Patinopecten yessoensis), with extremely high capture rate (98%-99%) and genotyping accuracy (97%-99%). We also demonstrate the potential of HD-Marker for high-throughput targeted genotyping of alternative marker types (e.g., microsatellites and indels). With its remarkable cost-effectiveness (as low as $0.002 per genotype) and high flexibility in choice of multiplex levels and marker types, HD-Marker provides a highly attractive tool over array-based platforms for fulfilling genome/transcriptome-wide targeted genotyping applications, especially in nonmodel organisms.

Sea cucumber genome provides insights into saponin biosynthesis and aestivation regulation.

Echinoderms exhibit several fascinating evolutionary innovations that are rarely seen in the animal kingdom, but how these animals attained such features is not well understood. Here we report the sequencing and analysis of the genome and extensive transcriptomes of the sea cucumber Apostichopus japonicus, a species from a special echinoderm group with extraordinary potential for saponin synthesis, aestivation and organ regeneration. The sea cucumber does not possess a reorganized Hox cluster as previously assumed for all echinoderms, and the spatial expression of Hox7 and Hox11/13b potentially guides the embryo-to-larva axial transformation. Contrary to the typical production of lanosterol in animal cholesterol synthesis, the oxidosqualene cyclase of sea cucumber produces parkeol for saponin synthesis and has "plant-like" motifs suggestive of convergent evolution. The transcriptional factors Klf2 and Egr1 are identified as key regulators of aestivation, probably exerting their effects through a clock gene-controlled process. Intestinal hypometabolism during aestivation is driven by the DNA hypermethylation of various metabolic gene pathways, whereas the transcriptional network of intestine regeneration involves diverse signaling pathways, including Wnt, Hippo and FGF. Decoding the sea cucumber genome provides a new avenue for an in-depth understanding of the extraordinary features of sea cucumbers and other echinoderms.

Scallop genome reveals molecular adaptations to semi-sessile life and neurotoxins.

Bivalve molluscs are descendants of an early-Cambrian lineage superbly adapted to benthic filter feeding. Adaptations in form and behavior are well recognized, but the underlying molecular mechanisms are largely unknown. Here, we investigate the genome, various transcriptomes, and proteomes of the scallop Chlamys farreri, a semi-sessile bivalve with well-developed adductor muscle, sophisticated eyes, and remarkable neurotoxin resistance. The scallop's large striated muscle is energy-dynamic but not fully differentiated from smooth muscle. Its eyes are supported by highly diverse, intronless opsins expanded by retroposition for broadened spectral sensitivity. Rapid byssal secretion is enabled by a specialized foot and multiple proteins including expanded tyrosinases. The scallop uses hepatopancreas to accumulate neurotoxins and kidney to transform to high-toxicity forms through expanded sulfotransferases, probably as deterrence against predation, while it achieves neurotoxin resistance through point mutations in sodium channels. These findings suggest that expansion and mutation of those genes may have profound effects on scallop's phenotype and adaptation.

Scallop genome provides insights into evolution of bilaterian karyotype and development.

Reconstructing the genomes of bilaterian ancestors is central to our understanding of animal evolution, where knowledge from ancient and/or slow-evolving bilaterian lineages is critical. Here we report a high-quality, chromosome-anchored reference genome for the scallop Patinopecten yessoensis, a bivalve mollusc that has a slow-evolving genome with many ancestral features. Chromosome-based macrosynteny analysis reveals a striking correspondence between the 19 scallop chromosomes and the 17 presumed ancestral bilaterian linkage groups at a level of conservation previously unseen, suggesting that the scallop may have a karyotype close to that of the bilaterian ancestor. Scallop Hox gene expression follows a new mode of subcluster temporal co-linearity that is possibly ancestral and may provide great potential in supporting diverse bilaterian body plans. Transcriptome analysis of scallop mantle eyes finds unexpected diversity in phototransduction cascades and a potentially ancient Pax2/5/8-dependent pathway for noncephalic eyes. The outstanding preservation of ancestral karyotype and developmental control makes the scallop genome a valuable resource for understanding early bilaterian evolution and biology.

Proteomic analysis of scallop hepatopancreatic extract provides insights into marine polysaccharide digestion.

Marine polysaccharides are used in a variety of applications, and the enzymes that degrade these polysaccharides are of increasing interest. The main food source of herbivorous marine mollusks is seaweed, and several polysaccharide-degrading enzymes have been extracted from mollusk digestive glands (hepatopancreases). Here, we used a comprehensive proteomic approach to examine the hepatopancreatic proteins of the Zhikong scallop (Chlamys farreri). We identified 435 proteins, the majority of which were lysosomal enzymes and carbohydrate and protein metabolism enzymes. However, several new enzymes related to polysaccharide metabolism were also identified. Phylogenetic and structural analyses of these enzymes suggest that these polysaccharide-degrading enzymes may have a variety of potential substrate specificities. Taken together, our study characterizes several novel polysaccharide-degrading enzymes in the scallop hepatopancreas and provides an enhanced view of these enzymes and a greater understanding of marine polysaccharide digestion.

Long Non-Coding RNAs (lncRNAs) of Sea Cucumber: Large-Scale Prediction, Expression Profiling, Non-Coding Network Construction, and lncRNA-microRNA-Gene Interaction Analysis of lncRNAs in Apostichopus japonicus and Holothuria glaberrima During LPS Challenge and Radial Organ Complex Regeneration.

Long non-coding RNA (lncRNA) structurally resembles mRNA but cannot be translated into protein. Although the systematic identification and characterization of lncRNAs have been increasingly reported in model species, information concerning non-model species is still lacking. Here, we report the first systematic identification and characterization of lncRNAs in two sea cucumber species: (1) Apostichopus japonicus during lipopolysaccharide (LPS) challenge and in heathy tissues and (2) Holothuria glaberrima during radial organ complex regeneration, using RNA-seq datasets and bioinformatics analysis. We identified A. japonicus and H. glaberrima lncRNAs that were differentially expressed during LPS challenge and radial organ complex regeneration, respectively. Notably, the predicted lncRNA-microRNA-gene trinities revealed that, in addition to targeting protein-coding transcripts, miRNAs might also target lncRNAs, thereby participating in a potential novel layer of regulatory interactions among non-coding RNA classes in echinoderms. Furthermore, the constructed coding-non-coding network implied the potential involvement of lncRNA-gene interactions during the regulation of several important genes (e.g., Toll-like receptor 1 [TLR1] and transglutaminase-1 [TGM1]) in response to LPS challenge and radial organ complex regeneration in sea cucumbers. Overall, this pioneer systematic identification, annotation, and characterization of lncRNAs in echinoderm pave the way for similar studies and future genetic, genomic, and evolutionary research in non-model species.

Evaluation of the 2b-RAD method for genomic selection in scallop breeding.

The recently developed 2b-restriction site-associated DNA (2b-RAD) sequencing method provides a cost-effective and flexible genotyping platform for aquaculture species lacking sufficient genomic resources. Here, we evaluated the performance of this method in the genomic selection (GS) of Yesso scallop (Patinopecten yessoensis) through simulation and real data analyses using six statistical models. Our simulation analysis revealed that the prediction accuracies obtained using the 2b-RAD markers were slightly lower than those obtained using all polymorphic loci in the genome. Furthermore, a small subset of markers obtained from a reduced tag representation (RTR) library presented comparable performance to that obtained using all markers, making RTR be an attractive approach for GS purpose. Six GS models exhibited variable performance in prediction accuracy depending on the scenarios (e.g., heritability, sample size, population structure), but Bayes-alphabet and BLUP-based models generally outperformed other models. Finally, we performed the evaluation using an empirical dataset composed of 349 Yesso scallops that were derived from five families. The prediction accuracy for this empirical dataset could reach 0.4 based on optimal GS models. In summary, the genotyping flexibility and cost-effectiveness make 2b-RAD be an ideal genotyping platform for genomic selection in aquaculture breeding programs.

The Rho GTPase Family Genes in Bivalvia Genomes: Sequence, Evolution and Expression Analysis.

BACKGROUND: Rho GTPases are important members of the Ras superfamily, which represents the largest signaling protein family in eukaryotes, and function as key molecular switches in converting and amplifying external signals into cellular responses. Although numerous analyses of Rho family genes have been reported, including their functions and evolution, a systematic analysis of this family has not been performed in Mollusca or in Bivalvia, one of the most important classes of Mollusca. RESULTS: In this study, we systematically identified and characterized a total set (Rho, Rac, Mig, Cdc42, Tc10, Rnd, RhoU, RhoBTB and Miro) of thirty Rho GTPase genes in three bivalve species, including nine in the Yesso scallop Patinopecten yessoensis, nine in the Zhikong scallop Chlamys farreri, and twelve in the Pacific oyster Crassostrea gigas. Phylogenetic analysis and interspecies comparison indicated that bivalves might possess the most complete types of Rho genes in invertebrates. A multiple RNA-seq dataset was used to investigate the expression profiles of bivalve Rho genes, revealing that the examined scallops share more similar Rho expression patterns than the oyster, whereas more Rho mRNAs are expressed in C. farreri and C. gigas than in P. yessoensis. Additionally, Rho, Rac and Cdc42 were found to be duplicated in the oyster but not in the scallops. Among the expanded Rho genes of C. gigas, duplication pairs with high synonymous substitution rates (Ks) displayed greater differences in expression. CONCLUSION: A comprehensive analysis of bivalve Rho GTPase family genes was performed in scallop and oyster species, and Rho genes in bivalves exhibit greater conservation than those in any other invertebrate. This is the first study focusing on a genome-wide characterization of Rho GTPase genes in bivalves, and the findings will provide a valuable resource for a better understanding of Rho evolution and Rho GTPase function in Bivalvia.

Construction of a High-Density Genetic Map and Quantitative Trait Locus Mapping in the Sea Cucumber Apostichopus japonicus.

Genetic linkage maps are critical and indispensable tools in a wide range of genetic and genomic research. With the advancement of genotyping-by-sequencing (GBS) methods, the construction of a high-density and high-resolution linkage maps has become achievable in marine organisms lacking sufficient genomic resources, such as echinoderms. In this study, high-density, high-resolution genetic map was constructed for a sea cucumber species, Apostichopus japonicus, utilizing the 2b-restriction site-associated DNA (2b-RAD) method. A total of 7839 markers were anchored to the linkage map with the map coverage of 99.57%, to our knowledge, this is the highest marker density among echinoderm species. QTL mapping and association analysis consistently captured one growth-related QTL located in a 5 cM region of linkage group (LG) 5. An annotated candidate gene, retinoblastoma-binding protein 5 (RbBP5), which has been reported to be an important regulator of cell proliferation, was recognized in the QTL region. This linkage map represents a powerful tool for research involving both fine-scale QTL mapping and marker assisted selection (MAS), and will facilitate chromosome assignment and improve the whole-genome assembly of sea cucumber in the future.

Integration of Transcriptomic and Proteomic Approaches Provides a Core Set of Genes for Understanding of Scallop Attachment.

Attachment is an essential physiological process in life histories of many marine organisms. Using a combination of transcriptomic and proteomic approach, scallop byssal proteins (Sbps) and their associated regulatory network genes were investigated for the first time. We built the first scallop foot transcriptome library, and 75 foot-specific genes were identified. Through integration of transcriptomic-proteomic approach, seven unique Sbps were identified. Of them, three showed significant amino acid sequence homology to known proteins. In contrast, the rest did not show significant protein matches, indicating they are possibly novel proteins. Our transcriptomic and proteomic analyses also suggest that post-translational modification may be one of the significant features for Sbps as well. Taken together, our study provides the first multidimensional collection of a core set of genes that may be potentially involved in scallop byssal attachment.

Mesostructured zeolites: bridging the gap between zeolites and MCM-41.

Surfactant-templating is one of the most versatile and useful techniques to implement mesoporous systems into solid materials. Various strategies based on various interactions between surfactants and solid precursors have been explored to produce new structures. Zeolites are invaluable as size- and shape-selective solid acid catalysts. Nevertheless, their micropores impose limitations on the mass transport of bulky feed and/or product molecules. Many studies have attempted to address this by utilizing surfactant-assisting technology to alleviate the diffusion constraints. However, most efforts have failed due to micro/mesopore phase separation. Recently, a new technique combining the uses of cationic surfactants and mild basic solutions was introduced to synthesise mesostructured zeolites. These materials sustain the unique characteristics of zeolites (i.e., strong acidity, crystallinity, microporosity, and hydrothermal stability), including tunable mesopore sizes and degrees of mesoporosity. The mesostructured zeolites are now commercially available through Rive Technology, and show superior performance in VGO cracking. This feature article provides an overview of recent explorations in the introduction of mesoporosity into zeolites using surfactant-templating techniques. Various porous materials, preparation methods, physical and catalytic properties of mesostructured zeolites will be discussed.

Large-scale development of gene-associated single-nucleotide polymorphism markers for molluscan population genomic, comparative genomic, and genome-wide association studies.

Mollusca is the second most diverse group of animals in the world. Despite their perceived importance, omics-level studies have seldom been applied to this group of animals largely due to a paucity of genomic resources. Here, we report the first large-scale gene-associated marker development and evaluation for a bivalve mollusc, Chlamys farreri. More than 21,000 putative single-nucleotide polymorphisms (SNPs) were identified from the C. farreri transcriptome. Primers and probes were designed and synthesized for 4500 SNPs, and 1492 polymorphic markers were successfully developed using a high-resolution melting genotyping platform. These markers are particularly suitable for population genomic analysis due to high polymorphism within and across populations, a low frequency of null alleles, and conformation to neutral expectations. Unexpectedly, high cross-species transferability was observed, suggesting that the transferable SNPs may largely represent ancestral genetic variations that have been preserved differentially among subfamilies of Pectinidae. Gene annotations were available for 73% of the markers, and 65% could be anchored to the recently released Pacific oyster genome. Large-scale association analysis revealed key candidate genes responsible for scallop growth regulation, and provided markers for further genetic improvement of C. farreri in breeding programmes.

High-resolution linkage and quantitative trait locus mapping aided by genome survey sequencing: building up an integrative genomic framework for a bivalve mollusc.

Genetic linkage maps are indispensable tools in genetic and genomic studies. Recent development of genotyping-by-sequencing (GBS) methods holds great promise for constructing high-resolution linkage maps in organisms lacking extensive genomic resources. In the present study, linkage mapping was conducted for a bivalve mollusc (Chlamys farreri) using a newly developed GBS method-2b-restriction site-associated DNA (2b-RAD). Genome survey sequencing was performed to generate a preliminary reference genome that was utilized to facilitate linkage and quantitative trait locus (QTL) mapping in C. farreri. A high-resolution linkage map was constructed with a marker density (3806) that has, to our knowledge, never been achieved in any other molluscs. The linkage map covered nearly the whole genome (99.5%) with a resolution of 0.41 cM. QTL mapping and association analysis congruously revealed two growth-related QTLs and one potential sex-determination region. An important candidate QTL gene named PROP1, which functions in the regulation of growth hormone production in vertebrates, was identified from the growth-related QTL region detected on the linkage group LG3. We demonstrate that this linkage map can serve as an important platform for improving genome assembly and unifying multiple genomic resources. Our study, therefore, exemplifies how to build up an integrative genomic framework in a non-model organism.

RADtyping: an integrated package for accurate de novo codominant and dominant RAD genotyping in mapping populations.

Genetic linkage maps are indispensable tools in genetic, genomic and breeding studies. As one of genotyping-by-sequencing methods, RAD-Seq (restriction-site associated DNA sequencing) has gained particular popularity for construction of high-density linkage maps. Current RAD analytical tools are being predominantly used for typing codominant markers. However, no genotyping algorithm has been developed for dominant markers (resulting from recognition site disruption). Given their abundance in eukaryotic genomes, utilization of dominant markers would greatly diminish the extensive sequencing effort required for large-scale marker development. In this study, we established, for the first time, a novel statistical framework for de novo dominant genotyping in mapping populations. An integrated package called RADtyping was developed by incorporating both de novo codominant and dominant genotyping algorithms. We demonstrated the superb performance of RADtyping in achieving remarkably high genotyping accuracy based on simulated and real mapping datasets. The RADtyping package is freely available at http://www2.ouc.edu.cn/mollusk/ detailen.asp?id=727.

Molecular characterization of RXR (Retinoid X Receptor) gene isoforms from the bivalve species Chlamys farreri.

BACKGROUND: Bivalves are among the oldest classes of invertebrates, and they exhibit diverse types of sexual patterning. However, our current understanding of the mechanisms of sex determination and differentiation in bivalves remains very limited. The retinoid X receptors (RXRs), which are members of the nuclear receptor family, are involved in sex differentiation in many organisms. RESULTS: In the present study, four full-length RXR-encoding cDNAs (CfRXRs) named CfRXRa, CfRXRb, CfRXRc and CfRXRd were retrieved from Zhikong scallop (Chlamys farreri). The four RXRs exhibited the conserved five-domain structure of nuclear receptor superfamily members and differed from each other only in the T-box of the C domain. The three variants, designated T (+4), T (+20) and T (+24), contained insertions of 4, 20 and 24 amino acids, respectively. The entire CfRXR gene is composed of eight exons and seven introns, and the four isoforms are generated via alternative mRNA splicing. Expression analysis showed that all four isoforms were expressed in both the testis and the ovary during the differentiation stage, whereas no expression was detected in the growth, mature or resting stages. This result suggests that CfRXRs are involved in germ cell differentiation in both sexes. The expression of the four isoforms was also detected in other tissues examined, including mantle, gill, digestive gland, and adductor muscle of sexually mature male and female Zhikong scallops, implying the multiple biological functions of CfRXRs. CONCLUSION: Our study presents the first report of RXR isoforms in bivalves. Further investigation of the functional roles of different RXR isoforms may provide deep insights into the regulatory mechanism of sex differentiation in C. farreri.

Transcriptome sequencing of Zhikong scallop (Chlamys farreri) and comparative transcriptomic analysis with Yesso scallop (Patinopecten yessoensis).

BACKGROUND: Bivalves play an important role in the ecosystems they inhabit and represent an important food source all over the world. So far limited genetic research has focused on this group of animals largely due to the lack of sufficient genetic or genomic resources. Here, we performed de novo transcriptome sequencing to produce the most comprehensive expressed sequence tag resource for Zhikong scallop (Chlamys farreri), and conducted the first transcriptome comparison for scallops. RESULTS: In a single 454 sequencing run, 1,033,636 reads were produced and then assembled into 26,165 contigs. These contigs were then clustered into 24,437 isotigs and further grouped into 20,056 isogroups. About 47% of the isogroups showed significant matches to known proteins based on sequence similarity. Transcripts putatively involved in growth, reproduction and stress/immune-response were identified through Gene ontology (GO) and KEGG pathway analyses. Transcriptome comparison with Yesso scallop (Patinopecten yessoensis) revealed similar patterns of GO representation. Moreover, 38 putative fast-evolving genes were identified through analyzing the orthologous gene pairs between the two scallop species. More than 46,000 single nucleotide polymorphisms (SNPs) and 350 simple sequence repeats (SSRs) were also detected. CONCLUSION: Our study provides the most comprehensive transcriptomic resource currently available for C. farreri. Based on this resource, we performed the first large-scale transcriptome comparison between the two scallop species, C. farreri and P. yessoensis, and identified a number of putative fast-evolving genes, which may play an important role in scallop speciation and/or local adaptation. A large set of single nucleotide polymorphisms and simple sequence repeats were identified, which are ready for downstream marker development. This transcriptomic resource should lay an important foundation for future genetic or genomic studies on C. farreri.

Reference-free SNP calling: improved accuracy by preventing incorrect calls from repetitive genomic regions.

BACKGROUND: Single nucleotide polymorphisms (SNPs) are the most abundant type of genetic variation in eukaryotic genomes and have recently become the marker of choice in a wide variety of ecological and evolutionary studies. The advent of next-generation sequencing (NGS) technologies has made it possible to efficiently genotype a large number of SNPs in the non-model organisms with no or limited genomic resources. Most NGS-based genotyping methods require a reference genome to perform accurate SNP calling. Little effort, however, has yet been devoted to developing or improving algorithms for accurate SNP calling in the absence of a reference genome. RESULTS: Here we describe an improved maximum likelihood (ML) algorithm called iML, which can achieve high genotyping accuracy for SNP calling in the non-model organisms without a reference genome. The iML algorithm incorporates the mixed Poisson/normal model to detect composite read clusters and can efficiently prevent incorrect SNP calls resulting from repetitive genomic regions. Through analysis of simulation and real sequencing datasets, we demonstrate that in comparison with ML or a threshold approach, iML can remarkably improve the accuracy of de novo SNP genotyping and is especially powerful for the reference-free genotyping in diploid genomes with high repeat contents. CONCLUSIONS: The iML algorithm can efficiently prevent incorrect SNP calls resulting from repetitive genomic regions, and thus outperforms the original ML algorithm by achieving much higher genotyping accuracy. Our algorithm is therefore very useful for accurate de novo SNP genotyping in the non-model organisms without a reference genome.