PubMed Health⌕ Search

Biomedical subjects

Mary L Marazita

Publications and source records attributed to Mary L Marazita.

At least 19 recordsLinked to original sources

Trio-based GWAS reveals loci associated with different forms of isolated cleft lip.

Orofacial clefts (OFCs) are the most common craniofacial birth defect and comprise a diverse group of traits with complex and heterogeneous etiologies. Genetic studies of OFCs typically approach this diversity by stratifying cases into broad diagnostic classes, including cleft lip (CL), cleft palate (CP), and cleft lip with palate (CLP). Although this strategy has yielded important insights into OFC risk, it ignores the phenotypic heterogeneity within each subtype. CL exhibits marked phenotypic variability, involving differences in alveolar involvement, laterality, and sidedness that may reflect distinct etiologies. Given this phenotypic diversity within CL, we assembled a multi-ancestry cohort of 837 nonsyndromic CL case-parent trios with whole-genome sequencing and detailed phenotyping. We performed genome-wide association scans (GWAS) via transmission disequilibrium tests for CL overall and for 14 CL subtypes defined by involvement of the alveolus (with and without), laterality (uni- and bilateral), and sidedness (left and right). We identified four genome-wide significant loci. Two loci, IRF6 and 8q24.21, were both detected in the overall CL GWAS. PLCB1/PLCB4 and MAFB were detected in GWASs of alveolar cleft involvement and CL left sidedness, respectively. These subtype-specific associations were followed by case-only comparisons that reflect the presence or absence of alveolus cleft or left-sided bias of CL to confirm the specificity of the association signal to the particular subtype. Our results provide evidence of within-class CL subtype-specific genetic links for loci previously discussed in the context of primary OFC classes and demonstrate the value of granular OFC subtype characterization to capture trait-specific associations.

Alveolus Cleft↗

Comprehensive analysis of de novo variants across 2,497 orofacial cleft trios reveals novel genetic drivers of disease.

BACKGROUND: Orofacial clefts (OFCs) and other palate abnormalities (PAs) are among the most common birth defects worldwide and are characterized by the abnormal formation of the lip and/or palate. Genetic studies have traditionally classified OFC cases as either syndromic, involving OFCs alongside other congenital anomalies, or nonsyndromic, which represent the majority of cases and occur in isolation. Emerging genomic evidence indicates that genes traditionally associated with syndromic forms of OFC can also harbor variants contributing to isolated cases, challenging the notion of a strict dichotomy between these categories and supporting their integration for gene discovery. METHODS: In this study, we applied multiple analytic approaches to characterize the genetic architecture of OFC and PAs by integrating genomic data from 2,497 trios with probands diagnosed with an OFC (n=2,080) or PA (n=417). We compared these findings across OFC subtypes and syndromic status with those from 5,515 control trios to identify enriched biological pathways and mechanisms and to prioritize candidate genes using variant burden testing. RESULTS: We observed a significant enrichment of de novo protein-truncating and damaging missense variants in cases compared to controls (OR = 2.17, p = 1.21×10-32), with particularly strong signals in biologically relevant gene sets involving OFC-associated, constrained, Mendelian disorder, and mouse candidate genes. Variant burden testing identified 39 OFC risk genes at FDR ≤ 0.05, which we then integrated with 593 established OFC genes to interrogate the functional underpinnings of OFC via network analysis. This analysis revealed 309 high-order interactor genes not previously associated with OFC. Notably, this OFC network clustered into ten distinct biological pathways, with nucleosome-associated genes showing significant enrichment among cases in our cohort (OR = 14.8, p = 8.1×10-4). In a final integrative step, we combined evidence across all analyses to nominate 231 candidate genes, 32 of which contained at least two deleterious de novo variants in our cohort. CONCLUSIONS: These findings underscore the value of integrating diverse OFC and PA subtypes, syndromic status, and variant classes to elucidate the genetic architecture of these disorders, highlighting both phenotypic expansion of known disease genes and the emergence of novel gene-phenotype associations.

De Novo Variant Enrichment↗

