Application of a newly constructed NGS panel with 45 X-linked microhaplotypes demonstrates the unique value of X-MH for kinship testing and mixture analysis.

X-linked microhaplotypes (X-MHs) have the potential to be a valuable supplementary tool in complex kinship identification or the resolution of DNA mixtures, because they bring together the distinctive genetic pattern of X chromosomal markers and the benefits of microhaplotypes (MHs). In this study, we used the 1000 Genome database to screen and select 63 X-MHs; 18 MHs were filtered out though a batch sequencing assessment of the DNA samples collected from 112 unrelated Chinese Han individuals. The resulting 45-plex panel performed well in comprehensive assessments including repeatability, sensitivity, species specificity, resistance to PCR inhibitors or degradation, mutation rate, and accuracy in detecting DNA mixture samples. The minimum amount of DNA template that can be tested with this panel is 0.5 ng. Additionally, the alleles of the minor contributor can be accurately detected when the mixture rate is larger than 1:9 in female-male mixture or 1:19 in male-male mixture. Then, we calculated population parameters on each MH based on the allele frequency data obtained from the sequence results of the aforementioned 112 unrelated samples. Combining these parameters on each MH, it can be calculated that TDPm, TDPf, CPET, CPEDFM, CPEDFF and CNCEP3 of the 45-plex system were 1-8.99×10-13, 1-1.62×10-19, 0.9999999995, 0.9999981, 0.9955, 0.9999971 and 0.99940, respectively, indicating that the panel is capable in personal identification and parentage testing. To reveal the unique advantage of X-MHs in the analyses of complex kinship and male DNA mixture, further assessments were made. For complex kinship identification, 22 types of individual pairs with different second-degree kinship were simulated and different types of likelihood ratios (LR) were calculated for each. The results revealed that the panel can achieve accuracy of approximately 70 %∼80 % when dividing each of the three types of second-degree kinships into three or four groups. Theoretically, such sub-division cannot be done by using independent autosomal markers. For male DNA mixture analysis without suspects, the maximum likelihood ratio strategy was derived and employed in the estimation of the number of male contributors (NOMC). Simulations were conducted to verify the efficacy of the 45-plex panel in the field and to compare it with autosomal markers by assuming the 45 MHs as autosomal ones. The results showed that X-MHs can achieve higher accuracy in the estimation of NOMC than autosomal ones when the mixed males were unrelated. The results highlighted the unique value of X-linked MHs in complex kinship and male mixture analyses.

The association between DNA methylation and human height and a prospective model of DNA methylation-based height prediction.

As a vital anthropometric characteristic, human height information not only helps to understand overall developmental status and genetic risk factors, but is also important for forensic DNA phenotyping. We utilized linear regression analysis to test the association between each CpG probe and the height phenotype. Next, we designed a methylation sequencing panel targeting 959 CpGs and subsequent height inference models were constructed for the Chinese population. A total of 11,730 height-associated sites were identified. By employing KPCA and deep neural networks, a prediction model was developed, of which the cross-validation RMSE, MAE and R2 were 5.62 cm, 4.45 cm and 0.64, respectively. Genetic factors could explain 39.4% of the methylation level variance of sites used in the height inference models. Collectively, we demonstrated an association between height and DNA methylation status through an EWAS analysis. Targeted methylation sequencing of only 959 CpGs combined with deep learning techniques could provide a model to estimate human height with higher accuracy than SNP-based prediction models.

Differences of circular RNA expression profiles between monozygotic twins' blood, with the forensic application in bloodstain and saliva.

