A Rapid PCR-RFLP Method for Assessing Heterozygosity of Murraya paniculata Germplasm / 中国实验方剂学杂志

ObjectiveTo establish a rapid method for evaluating the heterozygosity of Murraya paniculata germplasm materials and provide as a foundation for developing germplasm breeding and innovation measures for M. paniculata. MethodSingle nucleotide polymorphisms （SNPs） were screened from the genome resequencing data of 65 plants of M. paniculata. A self-written script was used to transform 20 SNPs into restriction fragment length polymorphism （RFLP） markers. Polymerase chain reaction-restriction fragment length polymorphism （PCR-RFLP） was employed to detect the 20 RFLP markers in 12 M. paniculata germplasm accessions， and the heterozygosity of M. paniculata germplasm accessions was calculated based on the number of enzyme-cutting bands at the 20 RFLP marker sites. Plink was used to calculate the whole genome heterozygosity of 12 M. paniculata germplasm accessions， and the results obtained with different methods were compared. ResultThere was no significant difference in the heterozygosity calculated by the PCR-RFLP method and the genome resequencing method. The PCR-RFLP and genome resequencing methods identified 8 and 9 germplasm accessions， respectively， with a heterozygosity level less than 30%. Seven germplasm accessions with heterozygosity less than 30.00% were calculated by both methods. ConclusionThe PCR-RFLP method established in this study for evaluating the heterozygosity of M. paniculata germplasm demonstrates the precision of 87.5% and the accuracy of 77.8%. This method serves as a reference for developing heterozygosity evaluation methods in other medicinal plant germplasm resources.

Molecular characteristics of clinically isolated Yersinia in Jiangsu Province from 2017 to 2021 / 中华地方病学杂志

Objective:To investigate the distribution and molecular characteristics of Yersinia isolated from diarrhea patients in Jiangsu Province. Methods:From 2017 to 2021, the stool samples of diarrhea patients were collected in Tongshan District of Xuzhou City and Dongtai City of Yancheng City, Jiangsu Province, where the national active monitoring sites of Yersinia enterocolitica, then Yersinia was isolated; meanwhile, suspected Yersinia strains were collected from sentinel hospitals in the province. The DNA of isolated strains was extracted for whole genome resequencing, and the data were uploaded to the EnteroBase database for Yersinia species identification; the original data were cleaned and processed for 16S ribosomal RNA (16S rRNA) gene polymorphism analysis. Five virulence genes (ail, ystA, ystB, yadA, virF) were scanned through the National Center for Biotechnology Information (NCBI) and Pathogen Virulence Factor Database (VFDB), and K-mer Tree was constructed and genomic characteristics were analyzed. Results:From 2017 to 2021, a total of 2 058 stool samples from diarrhea patients were collected, and 57 strains of Yersinia were isolated and identified; meanwhile, two Yersinia strains were collected from the sentinel hospital. Compared with EnteroBase database, 51 strains were identified as Yersinia enterocolitica, 4 strains as Yersinia proxima, 1 strain each as Yersinia aleksiciae, Yersinia massiliensis, Yersinia intermedia and Yersinia canariae. The 16S rRNA gene polymorphism analysis showed that all strains were clustered into 3 groups, which could distinguish Yersinia enterocolitica from other Yersinia. Among the 51 strains of Yersinia enterocolitica, 49 strains were virulence genotype Ⅲ(ail-, ystA-, ystB+, yadA-, virF-), two strains were virulence genotype Ⅱ(ail+, ystA+, ystB-, yadA-, virF-); and 8 other Yersinia strains were virulence genotype Ⅳ (ail-, ystA-, ystB-, yadA-, virF-). K-mer analysis could distinguish Yersinia enterocolitica from other Yersinia, JS-XZ-2020001 strain was far away from other Yersinia enterocolitica isolates, and serotype O8 strains were more concentrated. Conclusions:The clinical isolates of Yersinia enterocolitica from diarrhea patients are mainly Yersinia and other Yersinia co-exist in a small amount in Jiangsu Province, two new Yersinia species ( Yersinia proxima and Yersinia canariae) are discovered. The virulence genotype of Yersinia enterocolitica is mainly type Ⅲ. The 16S rRNA gene polymorphism analysis and K-mer analysis can effectively distinguish Yersinia enterocolitica from other Yersinia.

Genetic diversity of protopine-6-hydroxylase in three medicinal Papaver plants / 中国中药杂志

