To characterize type I and II collagen in the Amur sturgeon at the molecular level, mRNAs encoding the proα chain of both types of collagen were cloned and sequenced. Full sequences of both were obtained, and the molecular phylogeny based on the deduced amino acid sequence indicated that the correct sequences of the target genes were obtained. Analyses of primary structure of the proα chains revealed that type I and II collagen share the basic structure of the proα chain of fibril collagen, but have different characteristics, especially in residues related to thermal stability. In the triple helical domain, Gly-Pro-Pro sequence stabilizing the tripeptide unit was more frequent in type II than in type I, and Gly-Gly, which likely decline in thermal stability, was more frequent in type I than in type II. These results suggested that the denaturation temperature of type II would be remarkably higher than type I. The spatial pattern of gene expression was analyzed by quantitative real-time PCR, which showed that relatively ubiquitous type I gene and strongly skewed distribution of type II gene, which highly expressed only in vertebra, snout cartilage, and notochord. This pattern was similar to the distribution pattern of each collagen protein detected by previous biochemical analyses using Amur and Bester sturgeons. The present study is the first report of the cloning of the full-length cDNAs for both of type I and type II collagen in the Amur sturgeon, and is the first comparative analysis of type I and II collagens in a sturgeon species at the molecular level. The results provide basic and general information on collagens in sturgeons.

The analysis of transcript levels of specific genes is important for understanding transcriptional regulation and for the characterization of gene function. Real-time quantitative reverse transcriptase PCR (RT-qPCR) has become a powerful tool to quantify gene expression. The objective of this study was to identify reliable housekeeping genes in Giardia lamblia. Twelve genes were selected for this purpose, and their expression was analyzed in the wild type WB strain and in two strains with resistance to nitazoxanide (NTZ) and metronidazole (MTZ), respectively. RefFinder software analysis showed that the expression of the genes is different in the three strains. The integrated data from the four analyses showed that the NADH oxidase (NADH) and aldolase (ALD) genes were the most steadily expressed genes, whereas the glyceraldehyde-3-phosphate dehydrogenase gene was the most unstable. Additionally, the relative expression of seven genes were quantified in the NTZ- and MTZ-resistant strains by RT-qPCR, using the aldolase gene as the internal control, and the results showed a consistent differential pattern of expression in both strains. The housekeeping genes found in this work will facilitate the analysis of mRNA expression levels of other genes of interest in G. lamblia.

The goats are widely kept as livestock throughout the world. Two excellent domestic breeds in China, the Laiwu Black and Jining Grey goats, have different fecundities and prolificacies. Although the goat genome sequences have been resolved recently, little is known about the gene regulations at the transcriptional level in goat. To understand the molecular and genetic mechanisms related to the fecundities and prolificacies, we performed genome-wide sequencing of the mRNAs from two breeds of goat using the next-generation RNA-Seq technology and used functional annotation to identify pathways of interest. Digital gene expression analysis showed 338 genes were up-regulated in the Jining Grey goats and 404 were up-regulated in the Laiwu Black goats. Quantitative real-time PCR verified the reliability of the RNA-Seq data. This study suggests that multiple genes responsible for various biological functions and signaling pathways are differentially expressed in the two different goat breeds, and these genes might be involved in the regulation of goat fecundity and prolificacy. Taken together, our study provides insight into the transcriptional regulation in the ovaries of 2 species of goats that might serve as a key resource for understanding goat fecundity, prolificacy and genetic diversity between species.