Monozygotic twins (MZTs) possess identical genomic DNA sequences and are usually indistinguishable through routine forensic DNA typing methods, which can be relevant in criminal and paternity cases. Recently, novel epigenetic methods involving DNA methylation and microRNA analysis have been introduced to differentiate MZTs. In this study, we explore the potential of using epigenetic markers, specifically circular RNAs (circRNAs), a type of non-coding RNA (ncRNA), to identify MZTs, and investigate the unique expression patterns of circRNAs within pairs of MZTs, enabling effective differentiation. Epigenetics regulates gene expression at the post-transcriptional level and plays a crucial role in cell growth and aging. CircRNAs, a recently characterized subclass of ncRNA, have a distinct covalent loop structure without the typical 5' cap or 3' tail. They have been reported to modulate various cellular processes and play roles in embryogenesis and eukaryotic development. To achieve this, we conducted a comprehensive circRNA sequencing analysis (circRNA-seq) using total RNA extracted from the blood samples of five pairs of MZTs. We identified a total of 15,257 circRNAs in all MZTs using circRNA-seq. Among them, 3, 21, 338, and 2967 differentially expressed circRNAs (DEcircRNAs) were shared among five, four, three, and two pairs of MZTs, respectively. Subsequently, we validated twelve selected DEcircRNAs using real-time quantitative polymerase chain reaction (RT-qPCR) assays, which included hsa_circ_0004724, hsa_circ_0054196, hsa_circ_004964, hsa_circ_0000591, hsa_circ_0005077, hsa_circ_0054853, hsa_circ_0054716, hsa_circ_0002302, hsa_circ_0004482, hsa_circ_0001103, novel_circ_0030288 and novel_circ_0056831. Among them, hsa_circ_0005077 and hsa_circ_0004482 exhibited the best performance, showing differences in 7 out of 10 pairs of MZTs. These twelve differentially expressed circRNAs also demonstrated strong discriminative power when tested on saliva samples from 10 pairs of MZTs. Notably, hsa_circ_0004724 displayed differential expression in 8 out of 10 pairs of MZTs in their saliva. Additionally, we evaluated the detection sensitivity, longitudinal temporal stability, and suitability for aged bloodstains of these twelve DEcircRNAs in forensic scenarios. Our findings highlight the potential of circRNAs as molecular markers for distinguishing MZTs, emphasizing their suitability for forensic application.

Feature selection with a genetic algorithm can help improve the distinguishing power of microbiota information in monozygotic twins' identification.

Introduction: Personal identification of monozygotic twins (MZT) has been challenging in forensic genetics. Previous research has demonstrated that microbial markers have potential value due to their specificity and long-term stability. However, those studies would use the complete information of detected microbial communities, and low-value species would limit the performance of previous models. Methods: To address this issue, we collected 80 saliva samples from 10 pairs of MZTs at four different time points and used 16s rRNA V3-V4 region sequencing to obtain microbiota information. The data formed 280 inner-individual (Self) or MZT sample pairs, divided into four groups based on the individual relationship and time interval, and then randomly divided into training and testing sets with an 8:2 ratio. We built 12 identification models based on the time interval ( ≤ 1 year or ≥ 2 months), data basis (Amplicon sequence variants, ASVs or Operational taxonomic unit, OTUs), and distance parameter selection (Jaccard distance, Bray-Curist distance, or Hellinger distance) and then improved their identification power through genetic algorithm processes. The best combination of databases with distance parameters was selected as the final model for the two types of time intervals. Bayes theory was introduced to provide a numerical indicator of the evidence's effectiveness in practical cases. Results: From the 80 saliva samples, 369 OTUs and 1130 ASVs were detected. After the feature selection process, ASV-Jaccard distance models were selected as the final models for the two types of time intervals. For short interval samples, the final model can completely distinguish MZT pairs from Self ones in both training and test sets. Discussion: Our findings support the microbiota solution to the challenging MZT identification problem and highlight the importance of feature selection in improving model performance.

Development and evaluation of a novel panel containing 188 microhaplotypes for 2nd-degree kinship testing in the Hebei Han population.

Distant kinship identification is one of the critical problems in forensic genetics. As a new type of genetic marker defined and discussed in the last decade, the microhaplotype (MH) has drawn much attention in such identification owing to its specific advantages to traditional short tandem repeat (STR) or single nucleotide polymorphism (SNP) markers. In this study, MH markers were screened step by step from the 1000 Genomes Project database, and a novel multiplex panel containing 188 MHs (in which 181 are reported the first time, while 1 was reported in a previous study and the other 6 have partial overlaps with known markers) was constructed for application in 2nd- and 3rd-degree kinship identification. Along with the construction, a novel MH nomenclature was proposed, in which the SNP position information they contained was taken into account to eliminate the possibility that the same locus was named differently interlaboratory. After a series of evaluations, the panel was shown to have good sequencing accuracy, high sensitivity, species specificity, and resistance to anti-PCR inhibitors or degradation. Population data of the 188 MHs were calculated based on the genetic information of 221 unrelated Hebei Han individuals, and the effective number of alleles (Ae) ranged from 2.0925 to 8.2634 (with an average of 2.9267). For the whole system, the cumulative matching probability (CMP), the cumulative power of exclusion in paternity testing of duos (CPEduo) and that of trios (CPEtrio) reached 2.8422 × 10-137, 1-1.3109 × 10-21, and 1-2.8975 × 10-39, respectively, indicating that this panel was satisfactory for individual identification and paternity testing. Then, the efficiency of the 188 MHs in 2nd- and 3rd-degree kinship testing was studied based on 30 extended families consisting of 179 2nd-degree and 121 3rd-degree relatives, as well as simulations of 0.5 million pairs of those two kinships. The results showed that clear opinions would be given in 83.36% of 2nd-degree identifications with a false rate less than 10-5, when the confirming and excluding thresholds of cumulative likelihood ratio (CLR) were set as 104 and 10-4, respectively. This panel is still not sufficient to solve the problem of 3rd-degree kinship identification alone, and approximately 300 or 870 MH loci would be needed in 2nd- or 3rd-degree kinship identification, respectively, to achieve a system efficiency not less than 0.99 with such a threshold set; such necessary numbers would be used only as a reference in further research.