Sanguinarine is the main active component of the Papaver plants, and protopine-6-hydroxylase(P6 H), involved in the sanguinarine biosynthetic pathway, can oxidize protopine to 6-hydroxyprotopine. The investigation on the diversity of P6 H genes in the medicinal Papaver plants contributes to the acquirement of P6 H with high activity to increase the biosynthesis of sanguinarine. Five P6 H genes in P. somniferum, P. orientale, and P. rhoeas were discovered based on the re-sequencing data of the Papaver species, followed by bioinformatics analysis. With the elongation factor 1α(EF-1α), which exhibits stable expression in the root and stem, as the internal reference gene, the transcription levels of P6H genes in roots and stems of the Papaver plants were detected by real-time fluorescent quantitative PCR. As indicated by the re-sequencing results, there were two genotypes of P6H in P. somniferum and P. orientale, respectively, and only one in P. rhoeas. The bioinformatics analysis showed that the P6 H proteins of the three Papaver plants contained the conserved domain cl12078, which is the characteristic of p450 supergene family, and transmembrane regions. The existence of signal peptide remained verification. Real-time fluorescent quantitative PCR results revealed that the transcription level of P6 H in roots of P. somniferum was about 1.44 times of that in stems(α=0.05). The present study confirmed genetic diversity of P6 H in the three medicinal Papaver plants, which lays a basis for the research on the biosynthesis pathway and mechanism of sanguinarine in Papaver species.

Whole-genome resequencing reveals genetic indels of feathered-leg traits in domestic chickens

Whole-genome resequencing provides the opportunity to explore the genomic variations and pave way for further functional assays to map the economical trait loci. In this study, we sequenced the genomes of mixed chicken samples from a full-sib family, with feathered and unfeathered legs at an average effective depth of 4.43×, using Illumina Hiseq 2000 instruments. Over 2.1 million nonredundant short indels (1–71 bp) were obtained. Among them, 16,375 common indels that were polymorphic between the comparison groups were revealed for further analysis. The majority of the common differential indels (76.52%) were novel. Follow-up validation assays confirmed that 80% randomly selected indels represented true variations. The indels were annotated based on the chicken genome sequence assembly. As a result, 16,375 indels were found to be located within 2756 annotated genes, with only 33 (0.202%) located in exons. By integrated analysis of the 2756 genes with gene function and known quantitative trait loci, we identified a total of 24 promising candidate genes potentially affecting feathered-leg trait, i.e. FGF1, FGF4, FGF10, FGFR1, FRZB, WNT1, WNT3A, WNT11, PCDH1, PCDH10, PCDH19, SOX3, BMP2, NOTCH2, TGF-β2, DLX5, REPS2, SCN3B, TCF20, FGF3, FSTL1, WNT7B, ELOVL2 and FGF8. Our findings provide a basis for further study and reveal key genes for feathered-leg trait in chickens.

Genomics: cracking the mysteries of walnuts

The Juglans plants are economically important as they provide nuts, wood and garden trees. They also play an important ecological role by supplying food for wild insects and animals. The decoding of genome sequences has fundamental values for understanding the evolution of Juglans plants and molecules, and is also a prerequisite for molecular breeding. During the last three years, the rapid development of sequencing technology has made walnut research into the genome era. Here, we reviewed the progress of genome sequencing of six Juglans species, the resequencing of four Juglans populations as well as the genome sequencing of the closely related species Pterocarya stenoptera. The analysis of the Juglans regia genome uncovers a whole genome duplication (WGD) event. Based on the molecular dating of the divergence time of six Juglans species, we proposed this WGD event was associated with Cretaceous–Palaeogene (K-Pg) boundary occurred ∼65 million years ago. Genomic sequences also provide clear details for understanding the evolution and development of GGT and PPO genes involved in fruit development. The decoding of these genomes has made it easier for us to understand and enhance the use of walnuts. We expect that the functional genomics research of walnut will also develop rapidly in the near future.

DNA-pool high-throughput whole genome resequencing for exploring essential hypertension single nucleotide polymorphism mutation / 国际检验医学杂志

