Toward clinical exomes in diagnostics and management of male infertility.

Infertility, affecting ∼10% of men, is predominantly caused by primary spermatogenic failure (SPGF). We screened likely pathogenic and pathogenic (LP/P) variants in 638 candidate genes for male infertility in 521 individuals presenting idiopathic SPGF and 323 normozoospermic men in the ESTAND cohort. Molecular diagnosis was reached for 64 men with SPGF (12%), with findings in 39 genes (6%). The yield did not differ significantly between the subgroups with azoospermia (20/185, 11%), oligozoospermia (18/181, 10%), and primary cryptorchidism with SPGF (26/155, 17%). Notably, 19 of 64 LP/P variants (30%) identified in 28 subjects represented recurrent findings in this study and/or with other male infertility cohorts. NR5A1 was the most frequently affected gene, with seven LP/P variants in six SPGF-affected men and two normozoospermic men. The link to SPGF was validated for recently proposed candidate genes ACTRT1, ASZ1, GLUD2, GREB1L, LEO1, RBM5, ROS1, and TGIF2LY. Heterozygous truncating variants in BNC1, reported in female infertility, emerged as plausible causes of severe oligozoospermia. Data suggested that several infertile men may present congenital conditions with less pronounced or pleiotropic phenotypes affecting the development and function of the reproductive system. Genes regulating the hypothalamic-pituitary-gonadal axis were affected in >30% of subjects with LP/P variants. Six individuals had more than one LP/P variant, including five with two findings from the gene panel. A 4-fold increased prevalence of cancer was observed in men with genetic infertility compared to the general male population (8% vs. 2%; p = 4.4 × 10-3). Expanding genetic testing in andrology will contribute to the multidisciplinary management of SPGF.

DDB1- and CUL4-associated factor 12-like protein 1 (Dcaf12l1) is not essential for male fertility in mice.

Male infertility is a common condition affecting at least 7% of men worldwide and is often genetic in origin. Using whole exome sequencing, we recently discovered three hemizygous, likely damaging variants in DDB1- and CUL4-associated factor 12-like protein 1 (DCAF12L1) in men with azoospermia. DCAF12L1 is located on the X-chromosome and as identified by single cell sequencing studies, its expression is enriched in human testes and specifically in Sertoli cells and spermatogonia. However, very little is known about the role of DCAF12L1 in spermatogenesis, thus we generated a knockout mouse model to further explore the role of DCAF12L1 in male fertility. Knockout mice were generated using CRISPR/Cas9 technology to remove the entire coding region of Dcaf12l1 and were assessed for fertility over a broad range of ages (2-8 months of age). Despite outstanding genetic evidence in men, loss of DCAF12L1 had no discernible impact on male fertility in mice, as highlighted by breeding trials, histological assessment of the testis and epididymis, daily sperm production and evaluation of sperm motility using computer assisted methods. This disparity is likely due to the parallel evolution, and subsequent divergence, of DCAF12 family members in mice and men or the presence of compounding environmental factors in men.

Actionable secondary findings following exome sequencing of 836 non-obstructive azoospermia cases and their value in patient management.

