Exome sequencing in pooled DNA samples to identify maternal pre-eclampsia risk variants

Protein Network Analysis of Whole Exome Sequencing of Severe Preeclampsia

Preeclampsia is a hypertensive disorder of pregnancy, which complicates up to 15% of US deliveries. It is an idiopathic disorder associated with several different phenotypes. We sought to determine if the genetic architecture of preeclampsia can be described by clusters of patients with variants in genes in shared protein interaction networks. We performed a case-control study using whole exome sequencing on early onset preeclamptic mothers with severe clinical features and control mothers with uncomplicated pregnancies between 2016 and 2020. A total of 143 patients were enrolled, 61 women with early onset preeclampsia with severe features based on ACOG criteria, and 82 control women at term, matched for race and ethnicity. A network analysis and visualization tool, Proteinarium, was used to confirm there are clusters of patients with shared gene networks associated with severe preeclampsia. The majority of the sequenced patients appear in two significant clusters. We identified one case dominant and one control dominant cluster. Thirteen genes were unique to the case dominated cluster. Among these genes, LAMB2, PTK2, RAC1, QSOX1, FN1, and VCAM1 have known associations with the pathogenic mechanisms of preeclampsia. Using bioinformatic analysis, we were able to identify subsets of patients with shared protein interaction networks, thus confirming our hypothesis about the genetic architecture of preeclampsia.

Variant identification in BARD1, PRDM9, RCC1, and RECQL in patients with ovarian cancer by targeted next-generation sequencing of DNA pools

Several ovarian cancer susceptibility genes have been discovered, but more are likely to exist. In this study, we aimed to analyze knowledge-based selected genes, i.e., BARD1, PRDM9, RCC1, and RECQL, in which pathogenic germline variants have been reported in patients with breast and/or ovarian cancer. As deep sequencing of DNA samples remains costly, targeted next-generation sequencing of DNA pools was utilized to screen the exons of BARD1, PRDM9, RCC1, and RECQL in ~400 Polish ovarian cancer cases. 25 pools of 16 samples (including several duplicated samples with known variants) were sequenced on the NovaSeq6000 and analyzed with SureCall (Agilent) application. The set of variants was filtrated to exclude spurious variants, and, subsequently, the identified rare genetic variants were validated using Sanger sequencing. No pathogenic mutation was found within the analyzed cohort of ovarian cancer patients. Validation genotyping of filtered rare silent and missense variants revealed that the majority of them were true alterations, especially those with a higher mutation quality value. The high concordance (R2=0.95) of population allele frequency for 44 common SNPs in the European control population (gnomAD) and our experiment confirmed the reliability of pooled sequencing. Mutations in BARD1, PRDM9, RCC1, and RECQL do not contribute substantially to the risk of ovarian cancer. Pooled DNA sequencing is a cost-effective and reliable method for the initial screening of candidate genes; however, it still requires validation of identified rare variants.

Coordinated Expressional Landscape of the Human Placental miRNome and Transcriptome

Placenta is a unique organ that serves its own function, and contributes to maternal gestational adaptation and fetal development. Coordination of its transcriptome to satisfy all the maternal-fetal needs across gestation is not fully understood. MicroRNAs are powerful transcriptome modulators capable to adjust rapidly the expression level and dynamics of large gene sets. This MiR-Seq based study presents a multi-omics investigation of the human placental miRNome and its synergy with the transcriptome. The analysis included 52 placentas representing three trimesters of normal pregnancy, and term cases of late-onset preeclampsia (LO-PE), gestational diabetes and affected fetal growth. Gestational-age dependent differential expression (FDR < 0.05) was detected for 319 of 417 tested miRNAs (76.5%). A shared list of target genes of dynamic miRNAs suggested their coordinated action. The most abundant miR-143-3p revealed as a marker for pregnancy progression. The data suggested critical, but distinct roles of placenta-specific imprinted C19MC and C14MC miRNA clusters. Paternally encoded primate-specific C19MC was highly transcribed during first trimester, potentially fine-tuning the early placental transcriptome in dosage-sensitive manner. Maternally encoded eutherian C14MC showed high expression until term, underlining its key contribution across gestation. A major shift in placental miRNome (16% miRNAs) was observed in LO-PE, but not in other term pregnancy complications. Notably, 13/38 upregulated miRNAs were transcribed from C19MC and only one from C14MC, whereas 11/28 downregulated miRNAs represented C14MC and none C19MC. miR-210-3p, miR-512-5p, miR-32-5p, miR-19a-3p, miR-590-3p, miR-379-5p were differentially expressed in LO-PE and cases of small-for-gestational-age newborns, supporting a shared etiology. Expression correlation analysis with the RNA-Seq data (16,567 genes) of the same samples clustered PE-linked miRNAs into five groups. Large notable clusters of miRNA–gene pairs showing directly and inversely correlated expression dynamics suggested potential functional relationships in both scenarios. The first genome-wide study of placental miR-eQTLs identified 66 placental SNVs associated with the expression of neighboring miRNAs, including PE-linked miRNAs miR-30a-5p, miR-210-3p, miR-490-3p and miR-518-5p. This study provided a rich catalog of miRNAs for further in-depth investigations of their individual and joint effect on placental transcriptome. Several highlighted miRNAs may serve as potential biomarkers for pregnancy monitoring and targets to prevent or treat gestational disorders.

