Developmental and validation of a novel small and high-efficient panel of microhaplotypes for forensic genetics by the next generation sequencing.

BACKGROUND: In the domain of forensic science, the application of kinship identification and mixture deconvolution techniques are of critical importance, providing robust scientific evidence for the resolution of complex cases. Microhaplotypes, as the emerging class of genetic markers, have been widely studied in forensics due to their high polymorphisms and excellent stability. RESULTS AND DISCUSSION: In this research, a novel and high-efficient panel integrating 33 microhaplotype loci along with a sex-determining locus was developed by the next generation sequencing technology. In addition, we also assessed its forensic utility and delved into its capacity for kinship analysis and mixture deconvolution. The average effective number of alleles (Ae) of the 33 microhaplotype loci in the Guizhou Han population was 6.06, and the Ae values of 30 loci were greater than 5. The cumulative power of discrimination and cumulative power of exclusion values of the novel panel in the Guizhou Han population were 1-5.6 × 10- 43 and 1-1.6 × 10- 15, respectively. In the simulated kinship analysis, the panel could effectively distinguish between parent-child, full-sibling, half-sibling, grandfather-grandson, aunt-nephew and unrelated individuals, but uncertainty rates clearly increased when distinguishing between first cousins and unrelated individuals. For the mixtures, the novel panel had demonstrated excellent performance in estimating the number of contributors of mixtures with 1 to 5 contributors in combination with the machine learning methods. CONCLUSIONS: In summary, we have developed a small and high-efficient panel for forensic genetics, which could provide novel insights into forensic complex kinships testing and mixture deconvolution.

Improving DNA mixtures analysis using compound markers composed of InDels and SNPs screened from the whole genome with next-generation sequencing.

Next-generation sequencing (NGS) allows for better identification of insertion and deletion polymorphisms (InDels) and their combination with adjacent single nucleotide polymorphisms (SNPs) to form compound markers. These markers can improve the polymorphism of microhaplotypes (MHs) within the same length range, and thus, boost the efficiency of DNA mixture analysis. In this study, we screened InDels and SNPs across the whole genome and selected highly polymorphic markers composed of InDels and/or SNPs within 300 bp. Further, we successfully developed and evaluated an NGS-based panel comprising 55 loci, of which 24 were composed of both SNPs and InDels. Analysis of 124 unrelated Southern Han Chinese revealed an average effective number of alleles (Ae ) of 7.52 for this panel. The cumulative power of discrimination and cumulative probability of exclusion values of the 55 loci were 1-2.37 × 10-73 and 1-1.19 × 10-28 , respectively. Additionally, this panel exhibited high allele detection rates of over 97% in each of the 21 artificial mixtures involving from two to six contributors at different mixing ratios. We used EuroForMix to calculate the likelihood ratio (LR) and evaluate the evidence strength provided by this panel, and it could assess evidence strength with LR, distinguishing real and noncontributors. In conclusion, our panel holds great potential for detecting and analyzing DNA mixtures in forensic applications, with the capability to enhance routine mixture analysis.

Solution to a case involving the interpretation of trace degraded DNA mixtures.

DNA mixture analysis poses a significant challenge in forensic genetics, particularly when dealing with degraded and trace amount DNA samples. Multi-SNPs (MNPs) are genetic markers similar to microhaplotypes but with smaller molecular sizes (< 75 bp), making them theoretically more suitable for analyzing degraded and trace amount samples. In this case report, we investigated a cold case involving a campstool stored for over a decade, aiming to detect and locate the suspect's DNA. We employed both conventional capillary electrophoresis-based short tandem repeat (CE-STR) analysis and next-generation sequencing-based multi-SNP (NGS-MNP) analysis. The typing results and deconvolution of the mixed CE-STR profiles were inconclusive regarding the presence of the suspect's DNA in the mixed samples. However, through NGS-MNP analysis and presence probability calculations, we determined that the suspect's DNA was present in the samples from Sect. 4-1 with a probability of 1-8.41 × 10- 6 (99.999159%). This evidence contradicted the suspect's statement and aided in resolving the case. Our findings demonstrate the significant potential of MNP analysis for examining degraded and trace amount DNA mixtures in forensic investigations.

