Transmission Disequilibrium Test Research Articles

An Expression Quantitative Trait Locus of Fc Gamma Receptor Genes is Associated with Anti-Malarial IgG Responses and Infection Levels in Burkinabe Families.

The interaction between antibodies and Fc gamma receptors (FcγRs) plays a critical role in regulating immune responses to Plasmodium falciparum. Polymorphisms in genes encoding FcγRs influence the host's capacity to control parasite infection. This study investigates whether non-coding variants influencing FcγR expression are associated with anti-malarial immunization and infection traits. We utilized eQTL databases and functional annotations to identify non-coding variants, specifically rs1771575, rs2099684, and rs6700241, within the FCGR gene cluster. In addition, we examined the coding variants rs1801274 (p.His167Arg) and rs1050501 (p.Ile231Thr), which affect the affinity of FcγRIIa and FcγRIIb for IgG. These variants were genotyped in 163 individuals from Burkinabe families. Family-based linear mixed regression and Quantitative Transmission Disequilibrium Tests (QTDT) analyses were performed to assess associations with IgG levels and malaria infection, accounting for relevant covariates. Linear mixed models identified rs1771575 as associated with total IgG levels, while both rs1771575 and rs1801274 were linked to IgG2, and rs1050501 to IgG1 levels. A haplotype combining rs2099684 and rs6700241 was positively associated with IgG1. The rs1771575-CC and rs1050501-TT genotypes correlated with higher infection levels in children. QTDT models confirmed the association of rs1771575 with IgG2 and infection in children. Our findings suggest that the intergenic variant rs1771575 serves as an independent marker for IgG levels and blood infection in children. This highlights the interplay between regulatory variants and coding mutations in FCGR, which may influence immune function and antibody production. These results underscore the potential for personalized strategies to monitor humoral responses in malaria-endemic regions.

Genetic Variants in METTL16 Affect the Risk of Non-Syndromic Orofacial Clefts.

N6-methyladenosine (m6A) is the most prevalent modification of RNA in eukaryotes which is associated with many cellular processes and diseases. Here, our objective is to explore whether genetic variants in m6A modification genes are associated with the risk of non-syndrome orofacial clefts (NSOCs). The transmission disequilibrium test (TDT) was performed to calculate the association between single nucleotide polymorphisms (SNPs) in m6A modification genes and NSOCs risk in 944 case-parent trios. The function of SNP was predicted by HaploReg, RegulomeDB and histone enrichment data. The expression quantitative trait locus (eQTL) analysis was examined using Genotype-Tissue Expression (GTEx) and eQTLGen. The role of gene in the development of NSOCs was assessed with correlation and enrichment analysis based on gene expression data in mice craniofacial tissue and zebrafish embryo. We identified that rs8078195 (A > C) in METTL16 was suggestively associated with the increased risk of NSOCs (OR = 1.32, p = 1.80E - 03). The region surrounding rs8078195 was subjected to deoxyribonuclease hypersensitivity and enriched with multiple histone modifications. In addition, it had a significant eQTL effect with METTL16 in skin tissue and human peripheral blood, which played an important role in NSOCs development. Bioinformatic analysis indicated that METTL16 contributed to the development of NSOCs probably by regulating cell cycle process. Rs8078195 in METTL16 was associated with the occurrence of NSOCs.

Statistics to prioritize rare variants in family-based sequencing studies with disease subtypes.

Family-based sequencing studies are increasingly used to find rare genetic variants of high risk for disease traits with familial clustering. In some studies, families with multiple disease subtypes are collected and the exomes of affected relatives are sequenced for shared rare variants (RVs). Since different families can harbor different causal variants and each family harbors many RVs, tests to detect causal variants can have low power in this study design. Our goal is rather to prioritize shared variants for further investigation by, for example, pathway analyses or functional studies. The transmission-disequilibrium test prioritizes variants based on departures from Mendelian transmission in parent-child trios. Extending this idea to families, we propose methods to prioritize RVs shared in affected relatives with two disease subtypes, with one subtype more heritable than the other. Global approaches condition on a variant being observed in the study and assume a known probability of carrying a causal variant. In contrast, local approaches condition on a variant being observed in specific families to eliminate the carrier probability. Our simulation results indicate that global approaches are robust to misspecification of the carrier probability and prioritize more effectively than local approaches even when the carrier probability is misspecified.

