Quality control and analytic best practices for testing genetic models of sex differences in large populations.

Phenotypic sex-based differences exist for many complex traits. In other cases, phenotypes may be similar, but underlying biology may vary. Thus, sex-aware genetic analyses are becoming increasingly important for understanding the mechanisms driving these differences. To this end, we provide a guide outlining the current best practices for testing various models of sex-dependent genetic effects in complex traits and disease conditions, noting that this is an evolving field. Insights from sex-aware analyses will not only teach us about the biology of complex traits but also aid in achieving the goals of precision medicine and health equity for all.

Long-term prognosis of fatty-acid oxidation disorders in adults: Optimism despite the limited effective therapies available.

INTRODUCTION: Fatty-acid oxidation disorders (FAODs) are recessive genetic diseases. MATERIALS AND METHODS: We report here clinical and paraclinical data from a retrospective study of 44 adults with muscular FAODs from six French reference centers for neuromuscular or metabolic diseases. RESULTS: The study cohort consisted of 44 adult patients: 14 with carnitine palmitoyl transferase 2 deficiency (32%), nine with multiple acyl-CoA deficiency (20%), 13 with very long-chain acyl-CoA dehydrogenase deficiency (30%), three with long-chain 3-hydroxyacyl-CoA dehydrogenase deficiency (7%), and five with short-chain acyl-CoA dehydrogenase deficiency (11%). Disease onset occurred during childhood in the majority of patients (59%), with a mean age at onset of 15 years (range = 0.5-35) and a mean of 12.6 years (range = 0-58) from disease onset to diagnosis. The principal symptoms were acute muscle manifestations (rhabdomyolysis, exercise intolerance, myalgia), sometimes associated with permanent muscle weakness. Episodes of rhabdomyolysis were frequent (84%), with a mean creatinine kinase level of 68,958 U/L (range = 660-300,000). General metabolic complications were observed in 58% of patients, respiratory manifestations in 18% of cases, and cardiological manifestations in 9% of cases. Fasting acylcarnitine profile was used to orient genetic explorations in 65% of cases. After a mean follow-up of 10 years, 33% of patients were asymptomatic and 56% continued to display symptoms after exercise. The frequency of rhabdomyolysis decreased after diagnosis in 64% of cases. CONCLUSION: A standardized register would complete this cohort description of muscular forms of FAODs with exhaustive data, making it possible to assess the efficacy of therapeutic protocols in real-life conditions and during the long-term follow-up of patients.

Chemical investigation of biological trace evidence; toxicological screening of waste residues obtained from DNA extraction processes.

In a forensic scenario, if biological stains are found in very small quantities, these are usually collected for DNA analyses, considered essential for the forensic investigation and thus excluding possible investigations by other forensic disciplines as forensic toxicology. We developed an experimental study to evaluate the feasibility of analyzing DNA extraction residues obtained from DNA extraction procedures to perform toxicological analysis, with the aim to extract both genetic and toxicological information without affecting or compromising the genetic sample and/or DNA extraction. DNA extraction from four blood samples (fortified with 5 molecules of interest with a final concentrations of 1 µg/mL, 100 ng/mL, 10 ng/mL and 5 ng/mL, respectively) were analyzed with QIAGEN QIAmp® DNA Mini kit. Three waste residues collected from the DNA extraction were analyzed for the toxicological investigation via Solid-Phase Extraction and High-Performance Liquid Chromatography-Tandem Mass Spectrometry analyses (Thermo Scientific™ TSQ Fortis™ II Triple-Quadrupole Mass Spectrometer). The analytical investigation revealed that our analytes of interest were detected in two different residues of the DNA extraction procedure, allowing both genetic and toxicological analyses without affecting the DNA identification. At last, the experimental protocol was applied to a hypothetical case, with encouraging results and allowing the identification of our molecules of interest.

Stroke and frailty index: a two-sample Mendelian randomisation study.