Integrative lncRNA, circRNA, and mRNA analysis reveals expression profiles of six forensic body fluids/tissue.

RNAs have attracted much attention in forensic body fluid/tissue identification (BFID) due to their tissue-specific expression characteristics. Among RNAs, long RNAs (e.g., mRNA) have a higher probability of containing more polymorphic sites that can be used to assign the specific donor of the body fluid/tissue. However, few studies have characterized their overall profiles in forensic science. In this study, we sequenced the transcriptomes of 30 samples from venous blood, menstrual blood, semen, saliva, vaginal secretion, and skin tissue, obtaining a comprehensive picture of mRNA, lncRNA, and circRNA profiles. A total of 90,305 mRNAs, 102,906 lncRNAs (including 19,549 novel lncRNAs), and 40,204 circRNAs were detected. RNA type distribution, length distribution, and expression distribution were presented according to their annotation and expression level, and many novel body fluid/tissue-specific RNA markers were identified. Furthermore, the cognate relations among the three RNAs were analyzed according to gene annotations. Finally, SNPs and InDels from RNA transcripts were genotyped, and 21,611 multi-SNP and 4,471 multi-InDel transcriptomic microhaplotypes (tMHs) were identified. These results provide a comprehensive understanding of transcriptome profiles, which could provide new avenues for tracing the origin of the body fluid/tissue and identifying an individual.

Establishment and Application of a 42-plex Microhaplotype Assay in Forensic Genetics.

OBJECTIVES: To establish and forensically verify a 42 microhaplotypes (mircohaps, MHs) multiplex assay system based on next-generation sequencing (NGS), and to explore the application value of this system in the practice of forensic genetics. METHODS: A total of 42 highly polymorphic MHs were selected from previous studies, and sequenced by the MiSeq FGxTM platform to verify the repeata-bility, sensitivity, specificity, stability, and mixture analysis ability of the detection system. Through population genetic investigation of 102 unrelated Chinese Han individuals in Liyang City, Jiangsu Province, China, the application value of this system in forensic genetics was evaluated. RESULTS: The sequencing repeatability of the 42-plex MHs assay was 100% and the sensitivity was as low as 0.062 5 ng. The system had the ability to withstand the interference of indigo (≤2 500 ng/µL), humic acid (≤9 ng/µL), hemoglobin(≤20 µmol), and urea (≤200 ng/µL) and to detect mixtures of 2 people (1∶19), 3 people (1∶1∶9) and 4 people (1∶1∶1∶9). Based on 102 individual data, the combined power of discrimination and the combined power of exclusion were 1-3.45×10-30 and 1-3.77×10-11, respectively, and the average effect value of alleles was 2.899. CONCLUSIONS: The 42-plex MHs assay was successfully established in this study and this system has high repeatability and sensitivity, good anti-jamming ability and mixture analysis ability. The 42 MHs are highly polymorphism and have good application value in individual identification and paternity testing.

Estimating microhaplotype allele frequencies from low-coverage or pooled sequencing data.

BACKGROUND: Microhaplotypes have the potential to be more cost-effective than SNPs for applications that require genetic panels of highly variable loci. However, development of microhaplotype panels is hindered by a lack of methods for estimating microhaplotype allele frequency from low-coverage whole genome sequencing or pooled sequencing (pool-seq) data. RESULTS: We developed new methods for estimating microhaplotype allele frequency from low-coverage whole genome sequence and pool-seq data. We validated these methods using datasets from three non-model organisms. These methods allowed estimation of allele frequency and expected heterozygosity at depths routinely achieved from pooled sequencing. CONCLUSIONS: These new methods will allow microhaplotype panels to be designed using low-coverage WGS and pool-seq data to discover and evaluate candidate loci. The python script implementing the two methods and documentation are available at https://www.github.com/delomast/mhFromLowDepSeq .