Genetic associations and parent-of-origin effects of PVRL1 in non-syndromic cleft lip with or without cleft palate across multiple ethnic populations.

This study investigated the associations of PVRL1 gene variants with non-syndromic cleft lip with or without cleft palate (NSCL/P) by evaluating transmission distortion and parent-of-origin (POO) effects in multiple ethnic populations. We conducted allelic and genotypic transmission disequilibrium tests (TDT) on 10 single-nucleotide variants (SNVs) in PVRL1 using data from 142 Korean families with an affected child. POO effects were analyzed using the POO likelihood ratio test, comparing transmission rates of maternally and paternally inherited alleles. To assess generalizability and ethnic heterogeneity, we compared results from Korean families with data from the Center for Craniofacial and Dental Genetics, which included 2,226 individuals from 497 European and 245 Asian trios. TDT analysis identified significant over-transmission of the rs7940667 (G361V) C allele in Korean families (p=0.007), a finding replicated in both Asian (p=6.5×10-7) and European families (p=1.6×10-10). Eight SNVs showed strong TDT evidence in larger Asian and European datasets after multiple comparison corrections (p<0.0073). Of these, 4 SNVs (rs7940667, rs7103685, rs7129848, and rs4409845) showed particularly robust association (p<5×10-8). POO analysis revealed significant maternal over-transmission of the rs10790330-A allele in Korean families (p=0.044). This finding was replicated in European families (p=9.0×10-4). Additionally, 3 other SNVs, rs7129848 (p=0.001) and the linked SNVs rs3935406 and rs10892434 (p=0.025), exhibited maternal over-transmission in the validation datasets. Our findings provide robust evidence supporting the associations of PVRL1 variants with NSCL/P susceptibility. Further research is necessary to explore the potential clinical applications of these findings.

Combined maternal KIR2DL4 and fetal HLA-G polymorphisms were associated with preeclampsia in a Han Chinese population.

Preeclampsia is the main cause of maternal and infant mortality and morbidity during pregnancy. Killer cell immunoglobulin-like receptor 2DL4 (KIR2DL4) and human leukocyte antigen G (HLA-G) play crucial roles in immune tolerance at the maternal-fetal interface. In this case‒control study, 154 maternal-fetal pairs were recruited, including 74 pairs with preeclampsia (56 of 74 pairs from family triads) and 80 pairs with a normal pregnancy (78 of 80 pairs from family triads). SNaPshot technology was used to detect genetic polymorphisms for 7 TagSNPs in the KIR2DL4 and HLA-G genes. Among the fetal HLA-G gene polymorphisms, rs9380142 (A vs. G: OR = 2.802, 95% CI = 1.761-4.458) and rs1063320 (G vs. C: OR = 1.807, 95% CI = 1.144-2.852) differed between the preeclampsia group and the control group. The transmission disequilibrium test (TDT) suggested that the differences in the rs9380142G/A polymorphism in foetuses between preeclampsia triads and control triads were due to differences in transmission from the parents (P = 0.001). There was no significant difference in the distribution of maternal KIR2DL4 alleles or genotype frequency between the preeclampsia group and the control group. Gene‒gene interaction analysis revealed that the combined genotypes of maternal rs649216-CC and fetal rs9380142-GG, maternal rs1051456-CG/GG and fetal rs9380142-GG, maternal rs34785252-CC and fetal rs9380142-AA/GA, and maternal rs34785252-CC/AA and fetal rs9380142-GG were associated with a significantly lower risk of preeclampsia. Therefore, this study suggested that the combination of maternal KIR2DL4 and fetal HLA-G polymorphisms was associated with preeclampsia in a Han Chinese population.

Whole genome sequencing of a family with autosomal dominant features within the oculoauriculovertebral spectrum.