INTRODUCTION: Previous observational studies have found an increased risk of frailty in patients with stroke. However, evidence of a causal relationship between stroke and frailty is scarce. The aim of this study was to investigate the potential causal relationship between stroke and frailty index (FI). METHODS: Pooled data on stroke and debility were obtained from genome-wide association studies (GWAS).The MEGASTROKE Consortium provided data on stroke (N = 40,585), ischemic stroke (IS,N = 34,217), large-vessel atherosclerotic stroke (LAS,N = 4373), and cardioembolic stroke (CES,N = 7 193).Summary statistics for the FI were obtained from the most recent GWAS meta-analysis of UK BioBank participants and Swedish TwinGene participants of European ancestry (N = 175,226).Two-sample Mendelian randomization (MR) analyses were performed by inverse variance weighting (IVW), weighted median, MR-Egger regression, Simple mode, and Weighted mode, and heterogeneity and horizontal multiplicity of results were assessed using Cochran's Q test and MR-Egger regression intercept term test. RESULTS: The results of the current MR study showed a significant correlation between stroke gene prediction and FI (odds ratio 1.104, 95% confidence interval 1.064 - 1.144, P < 0.001). In terms of stroke subtypes, IS (odds ratio 1.081, 95% confidence interval 1.044 - 1.120, P < 0.001) and LAS (odds ratio 1.037, 95% confidence interval 1.012 - 1.062, P = 0.005). There was no causal relationship between gene-predicted CES and FI. Horizontal multidimensionality was not found in the intercept test for MR Egger regression (P > 0.05), nor in the heterogeneity test (P > 0.05). CONCLUSIONS: This study provides evidence for a causal relationship between stroke and FI and offers new insights into the genetic study of FI.

Novel lambdapapillomavirus in northern sea otters Enhydra lutris kenyoni, associated with oral hyperplastic nodules.

A novel papillomavirus (PV) associated with hyperplastic nodules scattered over the muco-cutaneous border of the oral cavity of a dead, wild, subadult northern sea otter Enhydra lutris kenyoni (NSO) in 2004 in Homer, Alaska, USA, was genetically characterized. Primers for the amplification of 2 large overlapping DNA fragments that contained the complete genome of the NSO PV were designed. Sanger methodology generated sequences from which new specific primers were designed for the primer-walking approach. The NSO PV genome consists of 8085 nucleotides and contains an early region composed of E6, E7, E1, and E2 open reading frames (ORFs), an E4 ORF (contained within E2) lacking an in-frame proximal ATG start codon, an unusually long (907 nucleotide) stretch lacking any ORFs, a late region that contains the capsid genes L2 and L1, and a non-coding regulatory region (ncRR). This NSO PV has been tentatively named Enhydra lutris kenyoni PV2 (ElkPV2). Pairwise and multiple sequence alignments of the complete L1 ORF nucleotides and concatenated E1-E2-L1 amino acid sequences showed that the NSO PV is a novel PV, phylogenetically most closely related to southern sea otter PV1. The carboxy end of the E6 oncoprotein does not contain the PDZ-binding motif with a strong correlation with oncogenicity, suggesting a low-risk PV, which is in agreement with histopathological findings. However, the ElkPV2 E7 oncoprotein does contain the retinoblastoma (pRb) binding domain LXCXE (LQCYE in ElkPV2), associated with oncogenicity in some high-risk PVs. Further studies on the prevalence and clinical significance of ElkPV2 infections in NSO are needed.

Dipylidium caninum draft genome - a new resource for comparative genomic and genetic explorations of flatworms.

Here, we present a draft genome of the tapeworm Dipylidium caninum (family Dipylidiidae) and compare it with other cestode genomes. This draft genome of D. caninum is 110 Mb in size, has a repeat content of ~13.4% and is predicted to encode ~10,000 protein-coding genes. We inferred excretory/secretory molecules (representing the secretome), other key groups of proteins (including peptidases, kinases, phosphatases, GTPases, receptors, transporters and ion-channels) and predicted potential intervention targets for future evaluation. Using 144 shared single-copy orthologous sequences, we investigated the genetic relationships of cestodes for which nuclear genomes are available. This study provides first insights into the molecular biology of D. caninum and a new resource for comparative genomic and genetic explorations of this and other flatworms.

The Mouse Bile Duct Tapeworm, Hymenolepis Microstoma in Free-living Small Mammals in Slovakia: Occurrence and Genetic Analysis.