Objective To use the DNA-pool technology to sequence patients with essential hypertension(EH) for exploring the single nucleotide polymorphism(SNP) mutation situation in Chinese patients with EH.Methods One hundred EH outpatients in the Shenzhen Sun Yat-sen Cardiovascular Hospital from March to June 2014 were continuously collected.The genomic DNA was performed the fragmentation process to 400-800 bp for conducting the database creation and sequencing.The sequencing results were compared with hg19 in the human gene bank(National Center of Biotechnology Information).Results A total of 120.8 Gb original sequence data were generated.The sequencing depth was 36.13 times,the coverage rate reached 99.88%.A total of 4 305 668 SNP loci were detected by the bioinformatic analysis,in which the C:G→T:A motation types were miximal,reaching 12 314 variation sites.Conclusion This study verifies that the data obtained by using the DNA-pool whole genome resequencing method replenishes the Chinese gene database of EH and provides some help for EH gene reasearch in the future.

Genomic profiling in a homogeneous molecular subtype of non‑small cell lung cancer: An effort to explore new drug targets.

BACKGROUND: Patients’ who are positive for kinase domain activating mutations in epidermal growth factor receptor (EGFR) gene, constitute 30–40% of non‑small cell lung cancer (NSCLC), and are suitable candidates for Tyrosine Kinase Inhibitor based targeted/personalized therapy. In EGFR non‑mutated subset, 8–10% that show molecular abnormalities such as EML4‑ALK, ROS1‑ALK, KIP4‑ALK, may also derive the benefit of targeted therapy. However, 40% of NSCLC belong to a grey zone of tumours that are negative for the clinically approved biomarkers for personalized therapy. This pilot study aims to identify and classify molecular subtypes of this group to address the un‑met need for new drug targets in this category. Here we screened for known/novel oncogenic driver mutations using a 46 gene Ampliseq Panel V1.0 that includes Ser/Thr/ Tyr kinases, transcription factors and tumor suppressors. METHODS: NSCLC with tumor burden of at least 40% on histopathology were screened for 29 somatic mutations in the EGFR kinase domain by real‑time polymerase chain reaction methods. 20 cases which were EGFR non‑mutated for TK domain mutations were included in this study. DNA Quality was verified from each of the 20 cases by fluorimeter, pooled and subjected to targeted re‑sequencing in the Ion Torrent platform. Torrent Suite software was used for next generation sequencing raw data processing and variant calling. RESULTS: The clinical relevance and pathological role of all the mutations/variants that include SNPs and Indels was assessed using polyphen‑2/SIFT/PROVEAN/mutation assessor structure function prediction programs. There were 10 pathogenic mutations in six different oncogenes for which annotation was available in the COSMIC database; C420R mutation in PIK3CA, Q472H mutation in vascular endothelial growth factor receptor 2 (VEGFR2) (KDR), C630W and C634R in RET, K367M mutation in fibroblast growth factor receptor 2 (FGFR2), G12C in KRAS and 4 pathogenic mutations in TP53 in the DNA binding domain (E285K, R213L, R175H, V173G). CONCLUSION: Results suggest, a potential role for PIK3CA, VEGFR2, RET and FGFR2 as therapeutic targets in EGFR non‑mutated NSCLC that requires further clinical validation.

Prediction of Genes Related to Positive Selection Using Whole-Genome Resequencing in Three Commercial Pig Breeds

Selective sweep can cause genetic differentiation across populations, which allows for the identification of possible causative regions/genes underlying important traits. The pig has experienced a long history of allele frequency changes through artificial selection in the domestication process. We obtained an average of 329,482,871 sequence reads for 24 pigs from three pig breeds: Yorkshire (n = 5), Landrace (n = 13), and Duroc (n = 6). An average read depth of 11.7 was obtained using whole-genome resequencing on an Illumina HiSeq2000 platform. In this study, cross-population extended haplotype homozygosity and cross-population composite likelihood ratio tests were implemented to detect genes experiencing positive selection for the genome-wide resequencing data generated from three commercial pig breeds. In our results, 26, 7, and 14 genes from Yorkshire, Landrace, and Duroc, respectively were detected by two kinds of statistical tests. Significant evidence for positive selection was identified on genes ST6GALNAC2 and EPHX1 in Yorkshire, PARK2 in Landrace, and BMP6, SLA-DQA1, and PRKG1 in Duroc.These genes are reportedly relevant to lactation, reproduction, meat quality, and growth traits. To understand how these single nucleotide polymorphisms (SNPs) related positive selection affect protein function, we analyzed the effect of non-synonymous SNPs. Three SNPs (rs324509622, rs80931851, and rs80937718) in the SLA-DQA1 gene were significant in the enrichment tests, indicating strong evidence for positive selection in Duroc. Our analyses identified genes under positive selection for lactation, reproduction, and meat-quality and growth traits in Yorkshire, Landrace, and Duroc, respectively.

