Meta-analysis of Transcriptomic Data from Lung Autopsy and Cellular Models of SARS-CoV-2 Infection.

Severe COVID-19 is a systemic disorder involving excessive inflammatory response, metabolic dysfunction, multi-organ damage, and several clinical features. Here, we performed a transcriptome meta-analysis investigating genes and molecular mechanisms related to COVID-19 severity and outcomes. First, transcriptomic data of cellular models of SARS-CoV-2 infection were compiled to understand the first response to the infection. Then, transcriptomic data from lung autopsies of patients deceased due to COVID-19 were compiled to analyze altered genes of damaged lung tissue. These analyses were followed by functional enrichment analyses and gene-phenotype association. A biological network was constructed using the disturbed genes in the lung autopsy meta-analysis. Central genes were defined considering closeness and betweenness centrality degrees. A sub-network phenotype-gene interaction analysis was performed. The meta-analysis of cellular models found genes mainly associated with cytokine signaling and other pathogen response pathways. The meta-analysis of lung autopsy tissue found genes associated with coagulopathy, lung fibrosis, multi-organ damage, and long COVID-19. Only genes DNAH9 and FAM216B were found perturbed in both meta-analyses. BLNK, FABP4, GRIA1, ATF3, TREM2, TPPP, TPPP3, FOS, ALB, JUNB, LMNA, ADRB2, PPARG, TNNC1, and EGR1 were identified as central elements among perturbed genes in lung autopsy and were found associated with several clinical features of severe COVID-19. Central elements were suggested as interesting targets to investigate the relation with features of COVID-19 severity, such as coagulopathy, lung fibrosis, and organ damage.

Possible New Candidates Involved to Thalidomide-Related Limbs and Cardiac Defects: A Systems Biology Approach.

Thalidomide is a known teratogen that causes malformations especially in heart and limbs. Its mechanism of teratogenicity is still not fully elucidated. Recently, a new target of thalidomide was described, TBX5, and was observed a new interaction between HAND2 and TBX5 that is disrupted in the presence of thalidomide. Therefore, our study aimed to raise potential candidates for thalidomide teratogenesis, through systems biology, evaluating HAND2 and TBX5 interaction and heart and limbs malformations of thalidomide. Genes and proteins related to TBX5 and HAND2 were selected through TF2DNA, REACTOME, Human Phenotype Ontology, and InterPro databases. Networks were assembled using STRING © database. Network analysis were performed in Cytoscape © and R v3.6.2. Differential gene expression (DGE) analysis was performed through gene expression omnibus. We constructed a network for HAND2 and TBX5 interaction; a network for heart and limbs malformations of TE; and the two joined networks. We observed that EP300 protein seemed to be important in all networks. We also looked for proteins containing C2H2 domain in the assembled networks. ZIC3, GLI1, GLI3, ZNF148, and PRDM16 were the ones present in both heart and limbs malformations of TE networks. Furthermore, in the DGE analysis after treatment with thalidomide, we observed that FANCB, ESCO2, and XRCC2 were downregulated and present both in heart and limbs networks. Through systems biology, we were able to point to different new proteins and genes, and selected specially EP300, which was important in all the analyzed networks, to be further evaluated in the TE teratogenicity.

Brazilian women in Bioinformatics: Challenges and opportunities.

Bioinformatics is a growing research field that received great notoriety in the years of the COVID-19 pandemic. It is a very integrative area, comprising professionals from science, technology, engineering, and mathematics (STEM). In agreement with the other STEM areas, several women have greatly contributed to bioinformatics ascension; however, they had to surpass prejudice and stereotypes to achieve recognition and leadership positions, a path that studies have demonstrated to be more comfortable to their male colleagues. In this review, we discuss the several difficulties that women in STEM, including bioinformatics, surpass during their careers. First, we present a historical context on bioinformatics and the main applications for this area. Then, we discuss gender disparity in STEM and present the challenges that still contribute to women's inequality in STEM compared to their male colleagues. We also present the opportunities and the transformation that we can start, acting in academia, inside the family and school environments, and as a society, hence contributing to gender equality in STEM. Finally, we discuss specific challenges in the bioinformatics field and how we can act to overcome them, especially in low and middle-income countries, such as Brazil.

