Search

A male factor, commonly associated with poor semen quality, is revealed in about 50% of infertile couples. CFTR gene (Cystic Fibrosis Transmembrane Conduction Regulator) variants are one of the common genetic causes of azoospermia-related male infertility. Notably, the spectrum and frequency of pathogenic CFTR variants vary between populations and geographical regions. In this work, we made an attempt to evaluate the allele frequency (AF) of 12 common CFTR variants in infertile Russian men and healthy individuals from different districts of Russia. Because of the limited number of population-based studies on Russian individuals, we characterized the population AFs based on data from the Registry of Russian cystic fibrosis (CF) patients. In addition to the CF patient registry, we estimated the local frequencies of the same set of variants based on the results of genotyping of CF patients in local biocollections (from St. Petersburg and Yugra regions). AFs of common CFTR variants calculated based on registry and biocollection data showed good concordance with directly measured population AFs. The estimated region-specific frequencies of CFTR variants allowed us to uncover statistically significant regional differences in the frequencies of the F508del (c.1521_1523del; p.Phe508del) and CFTRdele2,3(21kb) (c.54-5940_273+10250del21kb; p.Ser18ArgfsX) variants. The data from population-based studies confirmed previous observations that F508del, CFTRdele2,3(21kb), and L138ins (c.413_415dup; p.Leu138dup)variants are the most abundant among infertile patients, and their frequencies are significantly lower in healthy individuals and should be taken into account during genetic monitoring of the reproductive health of Russian individuals. [ABSTRACT FROM AUTHOR]

Although high altitude training has been increasingly popular among endurance athletes, the molecular and cellular bases of this adaptation remain poorly understood. We aimed to define the underlying physiological changes and screen for potential biomarkers of adaptation using transcriptional profiling of whole blood. Seven elite female speed skaters were profiled on the 18th day of high-altitude adaptation. Whole blood RNA-seq before and after an intense 1 h skating bout was used to measure gene expression changes associated with exercise. In order to identify the genes specifically regulated at high altitudes, we have leveraged the data from eight previously published microarray datasets studying blood expression changes after exercise at sea level. Using cell type-specific signatures, we were able to deconvolute changes of cell type abundance from individual gene expression changes. Among these were PHOSPHO1, with a known role in erythropoiesis, and MARC1 with a role in endogenic NO metabolism. We find that platelet and erythrocyte counts uniquely respond to altitude exercise, while changes in neutrophils represent a more generic marker of intense exercise. Publicly available data from both single cell atlases and exercise-related blood profiling dramatically increases the value of whole blood RNA-seq for the dynamic evaluation of physiological changes in an athlete's body. [ABSTRACT FROM AUTHOR]

Background: Preeclampsia (PE) is the most common complication of pregnancy that remains to be a major cause of maternal and fetal mortality. Prediction and early diagnosis of PE would allow for timely initiation of preventive therapy. According to recent studies of ACVR2A gene polymorphism is associated with PE, but it is still unclear whether these findings reflect specific pathogenetic mechanisms of this disease. Methods: We performed targeted next-generation sequencing (NGS) sequencing of ACVR2A gene by means of Ion Torrent Personal Genome machine (PGM) Sequencer. A genetic analysis of patients with PE and control group was performed. Bioinformatics analysis using Polyphen2 (Boston, MA), SIFT (La Jolla, CA), and SnpSift software were used. To select genetic markers in PE patients two additive models and score analysis were applied. Results: Based on the score analysis, we detected two substitutions (rs145399059 and rs17692648) and one insertion insAA at position 148642724 that were associated with PE in our cohorts. We also detected a variant rs17742573 that can be considered as protective against preeclampsia. Conclusions: Our data suggest that some variants in ACVR2A gene are associated with PE. But more studies are required to reveal the role of ACVR2A gene in the pathogenesis of this disease during pregnancy. [ABSTRACT FROM AUTHOR]

Today, whole-exome sequencing (WES) is used to conduct the massive screening of structural and regulatory genes in order to identify the allele frequencies of disease-associated polymorphisms in various populations and thus detect pathogenic genetic changes (mutations or polymorphisms) conducive to malfunctional protein sequences. With its extensive capabilities, exome sequencing today allows both the diagnosis of monogenic diseases (MDs) and the examination of seemingly healthy populations to reveal a wide range of potential risks prior to disease manifestation (in the future, exome sequencing may outpace costly and less informative genome sequencing to become the first-line examination technique). This review establishes the human genetic passport as a new WES-based clinical concept for the identification of new candidate genes, gene variants, and molecular mechanisms in the diagnosis, prediction, and treatment of monogenic, oligogenic, and multifactorial diseases. Various diseases are addressed to demonstrate the extensive potential of WES and consider its advantages as well as disadvantages. Thus, WES can become a general test with a broad spectrum pf applications, including opportunistic screening. [ABSTRACT FROM AUTHOR]

Pre-eclampsia (PE) is a complication of pregnancy that affects 5-8% of women after 20 weeks of gestation. It is usually diagnosed based on the de novo onset of hypertension and proteinuria. Preexisting hypertension in women developing PE, also known as superimposed PE on chronic hypertension (SPE), leads to elevated risk of maternal and fetal mortality. PE is associated with an altered microRNA (miRNA) expression pattern in the placenta, suggesting that miRNA deregulation is involved in the pathogenesis of PE. Whether and how the miRNA expression pattern is changed in the SPE placenta remains unclear. The present study analyzed the placental miRNA expression profile in pregnancies complicated by SPE. miRNA expression profiles in SPE and normal placentas were investigated using an Ion Torrent sequencing system. Sequencing data were processed using a comprehensive analysis pipeline for deep miRNA sequencing (CAP-miRSeq). A total of 22 miRNAs were identified to be deregulated in placentas from patients with SPE. They included 16 miRNAs previously known to be associated with PE and 6 novel miRNAs. Among the 6 novel miRNAs, 4 were upregulated (miR-518a, miR-527, miR-518e and miR-4532) and 2 downregulated (miR-98 and miR-135b) in SPE placentas compared with controls. The present results suggest that SPE is associated with specific alterations in the placental miRNA expression pattern, which differ from alterations detected in PE placentas, and therefore, provide novel targets for further investigation of the molecular mechanisms underlying SPE pathogenesis. [ABSTRACT FROM AUTHOR]

