Whole-genome analysis of Escherichia coli isolated from wild Amur tiger (Panthera tigris altaica) and North China leopard (Panthera pardus japonensis).

Background: Escherichia coli is an important intestinal flora, of which pathogenic E. coli is capable of causing many enteric and extra-intestinal diseases. Antibiotics are essential for the treatment of bacterial infections caused by pathogenic E. coli; however, with the widespread use of antibiotics, drug resistance in E. coli has become particularly serious, posing a global threat to human, animal, and environmental health. While the drug resistance and pathogenicity of E. coli carried by tigers and leopards in captivity have been studied intensively in recent years, there is an extreme lack of information on E. coli in these top predators in the wild environment. Methods: Whole genome sequencing data of 32 E. coli strains collected from the feces of wild Amur tiger (Panthera tigris altaica, n = 24) and North China leopard (Panthera pardus japonensis, n = 8) were analyzed in this article. The multi-locus sequence types, serotypes, virulence and resistance genotypes, plasmid replicon types, and core genomic SNPs phylogeny of these isolates were studied. Additionally, antimicrobial susceptibility testing (AST) was performed on these E. coli isolates. Results: Among the E. coli isolates studied, 18 different sequence types were identified, with ST939 (21.9%), ST10 (15.6%), and ST3246 (9.4%) being the most prevalent. A total of 111 virulence genes were detected, averaging about 54 virulence genes per sample. They contribute to invasion, adherence, immune evasion, efflux pump, toxin, motility, stress adaption, and other virulence-related functions of E. coli. Sixty-eight AMR genes and point mutations were identified. Among the detected resistance genes, those belonging to the efflux pump family were the most abundant. Thirty-two E. coli isolates showed the highest rate of resistance to tetracycline (14/32; 43.8%), followed by imipenem (4/32; 12.5%), ciprofloxacin (3/32; 9.4%), doxycycline (2/32; 6.3%), and norfloxacin (1/32; 3.1%). Conclusions: Our results suggest that E. coli isolates carried by wild Amur tigers and North China leopards have potential pathogenicity and drug resistance.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

The effect of reflux ratio on sulfur disproportionation tendency in anaerobic baffled reactor with the heterotrophic combining sulfur autotrophic processes under high concentration perchlorate stress.

In a laboratory scale, an anaerobic baffled reactor (ABR) consisting of eight compartments, the heterotrophic combining sulfur autotrophic processes under different reflux ratios were constructed to achieve effective perchlorate removal and alleviate sulfur disproportionation reaction. Perchlorate was efficiently removed with effluent perchlorate concentration below 0.5 µg/L when the influent perchlorate concentration was 1030 mg/L during stages I ~ V, indicating that heterotrophic combining sulfur autotrophic perchlorate reduction processes can effectively achieve high concentration perchlorate removal. Furthermore, the 100% reflux ratio could reduce the contact time between sulfur particles and water; thus, the sulfur disproportionation reaction was inhibited. However, the inhibition effect of reflux on sulfur disproportionation was attenuated due to dilute perchlorate concentration when a reflux ratio of 150% and 200% was implemented. Meanwhile, the content of extracellular polymeric substances (EPS) in the heterotrophic unit (36.79 ~ 45.71 mg/g VSS) was higher than that in the sulfur autotrophic unit (22.19 ~ 25.77 mg/g VSS), indicating that high concentration perchlorate stress in the heterotrophic unit promoted EPS secretion. Thereinto, the PN content of sulfur autotrophic unit decreased in stage III and stage V due to decreasing perchlorate concentration in the autotrophic unit. Meanwhile, the PS content increased with increasing reflux in the autotrophic unit, which was conducive to the formation of biofilm. Furthermore, the high-throughput sequencing result showed that Proteobacteria, Chloroflexi, Firmicutes, and Bacteroidetes were the dominant phyla and Longilinea, Diaphorobacter, Acinetobacter, and Nitrobacter were the dominant genus in ABR, which were associated with heterotrophic or autotrophic perchlorate reduction and beneficial for effective perchlorate removal. The study indicated that reflux was a reasonable strategy for alleviating sulfur disproportionation in heterotrophic combining sulfur autotrophic perchlorate removal processes.

Y and Mitochondrial Chromosomes in the Heterogeneous Stock Rat Population.