Exploring the frequency of a TP53 polyadenylation signal variant in tumor DNA from patients diagnosed with lung adenocarcinomas, sarcomas and uterine leiomyomas.

The TP53 3'UTR variant rs78378222 A>C has been detected in different tumor types as a somatic alteration that reduces p53 expression through modification of polyadenylation and miRNA regulation. Its prevalence is not yet known in all tumors. Herein, we examine tumor tissue prevalence of rs7837822 in Brazilian cohorts of patients from south and southeast regions diagnosed with lung adenocarcinoma (LUAD, n=586), sarcoma (SARC, n=188) and uterine leiomyoma (ULM, n=41). The minor allele (C) was identified in heterozygosity in 6/586 LUAD tumors (prevalence = 1.02 %) and none of the SARC and ULM samples. Additionally, next generation sequencing analysis revealed that all variant-positive tumors (n=4) with sample availability had additional pathogenic or likely pathogenic somatic variants in the TP53 coding regions. Among them, 3/4 (75 %) had the same pathogenic or likely pathogenic sequence variant (allele frequency <0.05 in tumor DNA) namely c.751A>C (p.Ile251Leu). Our results indicate a low somatic prevalence of rs78378222 in LUAD, ULM and SARC tumors from Brazilian patients, which suggests that no further analysis of this variant in the specific studied regions of Brazil is warranted. However, these findings should not exclude tumor molecular testing of this TP53 3'UTR functional variant for different populations.

A New Strategy for the Old Challenge of Thalidomide: Systems Biology Prioritization of Potential Immunomodulatory Drug (IMiD)-Targeted Transcription Factors.

Several molecular mechanisms of thalidomide embryopathy (TE) have been investigated, from anti-angiogenesis to oxidative stress to cereblon binding. Recently, it was discovered that thalidomide and its analogs, named immunomodulatory drugs (IMiDs), induced the degradation of C2H2 transcription factors (TFs). This mechanism might impact the strict transcriptional regulation of the developing embryo. Hence, this study aims to evaluate the TFs altered by IMiDs, prioritizing the ones associated with embryogenesis through transcriptome and systems biology-allied analyses. This study comprises only the experimental data accessed through bioinformatics databases. First, proteins and genes reported in the literature as altered/affected by the IMiDs were annotated. A protein systems biology network was evaluated. TFs beta-catenin (CTNNB1) and SP1 play more central roles: beta-catenin is an essential protein in the network, while SP1 is a putative C2H2 candidate for IMiD-induced degradation. Separately, the differential expressions of the annotated genes were analyzed through 23 publicly available transcriptomes, presenting 8624 differentially expressed genes (2947 in two or more datasets). Seventeen C2H2 TFs were identified as related to embryonic development but not studied for IMiD exposure; these TFs are potential IMiDs degradation neosubstrates. This is the first study to suggest an integration of IMiD molecular mechanisms through C2H2 TF degradation.

Revealing the expression profile of genes that encode the Subcortical Maternal Complex in human reproductive failures.

The Subcortical Maternal Complex (SCMC) is composed of maternally encoded proteins required for the early stages of embryo development. Here we aimed to investigate the expression profile of the genes that encode the individual members of the SCMC in human reproductive failures. To accomplish that, we selected three datasets in the Gene Expression Omnibus repository for differential gene expression (DGE) analysis, comprising human endometrial and placental tissues of patients with recurrent implantation failure (RIF) or recurrent pregnancy loss (RPL). The SCMC genes KHDC3L, NLRP2, NLRP4, NLRP5, OOEP, PADI6, TLE6, and ZBED3 were included in the DGE analysis, as well as CFL1 and CFL2 that connect the SCMC with the actin cytoskeleton. Additionally, differential co-expression analysis and systems biology analysis of gene-gene co-expression were performed for KHDC3L, NLRP5, OOEP, and TLE6, demonstrating gene pairs differentially correlated under the two conditions, and the co-expression with genes involved in immune response, cell cycle, DNA damage repair, embryo development, and male reproduction. Compared to control groups, NLRP5 demonstrated upregulation in the endometrium of RIF patients, and KHDC3L was upregulated in the fetal placental tissue of RPL patients, shedding light on the importance of considering SCMC genes in reproductive failures.