Sensitive, Highly Multiplexed Sequencing of Microhaplotypes From the Plasmodium falciparum Heterozygome.

BACKGROUND: Targeted next-generation sequencing offers the potential for consistent, deep coverage of information-rich genomic regions to characterize polyclonal Plasmodium falciparum infections. However, methods to identify and sequence these genomic regions are currently limited. METHODS: A bioinformatic pipeline and multiplex methods were developed to identify and simultaneously sequence 100 targets and applied to dried blood spot (DBS) controls and field isolates from Mozambique. For comparison, whole-genome sequencing data were generated for the same controls. RESULTS: Using publicly available genomes, 4465 high-diversity genomic regions suited for targeted sequencing were identified, representing the P. falciparum heterozygome. For this study, 93 microhaplotypes with high diversity (median expected heterozygosity = 0.7) were selected along with 7 drug resistance loci. The sequencing method achieved very high coverage (median 99%), specificity (99.8%), and sensitivity (90% for haplotypes with 5% within sample frequency in dried blood spots with 100 parasites/µL). In silico analyses revealed that microhaplotypes provided much higher resolution to discriminate related from unrelated polyclonal infections than biallelic single-nucleotide polymorphism barcodes. CONCLUSIONS: The bioinformatic and laboratory methods outlined here provide a flexible tool for efficient, low-cost, high-throughput interrogation of the P. falciparum genome, and can be tailored to simultaneously address multiple questions of interest in various epidemiological settings.

Application of Microhaplotypes in Sibling Kinship Testing.

OBJECTIVES: To investigate the efficacy of different numbers of microhaplotype (MH) loci and the introduction of different reference samples on the identification of full sibling, half sibling and differentiation between full sibling and half sibling kinships, and to explore the effect of changing mutation rate on sibling testing. METHODS: First, a family map involving three generations was established, and four full sibling identification models, five half sibling identification models and five models distinguishing full and half siblings were constructed for different reference samples introduced. Based on the results of the previous study, two sets of nonbinary SNP-MH containing 34 and 54 loci were selected. Based on the above MH loci, 100 000 pairs of full sibling vs. unrelated individuals, 100 000 pairs of half sibling vs. unrelated individuals and 100 000 pairs of full sibling vs. half sibling were simulated based on the corresponding sibling kinship testing models, and the efficacy of each sibling kinship testing model was analyzed by the likelihood ratio algorithm under different thresholds. The mutant rate of 54 MH loci was changed to analyze the effect of mutation rate on sibling identification. RESULTS: In the same relationship testing model, the systematic efficacy of sibling testing was positively correlated with the number of MH loci detected. With the same number of MH loci, the efficacy of full sibling testing was better than that of uncle or grandfather when the reference sample introduced was a full sibling of A, but there was no significant difference in the identification efficacy of the four reference samples introduced for full sibling and half sibling differentiation testing. In addition, the mutation rate had a slight effect on the efficacy of sibling kinship testing. CONCLUSIONS: Increasing the number of MH loci and introducing reference samples of known relatives can increase the efficacy of full sibling testing, half sibling testing, and differentiation between full and half sibling kinships. The level of mutation rate in sibling testing by likelihood ratio method has a slight but insignificant effect on the efficacy.

Identification Strategy of Biological Half Sibling Relationship.