Strategy for the Customized Mass Screening of Genetic Sensorineural Hearing Loss in Koreans

Hearing loss is one of the most common sensorineural disorder. More than half of congenital bilateral profound deafness cases have been estimated to be attributed to genetic cause. Identification of genetic cause can provide valuable information. We developed new diagnostic strategy combining phenotype-driven candidate gene approach and targeted exome sequencing to find out the causative mutation of hearing loss. The causative mutation detection rates of this strategy were 78.1% and 54.8% in Korean multiplex families and sporadic severe to profound hearing loss families, respectively. The most frequent causative genes of Korean multiplex families were SLC26A4 and POU3F4. The other causative genes were MRNR1, WFS1, COCH, TECTA, MYO6, COL11A2, EYA4, GJB3, OTOF, STRC, MYO3A, and GJB2. The most frequent causative gene of Korean sporadic severe to profound hearing loss families was SLC26A4 followed by GJB2, CHD7, and CDH23. Based upon the results, the value of this strategy as a diagnostic tool seems to be promising. Although whole genome and exome sequencing have advanced as the development of next-generation sequencing, this new strategy could be a good screening and diagnostic tool to find the causative mutations.

Survey of the Applications of NGS to Whole-Genome Sequencing and Expression Profiling

Recently, the technologies of DNA sequence variation and gene expression profiling have been used widely as approaches in the expertise of genome biology and genetics. The application to genome study has been particularly developed with the introduction of the next-generation DNA sequencer (NGS) Roche/454 and Illumina/Solexa systems, along with bioinformation analysis technologies of whole-genome de novo assembly, expression profiling, DNA variation discovery, and genotyping. Both massive whole-genome shotgun paired-end sequencing and mate paired-end sequencing data are important steps for constructing de novo assembly of novel genome sequencing data. It is necessary to have DNA sequence information from a multiplatform NGS with at least 2x and 30x depth sequence of genome coverage using Roche/454 and Illumina/Solexa, respectively, for effective an way of de novo assembly. Massive short-length reading data from the Illumina/Solexa system is enough to discover DNA variation, resulting in reducing the cost of DNA sequencing. Whole-genome expression profile data are useful to approach genome system biology with quantification of expressed RNAs from a whole-genome transcriptome, depending on the tissue samples. The hybrid mRNA sequences from Rohce/454 and Illumina/Solexa are more powerful to find novel genes through de novo assembly in any whole-genome sequenced species. The 20x and 50x coverage of the estimated transcriptome sequences using Roche/454 and Illumina/Solexa, respectively, is effective to create novel expressed reference sequences. However, only an average 30x coverage of a transcriptome with short read sequences of Illumina/Solexa is enough to check expression quantification, compared to the reference expressed sequence tag sequence.

The Importance of Global Studies of the Genetics of Type 2 Diabetes

Genome wide association analyses have revealed large numbers of common variants influencing predisposition to type 2 diabetes and related phenotypes. These studies have predominantly featured European populations, but are now being extended to samples from a wider range of ethnic groups. The transethnic analysis of association data is already providing insights into the genetic, molecular and biological causes of diabetes, and the relevance of such studies will increase as human discovery genetics increasingly moves towards sequencing-based approaches and a focus on low frequency and rare variants.

Identification of Causal and/or Rare Genetic Variants for Complex Traits by Targeted Resequencing in Population-based Cohorts