Evaluation of the influence of genetic variants in Cereblon gene on the response to the treatment of erythema nodosum leprosum with thalidomide.

BACKGROUND: Erythema nodosum leprosum (ENL) is an acute and systemic inflammatory reaction of leprosy characterised by painful nodules and involvement of various organs. Thalidomide is an immunomodulatory and anti-inflammatory drug currently used to treat this condition. Cereblon (CRBN) protein is the primary target of thalidomide, and it has been pointed out as necessary for the efficacy of this drug in others therapeutics settings. OBJECTIVES: In this study, we aimed to evaluate the influence of CRBN gene variants on the dose of thalidomide as well as its adverse effects during treatment of ENL. METHODS: A total of 103 ENL patients in treatment with thalidomide were included in this study. DNA samples were obtained from saliva and molecular analysis of CRBN gene were performed to investigate the variants rs1620675, rs1672770 and rs4183. Different genotypes of CRBN variants were evaluated in relation to their influence on the dose of thalidomide and on the occurrence of adverse effects. FINDINGS: No association was found between CRBN variants and thalidomide dose variation. However, the genotypes of rs1672770 showed association with gastrointestinal effects (p = 0.040). Moreover, the haplotype DEL/C/T (rs4183/rs1672770/rs1620675) was also associated with gastrointestinal adverse effects (p = 0.015). MAIN CONCLUSIONS: Our results show that CRBN variants affect the treatment of ENH with thalidomide, especially on the adverse effects related to the drug.

Angiogenesis and oxidative stress-related gene variants in recurrent pregnancy loss.

Recurrent pregnancy loss (RPL) affects ~3-5% of couples attempting to conceive and in around 50% of cases the aetiology remains unknown. Adequate vascularisation and placental circulation are indispensable for the development of a normal pregnancy. Prostaglandin-endoperoxide synthase 2 (PTGS2), vascular endothelial growth factor (VEGF) and the nitric oxide (NO) systems play important roles in reproductive physiology, participating in several steps including implantation and apoptosis of trophoblast cells. In this study we evaluated genetic polymorphisms in the inducible nitric oxide synthase (NOS2), PTGS2 and VEGFA genes as susceptibility factors for RPL. A case-control study was conducted in 149 women having two or more miscarriages and 208 controls. Allele and genotype distributions of the polymorphisms studied in the two groups were not statistically different. However, the dominant model showed that the presence of variant T (TT/GT) of rs2779249 (-1290G>T) of NOS2 was significantly associated with RPL (OR=1.58, CI 95%=1.03-2.44; P=0.037). The increased risk remained significant when adjusted for number of pregnancies, alcohol consumption and ethnicity (OR=1.92, CI95%=1.18-3.11; P=0.008). These results suggest that the variant genotypes of the functional polymorphism rs2779249 in the NOS2 promoter are a potential risk for RPL, possibly due to oxidative stress mechanisms.

Music genetics research: Association with musicality of a polymorphism in the AVPR1A gene.

Musicality is defined as a natural tendency, sensibility, knowledge, or talent to create, perceive, and play music. Musical abilities involve a great range of social and cognitive behaviors, which are influenced by both environmental and genetic factors. Although a number of studies have yielded insights into music genetics research, genes and biological pathways related to these traits are not fully understood. Our hypothesis in the current study is that genes associated with different behaviors could also influence the musical phenotype. Our aim was to investigate whether polymorphisms in six genes (AVPR1A, SLC6A4, ITGB3, COMT, DRD2 and DRD4) related to social and cognitive traits are associated with musicality in a sample of children. Musicality was assessed through an individualized music therapy assessment profile (IMTAP) which has been validated in Brazil to measure musical ability. We show here that the RS1 microsatellite of the AVPR1A gene is nominally associated with musicality, corroborating previous results linking AVPR1A with musical activity. This study is one of the first to investigate musicality in a comprehensive way, and it contributes to better understand the genetic basis underlying musical ability.

Thalidomide embryopathy: Follow-up of cases born between 1959 and 2010.