OBJECTIVES: To compare the application value of the likelihood ratio (LR) method and identity by state (IBS) method in the identification involving half sibling relationships, and to provide a reference for the setting of relevant standards for identification of half sibling relationship. METHODS: (1) Based on the same genetic marker combinations, the reliability of computer simulation method was verified by comparing the distributions of cumulated identity by state score (CIBS) and combined full sibling index in actual cases with the distributions in simulated cases. (2) In different numbers of three genetic marker combinations, the simulation of full sibling, half sibling and unrelated individual pairs, each 1 million pairs, was obtained; the CIBS, as well as the corresponding types of cumulative LR parameters, were calculated. (3) The application value of LR method was compared with that of IBS method, by comparing the best system efficiency provided by LR method and IBS method when genetic markers in different amounts and of different types and accuracy were applied to distinguish the above three relational individual pairs. (4) According to the existing simulation data, the minimum number of genetic markers required to distinguish half siblings from the other two relationships using different types of genetic markers was estimated by curve fitting. RESULTS: (1) After the rank sum test, under the premise that the real relationship and the genetic marker combination tested were the same, there was no significant difference between the simulation method and the results obtained in the actual case. (2) In most cases, under the same conditions, the system effectiveness obtained by LR method was greater than that by IBS method. (3) According to the existing data, the number of genetic markers required for full-half siblings and half sibling identification could be obtained by curve fitting when the system effectiveness reached 0.95 or 0.99. CONCLUSIONS: When distinguishing half sibling from full sibling pairs or unrelated pairs, it is recommended to give preference to the LR method, and estimate the required number of markers according to the identification types and the population data, to ensure the identification effect.

Development a multiplex panel of AISNPs, multi-allelic InDels, microhaplotypes, and Y-SNP/InDel loci for multiple forensic purposes via the NGS.

Recently, next generation sequencing has shown the promising application value in forensic research. In this study, we constructed a multiplex amplification system of different molecular genetic markers based on the previous selected ancestry informative single nucleotide polymorphisms (SNPs), multi-allelic insertion/deletion (InDel) polymorphisms, microhaplotypes, and Y-chromosomal SNP/InDel loci, and evaluated forensic efficiencies of the system in Chinese Shaanxi Han, Chinese Hui, and Chinese Mongolian groups via the next generation sequencing platform. Ancestry information analyses of Shaanxi Han, Hui, and Mongolian groups revealed that most Mongolian individuals could be differentiated from Shaanxi Hans and Huis based on the selected ancestry informative SNPs. Multi-allelic InDels and microhaplotypes showed the multiple allele variations and possessed relatively high genetic polymorphisms in these three groups, indicating these loci could provide higher forensic efficiencies for individual identification and paternity testing. Based on Y-chromosomal SNPs, different haplogroup distributions were observed among Shaanxi Han, Hui, and Mongolian groups. In conclusion, the self-developed system could be used to simultaneously carry out the individual identification, paternity analysis, mixture deconvolution, forensic ancestry information analysis, and Y-chromosomal haplogroup inference, which could provide more valuable investigative clues in forensic practices.

An overview of SNP-SNP microhaplotypes in the 26 populations of the 1000 Genomes Project.

Microhaplotypes (MHs) are a promising new type of forensic markers that are defined by the combinations of two- or more single-nucleotide polymorphisms (SNPs) within 200 bp. Their advantages, such as low mutation rates, lack of stutter artifacts, and short amplicons, have improved human identification, kinship analysis, ancestry prediction, and mixture deconvolution capabilities. Information on published MHs, e.g., allele frequencies, is available in widely used public databases, ALlele FREquency Database, and MicroHapDB. However, there are abundant non-published MHs spread over the whole genome, and those databases do not incorporate other databases (e.g., the SNP Database) to provide users with more integrated information. Therefore, it is essential to establish a robust, responsive, and comprehensive MHs database. In this study, we thoroughly screened for SNP-SNP MHs among 26 populations from the 1000 Genomes Project (Phase 3). All genotype data of SNPs in each MH were converted to PHASE input files, and allele frequencies were estimated using PHASE. We compiled a detailed summary of SNP-SNPs at the global, continental, and population levels focused on haplotypes and the Ae value and supplemented our database using dbSNP data (last updated in 2015). We have successfully established a dual-SNP MH database (D-SNPsDB) of MHs within 50 bp for 26 populations in the integration of basic data such as physical positions in the human genome, mapping of variant identifiers (rsIDs), allele frequencies, and basic variant information. For public database queries, the D-SNPsDB web app was developed with the R Shiny package to get integrated information.