Oculoauriculovertebral Spectrum (OAVS) encompasses abnormalities on derivatives from the first and second pharyngeal arches including macrostomia, hemifacial microsomia, micrognathia, preauricular tags, ocular and vertebral anomalies. We present genetic findings on a three-generation family affected with macrostomia, preauricular tags and uni- or bilateral ptosis following an autosomal dominant pattern. We generated whole genome sequencing data for the proband, affected parent and unaffected paternal grandparent followed by Sanger sequencing on 23 family members for the top 10 candidate genes: KCND2, PDGFRA, CASP9, NCOA3, WNT10A, SIX1, MTF1, KDR/VEGFR2, LRRK1, and TRIM2 We performed parent and sibling-based transmission disequilibrium tests and burden analysis via a penalized linear mixed model, for segregation and mutation burden respectively. Next, via bioinformatic tools we predicted protein function, mutation pathogenicity and pathway enrichment to investigate the biological relevance of mutations identified. Rare missense mutations in SIX1, KDR/VEGFR2, and PDGFRA showed the best segregation with the OAV phenotypes in this family. When considering any of the 3 OAVS phenotypes as an outcome, SIX1 had the strongest associations in parent-TDTs and sib-TDTs (p=0.025, p=0.052) (unadjusted p-values). Burden analysis identified SIX1 (RC=0.87) and PDGFRA (RC=0.98) strongly associated with OAVS severity. Using phenotype-specific outcomes, sib-TDTs identified SIX1 with uni- or bilateral ptosis (p=0.049) and ear tags (p=0.01), and PDGFRA and KDR/VEGFR2 with ear tags (both p<0.01). SIX1, PDGFRA, and KDR/VEGFR2 are strongly associated to OAVS phenotypes. SIX1 has been previously associated with OAVS ear malformations and is co-expressed with EYA1 during ear development. Efforts to strengthen the genotype-phenotype co-relation underlying the OAVS are key to discover etiology, family counseling and prevention.

Family-based GWAS for dental class I malocclusion and clefts

BackgroundIndividuals born with cleft lip and/or palate who receive corrective surgery regularly have abnormal growth in the midface region such that they exhibit premaxillary hypoplasia. However, there are also genetic contributions to craniofacial morphology in the midface region, so although these individuals appear to have Class III skeletal discrepancy, their molar relationship may be Class I. Past genome-wide association studies (GWASs) on skeletal Class II and III malocclusion suggested that multiple genetic markers contribute to these phenotypes via a multifactorial inheritance model, but research has yet to examine the genetic markers associated with dental Class I malocclusion. Thus, our goal was to conduct a family based GWAS to identify genes across the genome that are associated with Class I malocclusion, as defined by molar relations, in humans with and without clefts.MethodsOur cohort consisted of 739 individuals from 47 Filipino families originally recruited in 2006 to investigate the genetic basis of orofacial clefts. All individuals supplied blood samples for DNA extraction and genotyping, and a 5,766 single nucleotide polymorphism (SNP) custom panel was used for the analyses. We performed a transmission disequilibrium test for participants with and without clefts to identify genetic contributors potentially involved with Class I malocclusion.ResultsIn the total cohort, 13 SNPs had associations that reached the genomic control threshold (p < 0.005), while five SNPs were associated with Class I in the cohort of participants without clefts, including four associations that were identified in the total cohort. The associations for the SNPs ABCA4 rs952499, SOX1-OT rs726455, and RORA rs877228 are of particular interest, as past research found associations between these genes and various craniofacial phenotypes, including cleft lip and/or palate.ConclusionsThese findings support the multifactorial inheritance model for dental Class I malocclusion and suggest a common genetic basis for different aspects of craniofacial development.

Attention-deficit/hyperactivity disorder and dopamine receptor D4 (DRD4) exon 3 variable number of tandem repeats (VNTR) 2-repeat allele.

To investigate the association of attention-deficit/hyperactivity disorder (ADHD) with the 48-base pair (bp) variable number of tandem repeats (VNTR) in exon 3 of the dopamine receptor D4 (DRD4) gene, we genotyped 240 ADHD patients and their parents from Hong Kong. The 4R allele was most common, followed by 2R. We examined association between the 2R allele (relative to 4R) and ADHD by Transmission Disequilibrium Test (TDT). The odds ratio (OR) (95% confidence interval) was 0.90 (0.64-1.3). The p-value was 0.6. Examining subgroups revealed nominally significant association of 2R with inattentive ADHD: OR=0.33 (0.12-0.92) and p=0.03. Because our study used TDT analysis, we meta-analyzed the association of 2R with ADHD in Asians (1329 patient alleles), revealing results similar to ours: OR=0.97 (0.80-1.2) and p=0.8. To examine the association of 2R with inattentive ADHD, we meta-analyzed all studies (regardless of analysis type or ethnicity, in order to increase statistical power): 702 patient alleles, 1420 control alleles, OR=0.81 (0.57-1.1) and p=0.2. Overall, there is no evidence of association between ADHD and the 2R allele, but the suggestive association with the inattentive type warrants further investigation.

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Family-based association of HLA-DRB1 and DQB1 alleles and haplotypes in a group of Iranian Type 1 diabetes children.