Lung cancer is a malignant tumor with high mortality in both women and men. To study the mechanisms of smoking-induced lung cancer, we analyzed microarray of GSE4115. GSE4115 was downloaded from Gene Expression Omnibus including 78 and 85 bronchial epithelium tissue samples separately from smokers with and without lung cancer. Limma package in R was used to screen differentially expressed genes (DEGs). Hierarchical cluster analysis for DEGs was conducted using orange software and visualized by distance map. Using DAVID software, functional and pathway enrichment analyses separately were conducted for the DEGs. And protein-protein interaction (PPI) network was constructed using Cytoscape software. Then, the pathscores of enriched pathways were calculated. Besides, functional features were screened and optimized using the recursive feature elimination (RFE) method. Additionally, the support vector machine (SVM) method was used to train model. Total 1923 DEGs were identified between the two groups. Hierarchical cluster analysis indicated that there were differences in gene level between the two groups. And SVM analysis indicated that the five features had potential diagnostic value. Importantly, MAPK1 (degree=30), SRC (degree=29), SMAD4 (degree=23), EEF1A1 (degree=21), TRAF2 (degree=21) and PLCG1 (degree=20) had higher degrees in the PPI network of the DEGs. They might be involved in smoking-induced lung cancer by interacting with each other (e.g. MAPK1-SMAD4, SMAD4-EEF1A1 and SRC-PLCG1). MAPK1, SRC, SMAD4, EEF1A1, TRAF2 and PLCG1 might be responsible for the development of smoking-induced lung cancer.

Nitrogen use efficiency (NUE) in rice crop is the need of the hour for reduction of nitrous oxide emission resulting from excess nitrogen (N) fertilizer application and also in reduction of cost of cultivation. Ten rice genotypes were grown under low and recommended dose of N application and characterized in terms of parameters related to yield, yield related components and NUE indicators. Wide genetic variability under low N conditions was observed with significant variation for 15 yield related parameters in interactions of genotypes and treatment. Limitation of N has led to the decrease of all yield and yield related parameters, but for grain filling % and 1000 grain weight. Two genotypes, Rasi and Varadhan have shown minimum differences between low and recommended N conditions. Correlation analysis of various yield components showed the importance of the secondary branches for the total grains under low N. Expression analysis of OsSPL14 (LOC_Os08g39890) gene reported to be associated with increased panicle branching and higher grain yield through real time PCR in leaf and three stages of panicle has shown differential temporal expression and its association with yield and yield related components across the genotypes. The expression of OsSPL14 at panicle stage 3, has shown correlation (P<0.05) with N% in grain. Since OsSPL14 is a functional transcription activator, its association of expression in leaf and three panicle stages with yield components as observed in the present study suggests the role of nitrogen metabolism related genes in plant growth and development and its conversion into yield components in rice.

Carbonic anhydrase (CA) is involved in ion transport, acid-base balance and pH regulation by catalyzing the interconversion of CO2 and HCO3(-). In this study, full-length cDNA sequences of two CA isoforms were identified from Portunus trituberculatus. One was Portunus trituberculatus cytoplasmic carbonic anydrase (PtCAc) and the other one was Portunus trituberculatus glycosyl-phosphatidylinositol-linked carbonic anhydrase (PtCAg). The sequence of PtCAc was formed by an ORF of 816 bp, encoding a protein of 30.18 kDa. The PtCAg was constituted by an ORF of 927 bp, encoding a protein of 34.09 kDa. The deduced amino acid sequences of the two CA isoforms were compared to other crustacean' CA sequences. Both of them reflected high conservation of the residues and domains essential to the function of the two enzymes. The tissue expression analysis of PtCAc and PtCAg were detected in gill, muscle, hepatopancreas, hemocytes and gonad. PtCAc and PtCAg gene expressions were studied under salinity and pH challenge. The results showed that when salinity decreased (30 to 20 ppt), the mRNA expression of PtCAc increased significantly at 24 and 48 h, and the highest value appeared at 24h. The mRNA expression of PtCAg had the same situation with PtCAc. However, when salinity increased (30 to 35 ppt), only the mRNA expression of PtCAc increased significantly at 48 h. When pH changed, only the mRNA expression of PtCAc increased significantly at 12h, which was under low pH situation. The mRNA expression of PtCAg increased significantly at 12-48 h, and there was no significant difference of the expression between the pH challenged group and the control group in other experimental time. The results provided the base of understanding CA' function and the underlying mechanism in response to environmental changes in crustaceans.