A novel set of short microhaplotypes based on non-binary SNPs for forensic challenging samples.

Short tandem repeats (STRs) are the most widely used genetic markers in forensic application, but they are not ideal genetic markers for the analysis of forensic challenging samples such as highly degraded or unbalanced mixed samples because of their relatively large amplicons and stutter peaks. In this study, we developed a set of short microhaplotypes based on non-binary SNPs with molecular extent sizes no longer than 60 bases and genotyped 100 unrelated individuals from northern Han groups. Our results showed this panel has similar discrimination power to STR kits, as the combined random match probability (CMP) reached 1.396 × 10-22 and mean effective number of alleles (Ae) was 3.59. The cumulative probability of exclusion for duos (CPE-duos) was 0.999919 and the cumulative probability of exclusion for trios (CPE-trios) was 0.9999999987, suggesting this panel could be applied for forensic personal identification and parentage testing independently. Population differentiation in 26 populations from the 1000 Genomes Project indicated this panel could distinguish populations from Africa, East Asia, South Asia, America, and Europe. These microhaplotypes based on non-binary SNPs have short amplicons, good discrimination power, no stutter artifacts, and have great potential in detection of highly degraded and unbalanced mixtures for personal identification, paternity testing, and ancestry inference.

Population structure in a continuously distributed coastal marine species, the harbor porpoise, based on microhaplotypes derived from poor-quality samples.

Harbor porpoise in the North Pacific are found in coastal waters from southern California to Japan, but population structure is poorly known outside of a few local areas. We used multiplexed amplicon sequencing of 292 loci and genotyped clusters of single nucleotide polymoirphisms as microhaplotypes (N = 271 samples) in addition to mitochondrial (mtDNA) sequence data (N = 413 samples) to examine the genetic structure from samples collected along the Pacific coast and inland waterways from California to southern British Columbia. We confirmed an overall pattern of strong isolation-by-distance, suggesting that individual dispersal is restricted. We also found evidence of regions where genetic differences are larger than expected based on geographical distance alone, implying current or historical barriers to gene flow. In particular, the southernmost population in California is genetically distinct (FST  = 0.02 [microhaplotypes]; 0.31 [mtDNA]), with both reduced genetic variability and high frequency of an otherwise rare mtDNA haplotype. At the northern end of our study range, we found significant genetic differentiation of samples from the Strait of Georgia, previously identified as a potential biogeographical boundary or secondary contact zone between harbor porpoise populations. Association of microhaplotypes with remotely sensed environmental variables indicated potential local adaptation, especially at the southern end of the species' range. These results inform conservation and management for this nearshore species, illustrate the value of genomic methods for detecting patterns of genetic structure within a continuously distributed marine species, and highlight the power of microhaplotype genotyping for detecting genetic structure in harbor porpoises despite reliance on poor-quality samples.

Considering the flanking region variants of nonbinary SNP and phenotype-informative SNP to constitute 30 microhaplotype loci for increasing the discriminative ability of forensic applications.

The flanking region variants of nonbinary SNPs and phenotype-informative SNPs (piSNPs) have been observed, which may greatly improve the discriminative ability after constituting microhaplotype. In this study, 30 microhaplotype loci based on the nonbinary SNPs and piSNPs (shown to be related to phenotypes such as hair and eye color) were selected. Genotyping were conducted on 100 unrelated northern Han Chinese, and the 26 populations from the 1000 Genome Project were also included for comparison of populations differentiation. The simulated study was conducted for evaluating the efficiency of kinship testing. These 30 microhaplotype loci we selected had good polymorphism, with a mean effective number of alleles (Ae) of 3.46. The average Ae increase was 1.27 compared with the target SNPs. The populations from the five regions worldwide could also be distinguished using these loci. The results of kinship testing showed that these microhaplotype loci had the similar ability as 15 STR loci of AmpFlSTRR IdentifilerR PCR Amplification Kit to identify the biological parent and a stronger ability to exclude the nonbiological parents. So, these 30 microhaplotype loci may be multifunctional for forensic application, including the ability of personal identification and kinship testing equivalent to 15 STR loci, and the power of ancestry inference for distinguishing the main intercontinental population. Moreover, our selected phenotypic microhaplotype loci may theoretically have phenotype prediction capabilities. But the phenotype prediction efficiency of these phenotypic microhaplotype loci may be worse than that of piSNPs and the detailed prediction accuracy of different populations needs to be further studied.