This family-based study was conducted in a group of Iranians with Type 1 diabetes (T1D) to investigate the transmission from parents of risk and non-risk HLA alleles and haplotypes, and to estimate the genetic risk score for this disease within this population. A total of 240 T1D subjects including 111 parent-child trios (111 children with T1D, 133 siblings, and 222 parents) and 330 ethnically matched healthy individuals were recruited. High-resolution HLA typing for DRB1/DQB1 loci was performed for all study subjects (n = 925) using polymerase chain reaction-sequence-specific oligonucleotide probe method. The highest predisposing effect on developing T1D was conferred by the following haplotypes both in all subjects and in probands compared to controls: DRB1*04:05-DQB1*03:02 (Pc = 2.97e-06 and Pc = 6.04e-10, respectively), DRB1*04:02-DQB1*03:02 (Pc = 5.94e-17 and Pc = 3.86e-09, respectively), and DRB1*03:01-DQB1*02:01 (Pc = 8.26e-29 and Pc = 6.56e-16, respectively). Conversely, the major protective haplotypes included DRB1*13:01-DQB1*06:03 (Pc = 6.99e-08), DRB1*15:01-DQB1*06:02 (Pc = 2.97e-06) in the cases versus controls. Also, DRB1*03:01-DQB1*02:01/DRB1*04:02|05-DQB1*03:02 and DRB1*03:01-DQB1*02:01/DRB1*03:01-DQB1*02:01 diplotypes conferred the highest predisposing effect in the cases (Pc = 8.65e-17 and Pc = 6.26e-08, respectively) and in probands (Pc = 5.4e-15 and Pc = 0.001, respectively) compared to controls. Transmission disequilibrium test showed that the highest risk was conferred by DRB1*04:02-DQB1*03:02 (Pc = 3.26e-05) and DRB1*03:01-DQB1*02:01 (Pc = 1.78e-12) haplotypes and the highest protection by DRB1*14:01-DQB1*05:03 (Pc = 8.66e-05), DRB1*15:01-DQB1*06:02 (Pc = 0.002), and DRB1*11:01-DQB1*03:01 (Pc = 0.0003) haplotypes. Based on logistic regression analysis, carriage of risk haplotypes increased the risk of T1D development 24.5 times in the Iranian population (p = 5.61e-13). Also, receiver operating characteristic curve analysis revealed a high predictive power of those risk haplotypes in discrimination of susceptible from healthy individuals (area under curve: 0.88, p = 5.5e-32). Our study highlights the potential utility of genetic risk assessment based on HLA diplotypes for predicting T1D risk in individuals, particularly among family members of affected children in our population.

Association of Vitamin D Deficiency and Vitamin D Receptor Gene Polymorphisms with Type 1 Diabetes Risk: A South Indian Familial Study