Attempts to elucidate the cellular function of MTCH1 (mitochondrial carrier homolog 1) have not yet rendered a clear insight into the function of this outer mitochondrial membrane protein. Classical biochemical and cell biology approaches have not produced the expected outcome. In vitro experiments have indicated a likely role in the regulation of cell death by apoptosis, and its reported interaction with presenilin 1 suggests a role in the cellular pathways in which this membrane protease participates, nevertheless in vivo data are missing. In an attempt to identify cellular pathways in which this protein might participate, we have studied its promoter looking for transcriptional regulators. We have identified several putative binding sites for EGR-1 (Early growth response 1; a protein involved in growth, proliferation and differentiation), in the proximal region of the MTCH1 promoter. Chromatin immunoprecipitation showed an enrichment of these sequences in genomic DNA bound to EGR-1 and transient overexpression of EGR-1 in cultured HEK293T cells induces an increase of endogenous MTCH1 levels. We also show that MTCH1 levels increase in response to treatment of cells with doxorubicin, an apoptosis inducer through DNA damage. The endogenous levels of MTCH1 decrease when EGR-1 levels are lowered by RNA interference. Our results indicate that EGR-1 is a transcriptional regulator of MTCH1 and give some clues about the cellular processes in which MTCH1 might participate.

Eye infections due to Pseudomonas aeruginosa is an important cause of ocular morbidity. We presents the whole genomic comparative analysis of two P. aeruginosa VRFPA03 and VRFPA05 isolated from alkaline chemical injury mediated keratitis and post-cataract surgery endophthalmitis patients, respectively. The blaDIM-1 gene in VRFPA03 and the blaGes-9 gene in VRFPA05 were identified and reported for the first time from an ocular isolate. The current study revealed novel integrons In1107 and In1108, comprised of multidrug-resistant genes. Ocular virulence factors mainly mediated by exoenzymes T, Y, and U and exotoxin A, elastase B, and phenazine-specific methyltransferase. Genomic analysis uncovered multiple known and unknown factors involved in P. aeruginosa mediated ocular infection, which may lead to drug discovery and diagnostic markers to improve human vision care.

The monoamine oxidase A (MAOA) gene is an important candidate gene for human behavior that encodes an enzyme regulating the metabolism of key neurotransmitters. The regulatory mechanisms of the MAOA gene in dogs are yet to be elucidated. We measured MAOA gene transcription and analyzed the VNTR genotype and methylation status of the gene promoter region in different dog breeds to determine whether MAOA expression is correlated with the MAOA genotype or epigenetic modification in dogs. We found brain-specific expression of the MAOA gene and different transcription levels in different dog breeds including Beagle, Sapsaree, and German shepherd, and also a robust association of the DNA methylation of the gene promoter with mRNA levels. However, the 90 bp tandem repeats that we observed near the transcription start site were not variable, indicating no correlation with canine MAOA activity. These results show that differential DNA methylation in the MAOA promoter region may affect gene expression by modulating promoter activity. Moreover, the distinctive patterns of MAOA expression and DNA methylation may be involved in breed-specific or individual behavioral characteristics, such as aggression, because behavioral phenotypes are related to different physiological and neuroendocrine responses.

In the current study, an artificial neural network (ANN)-based breast cancer prediction model was developed from the data of folate and xenobiotic pathway genetic polymorphisms along with the nutritional and demographic variables to investigate how micronutrients modulate susceptibility to breast cancer. The developed ANN model explained 94.2% variability in breast cancer prediction. Fixed effect models of folate (400 μg/day) and B12 (6 μg/day) showed 33.3% and 11.3% risk reduction, respectively. Multifactor dimensionality reduction analysis showed the following interactions in responders to folate: RFC1 G80A × MTHFR C677T (primary), COMT H108L × CYP1A1 m2 (secondary), MTR A2756G (tertiary). The interactions among responders to B12 were RFC1G80A × cSHMT C1420T and CYP1A1 m2 × CYP1A1 m4. ANN simulations revealed that increased folate might restore ER and PR expression and reduce the promoter CpG island methylation of extra cellular superoxide dismutase and BRCA1. Dietary intake of folate appears to confer protection against breast cancer through its modulating effects on ER and PR expression and methylation of EC-SOD and BRCA1.