BACKGROUND: Thalidomide is a known teratogen and it is estimated that more than ten thousand babies were affected by thalidomide embryopathy (TE), which is characterized mainly by limb defects, but can involve many organs and systems. Most people with TE were only evaluated at birth and it is not well established if thalidomide exposure during embryonic development leads to later effects. We analyzed the clinical history of adults with TE to better understand this gap in the clinical findings of TE. METHODS: Brazilian individuals with TE were invited to answer a clinical questionnaire which considered family history, social information, medical history, and current clinical and psychological health status. A clinical examination was also performed, including on the infant subjects to evaluate congenital anomalies. The characterization of the features was analyzed using descriptive statistics and Chi-square or Fisher's exact test. RESULTS: The congenital anomalies caused by thalidomide were reviewed in 28 Brazilian individuals, and the questionnaire was applied to the 23 adult subjects with TE (aged 19 to 55). Progressive deafness and dental loss were reported. From the comparison of TE individuals with the general Brazilian population, the early onset of cardiovascular diseases (p = 0.009) and a higher frequency of psychological disorders (p = 0.011) were observed. CONCLUSION: Although there is no sufficient evidence that thalidomide exposure caused or worsened the described events, this approach helps to better understand the TE phenotype, improves the clinical diagnosis, and can lead to adequate health support for these individuals.

Bioinformatics Methods for Transcriptome Analysis on Teratogenesis Testing.

Teratogenesis testing can be challenging due to the limitations of both in vitro and in vivo models. Test-systems, based especially on human embryonic cells, have been helping to overcome the difficulties when allied to omics strategies, such as transcriptomics. In these test-systems, cells exposed to different compounds are then analyzed in microarray or RNA-seq platforms regarding the impacts of the potential teratogens in the gene expression. Nevertheless, microarray and RNA-seq dataset processing requires computational resources and bioinformatics knowledge. Here, a pipeline for microarray and RNA-seq processing is presented, aiming to help researchers from any field to interpret the main transcriptome results, such as differential gene expression, enrichment analysis, and statistical interpretation. This chapter also discusses the main difficulties that can be encountered in a transcriptome analysis and the better alternatives to overcome these issues, describing both programming codes and user-friendly tools. Finally, specific issues in the teratogenesis field, such as time-course analysis, are also described, demonstrating how the pipeline can be applied in these studies.

Serum myostatin as a candidate disease severity and progression biomarker of spinal muscular atrophy.

The identification of biomarkers for spinal muscular atrophy is crucial for predicting disease progression, severity, and response to new disease-modifying therapies. This study aimed to investigate the role of serum levels of myostatin and follistatin as biomarkers for spinal muscular atrophy, considering muscle atrophy secondary to denervation as the main clinical manifestation of the disease. The study evaluated the differential gene expression of myostatin and follistatin in a lesional model of gastrocnemius denervation in mice, as well as in a meta-analysis of three datasets in transgenic mice models of spinal muscular atrophy, and in two studies involving humans with spinal muscular atrophy. Subsequently, a case-control study involving 27 spinal muscular atrophy patients and 27 controls was conducted, followed by a 12-month cohort study with 25 spinal muscular atrophy cases. Serum levels of myostatin and follistatin were analysed using enzyme-linked immunosorbent assay at a single centre in southern Brazil. Skeletal muscle gene expression of myostatin decreased and of follistatin increased following lesional muscle denervation in mice, consistent with findings in the spinal muscular atrophy transgenic mice meta-analysis and in the iliopsoas muscle of five patients with spinal muscular atrophy type 1. Median serum myostatin levels were significantly lower in spinal muscular atrophy patients (98 pg/mL; 5-157) compared to controls (412 pg/mL; 299-730) (P < 0.001). Lower myostatin levels were associated with greater disease severity based on clinician-rated outcomes (Rho = 0.493-0.812; P < 0.05). After 12 months, there was a further reduction in myostatin levels among spinal muscular atrophy cases (P = 0.021). Follistatin levels did not differ between cases and controls, and no significant changes were observed over time. The follistatin:myostatin ratio was significantly increased in spinal muscular atrophy subjects and inversely correlated with motor severity. Serum myostatin levels show promise as a novel biomarker for evaluating the severity and progression of spinal muscular atrophy. The decrease in myostatin levels and the subsequent favourable environment for muscle growth may be attributed to denervation caused by motor neuron dysfunction.

