The HuaBiao project: whole-exome sequencing of 5000 Han Chinese individuals

Natural selection shaped the protective effect of the mtDNA lineage against obesity in Han Chinese populations

Mitochondria play a key role in lipid metabolism, and mitochondrial DNA (mtDNA) mutations are thus considered to affect obesity susceptibility by altering oxidative phosphorylation and mitochondrial function. In this study, we investigated mtDNA variants that may affect obesity risk in 2,877 Han Chinese individuals from three independent populations. The association analysis of 16 basal mtDNA haplogroups with body mass index (BMI), waist circumference (WC) and waist-to-hip ratio (WHR) revealed that only haplogroup M7 was significantly negatively correlated with all three adiposity-related anthropometric traits in the overall cohort (P=0.003 for BMI, P=1×10-5 for WC, P=0.005 for WHR), which was verified by the analysis of a single population, i.e., the Zhengzhou population. Furthermore, subhaplogroup analysis suggested that M7b1a1 was the most likely haplogroup associated with a decreased obesity risk, and the variation T12811C (causing Y159H in ND5) harbored in M7b1a1 may be the most likely candidate for altering mitochondrial function. Specifically, we found that proportionally more nonsynonymous mutations accumulated in M7b1a1 carriers, indicating that M7b1a1 was either under positive selection or subject to a relaxation of selective constraints. We also found that nuclear variants, especially in DACT2 and PIEZO1, may functionally interact with M7b1a1.

Distinguished biological adaptation architecture aggravated population differentiation of Tibeto-Burman-speaking people

Tibeto-Burman (TB) people have endeavored to adapt to the hypoxic, cold, and high-UV high-altitude environments in the Tibetan Plateau and complex disease exposures in lowland rainforests since the late Paleolithic period. However, the full landscape of genetic history and biological adaptation of geographically diverse TB-speaking people, as well as their interaction mechanism, remain unknown. Here, we generate a whole-genome meta-database of 500 individuals from 39 TB-speaking populations and present a comprehensive landscape of genetic diversity, admixture history, and differentiated adaptative features of geographically different TB-speaking people. We identify genetic differentiation related to geography and language among TB-speaking people, consistent with their differentiated admixture process with incoming or indigenous ancestral source populations. A robust genetic connection between the Tibetan-Yi corridor and the ancient Yellow River people supports their Northern China origin hypothesis. We finally report substructure-related differentiated biological adaptative signatures between highland Tibetans and Loloish speakers. Adaptative signatures associated with the physical pigmentation (EDAR and SLC24A5) and metabolism (ALDH9A1) are identified in Loloish people, which differed from the high-altitude adaptative genetic architecture in Tibetan. TB-related genomic resources provide new insights into the genetic basis of biological adaptation and better reference for the anthropologically informed sampling design in biomedical and genomic cohort research.

Functional prediction of the potential NGLY1 mutations associated with rare disease CDG

Genetic diseases are currently diagnosed by functional mutations. However, only some mutations are associated with disease. It is necessary to establish a quick prediction model for clinical screening. Pathogenic mutations in NGLY1 cause a rare autosomal recessive disease known as congenital disorder of deglycosylation (NGLY1-CDDG). Although NGLY1-CDDG can be diagnosed through gene sequencing, clinical relevance of a detected mutation in NGLY1 needs to be further confirmed. In this study, taken NGLY1-CDDG as an example, a comprehensive and practical predictive model for pathogenic mutations on NGLY1 through an NGLY1/Glycopeptide complex model was constructed, the binding sites of NGLY1 and glycopeptides were simulated, and an in vitro enzymatic assay system was established to facilitate quick clinical decisions for NGLY1-CDDG patients. The docking model covers 42 % of reported NGLY1-CDDG missense mutations (5/12). All reported mutations were subjected to in vitro enzymatic assay in which 18 mutations were dysfunctional (18/30). In addition, a full spectrum of functional R328 mutations was assayed and 11 mutations were dysfunctional (11/19). In this study, a model of NGLY1 and glycopeptides was built for potential functional mutations in NGLY1. In addition, the effect of potential regulatory compounds, including N-acetyl-l-cysteine and dithiothreitol, on NGLY1 was examined. The established in vitro assay may serve as a standard protocol to facilitate rapid diagnosis of all mutations in NGLY1-CDDG. This method could also be applied as a comprehensive and practical predictive model for the other rare genetic diseases.