AMP-acetyl CoA synthetase (AMP-AceCS) is a key enzyme which catalyzes the activation of acetate to acetyl CoA, an important intermediate at the cross roads of various anabolic and catabolic pathways. Multiple sequence alignment of Leishmania donovani AceCS with other organisms revealed the presence of a highly conserved leucine residue at 684 position which is known to be crucial for acetylation by protein acetyl transferases in other organisms. In an attempt to understand the role of leucine residue at 684 position in L. donovani acetyl CoA synthetase (LdAceCS), it was mutated to proline (P) by site directed mutagenesis. Kinetic analysis of the L684P-LdAceCS mutant revealed approximately two fold increased binding affinity with acetate, whereas fivefold decreased affinity was observed with ATP. There was insignificant change in secondary structure as revealed by CD however, two fold decreased fluorescence intensity was observed at an emission maxima of 340 nm. Interestingly, L684P mutation abolished the acetylation of the mutant enzyme indicating the importance of L684 in acetylation of the enzyme. Changes in biochemical parameters of the mutant protein were validated by homology modeling of the wild type and mutant LdAceCS enzyme using Salmonella enterica AceCS crystal structure as template. Our data provides evidence for the role of leucine 684 residue in substrate recognition, catalysis and acetylation of the AceCS enzyme.

In contrast with the past, the water buffalo is now not only a draft animal, but also an important food source of milk and meat. It is increasingly apparent that the water buffalo have huge potential for meat production, but its breeding needs to be investigated. Regarding the molecular mechanisms involved in the meat quality difference between the buffalo (Bubalus bulabis) and yellow cattle (Bos taurus), 12 chemical-physical characteristics related to the meat quality of longissimus thoracis muscles (LTM) have been compared at the age of 36 months. Intramuscular lipid and b* (yellowness) were greater in cattle than the buffalo, whereas a* (redness) was greater in the buffalo. Gene expression profiles were constructed by bovine genome array. A total of 8884 and 10,960 probes were detected in buffalo and cattle, respectively, with 1580 genes being differentially expressed. Over 400 probes were upregulated and nearly 1200 were downregulated in LTM of the buffalo, most being involved in ribosomal RNA (rRNA) processing, cholesterol homeostasis, regulation of transcription, response to hypoxia, and glycolysis. Quantitative real-time PCR was used to validate the microarray data. Enriched GO analyses of highly expressed genes in LTM showed that protein biosynthesis, striated muscle contraction, iron homeostasis, iron transport, glycolysis and glucose metabolism were similar between the buffalo and cattle. High protein content, low fat content and deep meat color of buffalo LTM may be closely associated with the increased expression of genes involved in cholesterol and iron homeostasis, while also reducing the expression of genes involved in ubiquitin-mediated proteolysis and protein oxidative phosphorylation. These results establish the groundwork for further studies on buffalo meat quality and will be beneficial in improving water buffalo breeding by molecular biotechnology.

Galectins constitute a group of lectins with binding specificity for β-galactoside sugars. Galectin-1 is a prototype galectin and the multifunctionality of mammalian galectin-1s is well-known, but only a few of fish galectin-1s have been identified. In this study, we obtained the full-length cDNA and genomic sequence of the galectin-1 gene (designated as Pdlgals1) from large scale loach (Paramisgurnus dabryanus), performed phylogenetic analysis, and characterized the expression pattern and the transcriptional activity of its 5' flanking region. The Pdlgals1 gene contains 4 exons that encode a peptide of 132 amino acids with all the galectin signature motifs. Phylogenetic analysis and sequence alignment indicated that Pdlgals1 is a homologue of human LGALS1. RT-PCR and whole-mount in situ hybridization revealed that Pdlgals1 is mainly expressed in the skin, muscle, intestine and cavum oropharyngeum. Transcriptional activity assays demonstrated that the basal promoter of Pdlgals1 is located in a region from -500bp to its transcriptional start site. Potential binding sites for transcription factors including C/EBP, AP-1, GATA, Oct-1, δEF1, NF-κB, c-Myb, SP-1, AP-2, AML-1α, and AP-4 were identified in the basal promoter, suggesting that these factors are associated with the regulation of Pdlgals1. These results provided clues for further investigation of galectin-1 functions in loaches.