Background Hypertrophic cardiomyopathy is a common genetic cardiac disease. Prevention and early diagnosis of this disease are very important. Because of the large number of causative genes and the high rate of mutations involved in the pathogenesis of this disease, traditional methods of early diagnosis are ineffective. Methods We developed a custom AmpliSeq panel for NGS sequencing of the coding sequences of ACTC1 , MYBPC3 , MYH7 , MYL2 , MYL3 , TNNI3 , TNNT2 , TPM1 , and CASQ2 . A genetic analysis of student cohorts (with and without cardiomyopathy risk in their medical histories) and patients with cardiomyopathies was performed. For the statistical and bioinformatics analysis, Polyphen2, SIFT, SnpSift and PLINK software were used. To select genetic markers in the patients with cardiomyopathy and in the students of the high risk group, four additive models were applied. Results Our AmpliSeq custom panel allowed us to efficiently explore targeted sequences. Based on the score analysis, we detected three substitutions in the MYBPC3 and CASQ2 genes and six combinations between loci in the MYBPC3 , MYH7 and CASQ2 genes that were responsible for cardiomyopathy risk in our cohorts. We also detected substitutions in the TNNT2 gene that can be considered as protective against cardiomyopathy. Conclusion We used NGS with AmpliSeq libraries and Ion PGM sequencing to develop improved predictive information for patients at risk of cardiomyopathy. [ABSTRACT FROM AUTHOR]

Background: Pre-eclampsia is the most common complication occurring during pregnancy. In the majority of cases, it is concurrent with other pathologies in a comorbid manner (frequent co-occurrences in patients), such as diabetes mellitus, gestational diabetes and obesity. Providing bronchial asthma, pulmonary tuberculosis, certain neurodegenerative diseases and cancers as examples, we have shown previously that pairs of inversely comorbid pathologies (rare co-occurrences in patients) are more closely related to each other at the molecular genetic level compared with randomly generated pairs of diseases. Data in the literature concerning the causes of pre-eclampsia are abundant. However, the key mechanisms triggering this disease that are initiated by other pathological processes are thus far unknown. The aim of this work was to analyse the characteristic features of genetic networks that describe interactions between comorbid diseases, using pre-eclampsia as a case in point. Results: The use of ANDSystem, Pathway Studio and STRING computer tools based on text-mining and databasemining approaches allowed us to reconstruct associative networks, representing molecular genetic interactions between genes, associated concurrently with comorbid disease pairs, including pre-eclampsia, diabetes mellitus, gestational diabetes and obesity. It was found that these associative networks statistically differed in the number of genes and interactions between them from those built for randomly chosen pairs of diseases. The associative network connecting all four diseases was composed of 16 genes (PLAT, ADIPOQ, ADRB3, LEPR, HP, TGFB1, TNFA, INS, CRP, CSRP1, IGFBP1, MBL2, ACE, ESR1, SHBG, ADA). Such an analysis allowed us to reveal differential gene risk factors for these diseases, and to propose certain, most probable, theoretical mechanisms of pre-eclampsia development in pregnant women. The mechanisms may include the following pathways: [TGFB1 or TNFA]-[IL1B]-[pre-eclampsia]; [TNFA or INS]-[NOS3]-[pre-eclampsia]; [INS]-[HSPA4 or CLU]-[pre-eclampsia]; [ACE]-[MTHFR]-[pre-eclampsia]. Conclusions: For pre-eclampsia, diabetes mellitus, gestational diabetes and obesity, we showed that the size and connectivity of the associative molecular genetic networks, which describe interactions between comorbid diseases, statistically exceeded the size and connectivity of those built for randomly chosen pairs of diseases. Recently, we have shown a similar result for inversely comorbid diseases. This suggests that comorbid and inversely comorbid diseases have common features concerning structural organization of associative molecular genetic networks. [ABSTRACT FROM AUTHOR]

Adaptation of humans to challenging environmental conditions, such as extreme temperature, malnutrition, or hypoxia, is an interesting phenomenon for both basic and applied research. Identification of the genetic factors contributing to human adaptation to these conditions enhances our understanding of the underlying molecular and physiological mechanisms. In our study, we analyzed the exomes of 22 high altitude mountaineers to uncover genetic variants contributing to hypoxic adaptation. To our surprise, we identified two putative loss‐of‐function variants, rs1385101139 in RTEL1 and rs1002726737 in COL6A1 in two extremely high altitude (personal record of more than 8500 m) professional climbers. Both variants can be interpreted as pathogenic according to medical geneticists' guidelines, and are linked to inherited conditions involving respiratory failure (late‐onset pulmonary fibrosis and severe Ullrich muscular dystrophy for rs1385101139 and rs1002726737, respectively). Our results suggest that a loss of gene function may act as an important factor of human adaptation, which is corroborated by previous reports in other human subjects. [ABSTRACT FROM AUTHOR]

Next-generation sequencing (NGS) has revolutionized the field of rare disease diagnostics. Whole exome and whole genome sequencing are now routinely used for diagnostic purposes; however, the overall diagnosis rate remains lower than expected. In this work, we review current approaches used for calling and interpretation of germline genetic variants in the human genome, and discuss the most important challenges that persist in the bioinformatic analysis of NGS data in medical genetics. We describe and attempt to quantitatively assess the remaining problems, such as the quality of the reference genome sequence, reproducible coverage biases, or variant calling accuracy in complex regions of the genome. We also discuss the prospects of switching to the complete human genome assembly or the human pan-genome and important caveats associated with such a switch. We touch on arguably the hardest problem of NGS data analysis for medical genomics, namely, the annotation of genetic variants and their subsequent interpretation. We highlight the most challenging aspects of annotation and prioritization of both coding and non-coding variants. Finally, we demonstrate the persistent prevalence of pathogenic variants in the coding genome, and outline research directions that may enhance the efficiency of NGS-based disease diagnostics. [ABSTRACT FROM AUTHOR]