Identifying new potential genetic biomarkers for HELLP syndrome using massive parallel sequencing

BACKGROUND

Preeclampsia (PE) is a frequently occurring multisystemic disease affecting ~5% of pregnancies. PE patients may develop HELLP syndrome (haemolysis, elevated liver enzymes, and low platelet), a mother and foetus life-threatening condition. Research into HELLP's genetic origin has been relatively unsuccessful, mainly because normal placental function and blood pressure regulation involve the fine-regulation of hundreds of genes.

OBJECTIVE

To identify new genes and mutations constituting potential biomarkers for HELLP syndrome.

STUDY DESIGN

The present case-control study involved whole-exome sequencing of 79 unrelated HELLP women. Candidate variants were screened in a control population constituted by 176 individuals. Stringent bioinformatics filters were used for selecting potentially etiological sequence variants in a subset of 487 genes. We used robust in silico mutation modelling for predicting the potential effect on protein structure.

RESULTS

We identified numerous sequence variants in genes related to angiogenesis/coagulation/blood pressure regulation, cell differentiation/communication/adhesion, cell cycle and transcriptional gene regulation, extracellular matrix biology, lipid metabolism and immunological response. Five sequence variants generated premature stop codons in genes playing an essential role in placental physiology (STOX1, PDGFD, IGF2, MMP1 and DNAH11). Six variants (ERAP1- p.Ile915Thr, ERAP2- p.Leu837Ser, COMT-p.His192Gln, CSAD-p.Pro418Ser, CDH1- p.Ala298Thr and CCR2-p.Met249Lys) led to destabilisation of protein structure as they had significant energy and residue interaction-related changes. We identified at least two mutations in 57% of patients, arguing in favour of a polygenic origin for the HELLP syndrome.

CONCLUSION

Our results provide novel evidence regarding PE/HELLP's genetic origin, leading to new biomarkers, having potential clinical usefulness, being proposed.

Whole-exome sequencing in multiplex preeclampsia families identifies novel candidate susceptibility genes

OBJECTIVE

Preeclampsia is a common and serious heritable disorder of human pregnancy. Although there have been notable successes in identification of maternal susceptibility genes a large proportion of the heritability of preeclampsia remains unaccounted for. It is has been postulated that rare variation may account for some of this missing heritability. In this study, we performed whole-exome sequencing (WES) in multiplex families to identify rare exonic risk variants.

METHODS

We conducted WES in 244 individuals from 34 Australian/New Zealand multiplex preeclampsia families. Variants were tested for association with preeclampsia using a threshold model and logistic regression.

RESULTS

We found significant association for two moderately rare missense variants, rs145743393 (Padj = 0.0032, minor allele frequency = 0.016) in the chromosome 1 open reading frame 35 (C1orf35) gene, and rs34270076 (Padj = 0.0128, minor allele frequency = 0.024) in the pyroglutamylated RFamide peptide receptor (QRFPR) gene. To replicate these associations we performed imputation in our Australian genome wide association scan for preeclampsia and found no significant exonic variants in either C1orf35 or QRFPR. However, 11 variants demonstrating nominal significance (P < 0.05) in the genomic region between QRFPR and annexin A5 (ANXA5) were identified. We further leveraged publicly available genome-wide available summary data from the UK Biobank to investigate association of these two variants with the underlying clinical phenotypes of preeclampsia and detected nominal association of the QRFPR variant (rs34270076, P = 0.03) with protein levels in females.