In present work, we described the mitochondrial genome (mitogenome) of the tea tussock moth Euproctis pseudoconspersa (Lepidoptera: Lymantriidae). The complete mitogenome of E. pseudoconspersa is a circular genome 15,461 bp in size. It contains 37 genes and an A+T-rich region usually presented in lepidopteran mitogenomes, which genes share a lot of features with other known lepidopteran mitogenomes. Nucleotide composition of A+T in this mitogenome is 79.92%, and the AT skew is slightly positive. Both codon distribution and relative synonymous codon usage of the 13 protein-coding genes (PCGs) are consistent with those published lepidopteran sequences. All tRNA genes have typical cloverleaf secondary structures, except for the tRNA(Ser(AGN)), in which the dihydrouridine (DHU) arm is simplified down to a loop. The A+T-rich region of E. pseudoconspersa mitogenome possess the motif 'ATAGA' and poly-T stretch as the formerly identified conserved elements of Lepidoptera mitogenomes. The phylogenetic relationships were reconstructed by using maximum likelihood (ML) and Bayesian inference (BI) methods based on nucleotide sequences of 13 PCGs of 38 moths. The results were very consistent with the traditional relationships within Noctuoidea from morphological data, and showed that Lymantriidae is more closely related to Erebidae than to Noctuidae.

Herein, Nile tilapia thioredoxin-interacting protein (On-TXNIP) and selenoprotein P (On-SEPP) cDNAs were cloned and characterized. The full-length On-TXNIP cDNA contained 2 arrestin domains, 2 conserved cysteine residues that bind to thioredoxin to inhibit thioredoxin function, and 2 PPXY motifs, which negatively regulate the protein by stimulating binding to E3 ubiquitin ligase. The On-SEPP cDNA contained 17 selenocysteines (Sec) encoded by the TGA codon, which can be recognized as either a stop codon or a Sec codon. The On-SEPP cDNA also carried 2 typical SECIS elements located in the 3'UTR that are important for selenocysteine translation. Evolutionary analyses of both the On-TXNIP and On-SEPP genes revealed that these genes are closely related to the TXNIP and SEPP genes in zebrafish (Danio rerio), with amino acid similarities of 91.8% and 61.9%, respectively. A normal tissue distribution analysis indicated that the On-TXNIP and On-SEPP genes were ubiquitously expressed in all tissues examined, and the highest expression levels of these genes were observed in peripheral blood leukocytes (PBLs) and the trunk kidney, respectively. The expression levels of On-TXNIP and On-SEPP transcripts were acutely and chronically analyzed following the injection of fish with 1, 10 or 100mg/kg silver nanoparticles (Ag NPs). Significant up-regulation of On-TXNIP and On-SEPP transcripts was observed in the liver, spleen, and head kidney at the early phase of Ag NP exposure (hours 6 through 48). Down-regulation of On-SEPP transcripts was clearly observed in the liver at weeks 1 to 4. Histopathology analysis demonstrated that the fish livers exhibited a dramatic infiltration of Kupffer cells, elevated bi-nucleated cells, expanded sinusoidal blood congestion and severe necrosis in a dose-dependent manner. Based on these findings, coupling of the expression analysis of these two cellular stress response genes and histopathological observation of fish exposed to Ag NPs should be reliable for the assessment of Ag NP contamination in teleost fish.

