Disruptions in a cluster of computationally identified enhancers near FOXC1 and GMDS may influence brain development

Genetic Variants Associated with Suspected Neonatal Hypoxic Ischaemic Encephalopathy: A Study in a South African Context

Neonatal encephalopathy suspected to be due to hypoxic ischaemic encephalopathy (NESHIE) carries the risk of death or severe disability (cognitive defects and cerebral palsy). Previous genetic studies on NESHIE have predominantly focused on exomes or targeted genes. The objective of this study was to identify genetic variants associated with moderate-severe NESHIE through whole-genome, unbiased analysis. Variant filtering and prioritization were performed, followed by association testing both on a case-control basis and to compare the grades of severity and/or progression. Association testing on neonates with NESHIE (N = 172) and ancestry-matched controls (N = 288) produced 71 significant genetic variants (false discovery rate corrected p-value < 6.2 × 10-4), all located in non-coding regions and not previously implicated in NESHIE. Disease-associated variants in non-coding regions are considered to affect regulatory functions, possibly by modifying gene expression, promoters, enhancers, or DNA structure. The most significant variant was at position 6:162010973 in the Parkin RBR E3 ubiquitin protein ligase (PRKN) intron. Intronic variants were also identified in genes involved in inflammatory processes (SLCO3A1), DNA repair (ZGRF1), synaptogenesis (CNTN5), haematopoiesis (ASXL2), and the transcriptional response to hypoxia (PADI4). Ten variants were associated with a higher severity or lack of improvement in NESHIE, including one in ADAMTS3, which encodes a procollagen amino protease with a role in angiogenesis and lymphangiogenesis. This analysis represents one of the first efforts to analyze whole-genome data to investigate the genetic complexity of NESHIE in diverse ethnolinguistic groups of African origin and provides direction for further study.

Diagnostic Utility of Whole Genome Sequencing After Negative Karyotyping/Chromosomal Microarray in Infants Born With Multiple Congenital Anomalies

BACKGROUND

Achieving a definitive genetic diagnosis of unexplained multiple congenital anomalies (MCAs) in neonatal intensive care units (NICUs) infants is challenging because of the limited diagnostic capabilities of conventional genetic tests. Although the implementation of whole genome sequencing (WGS) has commenced for diagnosing MCAs, due to constraints in resources and faculty, many NICUs continue to utilize chromosomal microarray (CMA) and/or karyotyping as the initial diagnostic approach. We aimed to evaluate the diagnostic efficacy of WGS in infants with MCAs who have received negative results from karyotyping and/or CMA.

METHODS

In this prospective study, we enrolled 80 infants with MCAs who were admitted to a NICU at a single center and had received negative results from CMA and/or karyotyping. The phenotypic characteristics were classified according to the International Classification of Diseases and the Human Phenotype Ontology. We assessed the diagnostic yield of trio-WGS in infants with normal chromosomal result and explored the process of diagnosing by analyzing both phenotype and genotype. Also, we compared the phenotype and clinical outcomes between the groups diagnosed with WGS and the undiagnosed group.

RESULTS

The diagnostic yield of WGS was 26% (21/80), of which 76% were novel variants. There was a higher diagnostic yield in cases of craniofacial abnormalities, including those of the eye and ear, and a lower diagnostic yield in cases of gastrointestinal and genitourinary abnormalities. In addition, higher rates of rehabilitation therapy and gastrostomy were observed in WGS-diagnosed infants than in undiagnosed infants.

CONCLUSION

This prospective cohort study assessed the usefulness of trio-WGS following chromosomal analysis for diagnosing MCAs in the NICU and revealed improvements in the diagnostic yield and clinical utility of WGS.

Genetic diversity, demographic history, and selective signatures of Silkie chicken

BACKGROUND

Silkie is a traditional Chinese chicken breed characterized by its unique combination of specialized morphological traits. While previous studies have focused on the genetic basis of these traits, the overall genomic characteristics of the Silkie breed remain largely unexplored. In this study, we employed whole genome resequencing data to examine the genetic diversity, selective signals and demographic history of the Silkie breed through comparative analyses with seven other Chinese indigenous breeds (IDGBs), a commercial breed, and the wild ancestor Red Jungle Fowl.

RESULTS