Molecular diagnosis, clinical evaluation and phenotypic spectrum of Townes-Brocks syndrome: insights from a large Chinese hearing loss cohort

BACKGROUND

Townes-Brocks syndrome (TBS) is a rare genetic disorder characterised by multiple malformations. Due to its phenotypic heterogeneity and rarity, diagnosis and recognition of TBS can be challenging and there has been a lack of investigation of patients with atypical TBS in large cohorts and delineation of their phenotypic characteristics.

METHODS

We screened SALL1 and DACT1 variants using next-generation sequencing in the China Deafness Genetics Consortium (CDGC) cohort enrolling 20 666 unrelated hearing loss (HL) cases. Comprehensive clinical evaluations were conducted on seven members from a three-generation TBS family. Combining data from previously reported cases, we also provided a landscape of phenotypes and genotypes of patients with TBS.

RESULTS

We identified five novel and two reported pathogenic/likely pathogenic (P/LP) SALL1 variants from seven families. Audiological features in patients differed in severity and binaural asymmetry. Moreover, previously undocumented malformations in the middle and inner ear were detected in one patient. By comprehensive clinical evaluations, we further provide evidence for the causal relationship between SALL1 variation and certain endocrine abnormalities. Penetrance analysis within familial contexts revealed incomplete penetrance among first-generation patients with TBS and a higher disease burden among their affected offspring.

CONCLUSION

This study presents the first insight of genetic screening for patients with TBS in a large HL cohort. We broadened the phenotypic-genotypic spectrum of TBS and our results supported an underestimated prevalence of TBS. Due to the rarity and phenotypic heterogeneity of rare diseases, broader spectrum molecular tests, especially whole genome sequencing, can improve the situation of underdiagnosis and provide effective recommendations for clinical management.

Both gain- and loss-of-function variants of KCNA1 are associated with paroxysmal kinesigenic dyskinesia

KCNA1 is the coding gene for Kv1.1 voltage-gated potassium-channel α subunit. Three variants of KCNA1 have been reported to manifest as paroxysmal kinesigenic dyskinesia (PKD), but the correlation between them remains unclear due to the phenotypic complexity of KCNA1 variants as well as the rarity of PKD cases. Using the whole exome sequencing followed by Sanger sequencing, we screen for potential pathogenic KCNA1 variants in patients clinically diagnosed with paroxysmal movement disorders and identify three previously unreported missense variants of KCNA1 in three unrelated Chinese families. The proband of one family (c.496G>A, p.A166T) manifests as episodic ataxia type 1, and the other two (c.877G>A, p.V293I and c.1112C>A, p.T371A) manifest as PKD. The pathogenicity of these variants is confirmed by functional studies, suggesting that p.A166T and p.T371A cause a loss-of-function of the channel, while p.V293I leads to a gain-of-function with the property of voltage-dependent gating and activation kinetic affected. By reviewing the locations of PKD-manifested KCNA1 variants in Kv1.1 protein, we find that these variants tend to cluster around the pore domain, which is similar to epilepsy. Thus, our study strengthens the correlation between KCNA1 variants and PKD and provides more information on genotype-phenotype correlations of KCNA1 channelopathy.

The prognostic value of systematic genetic screening in amyotrophic lateral sclerosis patients

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease with complex genetic architecture. Emerging evidence has indicated comorbidity between ALS and autoimmune conditions, suggesting a potential shared genetic basis. The objective of this study is to assess the prognostic value of systematic screening for rare deleterious mutations in genes associated with ALS and aberrant inflammatory responses.

METHODS

A discovery cohort of 494 patients and a validation cohort of 69 patients were analyzed in this study, with population-matched healthy subjects (n = 4961) served as controls. Whole exome sequencing (WES) was performed to identify rare deleterious variants in 50 ALS genes and 1177 genes associated with abnormal inflammatory responses. Genotype-phenotype correlation was assessed, and an integrative prognostic model incorporating genetic and clinical factors was constructed.

RESULTS