Genome-wide association studies typically evaluate the autosomes and sometimes the X Chromosome, but seldom consider the Y or mitochondrial Chromosomes. We genotyped the Y and mitochondrial chromosomes in heterogeneous stock rats (Rattus norvegicus), which were created in 1984 by intercrossing eight inbred strains and have subsequently been maintained as an outbred population for 100 generations. As the Y and mitochondrial Chromosomes do not recombine, we determined which founder had contributed these chromosomes for each rat, and then performed association analysis for all complex traits (n=12,055; intersection of 12,116 phenotyped and 15,042 haplotyped rats). We found the eight founders had 8 distinct Y and 4 distinct mitochondrial Chromosomes, however only two of each were observed in our modern heterogeneous stock rat population (Generations 81-97). Despite the unusually large sample size, the p-value distribution did not deviate from expectations; there were no significant associations for behavioral, physiological, metabolome, or microbiome traits after correcting for multiple comparisons. However, both Y and mitochondrial Chromosomes were strongly associated with expression of a few genes located on those chromosomes, which provided a positive control. Our results suggest that within modern heterogeneous stock rats there are no Y and mitochondrial Chromosomes differences that strongly influence behavioral or physiological traits. These results do not address other ancestral Y and mitochondrial Chromosomes that do not appear in modern heterogeneous stock rats, nor do they address effects that may exist in other rat populations, or in other species.

Exploration of Parascaris species in three different Equus populations in China.

BACKGROUND: The roundworms, Parascaris spp., are important nematode parasites of foals and were historically model organisms in the field of cell biology, leading to many important discoveries. According to karyotype, ascarids in Equus are commonly divided into Parascaris univalens (2n = 2) and Parascaris equorum (2n = 4). METHODS: Here, we performed morphological identification, karyotyping and sequencing of roundworms from three different hosts (horses, zebras and donkeys). Phylogenetic analysis was performed to study the divergence of these ascarids based on cytochrome c oxidase subunit I (COI) and internal transcribed spacer (ITS) sequences. RESULTS: Karyotyping, performed on eggs recovered from worms of three different Equus hosts in China, showed two different karyotypes (2n = 2 in P. univalens collected from horses and zebras; 2n = 6 in Parascaris sp. collected from donkeys). There are some differences in the terminal part of the spicula between P. univalens (concave) and Parascaris sp. (rounded). Additionally, it was found that the egg's chitinous layer was significantly thicker in Parascaris sp. (> 5 µm) than P. univalens (< 5 µm) (F(2537) = 1967, P < 0.01). Phylogenetic trees showed that the sequences of Parascaris from Equus hosts were divided into two distinct lineages based on sequences of the COI and ITS. CONCLUSIONS: Comparing the differences in roundworms collected from three different Equus hosts, this study describes a Parascaris species (Parascaris sp.) with six chromosomes in donkeys. It is worth noting that the thickness of the chitinous layer in the Parascaris egg may serve as a diagnostic indicator to distinguish the two roundworms (P. univalens and Parascaris sp.). The Parascaris sp. with six chromosomes in donkeys in the present study may be a species of P. trivalens described in 1934, but the possibility that it is a new Parascaris species cannot be ruled out. Both karyotyping and molecular analysis are necessary to solve the taxonomic problems in Parascaris species.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Genome-wide association study of delay discounting in Heterogenous Stock rats.

Delay discounting refers to the behavioral tendency to devalue rewards as a function of their delay in receipt. Heightened delay discounting has been associated with substance use disorders, as well as multiple co-occurring psychopathologies. Genetic studies in humans and animal models have established that delay discounting is a heritable trait, but only a few specific genes have been associated with delay discounting. Here, we aimed to identify novel genetic loci associated with delay discounting through a genome-wide association study (GWAS) using Heterogenous Stock rats, a genetically diverse outbred population derived from eight inbred founder strains. We assessed delay discounting in 650 male and female rats using an adjusting amount procedure in which rats chose between smaller immediate sucrose rewards or a larger reward at variable delays. Preference switch points were calculated for each rat and both exponential and hyperbolic functions were fitted to these indifference points. Area under the curve (AUC) and the discounting parameter k of both functions were used as delay discounting measures. GWAS for AUC, exponential k, and indifference points for a short delay identified significant loci on chromosomes 20 and 14. The gene Slc35f1, which encodes a member of the solute carrier family of nucleoside sugar transporters, was the only gene within the chromosome 20 locus. That locus also contained an eQTL for Slc35f1, suggesting that heritable differences in the expression of that gene might be responsible for the association with behavior. The gene Adgrl3, which encodes a member of the latrophilin family of G-protein coupled receptors, was the only gene within the chromosome 14 locus. These findings implicate novel genes in delay discounting and highlight the need for further exploration.