CONCLUSION

The study represents the first to use WES in multiplex families for preeclampsia and identifies two novel genes (QRFPR and C1orf35) not previously associated with preeclampsia and find nominal association of rs34270076 with protein levels, a key clinical feature of preeclampsia. We find further support for ANXA5 previously associated with pregnancy complications, including preeclampsia.

Exploring the Molecular Aetiology of Preeclampsia by Massive Parallel Sequencing of DNA

PURPOSE OF REVIEW

This manuscript aims to review (for the first time) studies describing NGS sequencing of preeclampsia (PE) women's DNA.

RECENT FINDINGS

Describing markers for the early detection of PE is an essential task because, although associated molecular dysfunction begins early on during pregnancy, the disease's clinical signs usually appear late in pregnancy. Although several biochemical biomarkers have been proposed, their use in clinical environments is still limited, thereby encouraging research into PE's genetic origin. Hundreds of genes involved in numerous implantation- and placentation-related biological processes may be coherent candidates for PE aetiology. Next-generation sequencing (NGS) offers new technical possibilities for PE studying, as it enables large genomic regions to be analysed at affordable cost. This technique has facilitated the description of genes contributing to the molecular origin of a significant amount of monogenic and complex diseases. Regarding PE, NGS of DNA has been used in familial and isolated cases, thereby enabling new genes potentially related to the phenotype to be proposed. For a better understanding of NGS, technical aspects, applications and limitations are presented initially. Thereafter, NGS studies of DNA in familial and non-familial cases are described, including pitfalls and positive findings. The information given here should enable scientists and clinicians to analyse and design new studies permitting the identification of novel clinically useful molecular PE markers.

Identification of genetic polymorphisms modulating nausea and vomiting in two series of opioid-treated cancer patients

Nausea and vomiting are often associated with opioid analgesia in cancer patients; however, only a subset of patients develop such side effects. Here, we tested the hypothesis that the occurrence of nausea and vomiting is modulated by the genetic background of the patients. Whole exome sequencing of DNA pools from patients with either low (n = 937) or high (n = 557) nausea and vomiting intensity, recruited in the European Pharmacogenetic Opioid Study, revealed a preliminary association of 53 polymorphisms. PCR-based genotyping of 45 of these polymorphisms in the individual patients of the same series confirmed the association for six SNPs in AIM1L, CLCC1, MUC16, PDE3A, POM121L2, and ZNF165 genes. Genotyping of the same 45 polymorphisms in 264 patients of the Italian CERP study, also treated with opioids for cancer pain, instead confirmed the association for two SNPs in ZNF568 and PDE3A genes. Only one SNP, rs12305038 in PDE3A, was confirmed in both series, although with opposite effects of the minor allele on the investigated phenotype. Overall, our findings suggest that genetic factors are indeed associated with nausea and vomiting in opioid-treated cancer patients, but the role of individual polymorphisms may be weak.

The genetic component of preeclampsia: A whole-exome sequencing study

Preeclampsia is a major cause of maternal and perinatal deaths. The aetiology of preeclampsia is largely unknown but a polygenetic component is assumed. To explore this hypothesis, we performed an in-depth whole-exome sequencing study in women with (cases, N = 50) and without (controls, N = 50) preeclampsia. The women were identified in an unselected cohort of 2,545 pregnant women based on data from the Danish National Patient Registry and the Medical Birth Registry. Matching DNA was obtained from a biobank containing excess blood from routine antenatal care visits. Novogene performed the whole-exome sequencing blinded to preeclampsia status. Variants for comparison between cases and controls were filtered in the Ingenuity Variant Analysis software. We applied two different strategies; a disease association panel approach, which included variants in single genes associated with established clinical risk factors for preeclampsia, and a gene panel approach, which included biological pathways harbouring genes previously reported to be associated with preeclampsia. Variant variability was compared in cases and controls at the level of biological processes, signalling pathways, and in single genes. Regardless of the applied strategy and the level of variability examined, we consistently found positive correlations between variant numbers in cases and controls (all R²s>0.88). Contrary to what was expected, cases carried fewer variants in biological processes and signalling pathways than controls (all p-values ≤0.02). In conclusion, our findings challenge the hypothesis of a polygenetic aetiology for preeclampsia with a common network of susceptibility genes. The greater genetic diversity among controls may suggest a protective role of genetic diversity against the development of preeclampsia.