New candidate genes potentially involved in Zika virus teratogenesis.

Despite efforts to elucidate Zika virus (ZIKV) teratogenesis, still several issues remain unresolved, particularly on the molecular mechanisms behind the pathogenesis of Congenital Zika Syndrome (CZS). To answer this question, we used bioinformatics tools, animal experiments and human gene expression analysis to investigate genes related to brain development potentially involved in CZS. Searches in databases for genes related to brain development and CZS were performed, and a protein interaction network was created. The expression of these genes was analyzed in a CZS animal model and secondary gene expression analysis (DGE) was performed in human cells exposed to ZIKV. A total of 2610 genes were identified in the databases, of which 1013 were connected. By applying centrality statistics of the global network, 36 candidate genes were identified, which, after selection resulted in nine genes. Gene expression analysis revealed distinctive expression patterns for PRKDC, PCNA, ATM, SMC3 as well as for FGF8 and SHH in the CZS model. Furthermore, DGE analysis altered expression of ATM, PRKDC, PCNA. In conclusion, systems biology are helpful tools to identify candidate genes to be validated in vitro and in vivo. PRKDC, PCNA, ATM, SMC3, FGF8 and SHH have altered expression in ZIKV-induced brain malformations.

Expression profiles of meiotic genes in male vs. female gonads and gametes: Insights into fertility issues.

Gametes are specialized cells that, at fertilization, give rise to a totipotent zygote capable of generating an entire organism. Female and male germ cells undergo meiosis to produce mature gametes; however, sex-specific events of oogenesis and spermatogenesis contribute to specific roles of gametes in reproductive issues. We investigate the differential gene expression (DGE) of meiosis-related genes in human female and male gonads and gametes in normal and pathological conditions. The transcriptome data for the DGE analysis was obtained through the Gene Expression Omnibus repository, comprising human ovary and testicle samples of the prenatal period and adulthood, additionally to male (non-obstructive azoospermia (NOA) and teratozoospermia), and female (polycystic ovary syndrome (PCOS) and advanced maternal age) reproductive conditions. Gene ontology terms related to meiosis were associated with 678 genes, of which 17 genes in common were differentially expressed between the testicle and ovary during the prenatal period and adulthood. Except for SERPINA5 and SOX9, the 17 meiosis-related genes were downregulated in the testicle during the prenatal period and upregulated in adulthood compared to the ovary. No differences were observed in the oocytes of PCOS patients; however, meiosis-related genes were differentially expressed according to the patient's age and maturity of the oocyte. In NOA and teratozoospermia, 145 meiosis-related genes were differentially expressed in comparison to the control, including OOEP; despite no recognized role in male reproduction, OOEP was co-expressed with genes related to male fertility. Taking together, these results shed light on potential genes that might be relevant to comprehend human fertility disorders.

Transcriptome meta-analysis of valproic acid exposure in human embryonic stem cells.

Valproic acid (VPA) is a widely used antiepileptic drug not recommended in pregnancy because it is teratogenic. Many assays have assessed the impact of the VPA exposure on the transcriptome of human embryonic stem-cells (hESC), but the molecular perturbations that VPA exerts in neurodevelopment are not completely understood. This study aimed to perform a transcriptome meta-analysis of VPA-exposed hESC to elucidate the main biological mechanisms altered by VPA effects on the gene expression. Publicly available microarray and RNA-seq transcriptomes were selected in the Gene Expression Omnibus (GEO) repository. Samples were processed according to the standard pipelines for each technology in the Galaxy server and R. Meta-analysis was performed using the Fisher-P method. Overrepresented genes were obtained by evaluating ontologies, pathways, and phenotypes' databases. The meta-analysis performed in seven datasets resulted in 61 perturbed genes, 54 upregulated. Ontology and pathway enrichments suggested neurodevelopment and neuroinflammatory effects; phenotype overrepresentation included epilepsy-related genes, such as SCN1A and GABRB2. The NDNF gene upregulation was also identified; this gene is involved in neuron migration and survival during development. Sub-network analysis proposed TGFß and BMP pathways activation. These results suggest VPA exerts effects in epilepsy-related genes even in embryonic cells. Neurodevelopmental genes, such as NDNF were upregulated and VPA might also disturb several development pathways. These mechanisms might help to explain the spectrum of VPA-induced congenital anomalies and the molecular effects on neurodevelopment.