The genetic characteristics of Sarcoptes scabiei from Chinese serow (Capricornis milneedwardsii) and goral (Naemorhedus goral arnouxianus) compared with other mites from different hosts and geographic locations using ITS2 and cox1 sequences.

Scabies is a common parasitic disease in many mammalian species, caused by the infestation of Sarcoptes scabiei. There is no consistent conclusion on whether Sarcoptes mites from different hosts or geographic locations have apparent genetic divergence. In this study, we collected and morphologically identified S. scabiei from Chinese serow and goral, and we described the genetic diversity of S. scabiei and other mites based on phylogenetic analyses of the ITS2 and cox1 sequence fragments, including data available in GenBank. The mites isolated from Chinese serow and goral were S. scabiei, and they were morphologically similar. The phylogenetic trees and haplotype networks showed that S. scabiei from other locations worldwide did not cluster according to host divergence or geographical distribution. Additionally, the Fst values were - 0.224 to 0.136 and - 0.045 to 1 between S. scabiei from different hosts, including humans and domestic and wild animals, based on partial ITS and cox1 sequences. Worldwide S. scabiei samples formed three clusters (with H2, H5, and H12 at their centers) in the ITS and one cluster (with C9 at the center) in the cox1 haplotype phylogenetic network. The S. scabiei collected from Chinese serow and goral were morphologically similar and had the same genotype. A study on the genetic characteristics of S. scabiei from Chinese serow and goral together with other mites from different hosts and geographic locations around the world showed no obvious divergence. These findings indicated that scabies likely is a zoonotic disease and that the global prevalence of scabies is probably related to the worldwide trade of domestic animals.

Comparative Genomic Analysis of Xanthomonas campestris pv. campestris Isolates BJSJQ20200612 and GSXT20191014 Provides Novel Insights Into Their Genetic Variability and Virulence.

Black rot is a disease that has a severe impact on cabbage yield and quality in China. Xanthomonas campestris pv. campestris (Xcc) is the causal agent of black rot of Brassicaceae crops. So far, the whole genomic sequences of more than 30 Xcc isolates have been sequenced; however, little information about genomic variability and virulence has been reported. In this study, 12 Xcc isolates were isolated from diseased cabbage leaves in seven Chinese provinces and two municipalities from July 2019 to November 2020. Pathogenicity analysis showed that isolate GSXT20191014 was more aggressive than BJSJQ20200612 and HRIW 3811 on cabbage inbred line 1371. Both BJSJQ20200612 and GSXT20191014 were sequenced and comparatively analyzed. The results showed that BJSJQ20200612 and GSXT20191014 have a single circular chromosome comprising 5,115,975 and 4,975,682 bp, respectively. Compared to the other six sequenced strains, 26 and 47 variable genomic regions were found in BJSJQ2020061 and GSXT20191014 genomic sequences, respectively. The variable genomic regions could be responsible for the genetic variation in Xcc strains and have led to the differences in type III secreted effector repertoires, virulence factors and secreted proteins between these two strains. Among the identified secreted proteins, two copies of peptidase S8/S53 were found in GSXT20191014-specific chromosomal segments. The common effectors xopR, xopH, avrBs1, and xopAH are found in most Xcc genomes, but they are absent in the GSXT20191014 genome. Variations in the composition of exopolysaccharides (EPS) and lipopolysaccharides (LPS) may aid GSXT20191014 isolate infections to evade recognition by the host immune system. Our results revealed a direct correlation between genomic variability and Xcc virulence. We also developed several markers for detecting BJSJQ20200612 and GSXT20191014 isolates and further tested the rest of our other 10 isolates. Finally, the isolated Xcc strains were classified into three genetic subgroups by specific molecular markers and multilocus sequence typing (MLST) approach. BJSJQ20200612 and GSXT20191014 isolates were also classified into two subgroups of Xcc according to the core-genome-based phylogenetic tree. This study extended our understanding of Xcc genomic features and provided the foundation to further characterize the mechanisms for Xcc virulence and a clue for black rot resistance breeding.

First isolation of Trichuris from Chinese serow, Capricornis milneedwardsii (Cetartiodactyla: Caprinae).