In total, 20.8 million high-quality single nucleotide polymorphisms and 86 large structural variations were obtained. We discovered that Silkie exhibits a relatively high level of inbreeding and is genetically distinct from other IDGBs. Furthermore, our analysis indicated that Silkie has experienced a stronger historical population bottleneck and has a smaller effective population size compared with other IDGBs. We identified 45 putatively selected genes that are enriched in the melanogenesis pathway, which probably is related to the feather color. Among these genes, LMBR1 and PDSS2 have been previously associated with the extra toe and the hookless feathers, respectively. Six of the selected genes (KITLG, GSK3B, SOBP, CTBP1, ELMO2, SNRPN) are known to be associated with neurodevelopment and mental diseases in human, and are possibly related to the distinct behavior of Silkie. We further identified structural variants in Silkie and found previously reported variants linked to hyperpigmentation (END3), muff and beard (HOXB8), and Rose-comb phenotype (MNR2). Additionally, we found a 0.61 Mb inversion overlapping with the GMDS gene, which was previously linked to neurodevelopmental defects in zebrafish and humans. This may also be related to the behavior distinctiveness of Silkie.

CONCLUSIONS

Our study revealed that Silkie is genetically distinct and relatively highly inbred compared to other IDGB chicken populations, possibly attributed to more prolong population bottlenecks and selective breeding practice. These results enhance our understanding of how domestication and selective breeding have shaped the genome of Silkie. These findings contribute to the broader field of domestication and avian genomics, and have implications for the future conservation and breeding efforts.

Novel maternal duplication of 6p22.3-p25.3 with subtelomeric 6p25.3 deletion: new clinical findings and genotype-phenotype correlations

BACKGROUND

Copy-number variants (CNVs) drive many neurodevelopmental-related disorders. Although many neurodevelopmental-related CNVs can give rise to widespread phenotypes, it is necessary to identify the major genes contributing to phenotypic presentation. Copy-number variations in chromosome 6, such as independent 6p deletion and 6p duplication, have been reported in several live-born infants and present widespread abnormalities such as intellectual disability, growth deficiency, developmental delay, and multiple dysmorphic facial features. However, a contiguous deletion and duplication in chromosome 6p regions have been reported in only a few cases.

CASE PRESENTATION

In this study, we reported the first duplication of chromosome band 6p25.3-p22.3 with deletion of 6p25.3 in a pedigree. This is the first case reported involving CNVs in these chromosomal regions. In this pedigree, we reported a 1-year-old boy with maternal 6p25-pter duplication characterized by chromosome karyotype. Further analysis using CNV-seq revealed a 20.88-Mb duplication at 6p25.3-p22.3 associated with a contiguous 0.66-Mb 6p25.3 deletion. Whole exome sequencing confirmed the deletion/duplication and identified no pathogenic or likely pathogenic variants related with the patient´s phenotype. The proband presented abnormal growth, developmental delay, skeletal dysplasia, hearing loss, and dysmorphic facial features. Additionally, he presented recurrent infection after birth. CNV-seq using the proband´s parental samples showed that the deletion/duplication was inherited from the proband´s mother, who exhibited a similar phenotype to the proband. When compared with other cases, this proband and his mother presented a new clinical finding: forearm bone dysplasia. The major candidate genes contributing to recurrent infection, eye development, hearing loss features, neurodevelopmental development, and congenital bone dysplasia were further discussed.

CONCLUSIONS

Our results showed a new clinical finding of a contiguous deletion and duplication in chromosome 6p regions and suggested candidate genes associated with phenotypic features, such as FOXC1, SERPINB6, NRN1, TUBB2A, IRF4, and RIPK1.

CRISPR-Cas9-mediated functional dissection of the foxc1 genomic region in zebrafish identifies critical conserved cis-regulatory elements

FOXC1 encodes a forkhead-domain transcription factor associated with several ocular disorders. Correct FOXC1 dosage is critical to normal development, yet the mechanisms controlling its expression remain unknown. Together with FOXQ1 and FOXF2, FOXC1 is part of a cluster of FOX genes conserved in vertebrates. CRISPR-Cas9-mediated dissection of genomic sequences surrounding two zebrafish orthologs of FOXC1 was performed. This included five zebrafish-human conserved regions, three downstream of foxc1a and two remotely upstream of foxf2a/foxc1a or foxf2b/foxc1b clusters, as well as two intergenic regions between foxc1a/b and foxf2a/b lacking sequence conservation but positionally corresponding to the area encompassing a previously reported glaucoma-associated SNP in humans. Removal of downstream sequences altered foxc1a expression; moreover, zebrafish carrying deletions of two or three downstream elements demonstrated abnormal phenotypes including enlargement of the anterior chamber of the eye reminiscent of human congenital glaucoma. Deletions of distant upstream conserved elements influenced the expression of foxf2a/b or foxq1a/b but not foxc1a/b within each cluster. Removal of either intergenic sequence reduced foxc1a or foxc1b expression during late development, suggesting a role in transcriptional regulation despite the lack of conservation at the nucleotide level. Further studies of the identified regions in human patients may explain additional individuals with developmental ocular disorders.

Copy number variants in genomes of local sheep breeds from Russia