Research progress on application of microhaplotype in forensic genetics.

As a novel genetic marker, microhaplotype can be applied in the field of forensic genetics. Microhaplotype has the advantages of high polymorphism, low mutation rate, no stutter products and short amplification fragments. Microhaplotype can effectively detect mixture, and quantitatively analyze the contributors of mixture. DNA with severe fragmentation can be successfully genotyped by microhaplotype. It can be used as ancestry informative marker to effectively divide the global continental population according to genetic structure. Microhaplotype system can provide more information than traditional short tandem repeat and help to identify complex relationships. It can provide new ideas for tumor source identification, cell line identification and prenatal paternity testing. Here we review the applications of microhaplotype, intending to provide references for forensic practice.

Chromosome 6p SNP microhaplotypes and IgG3 levels in hemochromatosis probands with HFE p.C282Y homozygosity.

Subnormal IgG1 or IgG3 levels occurred in 30% of hemochromatosis probands with HFE p.C282Y homozygosity and were concordant in HLA-identical siblings. We sought to identify factors associated with IgG subclasses in Alabama probands with p.C282Y homozygosity evaluated for 500 kb microhaplotypes AAT and GGG defined by SNPs in chromosome 6p genes PGBD1, ZNF193, and ZNF165. In regressions on IgG subclasses, we used: age; sex; GGG (dichotomous); iron removed to achieve depletion; CD8+ T-lymphocytes; and other IgG subclasses. Among 49 probands, AAT and GGG occurred in 95.9% and 16.3%, respectively. Thirteen probands (26.5%) had subnormal IgG1; 11 probands (22.4%) had subnormal IgG3. Mean IgG3 was higher in probands with than without GGG (75 mg/dL [95% confidence interval 63, 89] vs. 58 mg/dL [49, 71], respectively; p = 0.0321). Regression on IgG3 revealed: GGG positivity (p = 0.0106); and IgG1 (p = 0.0015). In a replication cohort of 22 Portugal probands with p.C282Y homozygosity, mean IgG3 was higher in probands with than without GGG (46 ±â€¯16 vs. 31 ±â€¯12 mg/dL, respectively; p = 0.0410). We conclude that mean IgG3 levels are higher in hemochromatosis probands with p.C282Y homozygosity with chromosome 6p microhaplotype GGG than in probands homozygous for microhaplotype AAT.

A new approach to detect a set of SNP-SNP markers: Combining ARMS-PCR with SNaPshot technology.

Microhaplotypes are a new promising type of forensic genetic marker. Without the interference of stutter and high mutation rates as for STRs, and with short amplification lengths and a higher degree of polymorphism than single SNP, microhaplotypes composed of two SNPs, SNP-SNP, have a strong application potential. Currently, the most common method to detect microhaplotypes is massive parallel sequencing. However, the cost and extensive use of instruments limit its wide application in forensic laboratories. In this study, we screened 23 new SNP-SNP loci and established a new detection method by combining a multiplex amplification refractory mutation system-based PCR (ARMS-PCR) and SNaPshot technology based on CE. First, we introduced an additional deliberate mismatch at the antepenultimate base from the 3' end of primers when designing ARMS-PCR for SNP 1 (the first SNP of the SNP-SNP). Then, single base extension primers for SNaPshot assay were designed next to the position of SNP 2 (the second SNP). Finally, 15 loci were successfully built into four panels and these loci showed a relatively high level of polymorphism in the Southwest Chinese Han population. All the loci had an average probability of informative genotypes (I value) of 0.319 and a combined discrimination power of 0.999999999. Therefore, this new detection system will provide a valuable supplement to current methods.