Preeclampsia is a disorder that can occur during pregnancy and is one of the leading causes of death among pregnant women. This disorder occurs after the 20th week of pregnancy and is characterized by arterial hypertension, proteinuria, fetoplacental, and multiple organ dysfunctions. Despite the long history of studying preeclampsia, its etiology and pathogenesis remain poorly understood, and therapy is symptomatic. One of the factors of the disorder is believed to be misfolded proteins that are prone to form amyloid aggregates. The CRD tests, utilizing the binding of the amyloid-specific dye Congo red to urine components, demonstrate high efficiency in diagnosing preeclampsia. However, these tests have also been found to be positive in other disorders with proteinuria, presumably associated with concomitant amyloidosis. To assess the limitations of the CRD tests, we examined urine congophilia and protein components mediating Congo red positivity in patients with proteinuria, including preeclampsia, amyloid and non-amyloid nephropathies. We stained the urine samples and calculated congophilia levels. We also assessed the contribution of large protein aggregates to congophilia values using ultracentrifugation and determined the molecular weights of congophilic urinary proteins using centrifugal concentrators. All proteinuric groups demonstrate positive results in the CRD tests and congophilia levels were more than two times higher compared with the control non-proteinuric groups (p <0.01). There was a strong correlation between urine protein excretion and congophilia in amyloid nephropathy (rs = 0.76), non-amyloid nephropathies (rs = 0.90), and preeclampsia (rs = 0.90). Removal of large aggregates from urine did not affect the congophilia levels. Separation of urine protein fractions revealed congophilic components in the range of 30–100 kDa, including monomeric serum albumin. Our results indicate limitations of CRD tests in preeclampsia diagnostics in women with renal disorders and underscore the need for further research on the mechanisms of Congo red binding with urine components. [ABSTRACT FROM AUTHOR]

Simple Summary: Genome-wide association studies (GWAS) are a powerful tool for the identification of genes affecting human traits. Still, the interpretation of GWAS results is complicated, and new tools are actively being developed. Due to the scarcity of available datasets, simulation of GWAS data with known genetic effects is important as it enables accurate evaluation of such tools. In this study, we developed a flexible tool, bioGWAS, that provides a set of important functionalities for simulating GWAS results. We demonstrate that bioGWAS can efficiently generate GWAS results with predefined causal genes and biological processes and is capable of recapitulating the results of published GWAS studies. We thus believe that bioGWAS is an excellent method for testing bioinformatics software for GWAS results processing, as well as for the generation of datasets for educational purposes. Genome-wide association studies (GWAS) have proven to be a powerful tool for the identification of genetic susceptibility loci affecting human complex traits. In addition to pinpointing individual genes involved in a particular trait, GWAS results can be used to discover relevant biological processes for these traits. The development of new tools for extracting such information from GWAS results requires large-scale datasets with known biological ground truth. Simulation of GWAS results is a powerful method that may provide such datasets and facilitate the development of new methods. In this work, we developed bioGWAS, a simple and flexible pipeline for the simulation of genotypes, phenotypes, and GWAS summary statistics. Unlike existing methods, bioGWAS can be used to generate GWAS results for simulated quantitative and binary traits with a predefined set of causal genetic variants and/or molecular pathways. We demonstrate that the proposed method can recapitulate complete GWAS datasets using a set of reported genome-wide associations. We also used our method to benchmark several tools for gene set enrichment analysis for GWAS data. Taken together, our results suggest that bioGWAS provides an important set of functionalities that would aid the development of new methods for downstream processing of GWAS results. [ABSTRACT FROM AUTHOR]

Pregnancy loss is the most frequent complication of a pregnancy which is devastating for affected families and poses a significant challenge for the health care system. Genetic factors are known to play an important role in the etiology of pregnancy loss; however, despite advances in diagnostics, the causes remain unexplained in more than 30% of cases. In this review, we aggregated the results of the decade-long studies into the genetic risk factors of pregnancy loss (including miscarriage, termination for fetal abnormality, and recurrent pregnancy loss) in euploid pregnancies, focusing on the spectrum of point mutations associated with these conditions. We reviewed the evolution of molecular genetics methods used for the genetic research into causes of pregnancy loss, and collected information about 270 individual genetic variants in 196 unique genes reported as genetic cause of pregnancy loss. Among these, variants in 18 genes have been reported by multiple studies, and two or more variants were reported as causing pregnancy loss for 57 genes. Further analysis of the properties of all known pregnancy loss genes showed that they correspond to broadly expressed, highly evolutionary conserved genes involved in crucial cell differentiation and developmental processes and related signaling pathways. Given the features of known genes, we made an effort to construct a list of candidate genes, variants in which may be expected to contribute to pregnancy loss. We believe that our results may be useful for prediction of pregnancy loss risk in couples, as well as for further investigation and revealing genetic etiology of pregnancy loss. [ABSTRACT FROM AUTHOR]

Type 2 diabetes mellitus (T2D) is a chronic metabolic disease characterized by insulin resistance and β-cell dysfunction and leading to many micro- and macrovascular complications. In this study we analyzed the circulating miRNA expression profiles in plasma samples from 44 patients with T2D and 22 healthy individuals using next generation sequencing and detected 229 differentially expressed miRNAs. An increased level of miR-5588-5p, miR-125b-2-3p, miR-1284, and a reduced level of miR-496 in T2D patients was verified. We also compared the expression landscapes in the same group of patients depending on body mass index and identified differential expression of miR-144-3p and miR-99a-5p in obese individuals. Identification and functional analysis of putative target genes was performed for miR-5588-5p, miR-125b-2-3p, miR-1284, and miR-496, showing chromatin modifying enzymes and apoptotic genes being among the significantly enriched pathways. [ABSTRACT FROM AUTHOR]

Phenotypicheterogeneity is a phenomenon in which distinct phenotypes can develop in individuals bearing pathogenic variants in the same gene. Genetic factors, gene interactions, and environmental factors are usually considered the key mechanisms of this phenomenon. Phenotypic heterogeneity may impact the prognosis of the disease severity and symptoms. In our work, we used publicly available data on the association between genetic variants and Mendelian disease to investigate the genetic factors (such as the intragenic localization and type of a variant) driving the heterogeneity of gene–disease relationships. First, we showed that genes linked to multiple rare diseases (GMDs) are more constrained and tend to encode more transcripts with high levels of expression across tissues. Next, we assessed the role of variant localization and variant types in specifying the exact phenotype for GMD variants. We discovered that none of these factors is sufficient to explain the phenomenon of such heterogeneous gene–disease relationships. In total, we identified only 38 genes with a weak trend towards significant differences in variant localization and 30 genes with nominal significant differences in variant type for the two associated disorders. Remarkably, four of these genes showed significant differences in both tests. At the same time, our analysis suggests that variant localization and type are more important for genes linked to autosomal dominant disease. Taken together, our results emphasize the gene-level factors dissecting distinct Mendelian diseases linked to one common gene based on open-access genetic data and highlight the importance of exploring other factors that contributed to phenotypic heterogeneity. [ABSTRACT FROM AUTHOR]

