A small number of de novo assembled human genomes have been reported to date, and few have been complemented with population-based genetic variation, which is particularly important for North Africa, a region underrepresented in current genome-wide references. Here, we combine long- and short-read whole-genome sequencing data with recent assembly approaches into a de novo assembly of an Egyptian genome. The assembly demonstrates well-balanced quality metrics and is complemented with variant phasing via linked reads into haploblocks, which we associate with gene expression changes in blood. To construct an Egyptian genome reference, we identify genome-wide genetic variation within a cohort of 110 Egyptian individuals. We show that differences in allele frequencies and linkage disequilibrium between Egyptians and Europeans may compromise the transferability of European ancestry-based genetic disease risk and polygenic scores, substantiating the need for multi-ethnic genome references. Thus, the Egyptian genome reference will be a valuable resource for precision medicine.

Even though immunoglobulins are critical for immune responses and human survival, the diversity of the immunoglobulin heavy chain gene (IGH) is poorly known and mostly characterized only by serological methods. Moreover, this genomic region is not well-covered in genomic databases and genome-wide association studies due to particularities that impose technical difficulties for its analysis. Therefore, the IGH gene has never been systematically sequenced across populations. Here, we deliver an unprecedented and comprehensive characterization of the diversity of the IGHG1, IGHG2, and IGHG3 gene segments, which encode the constant region of the most abundant circulating immunoglobulins: IgG1, IgG2, and IgG3, respectively. We used Sanger sequencing to analyze 357 individuals from seven different Brazilian populations, including five Amerindian, one Japanese-descendant and one Euro-descendant population samples. We discovered 28 novel IGHG alleles and provided evidence that some of them may have been originated by gene conversion between common alleles of different gene segments. The rate of synonymous substitutions was significantly higher than the rate of the non-synonymous substitutions for IGHG1 and IGHG2 (p = 0.01 and 0.03, respectively), consistent with purifying selection. Fay and Wu's test showed significant negative values for most populations (p < 0.001), which indicates that positive selection in an adjacent position may be shaping IGHG variation by hitchhiking of variants in the vicinity, possibly the regions that encode the Ig variable regions. This study shows that the variation in the IGH gene is largely underestimated. Therefore, exploring its nucleotide diversity in populations may provide valuable information for comprehension of its evolution, its impact on diseases and vaccine research.

De novo variants in the myelin regulatory factor (MYRF), a transcription factor involved in the differentiation of oligodendrocytes, have been linked recently to the cardiac and urogenital syndrome, while familiar variants are associated with nanophthalmos. Here, we report for the first time on a patient with a de novo stop-gain variant in MYRF (p.Q838*) associated with Scimitar syndrome, 46,XY partial gonadal dysgenesis (GD) and severe hyperopia. Since variants in MYRF have been described in both 46,XX and 46,XY GD, we assumed a role of MYRF in the early development of the bipotential gonad. We used publicly available single cell sequencing data of human testis and ovary from different developmental stages and analysed them for MYRF expression. We identified MYRF expression in the subset of coelomic epithelial cells at stages of gonadal ridge development in 46,XX and 46,XY individuals. Differential gene expression analysis revealed significantly upregulated genes. Within these, we identified CITED2 as a gene containing a MYRF binding site. It has been shown that Cited2^−/− mice have gonadal defects in both testis and ovary differentiation, as well as defects in heart development and establishment of the left–right axis. This makes MYRF a potential candidate as an early regulator of gonadal and heart development via upregulation of the transcriptional cofactor CITED2.

Congenital defects of androgen synthesis or action in 46, XY individuals can result in impaired virilisation, despite the apparent testicular development. In a recent case report of a young adult with complete androgen insensitivity syndrome (CAIS), tumorous gonadal tissue was shown to express HSD17B3 in Sertoli cells (SCs) and not in Leydig cells (LCs). This expression pattern differs from the typical adult human testis and resembles foetal mouse testis, suggesting an underlying testicular development and function defect. Here, we investigate the effect of altered androgen signalling in gonads from five 46, XY individuals with defects in androgen synthesis or action. Gonadal tissue sections from four patients with CAIS, one with CYP17A1 deficiency, and one control were immunostained for LC developmental and steroidogenic markers. The expression of some of these markers during development was investigated by reanalysing previously published single-cell RNA sequencing (scRNA-seq) data from normal human testicular tissues. All gonadal tissues from the patients show exclusive expression of HSD17B3 in SCs and expression of the foetal/immature LC marker DLK1 in a subset of LCs, suggesting an androgen-dependent differentiation defect of adult SCs and LCs. Furthermore, reanalysis of scRNA-seq data reveals an expression of HSD17B3 in foetal and neonatal SCs that is downregulated in adult SCs. Androgen signalling may affect the differentiation of adult but possibly not foetal SCs or LCs and may induce a shift of testosterone production from the tubular compartment in the foetal phase to the interstitial compartment in the adult phase.