Genetic evaluation of HAND2 gene and its effects on thalidomide embryopathy.

BACKGROUND: HAND2 is a transcription factor important for embryonic development, required for limbs and cardiovascular development. Thalidomide is a drug responsible to a spectrum of congenital anomalies known as Thalidomide Embryopathy (TE), which includes mainly limb and heart defects. It is known that HAND2 interaction with TBX5, an important protein for limbs and heart development, is inhibited by Thalidomide. The aim of this study was to evaluate and characterize HAND2 in the context of TE, and to evaluate its variability in TE individuals. METHODS: DNA from 35 TE subjects was extracted from saliva samples and PCR was performed for amplification and Sanger sequencing of HAND2 coding sequence. RESULTS: The analysis showed only one variant; a synonymous variant p.P51 (rs59621536) in exon 1 found in three individuals. Further in silico evaluation confirmed highly HAND2 conservation, being the 3'UTR the most polymorphic region of the gene. Additional computational analyses classified the variant as neutral, without alteration in splicing and miRNA sites. In silico predictions pointed to alteration of two CpG islands adjacent to the variant; however, we did not observe any alterations on the methylation pattern of HAND2 gene in our sample. Moreover, alteration of the binding site of MeCP2, a nuclear protein involved in DNA methylation, was predicted along with alteration in HAND2 mRNA structure. CONCLUSIONS: Considering HAND2 being a well conserved gene, further studies with a larger sample should be performed to evaluate the role this gene on genetic susceptibility to TE.

Neurodevelopment in Children Exposed to Zika in utero: Clinical and Molecular Aspects.

Five years after the identification of Zika virus as a human teratogen, we reviewed the early clinical manifestations, collectively called congenital Zika syndrome (CZS). Children with CZS have a very poor prognosis with extremely low performance in motor, cognitive, and language development domains, and practically all feature severe forms of cerebral palsy. However, these manifestations are the tip of the iceberg, with some children presenting milder forms of deficits. Additionally, neurodevelopment can be in the normal range in the majority of the non-microcephalic children born without brain or eye abnormalities. Vertical transmission and the resulting disruption in development of the brain are much less frequent when maternal infection occurs in the second half of the pregnancy. Experimental studies have alerted to the possibility of other behavioral outcomes both in prenatally infected children and in postnatal and adult infections. Cofactors play a vital role in the development of CZS and involve genetic, environmental, nutritional, and social determinants leading to the asymmetric distribution of cases. Some of these social variables also limit access to multidisciplinary professional treatment.

Comparative Genomics Identifies Putative Interspecies Mechanisms Underlying Crbn-Sall4-Linked Thalidomide Embryopathy.

The identification of thalidomide-Cereblon-induced SALL4 degradation has brought new understanding for thalidomide embryopathy (TE) differences across species. Some questions, however, regarding species variability, still remain. The aim of this study was to detect sequence divergences between species, affected or not by TE, and to evaluate the regulated gene co-expression in a murine model. Here, we performed a comparative analysis of proteins experimentally established as affected by thalidomide exposure, evaluating 14 species. The comparative analysis, regarding synteny, neighborhood, and protein conservation, was performed in 42 selected genes. Differential co-expression analysis was performed, using a publicly available assay, GSE61306, which evaluated mouse embryonic stem cells (mESC) exposed to thalidomide. The comparative analyses evidenced 20 genes in the upstream neighborhood of NOS3, which are different between the species who develop, or not, the classic TE phenotype. Considering protein sequence alignments, RECQL4, SALL4, CDH5, KDR, and NOS2 proteins had the biggest number of variants reported in unaffected species. In co-expression analysis, Crbn was a gene identified as a driver of the co-expression of other genes implicated in genetic, non-teratogenic, limb reduction defects (LRD), such as Tbx5, Esco2, Recql4, and Sall4; Crbn and Sall4 were shown to have a moderate co-expression correlation, which is affected after thalidomide exposure. Hence, even though the classic TE phenotype is not identified in mice, a deregulatory Crbn-induced mechanism is suggested in this animal. Functional studies are necessary, especially evaluating the genes responsible for LRD syndromes and their interaction with thalidomide-Cereblon.