We report on the case of prenatal detection of trisomy 2 in placental biopsy and further algorithm of genetic counseling and testing. A 29-year-old woman with first-trimester biochemical markers refused chorionic villus sampling and preferred targeted non-invasive prenatal testing (NIPT), which showed low risk for aneuploidies 13, 18, 21, and X. A series of ultrasound examinations revealed increased chorion thickness at 13/14 weeks of gestation and fetal growth retardation, a hyperechoic bowel, challenging visualization of the kidneys, dolichocephaly, ventriculomegaly, increase in placental thickness, and pronounced oligohydramnios at 16/17 weeks of gestation. The patient was referred to our center for an invasive prenatal diagnosis. The patient's blood and placenta were sampled for whole-genome sequencing-based NIPT and array comparative genomic hybridization (aCGH), respectively. Both investigations revealed trisomy 2. Further prenatal genetic testing in order to confirm trisomy 2 in amniocytes and/or fetal blood was highly questionable because oligohydramnios and fetal growth retardation made amniocentesis and cordocentesis technically unfeasible. The patient opted to terminate the pregnancy. Pathological examination of the fetus revealed internal hydrocephalus, atrophy of brain structure, and craniofacial dysmorphism. Conventional cytogenetic analysis and fluorescence in situ hybridization revealed chromosome 2 mosaicism with a prevalence of trisomic clone in the placenta (83.2% vs. 16.8%) and a low frequency of trisomy 2, which did not exceed 0.6% in fetal tissues, advocating for low-level true fetal mosaicism. To conclude, in pregnancies at risk of fetal chromosomal abnormalities that refuse invasive prenatal diagnosis, whole-genome sequencing-based NIPT, but not targeted NIPT, should be considered. In prenatal cases of trisomy 2, true mosaicism should be distinguished from placental-confined mosaicism using cytogenetic analysis of amniotic fluid cells or fetal blood cells. However, if material sampling is impossible due to oligohydramnios and/or fetal growth retardation, further decisions should be based on a series of high-resolution fetal ultrasound examinations. Genetic counseling for the risk of uniparental disomy in a fetus is also required. [ABSTRACT FROM AUTHOR]

Single-cell RNA sequencing (scRNA-seq) is a method that focuses on the analysis of gene expression profile in individual cells. This method has been successfully applied to answer the challenging questions of the pathogenesis of multifactorial diseases and open up new possibilities in the prognosis and prevention of reproductive diseases. In this article, we have reviewed the application of scRNA-seq to the analysis of the various cell types and their gene expression changes in normal pregnancy and pregnancy complications. The main principle, advantages, and limitations of single-cell technologies and data analysis methods are described. We discuss the possibilities of using the scRNA-seq method for solving the fundamental and applied tasks related to various pregnancy-associated disorders. Finally, we provide an overview of the scRNA-seq findings for the common pregnancy-associated conditions, such as hyperglycemia in pregnancy, recurrent pregnancy loss, preterm labor, polycystic ovary syndrome, and pre-eclampsia. [ABSTRACT FROM AUTHOR]

Complications endangering mother or fetus affect around one in seven pregnant women. Investigation of the genetic susceptibility to such diseases is of high importance for better understanding of the disease biology as well as for prediction of individual risk. In this study, we collected and analyzed GWAS summary statistics from the FinnGen cohort and UK Biobank for 24 pregnancy complications. In FinnGen, we identified 11 loci associated with pregnancy hypertension, excessive vomiting, and gestational diabetes. When UK Biobank and FinnGen data were combined, we discovered six loci reaching genome-wide significance in the meta-analysis. These include rs35954793 in FGF5 ( p = 6.1 × 10 − 9 ), rs10882398 in PLCE1 ( p = 8.9 × 10 − 9 ), and rs167479 in RGL3 ( p = 5.2 × 10 − 9 ) for pregnancy hypertension, rs10830963 in MTNR1B ( p = 4.5 × 10 − 41 ) and rs36090025 in TCF7L2 ( p = 3.4 × 10 − 15 ) for gestational diabetes, and rs2963457 in the EBF1 locus ( p = 6.5 × 10 − 9 ) for preterm birth. In addition to the identified genome-wide associations, we also replicated 14 out of 40 previously reported GWAS markers for pregnancy complications, including four more preeclampsia-related variants. Finally, annotation of the GWAS results identified a causal relationship between gene expression in the cervix and gestational hypertension, as well as both known and previously uncharacterized genetic correlations between pregnancy complications and other traits. These results suggest new prospects for research into the etiology and pathogenesis of pregnancy complications, as well as early risk prediction for these disorders. [ABSTRACT FROM AUTHOR]

In recent years, great advances have been made in the field of collection, storage, and analysis of biological samples. Large collections of samples, biobanks, have been established in many countries. Biobanks typically collect large amounts of biological samples and associated clinical information; the largest collections include over a million samples. In this review, we summarize the main directions in which biobanks aid medical genetics and genomic research, from providing reference allele frequency information to allowing large-scale cross-ancestry meta-analyses. The largest biobanks greatly vary in the size of the collection, and the amount of available phenotype and genotype data. Nevertheless, all of them are extensively used in genomics, providing a rich resource for genome-wide association analysis, genetic epidemiology, and statistical research into the structure, function, and evolution of the human genome. Recently, multiple research efforts were based on trans-biobank data integration, which increases sample size and allows for the identification of robust genetic associations. We provide prominent examples of such data integration and discuss important caveats which have to be taken into account in trans-biobank research. [ABSTRACT FROM AUTHOR]