In the discovery cohort, 8.1% of patients carried confirmed ALS variants, and an additional 15.2% of patients carried novel ALS variants. Gene burden analysis revealed 303 immune-implicated genes with enriched rare variants, and 13.4% of patients harbored rare deleterious variants in these genes. Patients with ALS variants exhibited a more rapid disease progression (HR 2.87 [95% CI 2.03-4.07], p < 0.0001), while no significant effect was observed for immune-implicated variants. The nomogram model incorporating genetic and clinical information demonstrated improved accuracy in predicting disease outcomes (C-index, 0.749).

CONCLUSION

Our findings enhance the comprehension of the genetic basis of ALS within the Chinese population. It also appears that rare deleterious mutations occurring in immune-implicated genes exert minimal influence on the clinical trajectories of ALS patients.

NLRP14 deficiency causes female infertility with oocyte maturation defects and early embryonic arrest by impairing cytoplasmic UHRF1 abundance

Oocyte maturation is vital to attain full competence required for fertilization and embryogenesis. NLRP14 is preferentially expressed in mammalian oocytes and early embryos. Yet, the role and molecular mechanism of NLRP14 in oocyte maturation and early embryogenesis are poorly understood, and whether NLRP14 deficiency accounts for human infertility is unknown. Here, we found that maternal loss of Nlrp14 resulted in sterility with oocyte maturation defects and early embryonic arrest (EEA). Nlrp14 ablation compromised oocyte competence due to impaired cytoplasmic and nuclear maturation. Importantly, we revealed that NLRP14 maintained cytoplasmic UHRF1 abundance by protecting it from proteasome-dependent degradation and anchoring it from nuclear translocation in the oocyte. Furthermore, we identified compound heterozygous NLRP14 variants in women affected by infertility with EEA, which interrupted the NLRP14-UHRF1 interaction and decreased UHRF1 levels. Our data demonstrate NLRP14 as a cytoplasm-specific regulator of UHRF1 during oocyte maturation, providing insights into genetic diagnosis for female infertility.

Quantitative thresholds for variant enrichment in 13,845 cases: improving pathogenicity classification in genetic hearing loss

BACKGROUND

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) guidelines recommend using variant enrichment among cases as "strong" evidence for pathogenicity per the PS4 criterion. However, quantitative support for PS4 thresholds from real-world Mendelian case-control cohorts is lacking.

METHODS

To address this gap, we evaluated and established PS4 thresholds using data from the Chinese Deafness Genetics Consortium. A total of 9,050 variants from 13,845 patients with hearing loss (HL) and 6,570 ancestry-matched controls were analyzed. Positive likelihood ratio and local positive likelihood ratio values were calculated to determine the thresholds corresponding to each strength of evidence across three variant subsets.

RESULTS

In subset 1, consisting of variants present in both cases and controls with an allele frequency (AF) in cases ≥ 0.0005, an odds ratio (OR) ≥ 6 achieved strong evidence, while OR ≥ 3 represented moderate evidence. For subset 2, which encompassed variants present in both cases and controls with a case AF < 0.0005, and subset 3, comprising variants found only in cases and absent from controls, we defined the PS4_Supporting threshold (OR > 2.27 or allele count ≥ 3) and the PS4_Moderate threshold (allele count ≥ 6), respectively. Reanalysis applying the adjusted PS4 criteria changed the classification of 15 variants and enabled diagnosis of an additional four patients.

CONCLUSIONS

Our study quantified evidence strength thresholds for variant enrichment in genetic HL cases, highlighting the importance of defining disease/gene-specific thresholds to improve the precision and accuracy of clinical genetic testing.

The complete and fully-phased diploid genome of a male Han Chinese

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

Large-scale analysis of de novo mutations identifies risk genes for female infertility characterized by oocyte and early embryo defects

BACKGROUND

Oocyte maturation arrest and early embryonic arrest are important reproductive phenotypes resulting in female infertility and cause the recurrent failure of assisted reproductive technology (ART). However, the genetic etiologies of these female infertility-related phenotypes are poorly understood. Previous studies have mainly focused on inherited mutations based on large pedigrees or consanguineous patients. However, the role of de novo mutations (DNMs) in these phenotypes remains to be elucidated.

RESULTS