The mouse bile duct tapeworm Hymenolepis microstoma, is a potentially zoonotic species with a wide variety of reported definitive hosts of rodent genera. In the present study the occurrence of H. microstoma in free-living small mammals in selected areas of Slovakia and the retrospective analysis of epidemiological data published in Slovakia were performed. Hymenolepis microstoma was detected in two animal species, the common shrew (Sorex araneus) and the European hamster (Cricetus cricetus) of 186 small mammals examined from two ecosystems, urban and natural ecosystem of national park. No mention about the presence of this parasite in Slovakia in the past was found following a bibliographical search. Partial sequences of the nuclear paramyosin gene showed the shrew isolate placed in a subclade together with H. microstoma from Portugal, with high bootstrap value for its differentiation from the sister species Hymenolepis nana. Similarly, the analysis of the nuclear ribosomal ITS region placed the hamster isolate in the cluster composed of H. microstoma from Australia, Spain and Portugal. The Slovak isolate was the most distinctive sample among available H. microstoma, differing in 1.4 - 1.9% of nucleotides from the remaining isolates. The difference (seven of 17 nucleotide positions) was partially due to indel polymorphisms associated with two and five nucleotides. To our knowledge, these are the first reports of H. microstoma in Central Europe and also the first record of infection in the common shrew. A recently indicated zoonotic potential of H. microstoma along with a possibility of its direct transmission between animals and/or humans without the need of intermediate hosts pose a public health concern in contaminated areas of Slovakia. The use of molecular techniques may substantially facilitate more thorough understanding of the epidemiological situation of H. microstoma and related tapeworms in various ecosystems of the country.

Genetic analysis of Rhipicephalus sanguineus sensu lato ticks, parasites of dogs in the Canary Islands, Cyprus, and Croatia, based on mitochondrial 16S rRNA gene sequences.

The aim of this work was to perform an analysis based on mtDNA sequences of the 16S rRNA gene in order to determine the phylogenetic position of ticks belonging to the Rhipicephalus sanguineus group from the Canary Islands, Cyprus, and Croatia. All the haplotypes obtained from ticks collected in the Canary Islands and Croatia grouped with R. sanguineus sensu stricto from France, Portugal, Italy, Switzerland, Argentina, Chile, Uruguay, and the USA. The sequences of R. sanguineus sensu lato from Cyprus formed a clade with R. sanguineus s.l. from Egypt, Turkey, and Romania, which belongs to the "Rhipicephalus sp. morphotype I" or "southeastern European lineage." Ticks determined as R. turanicus s.l. from Cyprus clustered separately from the remaining clades of the R. sanguineus group, including R. turanicus s.s. The data show that R. sanguineus s.s. is present in the Canary Islands and Croatia, while R. sanguineus "southeastern lineage" is found in Cyprus.

Comparison of 2016-17 and Previous Epizootics of Highly Pathogenic Avian Influenza H5 Guangdong Lineage in Europe.

We analyzed the highly pathogenic avian influenza (HPAI) H5 epizootic of 2016-17 in Europe by epidemiologic and genetic characteristics and compared it with 2 previous epizootics caused by the same H5 Guangdong lineage. The 2016-17 epizootic was the largest in Europe by number of countries and farms affected and greatest diversity of wild birds infected. We observed significant differences among the 3 epizootics regarding region affected, epidemic curve, seasonality, and outbreak duration, making it difficult to predict future HPAI epizootics. However, we know that in 2005-06 and 2016-17 the initial peak of wild bird detections preceded the peak of poultry outbreaks within Europe. Phylogenetic analysis of 2016-17 viruses indicates 2 main pathways into Europe. Our findings highlight the need for global surveillance of viral changes to inform disease preparedness, detection, and control.

Genetic analyses in a cohort of 191 pulmonary arterial hypertension patients.

BACKGROUND: Pulmonary arterial hypertension (PAH) is a progressive and fatal disorder associated with high pulmonary artery pressure. Genetic testing enables early diagnosis and offers an opportunity for family screening. To identify genetic mutations and help make a precise diagnosis, we performed genetic testing in 191 probands with PAH and tried to analyze the genotype-phenotype correlation. METHODS: Initially, PAH samples (n = 119) were submitted to BMPR2 screening using Sanger sequencing. Later, we developed a PAH panel test to identify causal mutations in 13 genes related to PAH and tried to call BMPR2 copy number variations (CNVs) with the panel data. Multiplex ligation-dependent probe amplification (MLPA) was used to search for CNVs in BMPR2, ACVRL1 and ENG. Notably, EIF2AK4 gene was also involved in the panel, which allowed to distinguish pulmonary veno-occlusive disease (PVOD)/pulmonary capillary hemangiomatosis (PCH) patients from idiopathic PAH (IPAH). Characteristics of patients were compared using t test for continuous variables. RESULTS: Pathogenic BMPR2 mutations were detected most frequently in 32 (17.9%) IPAH and 5 (41.7%) heritable PAH (HPAH) patients by sequencing, and 12 BMPR2 CNVs called from the panel data were all successfully confirmed by MLPA analysis. In addition, homozygous or compound heterozygous EIF2AK4 mutations were identified in 6 patients, who should be corrected to a diagnosis of PVOD/PCH. Genotype-phenotype correlation analysis revealed that PAH patients with BMPR2 mutations were younger at diagnosis (27.2y vs. 31.6y, p = 0.0003) and exhibited more severe pulmonary hemodynamic impairment and a worse cardiac index compared with those without BMPR2 mutations. CONCLUSIONS: The panel assay represented a highly valuable tool in PAH genetic testing, not only for the detection of small sequence alterations, but also for an indication of BMPR2 CNVs, which had implications for the specific samples to perform further MLPA assay. Analyses of PAH causal genes have a great help to clinical diagnosis and deep implications in disease treatment.