In recent years evidence has been accumulated showing that miRNAs can act as potential biomarkers or targets for therapy of preterm birth, one of the most important problems in modern obstetrics. We have performed a prospective study of the miRNA profile in the plasma during the first and second trimesters in pregnant women with high risk of preterm birth (n = 13 cases and n = 11 controls). For the study group plasma blood samples at 9–13 weeks before diagnosis and at 22–24 weeks after start of therapy were selected. Using high-throughput sequencing technology we detected differences in the levels of 15 miRNAs (3 upregulated—hsa-miR-122-5p, hsa-miR-34a-5p, hsa-miR-34c-5p; 12 downregulated—hsa-miR-487b-3p, hsa-miR-493-3p, hsa-miR-432-5p, hsa-miR-323b-3p, hsa-miR-369-3p, hsa-miR-134-5p, hsa-miR-431-5p, hsa-miR-485-5p, hsa-miR-382-5p, hsa-miR-369-5p, hsa-miR-485-3p, hsa-miR-127-3p) (log2(FC) ≥ 1.5; FDR ≤ 0.05) during the first trimester compared with the control (non-high-risk of preterm birth pregnant women). All downregulated miRNAs in the first trimester from the placenta-specific C14MC cluster. During the second trimester no differentially expressed miRNAs were found. Our results suggest that the miRNA profile in plasma during early pregnancy may predict a high risk of preterm birth, which is important in preventing gestational problems as early as possible. [ABSTRACT FROM AUTHOR]

Lopes–Maciel–Rodan syndrome (LOMARS) is an extremely rare disorder, with only a few cases reported worldwide. LOMARS is caused by a compound heterozygous mutation in the HTT gene. Little is known about LOMARS pathogenesis and clinical manifestations. Whole exome sequencing (WES) was performed to achieve a definitive molecular diagnosis of the disorder. All NGS-identified variants underwent the Sanger confirmation. In addition, a literature review on genetic variations in the HTT gene was conducted. The paper reports a case of LOMARS in a pediatric patient in Russia. A preterm girl of non-consanguineous parents demonstrated severe psychomotor developmental delays in her first 12 months. By the age of 6 years, she failed to develop speech but was able to understand everyday phrases and perform simple commands. Autism-like behaviors, stereotypies, and bruxism were noted during the examination. WES revealed two undescribed variants of unknown clinical significance in the HTT gene, presumably associated with the patient's phenotype (c.2350C>T and c.8440C>A). Medical re-examination of parents revealed that the patient inherited these variants from her father and mother. Lopes–Maciel–Rodan syndrome was diagnosed based on overlapping clinical findings and the follow-up genetic examination of parents. Our finding expands the number of reported LOMARS cases and provides new insights into the genetic basis of the disease. [ABSTRACT FROM AUTHOR]

Metformin is an oral hypoglycemic agent widely used in clinical practice for treatment of patients with type 2 diabetes mellitus (T2DM). The wide interindividual variability of response to metformin therapy was shown, and recently the impact of several genetic variants was reported. To assess the independent and combined effect of the genetic polymorphism on glycemic response to metformin, we performed an association analysis of the variants in ATM, SLC22A1, SLC47A1, and SLC2A2 genes with metformin response in 299 patients with T2DM. Likewise, the distribution of allele and genotype frequencies of the studied gene variants was analyzed in an extended group of patients with T2DM (n = 464) and a population group (n = 129). According to our results, one variant, rs12208357 in the SLC22A1 gene, had a significant impact on response to metformin in T2DM patients. Carriers of TT genotype and T allele had a lower response to metformin compared to carriers of CC/CT genotypes and C allele (p-value = 0.0246, p-value = 0.0059, respectively). To identify the parameters that had the greatest importance for the prediction of the therapy response to metformin, we next built a set of machine learning models, based on the various combinations of genetic and phenotypic characteristics. The model based on a set of four parameters, including gender, rs12208357 genotype, familial T2DM background, and waist–hip ratio (WHR) showed the highest prediction accuracy for the response to metformin therapy in patients with T2DM (AUC = 0.62 in cross-validation). Further pharmacogenetic studies may aid in the discovery of the fundamental mechanisms of type 2 diabetes, the identification of new drug targets, and finally, it could advance the development of personalized treatment. [ABSTRACT FROM AUTHOR]

Type 2 diabetes (T2D) is a common chronic disease whose etiology is known to have a strong genetic component. Standard genetic approaches, although allowing for the detection of a number of gene variants associated with the disease as well as differentially expressed genes, cannot fully explain the hereditary factor in T2D. The explosive growth in the genomic sequencing technologies over the last decades provided an exceptional impetus for transcriptomic studies and new approaches to gene expression measurement, such as RNA-sequencing (RNA-seq) and single-cell technologies. The transcriptomic analysis has the potential to find new biomarkers to identify risk groups for developing T2D and its microvascular and macrovascular complications, which will significantly affect the strategies for early diagnosis, treatment, and preventing the development of complications. In this article, we focused on transcriptomic studies conducted using expression arrays, RNA-seq, and single-cell sequencing to highlight recent findings related to T2D and challenges associated with transcriptome experiments. [ABSTRACT FROM AUTHOR]

The COVID-19 pandemic has drawn the attention of many researchers to the interaction between pathogen and host genomes. Over the last two years, numerous studies have been conducted to identify the genetic risk factors that predict COVID-19 severity and outcome. However, such an analysis might be complicated in cohorts of limited size and/or in case of limited breadth of genome coverage. In this work, we tried to circumvent these challenges by searching for candidate genes and genetic variants associated with a variety of quantitative and binary traits in a cohort of 840 COVID-19 patients from Russia. While we found no gene- or pathway-level associations with the disease severity and outcome, we discovered eleven independent candidate loci associated with quantitative traits in COVID-19 patients. Out of these, the most significant associations correspond to rs1651553 in MYH14p = 1.4 × 10−7), rs11243705 in SETX (p = 8.2 × 10−6), and rs16885 in ATXN1 (p = 1.3 × 10−5). One of the identified variants, rs33985936 in SCN11A, was successfully replicated in an independent study, and three of the variants were found to be associated with blood-related quantitative traits according to the UK Biobank data (rs33985936 in SCN11A, rs16885 in ATXN1, and rs4747194 in CDH23). Moreover, we show that a risk score based on these variants can predict the severity and outcome of hospitalization in our cohort of patients. Given these findings, we believe that our work may serve as proof-of-concept study demonstrating the utility of quantitative traits and extensive phenotyping for identification of genetic risk factors of severe COVID-19. [ABSTRACT FROM AUTHOR]