A study of strong nucleosomes in the human genome.

Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, nucleosomes are highly dynamic and susceptible to experimental conditions, resulting in extreme variability across nucleosome maps, which complicates the generation of accurate nucleosome organization data. We mapped nucleosomes from different individuals using improved MNase-seq. The improvements included setting different digestion levels (low, medium, high) and naked DNA correction to remove the noise caused by experimental manipulation and comparing maps to obtain the accurate position and occupancy of strong nucleosomes (SNs) in the whole genome. In addition, the characteristics of SNs were further excavated. SNs were enriched in Alu elements and near the centromere of Chr12. SNs contain some specific sequences, and the GC content of SNs is different from that of dynamic nucleosomes. The findings suggest that nucleosome location in the genome and the DNA sequence may affect nucleosome stability.

Response Surface of Speed-Loading Path to Grain Refinement during Current-Heating Compression at SAE 5137H Steel.

In thermal deformation of materials, grain refinement induced by dynamic recrystallization (DRX) is often pursued to obtain excellent mechanical properties. Here, the thermal deformation behaviors of SAE 5137H steel were investigated and characterized at temperature and strain rate range of 1123-1483 K and 0.01-10 s-1. Meanwhile, a design approach in speed-loading paths for grain refinement during current-heating compression was proposed, and these paths are linked to a typical three-dimensional (3D) response surface. Depending on the acquired stress-strain curves, the flow behaviors of this steel were analyzed and the typical 3D processing map was constructed to clarify the stable processing parameter domains during the continuous deformation process. Then, by the typical 3D processing map and microstructure observation, the 3D deformation mechanism map was constructed to connect the processing parameters and microstructural mechanisms. Subsequently, the 3D activation energy map was constructed to evaluate these deformation mechanisms, and the enhanced deformation mechanism map was constructed. Eventually, based on the enhanced deformation mechanism map, the speed-loading paths for SAE 5137H steel during current-heating compression were designed and they are mapped in a 3D response surface.

Circular RNA as a Potential Biomarker for Forensic Age Prediction.

In forensic science, accurate estimation of the age of a victim or suspect can facilitate the investigators to narrow a search and aid in solving a crime. Aging is a complex process associated with various molecular regulations on DNA or RNA levels. Recent studies have shown that circular RNAs (circRNAs) upregulate globally during aging in multiple organisms such as mice and C.elegans because of their ability to resist degradation by exoribonucleases. In the current study, we attempted to investigate circRNAs' potential capability of age prediction. Here, we identified more than 40,000 circRNAs in the blood of thirteen Chinese unrelated healthy individuals with ages of 20-62 years according to their circRNA-seq profiles. Three methods were applied to select age-related circRNA candidates including the false discovery rate, lasso regression, and support vector machine. The analysis uncovered a strong bias for circRNA upregulation during aging in human blood. A total of 28 circRNAs were chosen for further validation in 30 healthy unrelated subjects by RT-qPCR, and finally, 5 age-related circRNAs were chosen for final age prediction models using 100 samples of 19-73 years old. Several different algorithms including multivariate linear regression (MLR), regression tree, bagging regression, random forest regression (RFR), and support vector regression (SVR) were compared based on root mean square error (RMSE) and mean average error (MAE) values. Among five modeling methods, regression tree and RFR performed better than the others with MAE values of 8.767 years (S.rho = 0.6983) and 9.126 years (S.rho = 0.660), respectively. Sex effect analysis showed age prediction models significantly yielded smaller prediction MAE values for males than females (MAE = 6.133 years for males, while 10.923 years for females in the regression tree model). In the current study, we first used circRNAs as additional novel age-related biomarkers for developing forensic age estimation models. We propose that the use of circRNAs to obtain additional clues for forensic investigations and serve as aging indicators for age prediction would become a promising field of interest.