Copy number variants (CNVs) are genomic structural variations that contribute to many adaptive and economically important traits in livestock. In this study, we detected CNVs in 354 animals from 16 Russian indigenous sheep breeds and analysed their possible functional roles. Our analysis of the entire sample set resulted in 4527 CNVs forming 1450 CNV regions (CNVRs). When constructing CNVRs for individual breeds, a total of 2715 regions ranging from 88 in Groznensk to 337 in Osetin breeds were identified. To make interbreed CNVR frequency comparison possible, we also identified core CNVRs using CNVs with overlapping chromosomal locations found in different breeds. This resulted in 137 interbreed CNVRs with frequency >15% in at least one breed. Functional enrichment analysis of genes affected by CNVRs in individual breeds revealed 12 breeds with significant enrichments in olfactory perception, PRAME family proteins, and immune response. Function of genes affected by interbreed and breed-specific CNVRs revealed candidates related to domestication, adaptation to high altitudes and cold climates, reproduction, parasite resistance, milk and meat qualities, wool traits, fat storage, and fat metabolism. Our work is the first attempt to uncover and characterise the CNV makeup of Russian indigenous sheep breeds. Further experimental and functional validation of CNVRs would help in developing new and improving existing sheep breeds.

Cis-Repression of Foxq1 Expression Affects Foxf2-Mediated Gene Expression in Palate Development

Disruption of FOXF2, encoding a member of the Forkhead family transcription factors, has been associated with cleft palate in humans and mice. FOXF2 is located in a conserved gene cluster containing FOXQ1, FOXF2, and FOXC1. We found that expression of Foxq1 is dramatically upregulated in the embryonic palatal mesenchyme in Foxf2 ^-/- mouse embryos. We show here that the Foxf2 promoter-deletion mutation caused dramatically increased expression of the cis-linked Foxq1 allele but had little effect on the Foxq1 allele in trans. We analyzed effects of the Foxf2 mutation on the expression of other neighboring genes and compared those effects with the chromatin domain structure and recently identified enhancer-promoter associations as well as H3K27ac ChIP-seq data. We show that the Foxf2 mutation resulted in significantly increased expression of the Foxq1 and Exoc2 genes located in the same topologically associated domain with Foxf2 but not the expression of the Foxc1 and Gmds genes located in the adjacent chromatin domain. We inactivated the Foxq1 gene in mice homozygous for a Foxf2 conditional allele using CRISPR genome editing and generated (Foxf2/Foxq1)^+/- mice with loss-of-function mutations in Foxf2 and Foxq1 in cis. Whereas the (Foxf2/Foxq1)^-/- mice exhibited cleft palate at birth similar as in the Foxf2 ^-/- mice, systematic expression analyses of a large number of Foxf2-dependent genes revealed that the (Foxf2/Foxq1)^-/- embryos exhibited distinct effects on the domain-specific expression of several important genes, including Foxf1, Shox2, and Spon1, in the developing palatal shelves compared with Foxf2 ^-/- embryos. These results identify a novel cis-regulatory effect of the Foxf2 mutation and demonstrate that cis-regulation of Foxq1 contributed to alterations in palatal gene expression in Foxf2 ^-/- embryos. These results have important implications for interpretation of results and mechanisms from studies of promoter- or gene-deletion alleles. In addition, the unique mouse lines generated in this study provide a valuable resource for understanding the cross-regulation and combinatorial functions of the Foxf2 and Foxq1 genes in development and disease.

Extracellular Vesicles and the Oviduct Function

In mammals, the oviduct (or the Fallopian tube in humans) can be divided into the infundibulum (responsible for oocyte pick-up), ampulla (site of fertilization), isthmus (where preimplantation embryos develop), and uterotubal junction (where embryos transit to the uterus). The oviductal fluid, as well as extracellular vesicles produced from the oviduct epithelial cells, referred to as oEVs, have been shown to improve the fertilization process, prevent polyspermy, and aid in embryo development. oEVs contain molecular cargos (such as miRNAs, mRNAs, proteins, and lipids) that can be delivered and fuse to recipient cells. oEVs produced from the ampulla appear to be functionally distinct from those produced from the isthmus. In multiple species including mice, cats, dogs, pigs, and cows, oEVs can be incorporated into the oocytes, sperm, and embryos. In this review, we show the positive impact of oEVs on gamete function as well as blastocyst development and how they may improve embryo quality in in vitro conditions in an assisted reproductive technology setting for rodents, domestic animals, farm animals, and humans.

The distributions of protein coding genes within chromatin domains in relation to human disease

Background