Population genetic analyses and evaluation of 22 autosomal STRs in Indian populations.

The allele frequencies and forensically relevant parameters for the 22 autosomal short tandem repeats (STRs) present in PowerPlex® Fusion (Promega, Madison, WI) were determined for 357 unrelated individuals from 11 states across India. The combined power of discrimination and probability of exclusion were 0.99999999999999999999999999875 and 0.999999997200846, respectively. The panel was found to be informative for Indian populations and generated a total of 275 alleles. Further, analyses with these loci did not show any noticeable clustering among the Indian populations.

The Genetic Basis of Pheromone Evolution in Moths.

Moth sexual pheromones are widely studied as a fine-tuned system of intraspecific sexual communication that reinforces interspecific reproductive isolation. However, their evolution poses a dilemma: How can the female pheromone and male preference simultaneously change to create a new pattern of species-specific attraction? Solving this puzzle requires us to identify the genes underlying intraspecific variation in signals and responses and to understand the evolutionary mechanisms responsible for their interspecific divergence. Candidate gene approaches and functional analyses have yielded insights into large families of biosynthetic enzymes and pheromone receptors, although the factors controlling their expression remain largely unexplored. Intra- and interspecific crosses have provided tantalizing evidence of regulatory genes, although, to date, mapping resolution has been insufficient to identify them. Recent advances in high-throughput genome and transcriptome sequencing, together with established techniques, have great potential to help scientists identify the specific genetic changes underlying divergence and resolve the mystery of how moth sexual communication systems evolve.

Prenatal diagnosis of Simpson-Golabi-Behmel syndrome.

Simpson-Golabi-Behmel syndrome (SGBS) is an overgrowth syndrome and it is usually diagnosed postnatally, on the basis of phenotype. Prenatal ultrasonography may show fetal alterations, but they are not pathognomonic and most of them are frequently detectable only from the 20th week of gestation. Nevertheless, early diagnosis is important to avoid neonatal complications and make timely and informed decisions about the pregnancy. We report on four fetuses from two unrelated families, in whom the application of whole exome sequencing and array-CGH allowed the identification of GPC3 alterations causing SGBS. The careful follow up of pregnancies and more sophisticated analysis of ultrasound findings led to the identification of early prenatal alterations, which will improve the antenatal diagnosis of SGBS. © 2016 Wiley Periodicals, Inc.

Metalloproteins as risk factors for osteoarthritis: improving and understanding causal estimates using Mendelian randomization.

Osteoarthritis (OA) is one of the most prevalent musculoskeletal disorders and a primary cause of pain and disability among the elderly population. Research on the relationship between metalloproteins (MPs) and OA is limited, and causality remains unclear. Our objective is to utilize Mendelian randomization (MR) to explore the possible causal relationship between MPs and OA. The data on MPs were derived from a Genome-Wide Association Study (GWAS) analysis involving 3301 samples. The GWAS data for OA were obtained from an analysis involving 462,933 European individuals. In this study, a variety of two-sample Mendelian randomization methods (two-sample MR) to evaluate the causal effect of MPs on OA, including inverse variance weighted method (IVW), MR-Egger method, weighted median method (WM), simple mode, weight mode, and Wald ratio. The primary MR analysis using the IVW method reveals a significant negative correlation between Metallothionein-1F (MT-1F), zinc finger protein 134 (ZNF134), calcium/calmodulin-dependent protein kinase type 1D (CAMK1D), and EF-hand calcium-binding domain-containing protein 14 (EFCAB14) with the occurrence of osteoarthritis (OA) (p value < 0.05). However, no causal relationship was observed in the opposite direction between these MPs and OA. Notably, even in combined models accounting for confounding factors, the negative association between these four MPs and OA remained significant. Sensitivity analysis demonstrated no evidence of horizontal pleiotropy or heterogeneity, and leave-one-out analysis confirmed the robustness of the results. In this study, we have established a conspicuous association between four distinct MPs and OA. This discovery augments our understanding of potential avenues for the diagnosis and treatment of this condition. Key Points • The MR method was employed to assess the relationship between MPs and OA. • A total of four types of MPs have demonstrated inhibitory effects on the occurrence of OA.