Investigating the role of EGF-CFC gene family in recurrent pregnancy loss through bioinformatics and molecular approaches.

Recurrent pregnancy loss (RPL) is the most common reproductive failure, reaching 1-5% of women throughout their lives, and having unknown etiology in 50% of the cases. In humans, EGF-CFC1 (Epidermal Growth Factors & Cripto/FRL-1/Cryptic) gene family is composed by TDGF1 and CFC1, two developmental genes. The aim of this study was to investigate the role of EGF-CFC on RPL. To this, multiple approaches were performed; we conducted an expression analysis of TDGF1 and CFC1 using publicly available data from Gene Omnibus Expression (GEO), systems biology analyses and functional prediction; and a molecular analysis carried out in a case-control study. Our GEO analysis showed a decrease in TDGF1 expression in the endometrium (p=0.049) and CFC1 expression in placenta (p=0.015) of women with RPL. Network analysis, gene ontology and literature pointed to a strong connection between EGF-CFC1 gene family to pathways that play key roles during pregnancy, including TGF-ß, c-Src/MAPK/AKT, Notch, TNFα, IFNγ and IL-6. A pathogenicity score developed for this gene family showed that the c.-14+1429T>C (rs3806702) variant in the TDGF1 and the p.Arg47Gln (rs201431919) variant in CFC1 gene would be the ones with the highest deleterious effect for RPL. In the case-control study, which involved 149 women with RPL and 159 controls, no statistical difference was observed in the allele and genotype distributions of the variants studied in the two groups. In this study, we performed extensive bioinformatics analysis for biomarker prioritization followed by experimental validation of proposed selected markers. Although there is no statistical difference in the frequencies of these variants between RPL and controls, the expression analysis results suggest that TDGF1 and CFC1 genes might play a role in RPL. In addition, systems biology analyzes raise the hypothesis that genes in other signaling pathways that may be related to RPL as good candidates for future studies.Abbreviations RPL: recurrent pregnancy loss; EGF-CFC1: Epidermal Growth Factors - Cripto/FRL-1; GEO: Gene Omnibus Expression; KEGG: Kyoto Encyclopedia of Genes and Genomes.

Genetic Susceptibility to Drug Teratogenicity: A Systematic Literature Review.

Since the 1960s, drugs have been known to cause teratogenic effects in humans. Such teratogenicity has been postulated to be influenced by genetics. The aim of this review was to provide an overview of the current knowledge on genetic susceptibility to drug teratogenicity in humans and reflect on future directions within the field of genetic teratology. We focused on 12 drugs and drug classes with evidence of teratogenic action, as well as 29 drugs and drug classes with conflicting evidence of fetal safety in humans. An extensive literature search was performed in the PubMed and EMBASE databases using terms related to the drugs of interest, congenital anomalies and fetal development abnormalities, and genetic variation and susceptibility. A total of 29 studies were included in the final data extraction. The eligible studies were published between 1999 and 2020 in 10 different countries, and comprised 28 candidate gene and 1 whole-exome sequencing studies. The sample sizes ranged from 20 to 9,774 individuals. Several drugs were investigated, including antidepressants (nine studies), thalidomide (seven studies), antiepileptic drugs (five studies), glucocorticoids (four studies), acetaminophen (two studies), and sex hormones (estrogens, one study; 17-alpha hydroxyprogesterone caproate, one study). The main neonatal phenotypic outcomes included perinatal complications, cardiovascular congenital anomalies, and neurodevelopmental outcomes. The review demonstrated that studies on genetic teratology are generally small, heterogeneous, and exhibit inconsistent results. The most convincing findings were genetic variants in SLC6A4, MTHFR, and NR3C1, which were associated with drug teratogenicity by antidepressants, antiepileptics, and glucocorticoids, respectively. Notably, this review demonstrated the large knowledge gap regarding genetic susceptibility to drug teratogenicity, emphasizing the need for further efforts in the field. Future studies may be improved by increasing the sample size and applying genome-wide approaches to promote the interpretation of results. Such studies could support the clinical implementation of genetic screening to provide safer drug use in pregnant women in need of drugs.