Vitamin D is a potent immune modulator and is associated with autoimmune diseases, including type 1 diabetes (T1D). The vitamin D levels and its receptor gene polymorphisms together in T1D are not yet investigated in the South Indian population. The present study focused on exploring the significance of vitamin D levels and vitamin D receptor (VDR) gene polymorphisms with the risk of developing T1D in the South Indian population. Patients with T1D and unaffected first-degree relatives (FDRs) were included in this study. Genotyping of VDR polymorphisms at four different loci (FokI- F/f, BsmI- B/b, TaqI- T/t, and ApaI- A/a) was assessed through the amplification refractive mutation system-polymerase chain reaction method. Serum vitamin D levels were measured in 98 T1D patients and 75 age- and sex-matched siblings. A total of 120 patients with T1D and 214 FDRs were included. Vitamin D deficiency (VDD) was observed in a higher proportion of T1D patients than in controls (52% vs. 32%; p<0.03). The frequency of the FokI-FF genotype was significantly higher [odds ratio (OR)=1.66; p<0.03] in T1D patients conferring a susceptible association with the disease. Nevertheless, the increased frequency of heterozygous Ff genotype (OR=0.57; p<0.02) among controls may confer a protective association with T1D. Furthermore, the transmission disequilibrium test revealed over-transmission of ApaI-A (T: U=15/5; p<0.006) and BsmI-B alleles (T: U=17/5; p<0.01) and under-transmission of BsmI-b/ApaI-a/TaqI-T haplotype (T: U=5.4/14.4; p=0.04) from parents to T1D patients. The present study concludes that VDD is the major contributing risk factor to T1D development in the South Indian population. Furthermore, the FokI-FF genotype, BsmI-B, and ApaI-A alleles were positively associated with T1D. In contrast, the FokI-Ff genotype and BsmI-b/ApaI-a/TaqI-T haplotype were negatively associated with T1D.

Genome-wide analysis of spina bifida risk variants in a case-control study from Bangladesh.

Human studies of genetic risk factors for neural tube defects, severe birth defects associated with long-term health consequences in surviving children, have predominantly been restricted to a subset of candidate genes in specific biological pathways including folate metabolism. In this study, we investigated the association of genetic variants spanning the genome with risk of spina bifida (i.e., myelomeningocele and meningocele) in a subset of families enrolled from December 2016 through December 2022 in a case-control study in Bangladesh, a population often underrepresented in genetic studies. Saliva DNA samples were analyzed using the Illumina Global Screening Array. We performed genetic association analyses to compare allele frequencies between 112 case and 121 control children, 272 mothers, and 128 trios. In the transmission disequilibrium test analyses with trios only, we identified three novel exonic spina bifida risk loci, including rs140199800 (SULT1C2, p = 1.9 × 10-7), rs45580033 (ASB2, p = 4.2 × 10-10), and rs75426652 (LHPP, p = 7.2 × 10-14), after adjusting for multiple hypothesis testing. Association analyses comparing cases and controls, as well as models that included their mothers, did not identify genome-wide significant variants. This study identified three novel single nucleotide polymorphisms involved in biological pathways not previously associated with neural tube defects. The study warrants replication in larger groups to validate findings and to inform targeted prevention strategies.

Rare X-linked variants carry predominantly male risk in autism, Tourette syndrome, and ADHD

Autism spectrum disorder (ASD), Tourette syndrome (TS), and attention-deficit/hyperactivity disorder (ADHD) display strong male sex bias, due to a combination of genetic and biological factors, as well as selective ascertainment. While the hemizygous nature of chromosome X (Chr X) in males has long been postulated as a key point of “male vulnerability”, rare genetic variation on this chromosome has not been systematically characterized in large-scale whole exome sequencing studies of “idiopathic” ASD, TS, and ADHD. Here, we take advantage of informative recombinations in simplex ASD families to pinpoint risk-enriched regions on Chr X, within which rare maternally-inherited damaging variants carry substantial risk in males with ASD. We then apply a modified transmission disequilibrium test to 13,052 ASD probands and identify a novel high confidence ASD risk gene at exome-wide significance (MAGEC3). Finally, we observe that rare damaging variants within these risk regions carry similar effect sizes in males with TS or ADHD, further clarifying genetic mechanisms underlying male vulnerability in multiple neurodevelopmental disorders that can be exploited for systematic gene discovery.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate

Cleft palate (CP) is one of the most common craniofacial birth defects; however, there are relatively few established genetic risk factors associated with its occurrence despite high heritability. Historically, CP has been studied as a single phenotype, although it manifests across a spectrum of defects involving the hard and/or soft palate. We performed a genome-wide association study using transmission disequilibrium tests of 435 case-parent trios to evaluate broad risks for any cleft palate (ACP) (n= 435), and subtype-specific risks for any cleft soft palate (CSP), (n= 259) and any cleft hard palate (CHP) (n= 125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p= 4.24×10-8) associated with CHP. One gene at this locus, angiopoietin-like 2 (ANGPTL2), plays a role in osteoblast differentiation. It is expressed both in craniofacial tissue of human embryos and developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p<5×10-6), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios for the 20 loci were most similar across all 3 groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either subtype. We also found nominal evidence of replication (p<0.05) for 22 SNPs previously associated with orofacial clefts. Our study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