To decipher the role of DNMs in ART failure and female infertility with oocyte and embryo defects, we explore the landscape of DNMs in 473 infertile parent-child trios and identify a set of 481 confident DNMs distributed in 474 genes. Gene ontology analysis reveals that the identified genes with DNMs are enriched in signaling pathways associated with female reproductive processes such as meiosis, embryonic development, and reproductive structure development. We perform functional assays on the effects of DNMs in a representative gene Tubulin Alpha 4a (TUBA4A), which shows the most significant enrichment of DNMs in the infertile parent-child trios. DNMs in TUBA4A disrupt the normal assembly of the microtubule network in HeLa cells, and microinjection of DNM TUBA4A cRNAs causes abnormalities in mouse oocyte maturation or embryo development, suggesting the pathogenic role of these DNMs in TUBA4A.

CONCLUSIONS

Our findings suggest novel genetic insights that DNMs contribute to female infertility with oocyte and embryo defects. This study also provides potential genetic markers and facilitates the genetic diagnosis of recurrent ART failure and female infertility.

Presence of Rare Variants is Associated with Poorer Survival in Chinese Patients with Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder with phenotypic and genetic heterogeneity. Recent studies have suggested an oligogenic basis of ALS, in which the co-occurrence of two or more genetic variants has additive or synergistic deleterious effects. To assess the contribution of possible oligogenic inheritance, we profiled a panel of 43 relevant genes in 57 sporadic ALS (sALS) patients and eight familial ALS (fALS) patients from five pedigrees in east China. We filtered rare variants using the combination of the Exome Aggregation Consortium, the 1000 Genomes and the HuaBiao Project. We analyzed patients with multiple rare variants in 43 known ALS causative genes and the genotype-phenotype correlation. Overall, we detected 30 rare variants in 16 different genes and found that 16 of the sALS patients and all the fALS patients examined harbored at least one variant in the investigated genes, among which two sALS and four fALS patients harbored two or more variants. Of note, the sALS patients with one or more variants in ALS genes had worse survival than the patients with no variants. Typically, in one fALS pedigree with three variants, the family member with three variants (Superoxide dismutase 1 (SOD1) p.V48A,  Optineurin (OPTN) p.A433V and TANK binding kinase 1 (TBK1) p.R573H) exhibited much more severe disease phenotype than the member carrying one variant (TBK1 p.R573H). Our findings suggest that rare variants could exert a negative prognostic effect, thereby supporting the oligogenic inheritance of ALS.

Landscape of pathogenic mutations in premature ovarian insufficiency

Premature ovarian insufficiency (POI) is a major cause of female infertility due to early loss of ovarian function. POI is a heterogeneous condition, and its molecular etiology is unclear. To identify genetic variants associated with POI, here we performed whole-exome sequencing in a cohort of 1,030 patients with POI. We detected 195 pathogenic/likely pathogenic variants in 59 known POI-causative genes, accounting for 193 (18.7%) cases. Association analyses comparing the POI cohort with a control cohort of 5,000 individuals without POI identified 20 further POI-associated genes with a significantly higher burden of loss-of-function variants. Functional annotations of these novel 20 genes indicated their involvement in ovarian development and function, including gonadogenesis (LGR4 and PRDM1), meiosis (CPEB1, KASH5, MCMDC2, MEIOSIN, NUP43, RFWD3, SHOC1, SLX4 and STRA8) and folliculogenesis and ovulation (ALOX12, BMP6, H1-8, HMMR, HSD17B1, MST1R, PPM1B, ZAR1 and ZP3). Cumulatively, pathogenic and likely pathogenic variants in known POI-causative and novel POI-associated genes contributed to 242 (23.5%) cases. Further genotype-phenotype correlation analyses indicated that genetic contribution was higher in cases with primary amenorrhea compared to that in cases with secondary amenorrhea. This study expands understanding of the genetic landscape underlying POI and presents insights that have the potential to improve the utility of diagnostic genetic screenings.

PGG.MHC: toward understanding the diversity of major histocompatibility complexes in human populations