Causality between Ankylosing Spondylitis and osteoarthritis in European ancestry: a bidirectional Mendelian randomization study.

Objective: To explore the bidirectional causal relationship between Ankylosing Spondylitis (AS) and Osteoarthritis (OA) at the genetic level within the European ancestry. Methods: We implemented a series of quality control steps to select instrumental variables (IVs) related to the exposure. We conducted two-sample Mendelian randomization (MR) using the inverse-variance weighted method as the primary approach. We adjusted significance levels using Bonferroni correction, assessed heterogeneity using Cochrane's Q test. Sensitivity analysis was conducted through leave-one-out method. Additionally, external datasets and relaxed IV selection criteria were employed, and multivariate MR analyses were performed for validation purposes. Finally, Bayesian colocalization (COLOC) analysis identified common genes, validating the MR results. Results: The investigation focused on the correlation between OA and AS in knee, hip, and hand joints. MR results revealed that individuals with AS exhibit a decreased risk of knee OA (OR = 0.9882, 95% CI: 0.9804-0.9962) but no significant increase in the risk of hip OA (OR = 0.9901, 95% CI: 0.9786-1.0018). Conversely, AS emerged as a risk factor for hand OA (OR = 1.0026, 95% CI: 1.0015-1.0036). In reverse-direction MR analysis, OA did not significantly influence the occurrence of AS. Importantly, minimal heterogeneity was observed in our MR analysis results (p > 0.05), and the robustness of these findings was confirmed through sensitivity analysis and multivariate MR analysis. COLOC analysis identified four colocalized variants for AS and hand OA (rs74707996, rs75240935, rs181468789, and rs748670681). Conclusion: In European population, individuals with AS have a relatively lower risk of knee OA, whereas AS serves as a risk factor for hand OA. However, no significant causal relationship was found between AS and hip OA. Additionally, it offers novel insights into genetic research on AS and OA.

Mitochondrial DNA: Inherent Complexities Relevant to Genetic Analyses.

Mitochondrial DNA (mtDNA) exhibits distinct characteristics distinguishing it from the nuclear genome, necessitating specific analytical methods in genetic studies. This comprehensive review explores the complex role of mtDNA in a variety of genetic studies, including genome-wide, epigenome-wide, and phenome-wide association studies, with a focus on its implications for human traits and diseases. Here, we discuss the structure and gene-encoding properties of mtDNA, along with the influence of environmental factors and epigenetic modifications on its function and variability. Particularly significant are the challenges posed by mtDNA's high mutation rate, heteroplasmy, and copy number variations, and their impact on disease susceptibility and population genetic analyses. The review also highlights recent advances in methodological approaches that enhance our understanding of mtDNA associations, advocating for refined genetic research techniques that accommodate its complexities. By providing a comprehensive overview of the intricacies of mtDNA, this paper underscores the need for an integrated approach to genetic studies that considers the unique properties of mitochondrial genetics. Our findings aim to inform future research and encourage the development of innovative methodologies to better interpret the broad implications of mtDNA in human health and disease.

Clinical Characteristics and Genetic Analyses of Patients with Idiopathic Hypogonadotropic Hypogonadism

Objective: Idiopathic hypogonadotropic hypogonadism (IHH) is classified into two groups-Kalman syndrome and normosmic IHH (nIHH). Half of all cases can be explained by mutations in >50 genes. Targeted gene panel testing with nexrt generation sequencing (NGS) is required for patients without typical phenotypic findings. The aim was to determine the genetic etiologies of patients with IHH using NGS, including 54 IHH-associated genes, and to present protein homology modeling and protein stability analyzes of the detected variations. Methods: Clinical and demographic data of 16 patients (eight female), aged between 11.6-17.8 years, from different families were assessed. All patients were followed up for a diagnosis of nIHH, had normal cranial imaging, were without anterior pituitary hormone deficiency other than gonadotropins, had no sex chromosome anomaly, had no additional disease, and underwent genetic analysis with NGS between the years 2008-2021. Rare variants were classified according to the variant interpretation framework of the American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology. Changes in protein structure caused by variations were modeled using RoseTTAFold and changes in protein stability resulting from variation were analyzed. Results: Half of the 16 had no detectable variation. Three (18.75%) had a homozygous (pathogenic) variant in the GNRHR gene, one (6.25%) had a compound heterozygous [likely pathogenic-variants of uncertain significance (VUS)] variant in PROK2 and four (25%) each had a heterozygous (VUS) variant in HESX1, FGF8, FLRT3 and DMXL2. Protein models showed that variants interpreted as VUS according to ACMG could account for the clinical IHH. Conclusion: The frequency of variation detection was similar to the literature. Modelling showed that the variant in five different genes, interpreted as VUS according to ACMG, could explain the clinical IHH.