TGFA gene polymorphisms are not associated with the aetiology of non-syndromic cleft lip with or without cleft palate (NSCL/P) in south Indian population

IntroductionTransforming Growth Factor Alpha (TGFA) plays a role in the development of the facial structures, including the lip and palate. MethodologyAbout 148 case/parent trios were genotyped for rs1058213, rs2166975, rs3771494, rs3821261 and rs426081 polymorphisms in TGFA gene. Transmission/disequilibrium test (TDT) and the parent-of-origin effects (POE) statistics were implemented. FAMHAP and Haploview packages were used to test the effect of multi-SNP haplotypes and linkage disequilibrium respectively. ResultsAt individual SNP level, TDT did not reveal significant excess transmission of parental alleles to their affected child. None of the analysed markers showed significant POE. Significant LD between SNPs of rs1058213, rs2166975 and rs377149, rs3821261 were observed. Analysis multi –SNP haplotypes did not reveal significant excess transmission of haplotypes to the affected child. ConclusionAnalysis of TDT and POE of five SNPs in TGFA gene indicated that the TGFA is not associated with the aetiology of NSCL/P.

Transmission disequilibrium analysis of whole genome data in childhood-onset systemic lupus erythematosus.

Childhood-onset systemic lupus erythematosus (cSLE) patients are unique, with hallmarks of Mendelian disorders (early-onset and severe disease) and thus are an ideal population for genetic investigation of SLE. In this study, we use the transmission disequilibrium test (TDT), a family-based genetic association analysis that employs robust methodology, to analyze whole genome sequencing data. We aim to identify novel genetic associations in an ancestrally diverse, international cSLE cohort. Forty-two cSLE patients and 84 unaffected parents from 3 countries underwent whole genome sequencing. First, we performed TDT with single nucleotide variant (SNV)-based (common variants) using PLINK 1.9, and gene-based (rare variants) analyses using Efficient and Parallelizable Association Container Toolbox (EPACTS) and rare variant TDT (rvTDT), which applies multiple gene-based burden tests adapted for TDT, including the burden of rare variants test. Applying the GWAS standard threshold (5.0 × 10-8) to common variants, our SNV-based analysis did not return any genome-wide significant SNVs. The rare variant gene-based TDT analysis identified many novel genes significantly enriched in cSLE patients, including HNRNPUL2, a DNA repair protein, and DNAH11, a ciliary movement protein, among others. Our approach identifies several novel SLE susceptibility genes in an ancestrally diverse childhood-onset lupus cohort.

Behavioural deficits of autism spectrum disorder and associations with different gene clusters: a study with the whole-genome transmission disequilibrium test

BackgroundAutism spectrum disorder (ASD) is a diverse neurodevelopmental disease primarily distinguished by limited and stereotyped activities as well as impaired social interaction. Due to the high heritability of ASD, research...

Genetic variants in genes involved in creatine biosynthesis in patients with severe obesity or anorexia nervosa.

Increased thermogenesis in brown adipose tissue might have an obesity-reducing effect in humans. In transgenic mice, depletion of genes involved in creatine metabolism results in disrupted thermogenic capacity and altered effects of high-fat feeding on body weight. Data analyses of a sex-stratified genome-wide association study (GWAS) for body mass index (BMI) within the genomic regions of genes of this pathway (CKB, CKMT1B, and GATM) revealed one sex-dimorphic BMI-associated SNP in CKB (rs1136165). The effect size was larger in females than in males. A mutation screen of the coding regions of these three candidate genes in a screening group (192 children and adolescents with severe obesity, 192 female patients with anorexia nervosa, and 192 healthy-lean controls) identified five variants in each, CKB and GATM, and nine variants in the coding sequence of CKMT1B. Non-synonymous variants identified in CKB and CKMT1B were genotyped in an independent confirmation study group (781 families with severe obesity (trios), 320 children and adolescents with severe obesity, and 253 healthy-lean controls). In silico tools predicted mainly benign yet protein-destabilizing potentials. A transmission disequilibrium test in trios with severe obesity indicated an obesity-protective effect of the infrequent allele at rs149544188 located in CKMT1B. Subsequent correlation analyses in 1,479 individuals of the Leipzig Obesity BioBank revealed distinct correlations of CKB with the other two genes in omental visceral adipose tissue (VAT) and abdominal subcutaneous adipose tissue (SAT). Furthermore, between-subject comparisons of gene expression levels showed generally higher expressions of all three genes of interest in VAT than in SAT. Future in vitro analyses are needed to assess the functional implications of these findings.