Whipworms are soil-transmitted helminths that can infect a variety of animals. A Chinese serow possibly infected by whipworms was found during a wildlife disease surveillance project in Baima Snow Mountain National Nature Reserve, Yunnan, China, in 2021. As convergent evolution is common in coinhabiting parasites, a mitochondrial gene sequence (cox1) and ribosomal gene sequence (ITS1) were used to identify species similar to Trichuris from Chinese serow. The phylogenetic trees and genetic distances of ten Trichuris samples from a Chinese serow together with other Trichuris spp. that have been previously reported were analysed based on the cox1 and ITS1 sequences. The combined results of the phylogenetic tree and genetic distances based on cox1 gene showed that the whipworms in Chinese serow are T. skrjabini. However, the whipworms in the present study were divided into two apparent clades in the phylogenic trees constructed by the cox1 sequences (Clades A and B) and the ITS1 sequences (Clades C and D). In addition, the Fst and Nm values were 0.82 and 0.23 between Clade A and Clade B for the cox1 gene, and 0.30 and 0.45 between Clade C and Clade D for the ITS1 sequences; both indicators showed low gene flow among the clades. Therefore, the genetic population structure of T. skrjabini was illustrated.

Enhanced simultaneous partial nitrification and denitrification performance of aerobic granular sludge via tapered aeration in sequencing batch reactor for treating low strength and low COD/TN ratio municipal wastewater.

The aerobic granular sludge simultaneous partial nitrification, denitrification and phosphorus removal (AGS-SPNDPR) process was carried out via tapered aeration in sequencing batch reactor (SBR) for treating low strength and low COD/TN ratio municipal wastewater. The results showed that aerobic granular sludge was successfully cultivated with good sedimentation performance when treating the municipal wastewater. Meanwhile, the median granule size increased to 270 (R1) and 257 (R2) µm on day 80. The excellent removal performance of COD (92%) and NH4+-N (95%) were achieved under different aeration modes, while the higher TN removal efficiency (76%) was achieved by tapered aeration. The accumulation of NO2--N in R2 indicated that the tapered aeration was beneficial to achieve simultaneously partial nitrification and denitrification. Meanwhile, the high-efficiency phosphorus (95%) removal was realized via additional carbon source, and SPNDPR process was formed under tapered aeration. The bacterial community analysis indicated denitrifying glycogen-accumulating organisms (DGAOs) Candidatus_Competibacter and ammonia-oxidizing bacteria (AOB) Nitrosomonas were more effectively enriched via tapered aeration, while phosphorus-accumulating organisms (PAOs) Candidatus_Accumulibacter were effectively enriched under additional organic carbon. AOB, denitrifying bacteria and PAOs were simultaneously enriched by tapered aeration and additional carbon source, which was beneficial to nutrients removal. This study might be conducive to the application of AGS-SPNDPR system for treating low strength and low COD/TN ratio municipal wastewater under tapered aeration.

Probiotic consumption influences universal adaptive mutations in indigenous human and mouse gut microbiota.

The adaptive evolution in indigenous intestinal microbes derived from probiotics is critical to safety and efficacy evaluation of probiotics, yet it is still largely underexplored. Here, through 11 publicly accessible datasets, we demonstrated that probiotic consumption can lead to widespread single-nucleotide variants (SNVs) in the native microbiota. Interestingly, the same probiotic strains introduced far more SNVs in mouse gut than humans. Furthermore, the pattern of probiotics-induced SNVs was highly probiotic-strain specific, and 17 common SNVs in Faecalibacterium prausnitzii genome were identified cross studies, which might lead to changes in bacterial protein structure. Further, nearly 50% of F. prausnitzii SNVs can be inherited for six months in an independent human cohort, whereas the other half only transiently occurred. Collectively, our study substantially extended our understanding of co-evolution of the probiotics and the indigenous gut microbiota, highlighting the importance of assessment of probiotics efficacy and safety in an integrated manner.

Comparative analysis of the complete mitochondrial genome sequences and anther development cytology between maintainer and Ogura-type cytoplasm male-sterile cabbage (B. oleracea Var. capitata).