A complex CLCN1 variant associated with hereditary myotonia in a mixed-breed dog.

Hereditary myotonia (HM) is characterized by delayed muscle relaxation after contraction as a result of a mutation in the CLCN1 gene. We describe here a complex CLCN1 variant in a mixed-breed dog with clinical and electromyographic signs of HM. Blood samples from the myotonic dog, as well as from his male littermate and parents, were analyzed via amplification of the 23 exons encoding CLCN1. After sequencing the CLCN1 gene, a complex variant was found in exon 6 c.[705T>G; 708del; 712_732del], resulting in a premature stop codon in exon 7 and a protein that was 717 amino acids shorter than the normal CLC protein. The myotonic dog was identified as homozygous recessive for the complex CLCN1 variant; its parents were heterozygous, and its male littermate was homozygous wild-type. Knowledge of the CLCN1 mutations responsible for the development of hereditary myotonia allows greater clarification of this condition.

Genetic structure and characteristics of Tibetan chickens.

Tibetan chicken is one of the most common and widely distributed highland breeds, and is often used as a model organism for understanding genetic adaptation to extreme environments in Tibet. Despite its apparent geographical diversity and large variations in plumage patterns, the genetic differences within breed were not accounted for in most studies and have not been systematically investigated. In order to reveal and genetically differentiate the current existing TBC sub-populations that might have major implications for genomic research in TBCs, we systematically evaluated the population structure and demography of current TBC populations. Based on 344 whole-genome sequenced birds including 115 Tibetan chickens that were mostly sampled from family-farms across Tibet, we revealed a clear separation of Tibetan chickens into 4 sub-populations that broadly aligns with their geographical distribution. Moreover, population structure, population size dynamics, and the extent of admixture jointly suggest complex demographic histories of these sub-populations, including possible multiple origins, inbreeding, and introgressions. While most of the candidate selected regions found between the TBC sub-populations and Red Jungle fowls were nonoverlapping, 2 genes RYR2 and CAMK2D were revealed as strong selection candidates in all 4 sub-populations. These 2 previously identified high altitude associated genes indicated that the sub-populations responded to similar selection pressures in an independent but functionally similar fashion. Our results demonstrate robust population structure in Tibetan chickens that will help inform future genetic analyses on chickens and other domestic animals alike in Tibet, recommending thoughtful experimental design.

Insertion/deletion polymorphism for genetic background and forensic performance exploration of the Sui group from Guizhou.

Insertion/deletion polymorphisms (InDels) as ideal genetic markers for forensic genetics are appreciated by scholars both nationally and internationally because they integrated the favorable features of single nucleotide polymorphisms (SNPs) and short tandem repeats (STRs). Nevertheless, with the limited identification efficiency of InDels, the multiplex amplification systems of InDels might just be applied as the supplementary methods in paternity testing with respect to commonly used STRs. In the current research, we successfully genotyped 105 unrelated individuals from the Guizhou Sui population based on a six-color fluorescence multiplex panel that could simultaneously detect 64 genetic markers (59 autosomal InDels, two autosomal miniSTRs and three Y chromosomal genetic markers). In addition, frequency distributions and forensic statistical parameters of these loci in the Sui group were assessed using the STRAF software. Phylogenetic relationships among the Sui group and other reference populations were dissected by two methods (principal component analysis and phylogenetic trees) based on 59 InDels. The combined discrimination power and probability of exclusion values of 61 autosomal genetic markers in the Sui group were nearly equal to 1-1.90063 × 10-27 and 0.999998272, respectively. Furthermore, we observed that the Sui group from Guizhou had closer genetic affinities with East Asian populations with respect to other continental populations. In summary, we stated that the multiplex amplification system might be utilized as a prospective independent tool for human individual identification and parentage testing in the Sui group residing in Guizhou.