The human leukocyte antigen (HLA) system, or the human version of the major histocompatibility complex (MHC), is known for its extreme polymorphic nature and high heterogeneity. Taking advantage of whole-genome and whole-exome sequencing data, we developed PGG.MHC to provide a platform to explore the diversity of the MHC in Asia as well as in global populations. PGG.MHC currently archives high-resolution HLA alleles of 53 254 samples representing 190 populations spanning 66 countries. PGG.MHC provides: (i) high-quality allele frequencies for eight classical HLA loci (HLA-A, -B, -C, -DQA1, -DQB1, -DRB1, -DPA1 and -DPB1); (ii) visualization of population prevalence of HLA alleles on global, regional, and country-wide levels; (iii) haplotype structure of 134 populations; (iv) two online analysis tools including 'HLA imputation' for inferring HLA alleles from SNP genotyping data and 'HLA association' to perform case/control studies for HLA-related phenotypes and (v) East Asian-specific reference panels for HLA imputation. Equipped with high-quality frequency data and user-friendly computer tools, we expect that the PGG.MHC database can advance the understanding and facilitate applications of MHC genomic diversity in both evolutionary and medical studies. The PGG.MHC database is freely accessible via https://pog.fudan.edu.cn/pggmhc or https://www.pggmhc.org/pggmhc.

An additional whole-exome sequencing study in 102 panel-undiagnosed patients: A retrospective study in a Chinese craniosynostosis cohort

Craniosynostosis (CRS) is a disease with prematurely fused cranial sutures. In the last decade, the whole-exome sequencing (WES) was widely used in Caucasian populations. The WES largely contributed in genetic diagnosis and exploration on new genetic mechanisms of CRS. In this study, we enrolled 264 CRS patients in China. After a 17-gene-panel sequencing designed in the previous study, 139 patients were identified with pathogenic/likely pathogenic (P/LP) variants according to the ACMG guideline as positive genetic diagnosis. WES was then performed on 102 patients with negative genetic diagnosis by panel. Ten P/LP variants were additionally identified in ten patients, increasing the genetic diagnostic yield by 3.8% (10/264). The novel variants in ANKH, H1-4, EIF5A, SOX6, and ARID1B expanded the mutation spectra of CRS. Then we designed a compatible research pipeline (RP) for further exploration. The RP could detect all seven P/LP SNVs and InDels identified above, in addition to 15 candidate variants found in 13 patients with worthy of further study. In sum, the 17-gene panel and WES identified positive genetic diagnosis for 56.4% patients (149/264) in 16 genes. At last, in our estimation, the genetic testing strategy of “Panel-first” saves 24.3% of the cost compared with “WES only”, suggesting the “Panel-first” is an economical strategy.

Estimation of hereditary fructose intolerance prevalence in the Chinese population

Background

Hereditary fructose intolerance (HFI) caused by aldolase B reduction or deficiency that results in fructose metabolism disorder. The disease prevalence in the Chinese population is unknown, which impedes the formulation of HFI screening and diagnosis strategies.

Materials and methods

By searching a local cohort (Chinese Children’s Rare Disease Genetic Testing Clinical Collaboration System, CCGT) and public databases (ClinVar and Human Gene Mutation Database) and reviewing HFI-related literature, we manually curated ALDOB pathogenic or likely pathogenic (P/LP) variants according to ACMG guidelines. Allele frequency (AF) information from the local database CCGT and the public databases HuaBiao and gnomAD for ALDOB P/LP variants was used to estimate and the HFI prevalence in the Chinese population and other populations by the Bayesian framework. We collected the genotype and clinical characteristics of HFI patients from the CCGT database and published literature to study genotype–phenotype relationships.

Result

In total, 81 variants of ALDOB were curated as P/LP. The estimated Chinese HFI prevalence was approximately 1/504,678, which was much lower than that for non-Finland European (1/23,147), Finnish in Finland (1/55,539), admixed American (1/132,801) and Ashkenazi Jewish (1/263,150) populations. By analyzing the genetic characteristics of ALDOB in the Chinese population, two variants (A338V, A338G) had significantly higher AFs in the Chinese population than in the non-Finland European population from gnomAD (all P values < 0.05). Five variants (A150P, A175D, N335K, R60*, R304Q) had significantly lower AFs (all P values < 0.1). The genotype–phenotype association analyses were based on 68 reported HFI patients from a literature review and the CCGT database. The results showed that patients carrying homozygous variant sites (especially A150P) were more likely to present nausea, and patients carrying two missense variant sites were more likely to present aversion to sweets and fruit (all P values < 0.05). Our research reveals that some gastrointestinal symptoms seem to be associated with certain genotypes.

Conclusion