Fundamentals of FAIR biomedical data analyses in the cloud using custom pipelines.

As the biomedical data ecosystem increasingly embraces the findable, accessible, interoperable, and reusable (FAIR) data principles to publish multimodal datasets to the cloud, opportunities for cloud-based research continue to expand. Besides the potential for accelerated and diverse biomedical discovery that comes from a harmonized data ecosystem, the cloud also presents a shift away from the standard practice of duplicating data to computational clusters or local computers for analysis. However, despite these benefits, researcher migration to the cloud has lagged, in part due to insufficient educational resources to train biomedical scientists on cloud infrastructure. There exists a conceptual lack especially around the crafting of custom analytic pipelines that require software not pre-installed by cloud analysis platforms. We here present three fundamental concepts necessary for custom pipeline creation in the cloud. These overarching concepts are workflow and cloud provider agnostic, extending the utility of this education to serve as a foundation for any computational analysis running any dataset in any biomedical cloud platform. We illustrate these concepts using one of our own custom analyses, a study using the case-parent trio design to detect sex-specific genetic effects on orofacial cleft (OFC) risk, which we crafted in the biomedical cloud analysis platform CAVATICA.

Cloud Computing↗

Clustering individuals using INMTD: a novel versatile multi-view embedding framework integrating omics and imaging data.

MOTIVATION: Combining omics and images can lead to a more comprehensive clustering of individuals than classic single-view approaches. Among the various approaches for multi-view clustering, nonnegative matrix tri-factorization (NMTF) and nonnegative Tucker decomposition (NTD) are advantageous in learning low-rank embeddings with promising interpretability. Besides, there is a need to handle unwanted drivers of clusterings (i.e. confounders). RESULTS: In this work, we introduce a novel multi-view clustering method based on NMTF and NTD, named INMTD, which integrates omics and 3D imaging data to derive unconfounded subgroups of individuals. According to the adjusted Rand index, INMTD outperformed other clustering methods on a synthetic dataset with known clusters. In the application to real-life facial-genomic data, INMTD generated biologically relevant embeddings for individuals, genetics, and facial morphology. By removing confounded embedding vectors, we derived an unconfounded clustering with better internal and external quality; the genetic and facial annotations of each derived subgroup highlighted distinctive characteristics. In conclusion, INMTD can effectively integrate omics data and 3D images for unconfounded clustering with biologically meaningful interpretation. AVAILABILITY AND IMPLEMENTATION: INMTD is freely available at https://github.com/ZuqiLi/INMTD.

Cluster Analysis↗

Optimized phenotyping of complex morphological traits: enhancing discovery of common and rare genetic variants.

Genotype-phenotype (G-P) analyses for complex morphological traits typically utilize simple, predetermined anatomical measures or features derived via unsupervised dimension reduction techniques (e.g. principal component analysis (PCA) or eigen-shapes). Despite the popularity of these approaches, they do not necessarily reveal axes of phenotypic variation that are genetically relevant. Therefore, we introduce a framework to optimize phenotyping for G-P analyses, such as genome-wide association studies (GWAS) of common variants or rare variant association studies (RVAS) of rare variants. Our strategy is two-fold: (i) we construct a multidimensional feature space spanning a wide range of phenotypic variation, and (ii) within this feature space, we use an optimization algorithm to search for directions or feature combinations that are genetically enriched. To test our approach, we examine human facial shape in the context of GWAS and RVAS. In GWAS, we optimize for phenotypes exhibiting high heritability, estimated from either family data or genomic relatedness measured in unrelated individuals. In RVAS, we optimize for the skewness of phenotype distributions, aiming to detect commingled distributions that suggest single or few genomic loci with major effects. We compare our approach with eigen-shapes as baseline in GWAS involving 8246 individuals of European ancestry and in gene-based tests of rare variants with a subset of 1906 individuals. After applying linkage disequilibrium score regression to our GWAS results, heritability-enriched phenotypes yielded the highest SNP heritability, followed by eigen-shapes, while commingling-based traits displayed the lowest SNP heritability. Heritability-enriched phenotypes also exhibited higher discovery rates, identifying the same number of independent genomic loci as eigen-shapes with a smaller effective number of traits. For RVAS, commingling-based traits resulted in more genes passing the exome-wide significance threshold than eigen-shapes, while heritability-enriched phenotypes lead to only a few associations. Overall, our results demonstrate that optimized phenotyping allows for the extraction of genetically relevant traits that can specifically enhance discovery efforts of common and rare variants, as evidenced by their increased power in facial GWAS and RVAS.