Genome-wide association studies (GWASs) have greatly contributed to the identification of common variants responsible for numerous complex traits. There are, however, unavoidable limitations in detecting causal and/or rare variants for traits in this approach, which depends on an LD-based tagging SNP microarray chip. In an effort to detect potential casual and/or rare variants for complex traits, such as type 2 diabetes (T2D) and triglycerides (TGs), we conducted a targeted resequencing of loci identified by the Korea Association REsource (KARE) GWAS. The target regions for resequencing comprised whole exons, exon-intron boundaries, and regulatory regions of genes that appeared within 1 Mb of the GWA signal boundary. From 124 individuals selected in population-based cohorts, a total of 0.7 Mb target regions were captured by the NimbleGen sequence capture 385K array. Subsequent sequencing, carried out by the Roche 454 Genome Sequencer FLX, generated about 110,000 sequence reads per individual. Mapping of sequence reads to the human reference genome was performed using the SSAHA2 program. An average of 62.2% of total reads was mapped to targets with an average 22X-fold coverage. A total of 5,983 SNPs (average 846 SNPs per individual) were called and annotated by GATK software, with 96.5% accuracy that was estimated by comparison with Affymetrix 5.0 genotyped data in identical individuals. About 51% of total SNPs were singletons that can be considered possible rare variants in the population. Among SNPs that appeared in exons, which occupies about 20% of total SNPs, 304 nonsynonymous singletons were tested with Polyphen to predict the protein damage caused by mutation. In total, we were able to detect 9 and 6 potentially functional rare SNPs for T2D and triglycerides, respectively, evoking a further step of replication genotyping in independent populations to prove their bona fide relevance to traits.

Sequence Variations and Haplotypes of the GJB2 Gene Revealed by Resequencing of 192 Chromosomes from the General Population in Korea

OBJECTIVES: Hearing impairment (HI) is the most common sensory deficit in human. The Gap Junction Protein, Beta-2 (GJB2) gene encodes the protein connexin 26, and this gene accounts for up to half of the cases of autosomal recessive nonsyndromic HI. This study was conducted to obtain a set of sequence variations (SVs) of the GJB2 gene among Koreans from the general population for making molecular genetic diagnoses and performing genetic counseling. METHODS: We resequenced the GJB2 gene in 192 chromosomes from 96 adult individuals of Korean descent and who were without a history of hearing difficulty. The data of the SVs was obtained and the haplotypes were reconstructed from the data. RESULTS: Five SVs were observed, including a novel one (c.558G>A; p.T186T), with the allele frequencies ranging from 0.5% (1/192) to 41% (79/192). The linkage disequilibrium study and haplotype construction showed that some of the SVs are in tight linkage, resulting in a limited number of haplotypes. CONCLUSION: We observed SVs of the GJB2 gene with different allele frequencies, and a limited number of haplotypes were constructed. The data from this study can be used as reference data for GJB2-related hearing genetic studies, including studies on the founder effect and population genetics, and this data is particularly relevant to people of East Asian decent.

Application of re-sequencing microarray in the pathogenic detection of the firstly confirmed human case of avian influenza A in Beijing / 中国人兽共患病学报

To screen and identify the possible pathogen of the firstly confirmed human case of avian influenza A in Beijing, the throat swabs and tracheal aspirates of this case were collected and the H5N1 viral nucleotide was tested with real time RT-PCR. The certification of result, screening of other pathogens in respiratory tract and sub-typing of influenza viruses were made by using re-sequencing microarray. It was found that the H5N1 viral nucleic acid was positive in the tracheal aspirate of this case by means of detection with real time RT-PCR and the specific sequence of the non-structural protein (NS) gene of H5N1 virus was obtained through the detection with re-sequencing clip. Through the comparative study with the sequence in Genbank, it was proved to be the H5N1 nucleic acid of avian influenza viruses and excluded the possibility of infections with 30 subtypes of influenza viruses and 33 other respiratory tract pathogens. It is apparent that the pathogen detection with re-sequencing clip shows the high sensitivity and specificity and it plays an important role in the pathogen screening and identification for the firstly confirmed human case of avian influenza A in Beijing.

An Optimized Strategy for Genome Assembly of Sanger/pyrosequencing Hybrid Data using Available Software

During the last four years, the pyrosequencing-based 454 platform has rapidly displaced the traditional Sanger sequencing method due to its high throughput and cost effectiveness. Meanwhile, the Sanger sequencing methodology still provides the longest reads, and paired-end sequencing that is based on that chemistry offers an opportunity to ensure accurate assembly results. In this report, we describe an optimized approach for hybrid de novo genome assembly using pyrosequencing data and varying amounts of Sanger-type reads. 454 platformderived contigs can be used as single non-breakable virtual reads or converted to simpler contigs that consist of editable, overlapping pseudoreads. These modified contigs maintain their integrity at the first jumpstarting assembly stage and are edited by fragmenting and rejoining. Pre-existing assembly software then can be applied for mixed assembly with 454-derived data and Sanger reads. An effective method for identifying genomic differences between reference and sample sequences in whole-genome resequencing procedures also is suggested.