The present study was undertaken to characterize the genetic variation present in lymphoxin A gene (LTA gene) encoding for the lymphotoxin A protein also known as tumor necrosis factor beta, a cytokine produced by lymphocytes, known to be cytotoxic for a wide range of tumor cells both in vitro and in vivo, and, which is essential for normal immunological development; in 40 animals of 5 diverse Bos indicus Indian zebu cattle breeds. These breeds survive under the harsh and tough tropical climatic conditions of various parts of the Indian subcontinent. The LTA gene in the present study was observed to contain 33 SNPs and 3 small insertion/deletion polymorphisms. Four SNPs occurred in the coding regions of the gene viz. g.1327A>G and g.1400C>T in exon 2 and g.1840C>T and g.1942C>T in exon 3, of which the SNP g.1327A>G in exon 2 resulted in a non-synonymous amino acid change G38D. This amino acid change was however predicted not be affecting the protein function in any manner. The gene contained putative transcription factor binding sites for the c-Re1 and for Pax-4 transcription factors. A putative promoter region was also predicted on the reverse DNA strand from position 894 to 644. Several repeat elements and microsatellite repeats were detected to be occurring across the 3.2kb LTA gene sequence. The study showed the occurrence of 40 genotypes and 48 most probable haplotypes. The genotypes at the observed SNP positions in the LTA gene were in near Hardy-Weinberg equilibrium. A negative Tajima's D value that was not significant statistically at P>0.10 indicated that the neutral mutation hypothesis could not be excluded. The genetic variations observed in the LTA gene in the present study have not been reported earlier and these could possibly be used as molecular markers for further studies involving association of the gene variability with disease resistance/tolerance traits.

Advanced glycation end-products (AGEs) are a diverse group of molecules produced by the non-enzymatic addition of glucose to proteins, lipids, and nucleic acids. AGE levels have been associated with hyperglycemia and diabetic complications, especially in animal models, but less clearly in human studies. We measured total serum AGEs using an enzyme linked immunosorbant assay (ELISA) in 506 subjects from 246 families in the Diabetes Heart Study (DHS)/DHS MIND Study (n=399 type 2 diabetes (T2D)-affected). Single nucleotide polymorphisms (SNPs) in several candidate genes, including known AGE receptors, were tested for their influence on circulating AGE levels. The genetic analysis was expanded to include an exploratory genome-wide association study (GWAS) and exome chip analysis of AGEs (≈440,000 SNPs). AGEs were found to be highly heritable (h(2)=0.628, p=8.96 × 10(-10)). While no SNPs from candidate genes were significantly associated after Bonferroni correction, rs1035798 in the gene AGER was the most significantly associated (p=0.007). Additionally, rs7198427, in MT1A, showed a nominally significant p-value (p=0.0099). No SNPs from the GWAS or exome studies were identified after correction for multiple comparisons; however, rs17054480 in the PALLD2 gene on chromosome 4 showed the strongest association (p=7.77 × 10(-7)). Five SNPs at two loci (ISCA2/NPC2 and FBXO33) had p-values of less than 2.0 × 10(-5) and three additional SNPs (rs716326 in MACROD2, and rs6795197 and rs6765857 in ZBTB38) showed a nominal association with p-values of less than 1.0 × 10(-5).These findings provide a foundation for further investigation into the genetic component of circulating AGEs.