Highly accurate mtGenome haplotypes from long-read SMRT sequencing can distinguish between monozygotic twins.

Discriminating between monozygotic twins (MZ) remains a challenge in the field of forensics globally. It is very difficult to find sequence variants within MZ twins, despite using ultra-deep next generation sequencing (NGS) for nuclear DNA. However, mitochondrial DNA might be a potential marker owing to its higher mutation rate and easier sequencing via NGS. Here, we aimed to introduce a long-read single molecule real-time sequencing (SMRT) strategy, with better continuity and fewer alignment errors, to obtain more accurate mitochondrial genome (mtGenome) sequence on the Sequel platform. Compared to Ion Torrent Personal Genome Machine (PGM), the long-read SMRT sequencing strategy generated highly accurate and mapped circular consensus sequence (CCS) reads and exhibited robust performance in terms of reliable repeatability, consistent coverage pattern, and balanced strands in mtGenome recovery. Moreover, the long-read SMRT strategy exhibited superior ability to not only identify accurate haplotypes but also discover a total of 785 low-level variants within 16 MZ twin pairs with threshold of 2% and 20 CCS reads with Q30 quality. Taken together, our findings suggested the long-read SMRT technology as an appreciable strategy for obtaining accurate mitotypes and providing a promising solution for distinguishing between MZ twin pairs in forensic genetics.

Development of an NGS panel containing 42 autosomal STR loci and the evaluation focusing on secondary kinship analysis.

High-throughput next-generation sequencing (NGS) is a feasible technique to detect considerably more markers and simultaneously obtain length and sequence information in a single reaction. In this study, we developed an NGS panel including 42 commonly used autosomal short tandem repeats (STRs) and amelogenin on the Illumina MiSeq FGx™. Sequencing accuracy was validated by the consistency of 2800M Control DNA detected using the ForenSeq™ DNA Signature Prep Kit and Sanger sequencing. Nomenclature incompatibility was found between NGS-STR and CE-STR typing at 9 loci (D3S3045, D6S477, D7S3048, D9S925, D14S608, D17S1290, D18S535, D21S1270, GATA198B05), despite the correct sequence. The difference was caused by the two different methods of identifying motif sequence and a one-to-one correspondence can be found. We evaluated the panel by investigating consistency, sequencing sensitivity and the effectiveness of the 2nd-degree relationship identification. Herein, we present sequencing results from 58 unrelated individuals of the Hebei Han population. The total discrimination power (TDP) and cumulative probability of exclusion for trio paternity testing (CPEtrio) of the 42 NGS-STR panels reached 1-2.84 × 10-57 and 1-9.87 × 10-21, respectively. By family simulation and likelihood ratio (LR) calculation, this panel was shown to have effectiveness for the 2nd-degree kinship identification similar to the ForenSeq™ DNA Signature Prep Kit and certain advantages compared with it due to the relatively small number of loci. As expected, it provides new data for the development of NGS-STR typing technology.

Uncovering the relationship and mechanisms of Tartary buckwheat (Fagopyrum tataricum) and Type II diabetes, hypertension, and hyperlipidemia using a network pharmacology approach.

BACKGROUND: Tartary buckwheat (TB), a crop rich in protein, dietary fiber, and flavonoids, has been reported to have an effect on Type II diabetes (T2D), hypertension (HT), and hyperlipidemia (HL). However, limited information is available about the relationship between Tartary buckwheat and these three diseases. The mechanisms of how TB impacts these diseases are still unclear. METHODS: In this study, network pharmacology was used to investigate the relationship between the herb as well as the diseases and the mechanisms of how TB might impact these diseases. RESULTS: A total of 97 putative targets of 20 compounds found in TB were obtained. Then, an interaction network of 97 putative targets for these compounds and known therapeutic targets for the treatment of the three diseases was constructed. Based on the constructed network, 28 major nodes were identified as the key targets of TB due to their importance in network topology. The targets of ATK2, IKBKB, RAF1, CHUK, TNF, JUN, and PRKCA were mainly involved in fluid shear stress and the atherosclerosis and PI3K-Akt signaling pathways. Finally, molecular docking simulation showed that 174 pairs of chemical components and the corresponding key targets had strong binding efficiencies. CONCLUSION: For the first time, a comprehensive systemic approach integrating drug target prediction, network analysis, and molecular docking simulation was developed to reveal the relationships and mechanisms between the putative targets in TB and T2D, HT, and HL.