Our understanding of the nuclear chromatin structure has increased hugely during the last years mainly as a consequence of the advances in chromatin conformation capture methods like Hi-C. The unprecedented resolution of genome-wide interaction maps shows functional consequences that extend the initial thought of an efficient DNA packaging mechanism: gene regulation, DNA repair, chromosomal translocations and evolutionary rearrangements seem to be only the peak of the iceberg. One key concept emerging from this research is the topologically associating domains (TADs) whose functional role in gene regulation and their association with disease is not fully untangled.

Results

We report that the lower the number of protein coding genes inside TADs, the higher the tendency of those genes to be associated with disease (p-value = 4 × \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-54}$$\end{document}10-54). Moreover, housekeeping genes are less associated with disease than other genes. Accordingly, they are depleted in TADs containing less than three protein coding genes (p-value = 3.9 × \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-34}$$\end{document}10-34). We observed that TADs with higher ratios of enhancers versus genes contained higher numbers of disease-associated genes. We interpret these results as an indication that sharing enhancers among genes reduces their involvement in disease. Larger TADs would have more chances to accommodate many genes and select for enhancer sharing along evolution.

Conclusions

Genes associated with human disease do not distribute randomly over the TADs. Our observations suggest general rules that confer functional stability to TADs, adding more evidence to the role of TADs as regulatory units.

Axenfeld-Rieger syndrome

Axenfeld-Rieger syndrome (ARS) is a clinically and genetically heterogeneous group of developmental disorders affecting primarily the anterior segment of the eye, often leading to secondary glaucoma. Patients with ARS may also present with systemic changes, including dental defects, mild craniofacial dysmorphism, and umbilical anomalies. ARS is inherited in an autosomal-dominant fashion; the underlying defect in 40% of patients is mutations in PITX2 or FOXC1. Here, an overview of the clinical spectrum of ARS is provided. As well, the known underlying genetic defects, clinical diagnostic possibilities, genetic counseling and treatments of ARS are discussed in detail.

Embryo culture in presence of oviductal fluid induces DNA methylation changes in bovine blastocysts

During the transit through the oviduct, the early embryo initiates an extensive DNA methylation reprogramming of its genome. Given that these epigenetic modifications are susceptible to environmental factors, components present in the oviductal milieu could affect the DNA methylation marks of the developing embryo. The aim of this study was to examine if culture of bovine embryos with oviductal fluid (OF) can induce DNA methylation changes at specific genomic regions in the resulting blastocysts. In vitro produced zygotes were cultured in medium with 3 mg/mL bovine serum albumin (BSA) or 1.25% OF added at the one- to 16-cell stage (OF1–16), one- to 8-cell stage (OF1–8) or 8- to 16-cell stage (OF8–16), and then were cultured until Day 8 in medium with 3 mg/mL BSA. Genomic regions in four developmentally important genes (MTERF2, ABCA7, OLFM1, GMDS) and within LINE-1 retrotransposons were selected for methylation analysis by bisulfite sequencing on Day 7–8 blastocysts. Blastocysts derived from OF1–16 group showed lower CpG methylation levels in MTERF2 and ABCA7 compared with the BSA group. However, CpG sites within MTERF2, ABCA7 and OLFM1 showed higher methylation levels in groups OF1–8 and OF8–16 than in OF1–16. For LINE-1 elements, higher CpG methylation levels were observed in blastocysts from the OF1–16 group than in the other experimental groups. In correlation with the methylation changes observed, mRNA expression level of MTERF2 was increased, while LINE-1 showed a decreased expression in blastocysts from OF1–16 group. Our results suggest that embryos show transient sensitivity to OF at early stages, which is reflected by specific methylation changes at the blastocyst stage.

Genome-wide association analysis of actigraphic sleep phenotypes in the LIFE Adult Study

The genetic basis of sleep is still poorly understood. Despite the moderate to high heritability of sleep-related phenotypes, known genetic variants explain only a small proportion of the phenotypical variance. However, most previous studies were based solely upon self-report measures. The present study aimed to conduct the first genome-wide association (GWA) of actigraphic sleep phenotypes. The analyses included 956 middle- to older-aged subjects (40-79 years) from the LIFE Adult Study. The SenseWear Pro 3 Armband was used to collect 11 actigraphic parameters of night- and daytime sleep and three parameters of rest (lying down). The parameters comprised measures of sleep timing, quantity and quality. A total of 7 141 204 single nucleotide polymorphisms (SNPs) were analysed after imputation and quality control. We identified several variants below the significance threshold of P ≤ 5× 10^-8 (not corrected for analysis of multiple traits). The most significant was a hit near UFL1 associated with sleep efficiency on weekdays (P = 1.39 × 10^-8 ). Further SNPs were close to significance, including an association between sleep latency and a variant in CSNK2A1 (P = 8.20 × 10^-8 ), a gene known to be involved in the regulation of circadian rhythm. In summary, our GWAS identified novel candidate genes with biological plausibility being promising candidates for replication and further follow-up studies.