Background: Accurate variant detection in the coding regions of the human genome is a key requirement for molecular diagnostics of Mendelian disorders. Efficiency of variant discovery from next-generation sequencing (NGS) data depends on multiple factors, including reproducible coverage biases of NGS methods and the performance of read alignment and variant calling software. Although variant caller benchmarks are published constantly, no previous publications have leveraged the full extent of available gold standard whole-genome (WGS) and whole-exome (WES) sequencing datasets. Results: In this work, we systematically evaluated the performance of 4 popular short read aligners (Bowtie2, BWA, Isaac, and Novoalign) and 9 novel and well-established variant calling and filtering methods (Clair3, DeepVariant, Octopus, GATK, FreeBayes, and Strelka2) using a set of 14 "gold standard" WES and WGS datasets available from Genome In A Bottle (GIAB) consortium. Additionally, we have indirectly evaluated each pipeline's performance using a set of 6 non-GIAB samples of African and Russian ethnicity. In our benchmark, Bowtie2 performed significantly worse than other aligners, suggesting it should not be used for medical variant calling. When other aligners were considered, the accuracy of variant discovery mostly depended on the variant caller and not the read aligner. Among the tested variant callers, DeepVariant consistently showed the best performance and the highest robustness. Other actively developed tools, such as Clair3, Octopus, and Strelka2, also performed well, although their efficiency had greater dependence on the quality and type of the input data. We have also compared the consistency of variant calls in GIAB and non-GIAB samples. With few important caveats, best-performing tools have shown little evidence of overfitting. Conclusions: The results show surprisingly large differences in the performance of cutting-edge tools even in high confidence regions of the coding genome. This highlights the importance of regular benchmarking of quickly evolving tools and pipelines. We also discuss the need for a more diverse set of gold standard genomes that would include samples of African, Hispanic, or mixed ancestry. Additionally, there is also a need for better variant caller assessment in the repetitive regions of the coding genome. [ABSTRACT FROM AUTHOR]

Familial hypercholesterolemia (FH) is caused by mutations in various genes, including the LDLR, APOB and PSCK9 genes; however, the spectrum of these mutations in Russian individuals has not been fully investigated. In the present study, mutation screening was performed on the LDLR gene and other FH‑associated genes in patients with definite or possible FH, using next‑generation sequencing. In total, 59 unrelated patients were recruited and sorted into two separate groups depending on their age: Adult (n=31; median age, 49; age range, 23‑70) and children/adolescent (n=28; median age, 11; age range, 2‑21). FH‑associated variants were identified in 18 adults and 25 children, demonstrating mutation detec‑ tion rates of 58 and 89% for the adult and children/adolescent groups, respectively. In the adult group, 13 patients had FH‑associated mutations in the LDLR gene, including two novel variants [NM_000527.4: c.433_434dupG p.(Val145Glyfs*35) and c.1186G>C p.(Gly396Arg)], 3 patients had APOB mutations and two had ABCG5/G8 mutations. In the chil‑ dren/adolescent group, 21 patients had FH‑causing mutations in the LDLR gene, including five novel variants [NM_000527.4: c.325T>G p.(Cys109Gly), c.401G>C p.(Cys134Ser), c.616A>C p.(Ser206Arg), c.1684_1691delTGGCCCAA p.(Pro563Hisfs*14) and c.940+1_c.940+4delGTGA], and 2 patients had APOB mutations, as well as ABCG8 and LIPA mutations, being found in different patients. The present study reported seven novel LDLR variants considered to be patho‑ genic or likely pathogenic. Among them, four missense variants were located in the coding regions, which corresponded to functional protein domains, and two frameshifts were identified that produced truncated proteins. These variants were observed only once in different patients, whereas a splicing variant in intron 6 (c.940+1_c.940+4delGTGA) was detected in four unre‑ lated individuals. Previously reported variants in the LDLR, APOB, ABCG5/8 and LIPA genes were observed in 33 patients. The LDLR p.(Gly592Glu) variant was detected in 6 patients, representing 10% of the FH cases reported in the present study, thus it may be a major variant present in the Russian population. In conclusion, the present study identified seven novel variants of the LDLR gene and broadens the spectrum of mutations in FH‑related genes in the Russian Federation. [ABSTRACT FROM AUTHOR]

The present study reports on the frequency and the spectrum of genetic variants causative of monogenic diabetes in Russian children with non-type 1 diabetes mellitus. The present study included 60 unrelated Russian children with non-type 1 diabetes mellitus diagnosed before the age of 18 years. Genetic variants were screened using whole-exome sequencing (WES) in a panel of 35 genes causative of maturity onset diabetes of the young (MODY) and transient or permanent neonatal diabetes. Verification of the WES results was performed using PCR-direct sequencing. A total of 38 genetic variants were identified in 33 out of 60 patients (55%). The majority of patients (27/33, 81.8%) had variants in MODY-related genes: GCK (n=19), HNF1A (n=2), PAX4 (n=1), ABCC8 (n=1), KCNJ11 (n=1), GCK+HNF1A (n=1), GCK+BLK (n=1) and GCK+BLK+WFS1 (n=1). A total of 6 patients (6/33, 18.2%) had variants in MODY-unrelated genes: GATA6 (n=1), WFS1 (n=3), EIF2AK3 (n=1) and SLC19A2 (n=1). A total of 15 out of 38 variants were novel, including GCK, HNF1A, BLK, WFS1, EIF2AK3 and SLC19A2. To summarize, the present study demonstrates a high frequency and a wide spectrum of genetic variants causative of monogenic diabetes in Russian children with non-type 1 diabetes mellitus. The spectrum includes previously known and novel variants in MODY-related and unrelated genes, with multiple variants in a number of patients. The prevalence of GCK variants indicates that diagnostics of monogenic diabetes in Russian children may begin with testing for MODY2. However, the remaining variants are present at low frequencies in 9 different genes, altogether amounting to ~50% of the cases and highlighting the efficiency of using WES in non-GCK-MODY cases. [ABSTRACT FROM AUTHOR]

A comparative analysis of whole genome sequencing (WGS) and genotype calling was initiated for ten human genome samples sequenced by St. Petersburg State University Peterhof Sequencing Center and by three commercial sequencing centers outside of Russia. The sequence quality, efficiency of DNA variant and genotype calling were compared with each other and with DNA microarrays for each of ten study subjects. We assessed calling of SNPs, indels, copy number variation, and the speed of WGS throughput promised. Twenty separate QC analyses showed high similarities among the sequence quality and called genotypes. The ten genomes tested by the centers included eight American patients afflicted with autoimmune hepatitis (AIH), plus one case’s unaffected parents, in a prelude to discovering genetic influences in this rare disease of unknown etiology. The detailed internal replication and parallel analyses allowed the observation of two of eight AIH cases carrying a rare allele genotype for a previously described AIH-associated gene (FTCD), plus multiple occurrences of known HLA-DRB1 alleles associated with AIH (HLA-DRB1-03:01:01, 13:01:01 and 7:01:01). We also list putative SNVs in other genes as suggestive in AIH influence. [ABSTRACT FROM AUTHOR]