The complete chloroplast genome sequence of Chrysanthemum indicum.

Chrysanthemum indicum, an important medicinal plant of Asteraceae, had a long history in use for medicine in China. In this study, the complete chloroplast genome of C. indicum was sequenced by a 454 sequencing platform, and the structure of the obtained chloroplast genome was also analyzed. The complete chloroplast genome of C. indicum was 150 972 bp in length and had a pair of inverted repeats (IR, 24 956 bp) separated by a large (LSC, 82 741 bp) and small single copy (SSC, 18 319 bp) regions. Its total GC content was 37.48%. There were 126 chloroplast genes including 83 protein-coding genes, 35 tRNAs and eight rRNAs were successfully annotated. Sixteen genes contained one or two introns. Phylogenetic analyses declared that the chloroplast genome could distinguish C. indicum from its closely related species and might become a potential super barcode for the identification of these species.

The complete chloroplast genome sequence of Safflower (Carthamus tinctorius L.).

Safflower (Carthamus tinctorius L.) is a traditional medical plants of Asia. In this study, the complete chloroplast genome of safflower was presented. The total genome size was 153,675 bp in length, containing a pair of inverted repeats (IRs) of 25,407 bp, separated by large single copy (LSC) and small single copy (SSC) of 83,606 bp and 19,156 bp, respectively. Overall GC content of the genome was 37.4%. The chloroplast genome harbored 127 annotated genes, including 89 protein coding genes, 30 tRNA genes and 8 rRNA genes. A total of 7 of these genes were duplicated in the inverted repeat regions. Twelve genes contained one intron.

The complete chloroplast genome sequence of Perilla frutescens (L.).

Perilla frutescens (L.) is a traditionally medical herb of East Asia. The complete chloroplast genome of P. frutescens (L.) Britton var. frutescens was assembled in this study. Total chloroplast genome size of Perilla was 153,666 bp in length, containing a pair of inverted repeats (IRs) of 25,677 bp, separated by large single copy (LSC) and small single copy (SSC) of 84,288 bp and 17,925 bp, respectively. Overall GC contents of the genome were 37.6%. The chloroplast genome harbored 127 annotated genes, including 89 protein-coding genes, 29 tRNA genes and 8 rRNA genes. Eleven genes contained one or two introns.

The complete chloroplast genome sequence of Fagopyrum cymosum.

Fagopyrum cymosum is a traditional medicinal plant. In this study, the complete chloroplast genome of Fagopyrum cymosum is presented. The total genome size is 160,546 bp in length, containing a pair of inverted repeats (IRs) of 32,598 bp, separated by large single copy (LSC) and small single copy (SSC) of 84,237 bp and 11,014 bp, respectively. Overall GC contents of the genome were 36.9%. The chloroplast genome harbors 126 annotated genes, including 91 protein coding genes, 29 tRNA genes, and six rRNA genes. Eighteen genes contain one or two introns. Phylogenetic analyses indicated a clear evolutionary relationship among species of Caryophyllales.

The development dynamics of the maize root transcriptome responsive to heavy metal Pb pollution.

Lead (Pb), as a heavy metal element, has become the most important metal pollutant of the environment. With allocating a relatively higher proportion of its biomass in roots, maize could be a potential important model to study the phytoremediation of Pb-contaminated soil. Here we analyzed the maize root transcriptome of inbred lines 9782 under heavy metal lead (Pb) pollution, which was identified as a non-hyperaccumulator for Pb in roots. In the present study, more than 98 millions reads were mapped to define gene structure and detect polymorphism, thereby to qualify transcript abundance along roots development under Pb treatment. A total of 17,707, 17,440, 16,998 and 16,586 genes were identified in maize roots at four developmental stages (0, 12 h, 24 h and 48 h) respectively and 2,825, 2,626, 2161 and 2260 stage-specifically expressed genes were also identified respectively. In addition, based on our RNA-Seq data, transcriptomic changes during maize root development responsive to Pb were investigated. A total of 384 differentially expressed genes (DEGs) (log2Ratio ≥ 1, FDR ≤ 0.001) were identified, of which, 36 genes with significant alteration in expression were detected in four developmental stages; 12 DEGs were randomly selected and successful validated by qRT-PCR. Additionally, many transcription factor families might act as the important regulators at different developmental stages, such as bZIP, ERF and GARP et al. These results will expand our understanding of the complex molecular and cellular events in maize root development and provide a foundation for future study on root development in maize under heavy metal pollution and other cereal crops.