STUDY QUESTION: What is the load, distribution and added clinical value of secondary findings (SFs) identified in exome sequencing (ES) of patients with non-obstructive azoospermia (NOA)? SUMMARY ANSWER: One in 28 NOA cases carried an identifiable, medically actionable SF. WHAT IS KNOWN ALREADY: In addition to molecular diagnostics, ES allows assessment of clinically actionable disease-related gene variants that are not connected to the patient's primary diagnosis, but the knowledge of which may allow the prevention, delay or amelioration of late-onset monogenic conditions. Data on SFs in specific clinical patient groups, including reproductive failure, are currently limited. STUDY DESIGN, SIZE, DURATION: The study group was a retrospective cohort of patients with NOA recruited in 10 clinics across six countries and formed in the framework of the international GEMINI (The GEnetics of Male INfertility Initiative) study. PARTICIPANTS/MATERIALS, SETTING, METHODS: ES data of 836 patients with NOA were exploited to analyze SFs in 85 genes recommended by the American College of Medical Genetics and Genomics (ACMG), Geisinger's MyCode, and Clinical Genome Resource. The identified 6374 exonic variants were annotated with ANNOVAR and filtered for allele frequency, retaining 1381 rare or novel missense and loss-of-function variants. After automatic assessment of pathogenicity with ClinVar and InterVar, 87 variants were manually curated. The final list of confident disease-causing SFs was communicated to the corresponding GEMINI centers. When patient consent had been given, available family health history and non-andrological medical data were retrospectively assessed. MAIN RESULTS AND THE ROLE OF CHANCE: We found a 3.6% total frequency of SFs, 3.3% from the 59 ACMG SF v2.0 genes. One in 70 patients carried SFs in genes linked to familial cancer syndromes, whereas 1 in 60 cases was predisposed to congenital heart disease or other cardiovascular conditions. Retrospective assessment confirmed clinico-molecular diagnoses in several cases. Notably, 37% (11/30) of patients with SFs carried variants in genes linked to male infertility in mice, suggesting that some SFs may have a co-contributing role in spermatogenic impairment. Further studies are needed to determine whether these observations represent chance findings or the profile of SFs in NOA patients is indeed different from the general population. LIMITATIONS, REASONS FOR CAUTION: One limitation of our cohort was the low proportion of non-Caucasian ethnicities (9%). Additionally, as comprehensive clinical data were not available retrospectively for all men with SFs, we were not able to confirm a clinico-molecular diagnosis and assess the penetrance of the specific variants. WIDER IMPLICATIONS OF THE FINDINGS: For the first time, this study analyzed medically actionable SFs in men with spermatogenic failure. With the evolving process to incorporate ES into routine andrology practice for molecular diagnostic purposes, additional assessment of SFs can inform about future significant health concerns for infertility patients. Timely detection of SFs and respective genetic counseling will broaden options for disease prevention and early treatment, as well as inform choices and opportunities regarding family planning. A notable fraction of SFs was detected in genes implicated in maintaining genome integrity, essential in both mitosis and meiosis. Thus, potential genetic pleiotropy may exist between certain adult-onset monogenic diseases and NOA. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by the Estonian Research Council grants IUT34-12 and PRG1021 (M.L. and M.P.); National Institutes of Health of the United States of America grant R01HD078641 (D.F.C., K.I.A. and P.N.S.); National Institutes of Health of the United States of America grant P50HD096723 (D.F.C. and P.N.S.); National Health and Medical Research Council of Australia grant APP1120356 (M.K.O'B., D.F.C. and K.I.A.); Fundação para a Ciência e a Tecnologia (FCT)/Ministério da Ciência, Tecnologia e Inovação grant POCI-01-0145-FEDER-007274 (A.M.L., F.C. and J.G.) and FCT: IF/01262/2014 (A.M.L.). J.G. was partially funded by FCT/Ministério da Ciência, Tecnologia e Ensino Superior (MCTES), through the Centre for Toxicogenomics and Human Health-ToxOmics (grants UID/BIM/00009/2016 and UIDB/00009/2020). M.L.E. is a consultant for, and holds stock in, Roman, Sandstone, Dadi, Hannah, Underdog and has received funding from NIH/NICHD. Co-authors L.K., K.L., L.N., K.I.A., P.N.S., J.G., F.C., D.M.-M., K.A., K.A.J., M.K.O'B., A.M.L., D.F.C., M.P. and M.L. declare no conflict of interest. TRIAL REGISTRATION NUMBER: N/A.

Value of soluble fms-like tyrosine kinase-1/placental growth factor test in third trimester of pregnancy for predicting preeclampsia in asymptomatic women.

OBJECTIVES: To estimate the value of screening maternal serum soluble fms-like tyrosine kinase/placental growth factor (sFlt-1/PlGF) ratio in asymptomatic women during 3rd trimester to predict preeclampsia (PE) development. METHODS: The investigated group comprised of 178 pregnant women. During this gestation, 24 cases had developed PE and 12 isolated gestational hypertension (GH); whereas 142 remained normotensive. Blood samples were collected between 180 and 259 gestational days (g.d.) when the participants were asymptomatic. Serums were analyzed using the BRAHMS sFlt-1 Kryptor/BRAHMS PlGF plus Kryptor PE ratio test (Thermo Fisher Scientific, Henningdorf, Germany). High-risk pregnancies for the PE development were defined as sFlt-1/PlGF>38. RESULTS: The detection rate (DR) for manifestation of PE≤30 days after sampling was 83.3% and overall DR during pregnancy 58.3%. Ten of 15 women having false positive prediction of PE suffered from GH, preterm birth and/or delivery of a small-for-gestational-age-newborn. False positive rate was significantly higher at 239-253 g.d. compared to sampling at 210-224 g.d. and 225-238 g.d. (21.9% vs. 7.8% and 5.3%; p < 0.05). CONCLUSIONS: The sFlt-1/PlGF test during 180-259 g.d. detected approximately half of subsequent PE cases. An optimal time to use the test for screening purposes was estimated 225-238 g.d. (DR 66.7%). False positive test results were more common to cases with other adverse pregnancy outcomes and samples drawn at higher gestational age.

Evolutionary and functional analysis of RBMY1 gene copy number variation on the human Y chromosome.

Human RBMY1 genes are located in four variable-sized clusters on the Y chromosome, expressed in male germ cells and possibly associated with sperm motility. We have re-investigated the mutational background and evolutionary history of the RBMY1 copy number distribution in worldwide samples and its relevance to sperm parameters in an Estonian cohort of idiopathic male factor infertility subjects. We estimated approximate RBMY1 copy numbers in 1218 1000 Genomes Project phase 3 males from sequencing read-depth, then chose 14 for valid ation by multicolour fibre-FISH. These fibre-FISH samples provided accurate calibration standards for the entire panel and led to detailed insights into population variation and mutational mechanisms. RBMY1 copy number worldwide ranged from 3 to 13 with a mode of 8. The two larger proximal clusters were the most variable, and additional duplications, deletions and inversions were detected. Placing the copy number estimates onto the published Y-SNP-based phylogeny of the same samples suggested a minimum of 562 mutational changes, translating to a mutation rate of 2.20 × 10-3 (95% CI 1.94 × 10-3 to 2.48 × 10-3) per father-to-son Y-transmission, higher than many short tandem repeat (Y-STRs), and showed no evidence for selection for increased or decreased copy number, but possible copy number stabilizing selection. An analysis of RBMY1 copy numbers among 376 infertility subjects failed to replicate a previously reported association with sperm motility and showed no significant effect on sperm count and concentration, serum follicle stimulating hormone (FSH), luteinizing hormone (LH) and testosterone levels or testicular and semen volume. These results provide the first in-depth insights into the structural rearrangements underlying RBMY1 copy number variation across diverse human lineages.

Monogenic causes of non-obstructive azoospermia: challenges, established knowledge, limitations and perspectives.

It is estimated that one in 100 men have azoospermia, the complete lack of sperm in the ejaculate. Currently, ~ 20% of azoospermia cases remain idiopathic. Non-obstructive azoospermia (NOA) is mostly explained by congenital factors leading to spermatogenic failure, such as chromosome abnormalities. The knowledge of the monogenic causes of NOA is very limited. High genetic heterogeneity due to the complexity of spermatogenesis and testicular function, lack of non-consanguineous familial cases and confirmatory studies challenge the field. The reported monogenic defects cause syndromic NOA phenotypes presenting also additional congenital problems and isolated NOA cases, explained by spermatogenic defects. The established and recently reported NOA genes (n = 38) represent essential guardians of meiosis, transcriptional and endocrine regulators of reproduction. Despite the list being short, 92% of these loci are predicted to functionally interact with each other (STRING analysis: average 5.21 connections/gene, enrichment P < 10-16). Notably, ~ 50% of NOA genes have also been implicated in primary ovarian insufficiency, amenorrhea and female genital anomalies, referring to overlapping mechanisms. Considering the knowledge from respective female phenotypes and animal models, exploring the scenarios of di/oligogenic and de novo mutations represent perspective directions in the genetic research of NOA. Knowing the exact genetic cause in each patient improves the management of infertility and other health risks (e.g., cancer), and facilitates the counseling of family members about their reproductive health. Uncovering the loci and biological processes implicated in NOA will also broaden the understanding of etiologies behind spermatogenic failure and promote the development of novel non-invasive treatments for male infertility.

Variants in GCNA, X-linked germ-cell genome integrity gene, identified in men with primary spermatogenic failure.

Male infertility impacts millions of couples yet, the etiology of primary infertility remains largely unknown. A critical element of successful spermatogenesis is maintenance of genome integrity. Here, we present a genomic study of spermatogenic failure (SPGF). Our initial analysis (n = 176) did not reveal known gene-candidates but identified a potentially significant single-nucleotide variant (SNV) in X-linked germ-cell nuclear antigen (GCNA). Together with a larger follow-up study (n = 2049), 7 likely clinically relevant GCNA variants were identified. GCNA is critical for genome integrity in male meiosis and knockout models exhibit impaired spermatogenesis and infertility. Single-cell RNA-seq and immunohistochemistry confirm human GCNA expression from spermatogonia to elongated spermatids. Five identified SNVs were located in key functional regions, including N-terminal SUMO-interacting motif and C-terminal Spartan-like protease domain. Notably, variant p.Ala115ProfsTer7 results in an early frameshift, while Spartan-like domain missense variants p.Ser659Trp and p.Arg664Cys change conserved residues, likely affecting 3D structure. For variants within GCNA's intrinsically disordered region, we performed computational modeling for consensus motifs. Two SNVs were predicted to impact the structure of these consensus motifs. All identified variants have an extremely low minor allele frequency in the general population and 6 of 7 were not detected in > 5000 biological fathers. Considering evidence from animal models, germ-cell-specific expression, 3D modeling, and computational predictions for SNVs, we propose that identified GCNA variants disrupt structure and function of the respective protein domains, ultimately arresting germ-cell division. To our knowledge, this is the first study implicating GCNA, a key genome integrity factor, in human male infertility.

Bi-allelic Recessive Loss-of-Function Variants in FANCM Cause Non-obstructive Azoospermia.

Infertility affects around 7% of men worldwide. Idiopathic non-obstructive azoospermia (NOA) is defined as the absence of spermatozoa in the ejaculate due to failed spermatogenesis. There is a high probability that NOA is caused by rare genetic defects. In this study, whole-exome sequencing (WES) was applied to two Estonian brothers diagnosed with NOA and Sertoli cell-only syndrome (SCOS). Compound heterozygous loss-of-function (LoF) variants in FANCM (Fanconi anemia complementation group M) were detected as the most likely cause for their condition. A rare maternally inherited frameshift variant p.Gln498Thrfs∗7 (rs761250416) and a previously undescribed splicing variant (c.4387-10A>G) derived from the father introduce a premature STOP codon leading to a truncated protein. FANCM exhibits enhanced testicular expression. In control subjects, immunohistochemical staining localized FANCM to the Sertoli and spermatogenic cells of seminiferous tubules with increasing intensity through germ cell development. This is consistent with its role in maintaining genomic stability in meiosis and mitosis. In the individual with SCOS carrying bi-allelic FANCM LoF variants, none or only faint expression was detected in the Sertoli cells. As further evidence, we detected two additional NOA-affected case subjects with independent FANCM homozygous nonsense variants, one from Estonia (p.Gln1701∗; rs147021911) and another from Portugal (p.Arg1931∗; rs144567652). The study convincingly demonstrates that bi-allelic recessive LoF variants in FANCM cause azoospermia. FANCM pathogenic variants have also been linked with doubled risk of familial breast and ovarian cancer, providing an example mechanism for the association between infertility and cancer risk, supported by published data on Fancm mutant mouse models.

Translational aspects of novel findings in genetics of male infertility-status quo 2021.

INTRODUCTION: Male factor infertility concerns 7-10% of men and among these 40-60% remain unexplained. SOURCES OF DATA: This review is based on recent published literature regarding the genetic causes of male infertility. AREAS OF AGREEMENT: Screening for karyotype abnormalities, biallelic pathogenic variants in the CFTR gene and Y-chromosomal microdeletions have been routine in andrology practice for >20 years, explaining ~10% of infertility cases. Rare specific conditions, such as congenital hypogonadotropic hypogonadism, disorders of sex development and defects of sperm morphology and motility, are caused by pathogenic variants in recurrently affected genes, which facilitate high diagnostic yield (40-60%) of targeted gene panel-based testing. AREAS OF CONTROVERSY: Progress in mapping monogenic causes of quantitative spermatogenic failure, the major form of male infertility, has been slower. No 'recurrently' mutated key gene has been identified and worldwide, a few hundred patients in total have been assigned a possible monogenic cause. GROWING POINTS: Given the high genetic heterogeneity, an optimal approach to screen for heterogenous genetic causes of spermatogenic failure is sequencing exomes or in perspective, genomes. Clinical guidelines developed by multidisciplinary experts are needed for smooth integration of expanded molecular diagnostics in the routine management of infertile men. AREAS TIMELY FOR DEVELOPING RESEARCH: Di-/oligogenic causes, structural and common variants implicated in multifactorial inheritance may explain the 'hidden' genetic factors. It is also critical to understand how the recently identified diverse genetic factors of infertility link to general male health concerns across lifespan and how the clinical assessment could benefit from this knowledge.

NR5A1 c.991-1G > C splice-site variant causes familial 46,XY partial gonadal dysgenesis with incomplete penetrance.

OBJECTIVE: The study aimed to identify the genetic basis of partial gonadal dysgenesis (PGD) in a non-consanguineous family from Estonia. PATIENTS: Cousins P (proband) 1 (12 years; 46,XY) and P2 (18 years; 46,XY) presented bilateral cryptorchidism, severe penoscrotal hypospadias, low bitesticular volume and azoospermia in P2. Their distant relative, P3 (30 years; 46,XY), presented bilateral cryptorchidism and cryptozoospermia. DESIGN: Exome sequencing was targeted to P1-P3 and five unaffected family members. RESULTS: P1-P2 were identified as heterozygous carriers of NR5A1 c.991-1G > C. NR5A1 encodes the steroidogenic factor-1 essential in gonadal development and specifically expressed in adrenal, spleen, pituitary and testes. Together with a previous PGD case from Belgium (Robevska et al 2018), c.991-1G > C represents the first recurrent NR5A1 splice-site mutation identified in patients. The majority of previous reports on NR5A1 mutation carriers have not included phenotype-genotype data of the family members. Segregation analysis across three generations showed incomplete penetrance (<50%) and phenotypic variability among the carriers of NR5A1 c.991-1G > C. The variant pathogenicity was possibly modulated by rare heterozygous variants inherited from the other parent, OTX2 p.P134R (P1) or PROP1 c.301_302delAG (P2). For P3, the pedigree structure supported a distinct genetic cause. He carries a previously undescribed likely pathogenic variant SOS1 p.Y136H. SOS1, critical in Ras/MAPK signalling and foetal development, is a strong novel candidate gene for cryptorchidism. CONCLUSIONS: Detailed genetic profiling facilitates counselling and clinical management of the probands, and supports unaffected mutation carriers in the family for their reproductive decision making.

The decline of FANCM immunohistochemical expression in prostate cancer stroma correlates with the grade group.

Prostate adenocarcinoma (PCa) stromal markers have recently gained attention as complementary diagnostic tools. The DNA reparation complex protein FANCM has been shown to express in the normal prostate stroma and FANCM gene alterations to be associated with PCa susceptibility; this has led to the hypothesis that an insufficient level of FANCM expression may provide additional information for the evaluation of PCa. The study cohort comprised 60 radical prostatectomy specimens. The controls involved 11 autopsies (CTRL) and non-cancerous tissue (NCT) areas from the prostatectomy specimen. The samples were stained with the FANCM antibody. The quantification of the stromal staining index (SSI) was made using ImageJ and QuPath. Overall, 655 regions of interest (ROI) were analyzed. FANCM expression appeared equally intense and stroma specific in both CTRL and NCT, indicating the absence of underlying baseline alterations. Within the age span of the cohort 47-89 years, no significant effect of the age of the patients on the FANCM expression was seen. FANCM demonstrated Gleason grade (G) dependent decline in PCa, being statistically significant in controls versus G1 and G2 versus G3. In other adjacent International Society of Urological Pathology (ISUP) groups, it remained insignificant, still being meaningful between high and low-grade cancers.

FSHB -211 G>T is a major genetic modulator of reproductive physiology and health in childbearing age women.

STUDY QUESTION: Are the genetic variants FSHB -211 G>T (rs10835638), FSHR c.2039 A>G (Asn680Ser, rs6166) and FSHR -29 G>A (rs1394205) associated with serum FSH, LH and anti-Müllerian hormone (AMH) levels in reproductive age women, their menstrual cycle parameters and risk of infertility? SUMMARY ANSWER: Only the FSHB -211 G>T variant was a major genetic determinant of serum gonadotropin levels in both, eumenorrheic healthy women and female infertility patients, and the T-allele carrier status was enriched among idiopathic infertility cases. WHAT IS KNOWN ALREADY: There are accumulating data on common genetic variants modulating reproductive parameters and fertility potential. FSHB -211 G>T represents the strongest acknowledged genetic factor contributing to male circulating gonadotropins levels. Respective data in women are limited and the two previously published studies have reached conflicting results. In addition, previous studies have consistently associated FSHR c.2039 A>G (but not FSHR -29 G>A) with female serum FSH level. STUDY DESIGN, SIZE, DURATION: The study aimed to test robust and clinically meaningful genetic effects (if present) of the FSHB -211 G>T, FSHR c.2039 A>G and FSHR -29 G>A variants on female basal FSH, LH and AMH levels, and linked reproductive parameters. Genetic association testing was performed in two independent and clinically different study groups (i) eumenorrheic healthy women without known fertility problems (n = 169; 27.6 ± 6.1 years) and (ii) female partners of infertile couples (n = 186; 32.4 ± 4.7 years). The study groups were compared for allelic and genotypic distributions of the analysed variants. PARTICIPANTS/MATERIALS, SETTING, METHODS: All participants were recruited during the HAPPY PREGNANCY study (2013-2015) at the Women's Clinic, Tartu University Hospital, Estonia. Serum FSH, LH and AMH were measured in the follicular phase (Days 2-6) of the menstrual cycle. All three single nucleotide polymorphisms (SNPs) were genotyped by PCR and Taqman allelic discrimination assay. The effect of the analysed variants on hormonal measurements and menstrual cycle data was assessed using linear regression under additive and recessive models adjusted by age, BMI and smoking status. Results of the two subgroups were combined in a meta-analysis applying the fixed effects model. Restricted maximum likelihood analysis was applied to estimate the proportion of total phenotypic variance of analysed reproductive parameters, explainable by the tested genetic variants. In case-control analysis, genetic association with infertility status was tested using Fisher's exact test and logistic regression adjusted by age, BMI and smoking status. MAIN RESULTS AND THE ROLE OF CHANCE: In both study groups, T-allele of the FSHB -211 G>T was associated with significantly higher serum levels of FSH and LH. Results of the meta-analysis (additive genetic model) remained significant after Bonferroni correction for multiple testing: FSH, T-allele effect 0.80 IU/L, P = 1.2 × 10-3; LH, 1.58 IU/L, P = 1.8×10-8. A more pronounced effect of T-allele of the FSHB -211 G>T on circulating LH was identified as a driving factor to increased LH/FSH ratio (meta-analysis, P = 4.7 × 10-3). In healthy women, the FSHB -211 G>T variant was estimated to explain 3.5 and 7.1% of the total variance of the measured serum FSH and LH levels, respectively. The corresponding numbers for the infertility patients were 1.6 and 10.5%. Women with idiopathic infertility compared to controls exhibited a doubled T-allele frequency (23.6 versus 12.4%; P = 8.9 × 10-3) and a >3-fold excess of TT homozygotes (5.6 versus 1.8%; P = 3.5 × 10-2). The only association of the FSHR c.2039 A>G was detected with serum FSH levels in eumenorrheic healthy women, explaining 3.9% of the total parameter variance (G-allele effect 0.56 IU/L, P = 4.6 × 10-3). In the study group of healthy reproductive age women, the highest serum FSH levels were detected among the FSHB -211 T-allele carriers with the FSHR c.2039 GG-genotype (median 7.7 IU/L). In contrast, the lowest hormone concentrations were measured for the women carrying the combination of the FSHB -211 GG- and the FSHR c.2039 AA-homozygosity (median 5.8 IU/L, P = 9.6 × 10-3). None of the analysed reproductive parameters was associated with the FSHR -29 G>A variant. In our study groups, the tested polymorphisms did not reach significant associations with serum AMH measurements, menstrual cycle length or age at menarche. LIMITATIONS, REASONS FOR CAUTION: Small sample size and the design involving two clinical groups with different reproductive histories may have limited the capacity to replicate the associations with the age at menarche and length of menstrual cycle, initially reported in large genome-wide association studies. Small sample size may have also affected the accuracy in estimating the contribution of the tested variants to the total phenotypic variance of measured gonadotropin concentrations. The group of eumenorrheic healthy women had its limitations as a control to estimate the true effect of analysed genetic variants on individual's fertility potential as the recruitment strategy had been targeted mostly towards younger women, who may not yet have planned to conceive a child by this age. WIDER IMPLICATIONS OF THE FINDINGS: We propose that like in men, also in women the FSHB -211 G>T represents a key genetic modulator of circulating gonadotropin, leading to various possible downstream effects on reproductive physiology. This claim is strongly supported by the reports of genome-wide association studies on various female reproductive traits and diseases. In perspective, FSHB -211 G>T may have a diagnostic value in fertility clinics to detect female patients with genetically inherited elevated basal FSH and LH levels. STUDY FUNDING/COMPETING INTEREST(S): The study was supported by Estonian Science Foundation Grant (ETF9030 for M.L.); Institutional Research Grant (IUT34-12 for M.L.) and European Union through the European Regional Development Fund (project HAPPY PREGNANCY, 3.2.0701.12-0047; for M.L. and K.R.). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the article. We have no competing interests to declare. TRAIL REGISTRATION NUMBER: Not applicable.

Effects of RNA integrity on transcript quantification by total RNA sequencing of clinically collected human placental samples.

RNA degradation is a ubiquitous process that occurs in living and dead cells, as well as during handling and storage of extracted RNA. Reduced RNA quality caused by degradation is an established source of uncertainty for all RNA-based gene expression quantification techniques. RNA sequencing is an increasingly preferred method for transcriptome analyses, and dependence of its results on input RNA integrity is of significant practical importance. This study aimed to characterize the effects of varying input RNA integrity [estimated as RNA integrity number (RIN)] on transcript level estimates and delineate the characteristic differences between transcripts that differ in degradation rate. The study used ribodepleted total RNA sequencing data from a real-life clinically collected set (n = 32) of human solid tissue (placenta) samples. RIN-dependent alterations in gene expression profiles were quantified by using DESeq2 software. Our results indicate that small differences in RNA integrity affect gene expression quantification by introducing a moderate and pervasive bias in expression level estimates that significantly affected 8.1% of studied genes. The rapidly degrading transcript pool was enriched in pseudogenes, short noncoding RNAs, and transcripts with extended 3' untranslated regions. Typical slowly degrading transcripts (median length, 2389 nt) represented protein coding genes with 4-10 exons and high guanine-cytosine content.-Reiman, M., Laan, M., Rull, K., Sõber, S. Effects of RNA integrity on transcript quantification by total RNA sequencing of clinically collected human placental samples.

CDH13 promoter SNPs with pleiotropic effect on cardiometabolic parameters represent methylation QTLs.

CDH13 encodes T-cadherin, a receptor for high molecular weight (HMW) adiponectin and low-density lipoprotein, promoting proliferation and migration of endothelial cells. Genome-wide association studies have mapped multiple variants in CDH13 associated with cardiometabolic traits (CMT) with variable effects across studies. We hypothesized that this heterogeneity might reflect interplay with DNA methylation within the region. Resequencing and EpiTYPER™ assay were applied for the HYPertension in ESTonia/Coronary Artery Disease in Czech (HYPEST/CADCZ; n = 358) samples to identify CDH13 promoter SNPs acting as methylation Quantitative Trait Loci (meQTLs) and to investigate their associations with CMT. In silico data were extracted from genome-wide DNA methylation and genotype datasets of the population-based sample Estonian Genome Center of the University of Tartu (EGCUT; n = 165). HYPEST-CADCZ meta-analysis identified a rare variant rs113460564 as highly significant meQTL for a 134-bp distant CpG site (P = 5.90 × 10(-6); ß = 3.19%). Four common SNPs (rs12443878, rs12444338, rs62040565, rs8060301) exhibited effect on methylation level of up to 3 neighboring CpG sites in both datasets. The strongest association was detected in EGCUT between rs8060301 and cg09415485 (false discovery rate corrected P value = 1.89 × 10(-30)). Simultaneously, rs8060301 showed association with diastolic blood pressure, serum high-density lipoprotein and HMW adiponectin (P < 0.005). Novel strong associations were identified between rare CDH13 promoter meQTLs (minor allele frequency <5%) and HMW adiponectin: rs2239857 (P = 5.50 × 10(-5), ß = -1,841.9 ng/mL) and rs77068073 (P = 2.67 × 10(-4), ß = -2,484.4 ng/mL). Our study shows conclusively that CDH13 promoter harbors meQTLs associated with CMTs. It paves the way to deeper understanding of the interplay between DNA variation and methylation in susceptibility to common diseases.

Structural genomic variation as risk factor for idiopathic recurrent miscarriage.

Recurrent miscarriage (RM) is a multifactorial disorder with acknowledged genetic heritability that affects ∼3% of couples aiming at childbirth. As copy number variants (CNVs) have been shown to contribute to reproductive disease susceptibility, we aimed to describe genome-wide profile of CNVs and identify common rearrangements modulating risk to RM. Genome-wide screening of Estonian RM patients and fertile controls identified excessive cumulative burden of CNVs (5.4 and 6.1 Mb per genome) in two RM cases possibly increasing their individual disease risk. Functional profiling of all rearranged genes within RM study group revealed significant enrichment of loci related to innate immunity and immunoregulatory pathways essential for immune tolerance at fetomaternal interface. As a major finding, we report a multicopy duplication (61.6 kb) at 5p13.3 conferring increased maternal risk to RM in Estonia and Denmark (meta-analysis, n = 309/205, odds ratio = 4.82, P = 0.012). Comparison to Estonian population-based cohort (total, n = 1000) confirmed the risk for Estonian female cases (P = 7.9 × 10(-4) ). Datasets of four cohorts from the Database of Genomic Variants (total, n = 5,846 subjects) exhibited similar low duplication prevalence worldwide (0.7%-1.2%) compared to RM cases of this study (6.6%-7.5%). The CNV disrupts PDZD2 and GOLPH3 genes predominantly expressed in placenta and it may represent a novel risk factor for pregnancy complications.

Undiagnosed RASopathies in infertile men.

RASopathies are syndromes caused by congenital defects in the Ras/mitogen-activated protein kinase (MAPK) pathway genes, with a population prevalence of 1 in 1,000. Patients are typically identified in childhood based on diverse characteristic features, including cryptorchidism (CR) in >50% of affected men. As CR predisposes to spermatogenic failure (SPGF; total sperm count per ejaculate 0-39 million), we hypothesized that men seeking infertility management include cases with undiagnosed RASopathies. Likely pathogenic or pathogenic (LP/P) variants in 22 RASopathy-linked genes were screened in 521 idiopathic SPGF patients (including 155 CR cases) and 323 normozoospermic controls using exome sequencing. All 844 men were recruited to the ESTonian ANDrology (ESTAND) cohort and underwent identical andrological phenotyping. RASopathy-specific variant interpretation guidelines were used for pathogenicity assessment. LP/P variants were identified in PTPN11 (two), SOS1 (three), SOS2 (one), LZTR1 (one), SPRED1 (one), NF1 (one), and MAP2K1 (one). The findings affected six of 155 cases with CR and SPGF, three of 366 men with SPGF only, and one (of 323) normozoospermic subfertile man. The subgroup "CR and SPGF" had over 13-fold enrichment of findings compared to controls (3.9% vs. 0.3%; Fisher's exact test, p = 5.5 × 10-3). All ESTAND subjects with LP/P variants in the Ras/MAPK pathway genes presented congenital genitourinary anomalies, skeletal and joint conditions, and other RASopathy-linked health concerns. Rare forms of malignancies (schwannomatosis and pancreatic and testicular cancer) were reported on four occasions. The Genetics of Male Infertility Initiative (GEMINI) cohort (1,416 SPGF cases and 317 fertile men) was used to validate the outcome. LP/P variants in PTPN11 (three), LZTR1 (three), and MRAS (one) were identified in six SPGF cases (including 4/31 GEMINI cases with CR) and one normozoospermic man. Undiagnosed RASopathies were detected in total for 17 ESTAND and GEMINI subjects, 15 SPGF patients (10 with CR), and two fertile men. Affected RASopathy genes showed high expression in spermatogenic and testicular somatic cells. In conclusion, congenital defects in the Ras/MAPK pathway genes represent a new congenital etiology of syndromic male infertility. Undiagnosed RASopathies were especially enriched among patients with a history of cryptorchidism. Given the relationship between RASopathies and other conditions, infertile men found to have this molecular diagnosis should be evaluated for known RASopathy-linked health concerns, including specific rare malignancies.

Genome-wide association study of blood pressure extremes identifies variant near UMOD associated with hypertension.

Hypertension is a heritable and major contributor to the global burden of disease. The sum of rare and common genetic variants robustly identified so far explain only 1%-2% of the population variation in BP and hypertension. This suggests the existence of more undiscovered common variants. We conducted a genome-wide association study in 1,621 hypertensive cases and 1,699 controls and follow-up validation analyses in 19,845 cases and 16,541 controls using an extreme case-control design. We identified a locus on chromosome 16 in the 5' region of Uromodulin (UMOD; rs13333226, combined P value of 3.6 × 10⁻¹¹). The minor G allele is associated with a lower risk of hypertension (OR [95%CI]: 0.87 [0.84-0.91]), reduced urinary uromodulin excretion, better renal function; and each copy of the G allele is associated with a 7.7% reduction in risk of CVD events after adjusting for age, sex, BMI, and smoking status (H.R. = 0.923, 95% CI 0.860-0.991; p = 0.027). In a subset of 13,446 individuals with estimated glomerular filtration rate (eGFR) measurements, we show that rs13333226 is independently associated with hypertension (unadjusted for eGFR: 0.89 [0.83-0.96], p = 0.004; after eGFR adjustment: 0.89 [0.83-0.96], p = 0.003). In clinical functional studies, we also consistently show the minor G allele is associated with lower urinary uromodulin excretion. The exclusive expression of uromodulin in the thick portion of the ascending limb of Henle suggests a putative role of this variant in hypertension through an effect on sodium homeostasis. The newly discovered UMOD locus for hypertension has the potential to give new insights into the role of uromodulin in BP regulation and to identify novel drugable targets for reducing cardiovascular risk.

Urinary microRNA profiling in pregnant women identifies placental microRNAs that are candidate biomarkers for monitoring placental health.

MicroRNAs are among the key modulators of placental transcriptome dynamics. This study aimed at comparative profiling of urinary (sampling at 228-230 gestational days, g.days), serum (217-230 g.days), and placental (279-286 g.days) microRNAs in three healthy pregnant women using miRNome sequencing. Placenta showed significant enrichment of microRNAs compared with serum and urine (1174, 341, and 193, respectively; P < 10-5 ). All sample types shared 153 microRNAs, representing candidate biomarkers for placental health. Urine samples contained eight of 56 transcripts from the placenta-specific chromosome 19 microRNA cluster C19MC and one of 91 transcripts (miR-432-5p) from the chromosome 14 cluster C14MC. These data suggest an active filtering at the maternal-fetal interface in passing through only selected microRNAs. Urine represents a valid source to monitor the signature of placenta-expressed microRNAs that are differentially expressed in pregnancy complications.

Microdeletions and microduplications linked to severe congenital disorders in infertile men.

Data on the clinical validity of DNA copy number variants (CNVs) in spermatogenic failure (SPGF) is limited. This study analyzed the genome-wide CNV profile in 215 men with idiopathic SPGF and 62 normozoospermic fertile men, recruited at the Andrology Clinic, Tartu University Hospital, Estonia. A two-fold higher representation of > 1 Mb CNVs was observed in men with SPGF (13%, n = 28) compared to controls (6.5%, n = 4). Seven patients with SPGF were identified as carriers of microdeletions (1q21.1; 2.4 Mb) or microduplications (3p26.3, 1.1 Mb; 7p22.3-p22.2, 1.56 Mb; 10q11.22, 1.42 Mb, three cases; Xp22.33; 2.3 Mb) linked to severe congenital conditions. Large autosomal CNV carriers had oligozoospermia, reduced or low-normal bitesticular volume (22-28 ml). The 7p22.3-p22.2 microduplication carrier presented mild intellectual disability, neuropsychiatric problems, and short stature. The Xp22.33 duplication at the PAR1/non-PAR boundary, previously linked to uterine agenesis, was detected in a patient with non-obstructive azoospermia. A novel recurrent intragenic deletion in testis-specific LRRC69 was significantly overrepresented in patients with SPGF compared to the general population (3.3% vs. 0.85%; χ2 test, OR = 3.9 [95% CI 1.8-8.4], P = 0.0001). Assessment of clinically valid CNVs in patients with SPGF will improve their management and counselling for general and reproductive health, including risk of miscarriage and congenital disorders in future offspring.