Humans↗

Minor physical anomalies in schizophrenia: a meta-analysis.

Numerous studies report an increased frequency of minor physical anomalies (MPAs) in schizophrenic individuals compared with controls. However, these studies vary considerably regarding the magnitude of the case-control disparity and the topographical distribution of the anomalies. A meta-analysis was carried out on the existing MPA literature in an effort to better understand the relationship between MPAs and schizophrenia. Following a literature search, 13 studies were identified that met our inclusion criteria. Mean total MPA scores were available for 11 of these studies, whereas only seven studies provided regional MPA scores. For both the total MPA and regional MPA analyses, pooled effect sizes (Hedges' g and pooled odds ratios, respectively) were calculated along with tests of heterogeneity. For the total MPA analyses, a meta-regression approach was used to explore the relationship between possible moderator variables (e.g., number of MPA scale items) and effect size heterogeneity. The magnitude of the pooled effect size for the total MPA scores was high (1.131; p<0.001), indicating significantly more overall MPAs in schizophrenic individuals. Significant effect size heterogeneity was present (p<0.001); however, this heterogeneity could not be explained by any of the included moderator variables. The regional MPA analysis revealed significantly increased MPAs in all six anatomical regions (p<0.05), although the pooled odds ratios for these regions did not differ significantly from one another. These results suggest a lack of regional specificity for MPAs in schizophrenia.

Adult↗

Nicotine metabolizing genes GSTT1 and CYP1A1 in sudden infant death syndrome.

Exposure to tobacco, both to the developing fetus as well as in the postnatal period, has been identified as a key risk factor in the etiology of sudden infant death syndrome (SIDS). Polymorphisms in both the GSTT1 and CYP1A1 genes have been reported to impact the metabolic detoxification process for cigarette smoke and have been associated with low birth weight. Thus, expression of polymorphisms in these genes may account for the varying susceptibility to the adverse health consequences of tobacco exposure, including SIDS. We hypothesized that functional polymorphisms in GSTT1 (gene deletion) and CYP1A1 (m1, m2, and m3) might be associated with SIDS risk. DNA was prepared from 106 SIDS cases and 106 ethnicity- and gender-matched controls using standard methods. Regions of interest were amplified using PCR, subjected to enzyme digestion, and analyzed on agarose gel. No association was observed between the GSTT1 gene deletion or the CYP1A1 m1, m2, and m3 polymorphisms with SIDS risk when considered independently or in combination. These results indicate that the GSTT1 gene deletion and polymorphisms of CYP1A1 are not responsible for increased SIDS risk in our dataset. However, because SIDS cases with confirmed history of nicotine exposure were limited (7/106 cases), a relationship that might be apparent in a cohort with a large subset of SIDS cases with known history of nicotine exposure cannot be ruled out. A prospective study of SIDS cases with nicotine exposure history is necessary to resolve the relationship between nicotine metabolizing genes and SIDS.

Case-Control Studies↗

Anthropometric precision and accuracy of digital three-dimensional photogrammetry: comparing the Genex and 3dMD imaging systems with one another and with direct anthropometry.

A variety of commercially available three-dimensional (3D) surface imaging systems are currently in use by craniofacial specialists. Little is known, however, about how measurement data generated from alternative 3D systems compare, specifically in terms of accuracy and precision. The purpose of this study was to compare anthropometric measurements obtained by way of two different digital 3D photogrammetry systems (Genex and 3dMD) as well as direct anthropometry and to evaluate intraobserver precision across these three methods. On a sample of 18 mannequin heads, 12 linear distances were measured twice by each method. A two-factor repeated measures analysis of variance was used to test simultaneously for mean differences in precision across methods. Additional descriptive statistics (e.g., technical error of measurement [TEM]) were used to quantify measurement error magnitude. Statistically significant (P < 0.05) mean differences were observed across methods for nine anthropometric variables; however, the magnitude of these differences was consistently at the submillimeter level. No significant differences were noted for precision. Moreover, the magnitude of imprecision was determined to be very small, with TEM scores well under 1 mm, and intraclass correlation coefficients ranging from 0.98 to 1. Results indicate that overall mean differences across these three methods were small enough to be of little practical importance. In terms of intraobserver precision, all methods fared equally well. This study is the first attempt to simultaneously compare 3D surface imaging systems directly with one another and with traditional anthropometry. Results suggest that craniofacial surface data obtained by way of alternative 3D photogrammetric systems can be combined or compared statistically.

Anthropometry↗

Candidate genes for oral-facial clefts in Guatemalan families.

Nonsyndromic cleft lip +/- cleft palate (CL/P) is a complex trait of unknown etiology. Most genetic studies of CL/P define affection status in a way that ignores subtle subclinical manifestations, resulting in a potential loss of statistical power. This study investigated 10 candidate genes in 155 individuals from 25 Guatemalan CL/P families. High-resolution ultrasound images of the orbicularis oris (OO) muscle were obtained. CL/P was present in 28 family members; an additional 10 had subcutaneous OO muscle defects. Family-based association studies were performed for both narrow (CL/P only) and broad (CL/P plus OO muscle defects) definitions of affection status. PVRL1 was significantly associated under both definitions (P = 0.04, narrow; P = 0.02, broad). Association with JAG2 improved from P = 0.09 under the narrow definition to P = 0.04 under the broad definition. Broadening the oral-facial cleft phenotype to include subclinical variants may improve power in genetic studies.

Chromosome Aberrations↗

Characterization of dermatoglyphics in PHOX2B-confirmed congenital central hypoventilation syndrome.

OBJECTIVE: Individuals with congenital central hypoventilation syndrome have characteristic variants in the PHOX2B gene (primarily polyalanine expansion mutations). The PHOX2B gene acts as a transcriptional activator in the promotion of pan-neuronal differentiation in the autonomic nervous system during early embryologic development, with a primary role in the sympathetic noradrenergic phenotype in vertebrates. Because sympathetic innervation has been hypothesized to affect the development of dermatoglyphic pattern types, we hypothesized that individuals with PHOX2B-confirmed congenital central hypoventilation syndrome would have characteristic dermatoglyphic patterning and that the dermatoglyphic phenotype would be related to the disease-defining PHOX2B genotype. METHODS: Dermatoglyphic pattern type frequency, left/right symmetry, and genotype/phenotype correlation were assessed for 33 individuals with PHOX2B-confirmed congenital central hypoventilation syndrome and compared with published control data. RESULTS: Dermatoglyphic pattern type frequencies were altered in congenital central hypoventilation syndrome cases versus controls. In particular, there was an increase of arches in females and ulnar loops in males, with the largest differences for the left hand and for individuals with both congenital central hypoventilation syndrome and Hirschsprung disease. Dissimilarity scores between the congenital central hypoventilation syndrome and congenital central hypoventilation syndrome + Hirschsprung disease cases were not significantly different, nor were dissimilarity scores between all of the female and all of the male cases. No significant association was found between the number of polyalanine repeats in the PHOX2B genotypic category and dermatoglyphic pattern frequencies in the congenital central hypoventilation syndrome study groups. CONCLUSIONS: These results represent the first report describing specific dermatoglyphic patterning in congenital central hypoventilation syndrome and suggest a relationship between PHOX2B and the expression of dermatoglyphic pattern types. An expanded congenital central hypoventilation syndrome data set to include the full spectrum of PHOX2B mutations is necessary to further delineate the role of PHOX2B in dermatoglyphic patterning.

Adolescent↗

The Pittsburgh Oral-Facial Cleft study: expanding the cleft phenotype. Background and justification.

The Pittsburgh Oral-Facial Cleft study was begun in 1993 with the primary goal of identifying genes involved in nonsyndromic orofacial clefts in a variety of populations worldwide. Based on the results from a number of pilot studies and preliminary genetic analyses, a new research focus was added to the Pittsburgh Oral-Facial Cleft study in 1999: to elucidate the role that associated phenotypic features play in the familial transmission patterns of orofacial clefts in order to expand the definition of the nonsyndromic cleft phenotype. The purpose of this paper is to provide a comprehensive review of phenotypic features associated with nonsyndromic orofacial clefts. These features include fluctuating and directional asymmetry, non-right-handedness, dermatoglyphic patterns, craniofacial morphology, orbicularis oris muscle defects, dental anomalies, structural brain and vertebral anomalies, minor physical anomalies, and velopharyngeal incompetence.

Brain↗

Asian oral-facial cleft birth prevalence.

OBJECTIVE: To determine the clefting birth prevalence among Asian populations, specifically Chinese and Japanese, using raw counts from nonoverlapping published studies of Asian populations, and to investigate whether Asian clefting rates have been interpreted accurately as being up to twice the Caucasian rate. DESIGN: A literature review of articles giving raw counts of clefting in Asian populations, primarily Japanese and Chinese. MAIN OUTCOME MEASURES: Where possible, clefts were identified by the patients' ethnicity, country of origin, cleft type, syndromic status, and birth status. RESULTS: Prevalence rates of cleft lip with or without cleft palate per 1000 live births are reported. Syndromic plus nonsyndromic cleft lip with or without cleft palate: Chinese, 1.30; Japanese, 1.34; Other Asian, 1.47; and total, 1.33. Nonsyndromic cleft lip with or without cleft palate: Chinese, 1.20; Japanese, 1.18; Other Asian, 1.22; and total, 1.19. CONCLUSIONS: Overall, Chinese and Japanese live birth prevalence rates for nonsyndromic cleft lip with or without cleft palate, based on the published reports of birth prevalence, are significantly lower than the oft-quoted rate of 2 per 1000 for Asians. The apparent reason for the discrepancy is that many published prevalence rates included all pregnancies (live births plus pregnancy losses) and do not distinguish between syndromic and nonsyndromic clefts or between cleft palate alone and cleft lip with or without cleft palate. These results demonstrate that it is extremely important for current population-based studies of clefts to include careful delineation of population groups, syndromes, cleft type, and birth status.

Asian People↗

Methods for detecting gene x gene interaction in multiplex extended pedigrees.

Complex diseases are multifactorial in nature and can involve multiple loci with gene x gene and gene x environment interactions. Research on methods to uncover the interactions between those genes that confer susceptibility to disease has been extensive, but many of these methods have only been developed for sibling pairs or sibships. In this report, we assess the performance of two methods for finding gene x gene interactions that are applicable to arbitrarily sized pedigrees, one based on correlation in per-family nonparametric linkage scores and another that incorporates candidate loci genotypes as covariates into an affected relative pair linkage analysis. The power and type I error rate of both of these methods was addressed using the simulated Genetic Analysis Workshop 14 data. In general, we found detection of the interacting loci to be a difficult problem, and though we experienced some modest success there is a clear need to continue developing new methods and approaches to the problem.

Epistasis, Genetic↗

Identifying genomic regions for fine-mapping using genome scan meta-analysis (GSMA) to identify the minimum regions of maximum significance (MRMS) across populations.

In order to detect linkage of the simulated complex disease Kofendrerd Personality Disorder across studies from multiple populations, we performed a genome scan meta-analysis (GSMA). Using the 7-cM microsatellite map, nonparametric multipoint linkage analyses were performed separately on each of the four simulated populations independently to determine p-values. The genome of each population was divided into 20-cM bin regions, and each bin was rank-ordered based on the most significant linkage p-value for that population in that region. The bin ranks were then averaged across all four studies to determine the most significant 20-cM regions over all studies. Statistical significance of the averaged bin ranks was determined from a normal distribution of randomly assigned rank averages. To narrow the region of interest for fine-mapping, the meta-analysis was repeated two additional times, with each of the 20-cM bins offset by 7 cM and 13 cM, respectively, creating regions of overlap with the original method. The 6-7 cM shared regions, where the highest averaged 20-cM bins from each of the three offsets overlap, designated the minimum region of maximum significance (MRMS). Application of the GSMA-MRMS method revealed genome wide significance (p-values refer to the average rank assigned to the bin) at regions including or adjacent to all of the simulated disease loci: chromosome 1 (p < 0.0001 for 160-167 cM, including D1), chromosome 3 (p-value < 0.0000001 for 287-294 cM, including D2), chromosome 5 (p-value < 0.001 for 0-7 cM, including D3), and chromosome 9 (p-value < 0.05 for 7-14 cM, the region adjacent to D4). This GSMA analysis approach demonstrates the power of linkage meta-analysis to detect multiple genes simultaneously for a complex disorder. The MRMS method enhances this powerful tool to focus on more localized regions of linkage.

Chromosomes, Human↗

Facial phenotype in children and young adults with PHOX2B-determined congenital central hypoventilation syndrome: quantitative pattern of dysmorphology.

Congenital central hypoventilation syndrome (CCHS) is caused by mutations in PHOX2B, which is essential for maturation of the neural crest into the autonomic nervous system and is expressed in the dorsal rhombencephalon, a region that gives rise to facial structures. Digital photographs of 45 individuals with PHOX2B-confirmed CCHS, and 45 matched controls were analyzed for 17 linear and 6 angular measurements, and 9 derived indices. Paired t tests were used to compare group means, correlation was calculated between PHOX2B polyalanine expansion number and facial measures, and stepwise logistic regression was used to predict case-control and genotype status. CCHS cases differed significantly from controls on 13 variables (6 after p value correction: nasolabial angle, upper lip height, lateral lip height, facial index, upper facial index, and presence of inferior inflection of the lateral segment of the upper lip vermillion border). Five variables were able to predict correctly 85.7% of CCHS cases and 82.2% of controls: upper lip height, biocular width, upper facial height, nasal tip protrusion, and inferior inflection of the upper lip vermillion border. A negative relationship between number of repeats and four anthropometric measures was observed: mandible breadth, nasolabial angle, lateral lip height, and mandible-face width index. These results suggest a characteristic facial phenotype in children and young adults with CCHS, due to an expansion mutation in PHOX2B.

Adolescent↗

Medical sequencing of candidate genes for nonsyndromic cleft lip and palate.

Nonsyndromic or isolated cleft lip with or without cleft palate (CL/P) occurs in wide geographic distribution with an average birth prevalence of 1/700. We used direct sequencing as an approach to study candidate genes for CL/P. We report here the results of sequencing on 20 candidate genes for clefts in 184 cases with CL/P selected with an emphasis on severity and positive family history. Genes were selected based on expression patterns, animal models, and/or role in known human clefting syndromes. For seven genes with identified coding mutations that are potentially etiologic, we performed linkage disequilibrium studies as well in 501 family triads (affected child/mother/father). The recently reported MSX1 P147Q mutation was also studied in an additional 1,098 cleft cases. Selected missense mutations were screened in 1,064 controls from unrelated individuals on the Centre d'Etude du Polymorphisme Humain (CEPH) diversity cell line panel. Our aggregate data suggest that point mutations in these candidate genes are likely to contribute to 6% of isolated clefts, particularly those with more severe phenotypes (bilateral cleft of the lip with cleft palate). Additional cases, possibly due to microdeletions or isodisomy, were also detected and may contribute to clefts as well. Sequence analysis alone suggests that point mutations in FOXE1, GLI2, JAG2, LHX8, MSX1, MSX2, SATB2, SKI, SPRY2, and TBX10 may be rare causes of isolated cleft lip with or without cleft palate, and the linkage disequilibrium data support a larger, as yet unspecified, role for variants in or near MSX2, JAG2, and SKI. This study also illustrates the need to test large numbers of controls to distinguish rare polymorphic variants and prioritize functional studies for rare point mutations.

Cleft Lip↗