The purpose of this study was to investigate the genetic mechanisms of the defense vascular reactions in response to the diving reflex in humans with polymorphisms in the genes ADBR2, ACE, AGTR1, BDKRB2, and REN. We hypothesized that protective vascular reactions, in response to the diving reflex, are genetically determined and are distinguished in humans with gene polymorphisms of the renin-angiotensin and kinin-bradykinin system. A total of 80 subjects (19 ± 1.4 yr) participated in the study. The intensity of the vascular response was estimated using photoplethysmogram. The I/D polymorphism (rs4340) of ACE was analyzed by PCR. REN (G/A, rs2368564), AGTR1 (A/C, rs5186), BDKRB2 (T/C, rs1799722), and ADBR2 (A/G, rs1042713) polymorphisms were examined using the two-step multiplex PCR followed by carrying allele hybridization on the biochip. Subjects with the BDKRB2 (C/C), ACE (D/D), and ADBR2 (G/G, G/A) genotypes exhibited the strongest peripheral vasoconstriction in response to diving. In subjects with a combination of the BDKRB2 (C/C) plus ACE (D/D) genotypes, we observed the lowest pulse wave amplitude and pulse transit time values and the highest arterial blood pressure during face immersion compared with the heterozygous individuals, suggesting that these subjects are more susceptible to diving hypoxia. This study observed that humans with gene polymorphisms of the renin-angiotensin and kinin-bradykinin systems demonstrate various expressions of protective vascular reactions in response to the diving reflex. The obtained results might be used in estimation of resistance to hypoxia of any origin in human beings or in a medical practice. NEW & NOTEWORTHY Our study demonstrates that the vascular reactions in response to the diving reflex are genetically determined and depend on gene polymorphisms of the kinin-bradykinin and the renin-angiotensin systems. [ABSTRACT FROM AUTHOR]

New Findings What is the central question of this study? Variations in genes are considered to be molecular determinants maintaining the expression of the slow or fast myosin heavy chains of adult skeletal muscle. The role of polymorphisms of candidate genes involved in skeletal muscle development, energy homeostasis and thyroid and calcium metabolism in the determination of muscle fibre type has not previously been reported., What is the main finding and its importance? We show that the AGTR2 rs11091046 C allele is associated with an increased proportion of slow-twitch muscle fibres, endurance athlete status and aerobic performance. Such findings have important implications for our understanding of muscle function in both health and disease., Muscle fibre type is a heritable trait and can partly predict athletic success. It has been proposed that polymorphisms of genes involved in the regulation of muscle fibre characteristics may predispose the muscle precursor cells of a given individual to be predominantly fast or slow. In the present study, we examined the association between 15 candidate gene polymorphisms and muscle fibre type composition of the vastus lateralis muscle in 55 physically active, healthy men. We found that rs11091046 C allele carriers of the angiotensin II type 2 receptor gene ( AGTR2; involved in skeletal muscle development, metabolism and circulatory homeostasis) had a significantly higher percentage of slow-twitch fibres than A allele carriers [54.2 (11.1) versus 45.2 (10.2)%; P = 0.003]. These data indicate that 15.2% of the variation in muscle fibre composition of the vastus lateralis muscle can be explained by the AGTR2 genotype. Next, we investigated the frequencies of the AGTR2 alleles in 2178 Caucasian athletes and 1220 control subjects. The frequency of the AGTR2 C allele was significantly higher in male and female endurance athletes compared with power athletes (males, 62.7 versus 51.7%, P = 0.0038; females, 56.6 versus 48.1%, P = 0.0169) and control subjects (males, 62.7 versus 51.0%, P = 0.0006; elite female athletes, 65.1 versus 55.2%, P = 0.0488). Furthermore, the frequency of the AGTR2 A allele was significantly over-represented in female power athletes (51.9%) in comparison to control subjects (44.8%, P = 0.0069). We also found that relative maximal oxygen consumption was significantly greater in male endurance athletes with the AGTR2 C allele compared with AGTR2 A allele carriers [ n = 28; 62.3 (4.4) versus 57.4 (6.0) ml min−1 kg−1; P = 0.0197]. Taken together, these results demonstrate that the AGTR2 gene C allele is associated with an increased proportion of slow-twitch muscle fibres, endurance athlete status and aerobic performance, while the A allele is associated with a higher percentage of fast-twitch fibres and power-oriented disciplines. [ABSTRACT FROM AUTHOR]

Although circulating microRNAs (miRNAs) in maternal blood may play an important role in regulation of pregnancy progression and serve as non-invasive biomarkers for different gestation complications, little is known about their profile in blood during normally developing pregnancy. In this study we evaluated the miRNA profiles in paired plasma and serum samples from pregnant women without health or gestational abnormalities at three time points using high-throughput sequencing technology. Sequencing revealed that the percentage of miRNA reads in plasma and serum decreased by a third compared to first and second trimesters. We found two miRNAs in plasma (hsa-miR-7853-5p and hsa-miR-200c-3p) and 10 miRNAs in serum (hsa-miR-203a-5p, hsa-miR-495-3p, hsa-miR-4435, hsa-miR-340-5p, hsa-miR-4417, hsa-miR-1266-5p, hsa-miR-4494, hsa-miR-134-3p, hsa-miR-5008-5p, and hsa-miR-6756-5p), that exhibit level changes during pregnancy (p-value adjusted < 0.05). In addition, we observed differences for 36 miRNAs between plasma and serum (p-value adjusted < 0.05), which should be taken into consideration when comparing the results between studies performed using different biosample types. The results were verified by analysis of three miRNAs using qRT-PCR (p < 0.05). The present study confirms that the circulating miRNA profile in blood changes during gestation. Our results set the basis for further investigation of molecular mechanisms, involved in regulation of pregnancy, and the search for biomarkers of gestation abnormalities. [ABSTRACT FROM AUTHOR]