Analysis of common genetic variation across targets of microRNAs dysregulated both in ASD and epilepsy reveals negative correlation.

Genetic overlap involving rare disrupting mutations may contribute to high comorbidity rates between autism spectrum disorders and epilepsy. Despite their polygenic nature, genome-wide association studies have not reported a significant contribution of common genetic variation to comorbidity between both conditions. Analysis of common genetic variation affecting specific shared pathways such as miRNA dysregulation could help to elucidate the polygenic mechanisms underlying comorbidity between autism spectrum disorders and epilepsy. We evaluated here the role of common predisposing variation to autism spectrum disorders and epilepsy across target genes of 14 miRNAs selected through bibliographic research as being dysregulated in both disorders. We considered 4,581 target genes from various in silico sources. We described negative genetic correlation between autism spectrum disorders and epilepsy across variants located within target genes of the 14 miRNAs selected (p = 0.0228). Moreover, polygenic transmission disequilibrium test on an independent cohort of autism spectrum disorders trios (N = 233) revealed an under-transmission of autism spectrum disorders predisposing alleles within miRNAs' target genes across autism spectrum disorders trios without comorbid epilepsy, thus reinforcing the negative relationship at the common genetic variation between both traits. Our study provides evidence of a negative relationship between autism spectrum disorders and epilepsy at the common genetic variation level that becomes more evident when focusing on the miRNA regulatory networks, which contrasts with observed clinical comorbidity and results from rare variation studies. Our findings may help to conceptualize the genetic heterogeneity and the comorbidity with epilepsy in autism spectrum disorders.

Ancestry-dependent genetic structure of the Xq28 risk haplotype in the Mexican population and its association with childhood-onset systemic lupus erythematosus.

Here we aimed to investigate the association of the Xq28 risk haplotype (H1) with susceptibility to childhood-onset systemic lupus erythematosus (SLE), and to compare its frequency and genetic structure in the Mexican population with those in other continental populations. We genotyped 15 single-nucleotide variants (SNVs) that form the H1 haplotype, using TaqMan real-time PCR. The association analysis [case-control and transmission disequilibrium test (TDT)] included 376 cases and 400 adult controls, all of whom were mestizos (MEZ). To identify risk alleles in Mexican Indigenous individuals, SNVs were imputed from whole-exome sequencing data of 1,074 individuals. The allelic frequencies determined in MEZ and Indigenous individuals were compared with those of the continental populations from the 1,000 Genomes database phase 3. Linkage disequilibrium (LD) analysis of risk alleles was performed on all populations. Interleukin-1 receptor associated kinase 1 (IRAK1) and methyl CpG binding protein 2 (MECP2) mRNA levels were determined using real-time PCR. Case-control analysis revealed genetic association with childhood-onset SLE for all 15 SNVs (OR = 1.49-1.75; p = 0.0095 to 1.81 × 10-4) and for the Xq28 risk haplotype (OR = 1.97, p = 4 × 10-6). Comparing with individuals of European ancestry (0.14-0.16), the frequencies of the risk alleles were significantly higher in the MEZ individuals (0.55-0.68) and even higher in Indigenous individuals (0.57-0.83). LD analysis indicated a differential haplotype structure within the Indigenous groups, which was inherited to the MEZ population as a result of genetic admixture. Individuals homozygous for the Xq28 risk haplotype exhibited decreased levels of both MECP2A and B transcripts. We found that the H1 risk haplotype differs in its conformation in the Mexican population. This difference could be attributed to positive selection within the Indigenous population, with its inheritance now having an autoimmune health impact in both the Mexican Indigenous and MEZ populations.