The crustacean hyperglycemic hormone (CHH) family is an important group of neuropeptides involved in controlling growth, reproduction, and stress response in decapod species. In this study, a new gene containing 4 exons-3 introns flanked by canonical 5'-GT-AG-3' intron splice-site junctions was isolated from Litopenaeus vannamei. Two full length transcripts of this CHH were isolated from eyestalk and pericardial tissue of males and females using rapid amplification of cDNA ends (RACE). Transcripts sequences were 1578bp in length in males pericardial tissues and in males and females eyestalk with 100% identity, but the transcript isolated from females pericardial tissues was shorter (974bp). The differences in transcripts length is a result of two polyadenylation sites present in the 3'UTR resulting in two transcription termination signals. Transcript sequences encoded one unique protein that can be classified as type I CHH subfamily because of the 4 exons and 3 introns structure, although the CPRP region is not-well conserved and there is no amidation in the C-terminal of the deduced amino acid sequence. Furthermore, there is a glycine inserted in the mature peptide not at position 12 as in type II CHHs but after amino acid 31 and the phylogenetic analysis did not group the peptide within type I, but closer to type II CHHs. We demonstrated by endpoint-PCR, qPCR, and in situ hybridization (ISH), that this gene is expressed in neuroendocrine organs known to express CHHs in penaeid shrimp, including X-organ and optic nerve in eyestalk, supraesophageal ganglion (SoG), but it is also expressed in other organs as gill, gut, pericardial cavity, as well as in terminal ampoule or spermatophore and vas deferens of males.

Since October 2010, porcine epidemic diarrhea (PED) caused by variant porcine epidemic diarrhea virus (PEDV) has led great economic losses to the global pig industry, especially in China. To study the genetic characteristics of PEDV strains in Chinese mainland, a total of 603 clinical samples from nine provinces/districts of Chinese mainland from January 2014 to December 2015 were collected for RT-PCR detection and 1-1323bp of S gene of 91 isolates and ORF3 gene of 46 isolates were sequenced. The results showed that the variant PEDV were the dominant pathogens of viral diarrhea diseases in these areas. Six novel variant PEDV strains (FJAX1, FJAX2, HeNPDS1, HeNPDS2, HeNPY3, and HeNPY4) with two amino acids (aa) deletion at the 56-57 aa of S protein were identified. A total of 405 Chinese PEDV strains were subjected to phylogenetic and phylogeographic analysis. The results revealed that the subgroup Va in variant PEDV group were the dominant subgroup and the spread trend of variant PEDV strains seemed to be from the southeast coastal districts to other coastal districts and interior districts. The N-terminal of S gene (1-750bp), to some extent, could represent S1 or full length S gene for phylogenetic, similarity, antigen index, hydrophilicity plot, and differentiation analyses. The 404-472bp of S gene contained the three genetic markers, i.e., "TAA" insertion at 404-405bp, "ACAGGT" deletion at 430-435bp, and "ATA" deletion at 455-457bp can be used to differentiate the classical and variant virulent parental/attenuated PEDV strains and help us to learn the infectious and genetic characteristics of PEDV strains more convenient and cheaper. This study has important implication for understanding the infectious, genetic, and evolutionary aspects of PEDV strains in Chinese mainland.

To establish a three-step programmed method to find gene mutations related to maturity onset diabetes of the young (MODY). Target region capture and next-generation sequencing (NGS) were performed using customized oligonucleotide probes designed to capture suspected genes for MODY in 11 probands with clinically diagnosed MODY. The suspected associations of certain genes with MODY were then confirmed by Sanger sequencing in the probands and their family members. Finally, to validate variants of one of the genes of interest (glucokinase, GCK) as pathogenic mutations, protein function editing by the variant genes was assessed. In the target region capture and NGS phase, a total of nine variants of seven genes (GCK, WFS1, SLC19A2, SH2B1, SERPINB4, RFX6, and GATA6) were identified in eight probands. Two heterozygous GCK mutations located on the same allele (p.Leu77Arg and p.Val101Met) were identified in a MODY family. Sanger sequencing was used to confirm the variants identified by NGS to be present in probands and their diabetic family members, but not in non-diabetic family members. Finally, enzyme kinetic and thermal stability analyses revealed that the p.Leu77Arg mutation or the p.Leu77Arg mutation in combination with the p.Val101Met mutation inactivates GCK function and stability, while mutation of p.Val101Met alone does not. The p.Leu77Arg but not p.Val101Met GCK mutation is therefore considered a pathogenic mutation associated with MODY. Genetic screening coupled with gene-editing protein function testing is an effective and reliable method by which causative gene mutations of MODY can be identified.