Type 2 diabetes mellitus (T2D) is a chronic metabolic disease resulting from insulin resistance and progressively reduced insulin secretion, which leads to impaired glucose utilization, dyslipidemia and hyperinsulinemia and progressive pancreatic beta cell dysfunction. The incidence of type 2 diabetes mellitus is increasing worldwide and nowadays T2D already became a global epidemic. The well-known interindividual variability of T2D drug actions such as biguanides, sulfonylureas/meglitinides, DPP-4 inhibitors/GLP1R agonists and SGLT-2 inhibitors may be caused, among other things, by genetic factors. Pharmacogenetic findings may aid in identifying new drug targets and obtaining in-depth knowledge of the causes of disease and its physiological processes, thereby, providing an opportunity to elaborate an algorithm for tailor or precision treatment. The aim of this article is to summarize recent progress and discoveries for T2D pharmacogenetics and to discuss the factors which limit the furthering accumulation of genetic variability knowledge in patient response to therapy that will allow improvement the personalized treatment of T2D. [ABSTRACT FROM AUTHOR]

Over the recent years, many advances have been made in the research of the genetic factors of pregnancy complications. In this work, we use publicly available data repositories, such as the National Human Genome Research Institute GWAS Catalog, HuGE Navigator, and the UK Biobank genetic and phenotypic dataset to gain insights into molecular pathways and individual genes behind a set of pregnancy-related traits, including the most studied ones—preeclampsia, gestational diabetes, preterm birth, and placental abruption. Using both HuGE and GWAS Catalog data, we confirm that immune system and, in particular, T-cell related pathways are one of the most important drivers of pregnancy-related traits. Pathway analysis of the data reveals that cell adhesion and matrisome-related genes are also commonly involved in pregnancy pathologies. We also find a large role of metabolic factors that affect not only gestational diabetes, but also the other traits. These shared metabolic genes include IGF2, PPARG, and NOS3. We further discover that the published genetic associations are poorly replicated in the independent UK Biobank cohort. Nevertheless, we find novel genome-wide associations with pregnancy-related traits for the FBLN7, STK32B, and ACTR3B genes, and replicate the effects of the KAZN and TLE1 genes, with the latter being the only gene identified across all data resources. Overall, our analysis highlights central molecular pathways for pregnancy-related traits, and suggests a need to use more accurate and sophisticated association analysis strategies to robustly identify genetic risk factors for pregnancy complications. [ABSTRACT FROM AUTHOR]

Preeclampsia (PE) is a multisystem heterogeneous complication of pregnancy remaining a leading cause of maternal and perinatal morbidity and mortality over the world. PE has a large spectrum of clinical features and symptoms, which make diagnosis challenging. Despite a long period of studying, PE etiology is still unclear and there are no reliable rapid tests for early diagnosis of this disease. During the last decade, it was shown that proteins misfolding and aggregation are associated with PE. Several proteins, including amyloid beta peptide, transthyretin, alpha-1 antitrypsin, albumin, IgG k-free light chains, and ceruloplasmin are dysregulated in PE, resulting in toxic deposition of amyloid-like aggregates in the placenta and body fluids. It is also possible that aggregated proteins induce defective trophoblast invasion, placental ischemia, ER stress, and promote PE manifestation. The fact that protein aggregation is an emerging biomarker of PE provides an opportunity to develop new diagnostic approaches based on amyloids special features, such as Congo red (CR) staining and thioflavin T (ThT) enhanced fluorescence. [ABSTRACT FROM AUTHOR]

Background: Allele frequency data from large exome and genome aggregation projects such as the Genome Aggregation Database (gnomAD) are of ultimate importance to the interpretation of medical resequencing data. However, allele frequencies might significantly differ in poorly studied populations that are underrepresented in large‐scale projects, such as the Russian population. Methods: In this work, we leveraged our access to a large dataset of 694 exome samples to analyze genetic variation in the Northwest Russia. We compared the spectrum of genetic variants to the dbSNP build 151, and made estimates of ClinVar‐based autosomal recessive (AR) disease allele prevalence as compared to gnomAD r. 2.1. Results: An estimated 9.3% of discovered variants were not present in dbSNP. We report statistically significant overrepresentation of pathogenic variants for several Mendelian disorders, including phenylketonuria (PAH, rs5030858), Wilson's disease (ATP7B, rs76151636), factor VII deficiency (F7, rs36209567), kyphoscoliosis type of Ehlers‐Danlos syndrome (FKBP14, rs542489955), and several other recessive pathologies. We also make primary estimates of monogenic disease incidence in the population, with retinal dystrophy, cystic fibrosis, and phenylketonuria being the most frequent AR pathologies. Conclusion: Our observations demonstrate the utility of population‐specific allele frequency data to the diagnosis of monogenic disorders using high‐throughput technologies. [ABSTRACT FROM AUTHOR]

Type 2 diabetes (T2D) and obesity are common chronic disorders with multifactorial etiology. In our study, we performed an exome sequencing analysis of 110 patients of Russian ethnicity together with a multi-perspective approach based on biologically meaningful filtering criteria to detect novel candidate variants and loci for T2D and obesity. We have identified several known single nucleotide polymorphisms (SNPs) as markers for obesity (rs11960429), T2D (rs9379084, rs1126930), and body mass index (BMI) (rs11553746, rs1956549 and rs7195386) (p < 0.05). We show that a method based on scoring of case-specific variants together with selection of protein-altering variants can allow for the interrogation of novel and known candidate markers of T2D and obesity in small samples. Using this method, we identified rs328 in LPL (p = 0.023), rs11863726 in HBQ1 (p = 8 × 10−5), rs112984085 in VAV3 (p = 4.8 × 10−4) for T2D and obesity, rs6271 in DBH (p = 0.043), rs62618693 in QSER1 (p = 0.021), rs61758785 in RAD51B (p = 1.7 × 10−4), rs34042554 in PCDHA1 (p = 1 × 10−4), and rs144183813 in PLEKHA5 (p = 1.7 × 10−4) for obesity; and rs9379084 in RREB1 (p = 0.042), rs2233984 in C6orf15 (p = 0.030), rs61737764 in ITGB6 (p = 0.035), rs17801742 in COL2A1 (p = 8.5 × 10−5), and rs685523 in ADAMTS13 (p = 1 × 10−6) for T2D as important susceptibility loci in Russian population. Our results demonstrate the effectiveness of whole exome sequencing (WES) technologies for searching for novel markers of multifactorial diseases in cohorts of limited size in poorly studied populations. [ABSTRACT FROM AUTHOR]