The prevalence of HFI in the Chinese population is extremely low, and there is no need to add HFI testing to the current newborn screening programs if medical costs are considered. A genetic testing strategy is suggested for early diagnosis of HFI.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-022-02487-3.

Genetic evidence of tri-genealogy hypothesis on the origin of ethnic minorities in Yunnan

BACKGROUND

Yunnan is located in Southwest China and consists of great cultural, linguistic, and genetic diversity. However, the genomic diversity of ethnic minorities in Yunnan is largely under-investigated. To gain insights into population history and local adaptation of Yunnan minorities, we analyzed 242 whole-exome sequencing data with high coverage (~ 100-150 ×) of Yunnan minorities representing Achang, Jingpo, Dai, and Deang, who were linguistically assumed to be derived from three ancient lineages (the tri-genealogy hypothesis), i.e., Di-Qiang, Bai-Yue, and Bai-Pu.

RESULTS

Yunnan minorities show considerable genetic differences. Di-Qiang populations likely migrated from the Tibetan area about 6700 years ago. Genetic divergence between Bai-Yue and Di-Qiang was estimated to be 7000 years, and that between Bai-Yue and Bai-Pu was estimated to be 5500 years. Bai-Pu is relatively isolated, but gene flow from surrounding Di-Qiang and Bai-Yue populations was also found. Furthermore, we identified genetic variants that are differentiated within Yunnan minorities possibly due to the living circumstances and habits. Notably, we found that adaptive variants related to malaria and glucose metabolism suggest the adaptation to thalassemia and G6PD deficiency resulting from malaria resistance in the Dai population.

CONCLUSIONS

We provided genetic evidence of the tri-genealogy hypothesis as well as new insights into the genetic history and local adaptation of the Yunnan minorities.

MECOM-related disorder: Radioulnar synostosis without hematological aberration due to unique variants

PURPOSE

The etiology for a considerable proportion of patients with congenital radioulnar synostosis (RUS) remains unclear. This study aimed to investigate the genetic cause of RUS without a known cause.

METHODS

Patients with RUS were investigated. Exome sequencing and/or Sanger sequencing was performed. Bioinformatics analysis was also performed. Pathogenicity was evaluated for variants of interest.

RESULTS

We identified unique missense variants in MECOM (encodes EVI1) associated with RUS in 8 families. Of them, 6 families had variants in residue R781, including 3 families with R781C (c.2341C>T), 2 families with R781H (c.2342G>A), and 1 family with R781L (c.2342G>T). Another 2 variants included I783T (c.2348T>C) in 1 family and Q777E (c.2329C>G) in 1 family. All these variants were clustered within the ninth zinc finger motif of EVI1. Phenotype evaluation identified that most of these patients with RUS harboring mutant MECOM had finger malformations, but none of them had identifiable hematological abnormalities. Functional experiments showed that MECOM R781C led to alterations in TGF-β-mediated transcriptional responses.

CONCLUSION

This study examined MECOM variants by focusing on RUS instead of hematological abnormalities. The R781 residue in EVI1 is a hotspot for human RUS variants. Mutant MECOM is the second most common cause for familial RUS.

Deleterious Rare Mutations of GLI1 Dysregulate Sonic Hedgehog Signaling in Human Congenital Heart Disease

The Glioma-associated oncogene (Gli) family members of zinc finger DNA-binding proteins are core effectors of Sonic hedgehog (SHH) signaling pathway. Studies in model organisms have identified that the Gli genes play critical roles during organ development, including the heart, brain, kidneys, etc. Deleterious mutations in GLI genes have previously been revealed in several human developmental disorders, but few in congenital heart disease (CHD). In this study, the mutations in GLI1-3 genes were captured by next generation sequencing in human cohorts composed of 412 individuals with CHD and 213 ethnically matched normal controls. A total of 20 patient-specific nonsynonymous rare mutations in coding regions of human GLI1-3 genes were identified. Functional analyses showed that GLI1 c.820G> T (p.G274C) is a gain-of-function mutation, while GLI1 c.878G>A (p.R293H) and c.1442T>A (p.L481X) are loss-of-function mutations. Our findings suggested that deleterious rare mutations in GLI1 gene broke the balance of the SHH signaling pathway regulation and may constitute a great contribution to human CHD, which shed new light on understanding genetic mechanism of embryo cardiogenesis regulated by SHH signaling.