Genetic Approaches in Preeclampsia

Preeclampsia (PE) is a serious hypertensive disorder that affects up to 8% of all pregnancies annually. An established risk factor for PE is family history, clearly demonstrating an underlying genetic component to the disorder. To date, numerous genetic studies, using both the candidate gene and genome-wide approach, have been undertaken to tease out the genetic basis of PE and understand its origins. Such studies have identified some promising candidate genes such as STOX1 and ACVR2A. Nevertheless, researchers face ongoing challenges of replicating these genetic associations in different populations and performing the functional validation of identified genetic variants to determine their causality in the disorder. This chapter will review the genetic approaches used in the study of PE, discuss their limitations and possible confounders, and describe current strategies.

Drug-sensitive FGFR3 mutations in lung adenocarcinoma

BACKGROUND

Lung cancer is the leading cause of cancer-related deaths across the world. In this study, we present therapeutically relevant genetic alterations in lung adenocarcinoma of Indian origin.

MATERIALS AND METHODS

Forty-five primary lung adenocarcinoma tumors were sequenced for 676 amplicons using RainDance cancer panel at an average coverage of 1500 × (reads per million mapped reads). To validate the findings, 49 mutations across 23 genes were genotyped in an additional set of 363 primary lung adenocarcinoma tumors using mass spectrometry. NIH/3T3 cells over expressing mutant and wild-type FGFR3 constructs were characterized for anchorage independent growth, constitutive activation, tumor formation and sensitivity to FGFR inhibitors using in vitro and xenograft mouse models.

RESULTS

We present the first spectrum of actionable alterations in lung adenocarcinoma tumors of Indian origin, and shows that mutations of FGFR3 are present in 20 of 363 (5.5%) patients. These FGFR3 mutations are constitutively active and oncogenic when ectopically expressed in NIH/3T3 cells and using a xenograft model in NOD/SCID mice. Inhibition of FGFR3 kinase activity inhibits transformation of NIH/3T3 overexpressing FGFR3 constructs and growth of tumors driven by FGFR3 in the xenograft models. The reduction in tumor size in the mouse is paralleled by a reduction in the amounts of phospho-ERK, validating the in vitro findings. Interestingly, the FGFR3 mutations are significantly higher in a proportion of younger patients and show a trend toward better overall survival, compared with patients lacking actionable alterations or those harboring KRAS mutations.

CONCLUSION

We present the first actionable mutation spectrum in Indian lung cancer genome. These findings implicate FGFR3 as a novel therapeutic in lung adenocarcinoma.

Investigation of rare and low-frequency variants using high-throughput sequencing with pooled DNA samples

High-throughput sequencing using pooled DNA samples can facilitate genome-wide studies on rare and low-frequency variants in a large population. Some major questions concerning the pooling sequencing strategy are whether rare and low-frequency variants can be detected reliably, and whether estimated minor allele frequencies (MAFs) can represent the actual values obtained from individually genotyped samples. In this study, we evaluated MAF estimates using three variant detection tools with two sets of pooled whole exome sequencing (WES) and one set of pooled whole genome sequencing (WGS) data. Both GATK and Freebayes displayed high sensitivity, specificity and accuracy when detecting rare or low-frequency variants. For the WGS study, 56% of the low-frequency variants in Illumina array have identical MAFs and 26% have one allele difference between sequencing and individual genotyping data. The MAF estimates from WGS correlated well (r = 0.94) with those from Illumina arrays. The MAFs from the pooled WES data also showed high concordance (r = 0.88) with those from the individual genotyping data. In conclusion, the MAFs estimated from pooled DNA sequencing data reflect the MAFs in individually genotyped samples well. The pooling strategy can thus be a rapid and cost-effective approach for the initial screening in large-scale association studies.