BACKGROUND: Cytoplasmic male sterility (CMS) has been widely used for commercial F1 hybrid seeds production. CMS is primarily caused by chimeric genes in mitochondrial genomes. However, which specific stages of anther development in cabbage are affected by the chimeric genes remain unclear. RESULTS: In the present study, the complete mitochondrial genomes were sequenced and assembled for the maintainer and Ogura CMS cabbage lines. The genome size of the maintainer and Ogura CMS cabbage are 219,962 bp and 236,648 bp, respectively. There are 67 and 69 unknown function ORFs identified in the maintainer and Ogura CMS cabbage mitochondrial genomes, respectively. Four orfs, orf102a, orf122b, orf138a and orf154a were specifically identified in the Ogura CMS mitochondrial genome, which were likely generated by recombination with Ogura type radish during breeding process. Among them, ORF138a and ORF154a possessed a transmembrane structure, and orf138a was co-transcribed with the atp8 and trnfM genes. orf154a is partially homologous to the ATP synthase subunit 1 (atpA) gene. Both these genes were likely responsible for the CMS phenotype. In addition, cytological sections showed that the abnormal proliferation of tapetal cells might be the immediate cause of cytoplasmic male-sterility in Ogura CMS cabbage lines. RNA-seq results showed that orf138a and orf154a in Ogura CMS might influence transcript levels of genes in energy metabolic pathways. CONCLUSIONS: The presence of orf138a and orf154a lead to increased of ATPase activity and ATP content by affecting the transcript levels of genes in energy metabolic pathways, which could provide more energy for the abnormal proliferation of tapetal cells. Our data provides new insights into cytoplasmic male-sterility from whole mitochondrial genomes, cytology of anther development and transcriptome data.

Detection of Clinically Relevant Genetic Variants in Chinese Patients With Nanophthalmos by Trio-Based Whole-Genome Sequencing Study.

Purpose: Nanophthalmos is a rare genetic disorder commonly characterized by a short axial length (AL) and severe hyperopia. Mutations that have been identified through Mendelian genetic analysis can only explain a fraction of nanophthalmic cases. We investigate the clinically relevant genetic variants in nanophthalmos by whole-genome sequencing (WGS), including de novo mutations (DNMs) and inherited mutations. Methods: Clinically relevant genetic variants of 11 trios (11 nanophthalmic probands and their unaffected parents) from the Zhongshan Ophthalmic Center, China, were analyzed by WGS. We further screened three trios and 10 sporadic cases to identify the MYRF mutations. Results: In two of 11 trios, without evidence of the presence of deleterious inherited autosomal variants, two DNMs of MYRF (c.789delC, p.S264fs and c.789dupC, p.S264fs) were identified in the probands. These loss-of-function DNMs were predicted to result in premature stop codons and protein structure damage in both probands. In addition, deleterious inherited genetic variants in PRSS56 and MFRP were found in eight probands of the other nine trios. Expanded screening found an additional MYRF DNM (c.1433G>C, p.R478P) in one trio and a stop-gain MYRF mutation (c.2956C>T, p.R986X) in one sporadic case, suggesting the recurrence of MYRF mutations in nanophthalmic patients. Conclusions: This is the first trio-based WGS study for nanophthalmos, revealing the potential role of DNMs in MYRF and rare inherited genetic variants in PRSS56 and MFRP. The underlying mechanism of MYRF in the development of nanophthalmos needs to be further investigated.

Anterior Segment Features in Nanophthalmos With Secondary Chronic Angle Closure Glaucoma: An Ultrasound Biomicroscopy Study.

Purpose: We compare anterior segment parameters and the risk of malignant glaucoma between nanophthalmos with secondary chronic angle closure glaucoma (NSCACG) and chronic primary angle-closure glaucoma (CPACG). Methods: This retrospective case-control study included 32 NSCACG (32 eyes) and 36 CPACG (36 eyes) patients. Anterior segment parameters, including anterior chamber depth (ACD) and width (ACW), pupil diameter (PD), lens vault (LV), ciliary process-ciliary process distance (CCD), angle opening distance500 (AOD500), trabecular-iris angle (TIA), trabecular-ciliary processes distance (TCPD), ciliary body max thickness (CBMT), iris convexity (IC), peripheral iris thickness (PIT), iris-zonule distance (IZD), trabecular ciliary process angle (TCPA), and anterior vault (AV), were measured by ultrasound biomicroscopy. A-scan ultrasonography measurements, including lens thickness (LT) and axial length (AL), also were reviewed. Results: ACD, ACW, CCD, TCPD, and AL were smaller, whereas the LV, LT/AL, LV/LT, LV/ACD, LT/ACD, LV/AV, and AV/AL were larger in NSCACG compared to CPACG eyes (all Bonferroni-corrected P < 0.05). Eyes with NSCACG had a higher risk of malignant glaucoma postoperatively than eyes with CPACG (P = 0.018). Conclusions: Characterized by narrower anterior segment, forward movement of larger lens, more anteriorly rotated ciliary bodies, and smaller CCD, NSCACG eyes have a higher risk of malignant glaucoma than CPACG eyes.