Developing and population analysis of a new multiplex panel of 18 microhaplotypes and compound markers using next generation sequencing and its application in the Shaanxi Han population.

Compound marker consists of two different types of genetic markers, like deletion/insertion polymorphism and single nucleotide polymorphism in the genomic region of 200 bp, and microhaplotype consists of a series of closely linked single nucleotide polymorphisms in a small DNA segment (<300 bp), which show great potential for human identifications and mixture analyses. In this study, we initially selected 23 novel genetic markers comprising 10 microhaplotypes and 13 compound markers according to previously reported single nucleotide polymorphism or deletion/insertion polymorphism loci. Genetic distributions of these 23 loci in different continental populations showed that they could be used as valuable loci for forensic human identification purpose. Besides, high informativeness values (>0.1) were observed in six loci which could be further employed for forensic ancestry analyses. Finally, 18 loci were successfully developed into a multiplex panel and detected by the next generation sequencing (NGS) technology. Further analyses of these 18 loci in the studied Shaanxi Han population showed that 15 loci exhibited relatively high expected heterozygosities (>0.5). Cumulative power of discrimination (0.999 999 999 99 4835) of these 18 loci revealed that the multiplex panel could also be utilized for human identifications in the studied Shaanxi Han population.

Dispersal of a nearshore marine fish connects marine reserves and adjacent fished areas along an open coast.

Marine species with pelagic larvae typically exhibit little population structure, suggesting long-distance dispersal and high gene flow. Directly quantifying dispersal of marine fishes is challenging but important, particularly for the design of marine protected areas (MPAs). Here, we studied kelp rockfish (Sebastes atrovirens) sampled along ~25 km of coastline in a boundary current-dominated ecosystem and used genetic parentage analysis to identify dispersal events and characterize them, because the distance between sedentary parents and their settled offspring is the lifetime dispersal distance. Large sample sizes and intensive sampling are critical for increasing the likelihood of detecting parent-offspring matches in such systems and we sampled more than 6,000 kelp rockfish and analysed them with a powerful set of 96 microhaplotype markers. We identified eight parent-offspring pairs with high confidence, including two juvenile fish that were born inside MPAs and dispersed to areas outside MPAs, and four fish born in MPAs that dispersed to nearby MPAs. Additionally, we identified 25 full-sibling pairs, which occurred throughout the sampling area and included all possible combinations of inferred dispersal trajectories. Intriguingly, these included two pairs of young-of-the-year siblings with one member each sampled in consecutive years. These sibling pairs suggest monogamy, either intentional or accidental, which has not been previously demonstrated in rockfishes. This study provides the first direct observation of larval dispersal events in a current-dominated ecosystem and direct evidence that larvae produced within MPAs are exported both to neighbouring MPAs and to proximate areas where harvest is allowed.

Evaluation of the microhaplotype markers in kinship analysis.

Microhaplotype markers are emerging forensic genetic markers, which may supplement existing markers. Consisting of two to four SNPs with an extent of <200 bp, microhaplotype can be genotyped through massively parallel sequencing technology. Articles that have been published suggested that microhaplotype markers have good application prospect in forensics. Multiallelic haplotype loci are potentially important in certain forensic works, as the stutter and high mutation rate of short tandem repeats and the low polymorphism of single nucleotide polymorphisms may limit the power of these two kind of regular markers. In this study, we explored the potential of 11 new microhaplotype loci in kinship analysis. The results suggested that these loci have relatively high polymorphic information in different populations worldwide and relatively high system effectiveness in the kinship analysis. Microhaplotypes have potential for forensic kinshipg analysis.