IDGenetics: a comprehensive database for genes and mutations of intellectual disability related disorders.

Intellectual disability (ID) is one of the most prevalent chronic developmental brain disorders or phenotype of syndromic ID, affecting nearly 1-2% of the general population worldwide. Over recent decades, tremendous effort and high-throughput platforms have been devised to explore the complex heterogeneity, numerous genes and variants have been associated with the ID, especially de novo mutations and copy number variants. An organized resource containing the increasing genetic data is imperative to assist ID research. In this study, the integrative and annotated intellectual disability database has been developed, named 'IDGenetics', which contains known information about ID, including 815 genes and 17102 variants associated with 918 clinical diseases (3001 clinical phenotype) collected from 3822 publications and ID-related databases. Furthermore, in-depth data mining was performed to obtain an understanding of each entry, including functional annotation, gene/disease/phenotype network establishment and overlap analysis focusing on comorbidity. 1478 candidate genes (483 high-confidence and 995 low-confidence) were collected and prioritized by adopting the annotations of 12 functional prediction tools and algorithm. In addition, IDGenetics database provides concise search methods, convenient browsing functions, intuitive graphical displays and constantly updated features. IDGenetics will be a valuable and integrative resource for deciphering the genetic and functional architecture of ID and the improvement of clinical diagnosis, intervention and treatment.

mirDNMR: a gene-centered database of background de novo mutation rates in human.

De novo germline mutations (DNMs) are the rarest genetic variants proven to cause a considerable number of sporadic genetic diseases, such as autism spectrum disorders, epileptic encephalopathy, schizophrenia, congenital heart disease, type 1 diabetes, and hearing loss. However, it is difficult to accurately assess the cause of DNMs and identify disease-causing genes from the considerable number of DNMs in probands. A common method to this problem is to identify genes that harbor significantly more DNMs than expected by chance, with accurate background DNM rate (DNMR) required. Therefore, in this study, we developed a novel database named mirDNMR for the collection of gene-centered background DNMRs obtained from different methods and population variation data. The database has the following functions: (i) browse and search the background DNMRs of each gene predicted by four different methods, including GC content (DNMR-GC), sequence context (DNMR-SC), multiple factors (DNMR-MF) and local DNA methylation level (DNMR-DM); (ii) search variant frequencies in publicly available databases, including ExAC, ESP6500, UK10K, 1000G and dbSNP and (iii) investigate the DNM burden to prioritize candidate genes based on the four background DNMRs using three statistical methods (TADA, Binomial and Poisson test). As a case study, we successfully employed our database in candidate gene prioritization for a sporadic complex disease: intellectual disability. In conclusion, mirDNMR (https://www.wzgenomics.cn/mirdnmr/) can be widely used to identify the genetic basis of sporadic genetic diseases.

mirVAFC: A Web Server for Prioritizations of Pathogenic Sequence Variants from Exome Sequencing Data via Classifications.

Exome sequencing has been widely used to identify the genetic variants underlying human genetic disorders for clinical diagnoses, but the identification of pathogenic sequence variants among the huge amounts of benign ones is complicated and challenging. Here, we describe a new Web server named mirVAFC for pathogenic sequence variants prioritizations from clinical exome sequencing (CES) variant data of single individual or family. The mirVAFC is able to comprehensively annotate sequence variants, filter out most irrelevant variants using custom criteria, classify variants into different categories as for estimated pathogenicity, and lastly provide pathogenic variants prioritizations based on classifications and mutation effects. Case studies using different types of datasets for different diseases from publication and our in-house data have revealed that mirVAFC can efficiently identify the right pathogenic candidates as in original work in each case. Overall, the Web server mirVAFC is specifically developed for pathogenic sequence variant identifications from family-based CES variants using classification-based prioritizations. The mirVAFC Web server is freely accessible at https://www.wzgenomics.cn/mirVAFC/.