Search

The primary goal of the unrelated population HLA diversity (UPHD) component of the 17th International HLA and Immunogenetics Workshop was to characterize HLA alleles at maximum allelic-resolution in worldwide populations and re-evaluate patterns of HLA diversity across populations. The UPHD project included HLA genotype and sequence data, generated by various next-generation sequencing methods, from 4,240 individuals collated from 12 different countries. Population data included well-defined large datasets from the USA and smaller samples from Europe, Australia, and Western Asia. Allele and haplotype frequencies varied across populations from distant geographical regions. HLA genetic diversity estimated at 2- and 4-field allelic resolution revealed that diversity at the majority of loci, particularly for European-descent populations, was lower at the 2-field resolution. Several common alleles with identical protein sequences differing only by intronic substitutions were found in distinct haplotypes, revealing a more detailed characterization of linkage between variants within the HLA region. The examination of coding and non-coding nucleotide variation revealed many examples in which almost complete biunivocal relations between common alleles at different loci were observed resulting in higher linkage disequilibrium. Our reference data of HLA profiles characterized at maximum resolution from many populations is useful for anthropological studies, unrelated donor searches, transplantation, and disease association studies.

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10-31).

Plasma lipids are known heritable risk factors for cardiovascular disease, but increasing evidence also supports shared genetics with diseases of other organ systems. We devised a comprehensive three-phase framework to identify new lipid-associated genes and study the relationships among lipids, genotypes, gene expression and hundreds of complex human diseases from the Electronic Medical Records and Genomics (347 traits) and the UK Biobank (549 traits). Aside from 67 new lipid-associated genes with strong replication, we found evidence for pleiotropic SNPs/genes between lipids and diseases across the phenome. These include discordant pleiotropy in the HLA region between lipids and multiple sclerosis and putative causal paths between triglycerides and gout, among several others. Our findings give insights into the genetic basis of the relationship between plasma lipids and diseases on a phenome-wide scale and can provide context for future prevention and treatment strategies.

Abstract Background In studies of time-to-events, it is common to collect information about events that occurred before the inclusion in a prospective cohort. When the studied risk factors are independent of time, including both pre- and post-inclusion events in the analyses, generally referred to as relying on an ambispective design, increases the statistical power but may lead to a selection bias. In the field of venous thromboembolism (VT), ABO blood groups have been the subject of extensive research due to their substantial effect on VT risk. However, few studies have investigated their effect on the risk of VT recurrence. Motivated by the study of the association of genetically determined ABO blood groups with VT recurrence, we propose a methodology to include pre-inclusion events in the analysis of ambispective studies while avoiding the selection bias due to mortality. Methods This work relies on two independent cohorts of VT patients, the French MARTHA study built on an ambispective design and the Dutch MEGA study built on a standard prospective design. For the analysis of the MARTHA study, a weighted Cox model was developed where weights were defined by the inverse of the survival probability at the time of data collection about the events. Thanks to the collection of information on the vital status of patients, we could estimate the survival probabilities using a delayed-entry Cox model on the death risk. Finally, results obtained in both studies were then meta-analysed. Results In the combined sample totalling 2,752 patients including 993 recurrences, the A1 blood group has an increased risk (Hazard Ratio (HR) of 1.18, p = 4.2 × 10–3) compared with the O1 group, homogeneously in MARTHA and in MEGA. The same trend (HR = 1.19, p = 0.06) was observed for the less frequent A2 group. Conclusion The proposed methodology increases the power of studies relying on an ambispective design which is frequent in epidemiologic studies about recurrent events. This approach allowed to clarify the association of ABO blood groups with the risk of VT recurrence. Besides, this methodology has an immediate field of application in the context of genome wide association studies.

Cumulative meta-analysis follows a methodical rather than a statistical approach, to show evolutions over time. Cumulative meta-analysis includes repeated meta-analyses over time of newly published studies on the research question. Thus, a meta-analysis is conducted each time a new study is published, and the pooled estimates and confidence intervals from all the meta-analyses are graphically displayed in chronological sequence. In this way, retrospectively it is possible to determine the point in time from which a statistically significant effect becomes apparent. Recursive cumulative meta-analyses find changes in the common estimator by performing a meta-analysis of the chronologically arranged papers at the conclusion of each year. The strong assumption of a delay bias (“lag time bias”) is supported by an ever-declining effect. Although there are a number of software programs that perform meta-analysis, none of them have been expressly created for genetic association studies, and using them frequently necessitates advanced programming or scripting skills. Of course, they are unable to do cumulative and recursive cumulative meta-analyses. We developed Genetical Evidence SYNthesis (GESYN), a web-application for cumulative and recursive cumulative meta-analysis of Genetic association studies (GAS). GESYN can be run in an environment that is easy to use.11 The Web-application is available free of charge after registration at https://biomath.med.uth.gr/.

Determining genotype-phenotype correlations in patients with hypodontia is important for understanding disease pathogenesis, although only a few studies have elucidated it. We aimed to identify genetic variants linked to non-syndromic bilateral mandibular second premolar hypodontia in a Korean population for the first time by specifying the phenotype of hypodontia. Twenty unrelated individuals with non-syndromic bilateral mandibular second premolar hypodontia were enrolled for whole-exome sequencing. Using a tooth agenesis gene set panel consisting of 112 genes based on literature, potential candidate variants were screened through variant filtering and prioritization. We identified 13 candidate variants in 12 genes, including a stop-gain variant (c.4750C>T) in LAMA3. Through the functional enrichment analysis of the prioritized genes, several terms related to tooth development were enriched in a protein-protein interaction network of candidate genes for mandibular premolar hypodontia. The hypodontia group also had approximately 2-fold as many mutated variants in all four genes related to these key terms, which are CDH1, ITGB4, LAMA3, LAMB3, as those in the 100 healthy control group individuals. The relationship between enriched terms and pathways and mandibular premolar hypodontia was also investigated. In addition, we identified some known oligodontia variants in patients with hypodontia, strengthening the possibility of synergistic effects in other genes. This genetic investigation may be a worthwhile preliminary attempt to reveal the pathogenesis of tooth agenesis and sets a background for future studies. [ABSTRACT FROM AUTHOR]

In recent research, there has been a growing awareness of the role of genetic factors in the positioning and eruption of teeth in the maxilla and mandible. This study aimed to evaluate the potential of specific polymorphic markers of single nucleotide polymorphisms (SNPs) located within the PAX9, MSX1, AXIN2, and IRF6 genes to determine the predisposition to tooth impaction. The study participants were divided into two groups: the first group consisted of individuals with at least one impacted secondary tooth. In contrast, the second group (control group) had no impacted teeth in their jaws. To analyze the genes, real-time PCR (polymerase chain reaction) and TaqMan probes were utilized to detect the selected polymorphisms. The findings suggest that disruptions in the structure and function of the mentioned genetic factors such as polymorphic and haplotype variants of PAX9, MSX1, AXIN2, and IRF6 genes, which play a direct role in tooth and periodontal tissue development, might be significant factors in tooth impaction in individuals with genetic variations. Therefore, it is reasonable to hypothesize that tooth impaction may be influenced, at least in part, by the presence of specific genetic markers, including different allelic variants of the PAX9, AXIN2, and IRF6 genes, and especially MSX1. [ABSTRACT FROM AUTHOR]

Background: Observational studies have shown an association between higher bilirubin levels and improved respiratory health outcomes. Targeting higher bilirubin levels has been proposed as a novel therapeutic strategy in COPD. However, bilirubin levels are influenced by multiple intrinsic and extrinsic factors, and these observational studies are prone to confounding. Genetic analyses are one approach to overcoming residual confounding in observational studies. Objectives: To test associations between a genetic determinant of bilirubin levels and respiratory health outcomes. Methods: COPDGene participants underwent genotyping at the baseline visit. We confirmed established associations between homozygosity for rs6742078 and higher bilirubin, and between higher bilirubin and decreased risk of acute respiratory events within this cohort. For our primary analysis, we used negative binomial regression to test associations between homozygosity for rs6742078 and rate of acute respiratory events. Results: 8,727 participants (n = 6,228 non-Hispanic white and 2,499 African American) were included. Higher bilirubin was associated with decreased rate of acute respiratory events [incidence rate ratio (IRR) 0.85, 95% CI 0.75 to 0.96 per SD increase in bilirubin intensity]. We did not find significant associations between homozygosity for rs6742078 and acute respiratory events (IRR 0.94, 95% CI 0.70 to 1.25 for non-Hispanic white and 1.09, 95% CI 0.91 to 1.31 for African American participants). Conclusions: A genetic determinant of higher bilirubin levels was not associated with better respiratory health outcomes. These results do not support targeting higher bilirubin levels as a therapeutic strategy in COPD. [ABSTRACT FROM AUTHOR]

The single-nucleotide polymorphism (SNP) rs2235371 (IRF6 V274I) is associated with nonsyndromic cleft lip with or without cleft palate (NSCL/P) in Han Chinese and other populations but appears to be without a functional effect. To find the common etiologic variant or variants within the haplotype tagged by rs2235371, we carried out targeted sequencing of an interval containing IRF6 in 159 Han Chinese with NSCL/P. This study revealed that the SNP rs12403599, within the IRF6 promoter, is associated with all phenotypes of NSCL/P, especially nonsyndromic cleft lip (NSCLO) and a subphenotype of it, microform cleft lip (MCL). This association was replicated in 2 additional much larger cohorts of cases and controls from the Han Chinese. Conditional logistic analysis indicated that association of rs2235371 with NSCL/P was lost if rs12403599 was excluded. rs12403599 contributes the most risk to MCL: its G allele is responsible for 38.47% of the genetic contribution to MCL, and the odds ratios of G/C and G/G genotypes were 2.91 and 6.58, respectively, for MCL. To test if rs12403599 is functional, we carried out reporter assays in a fetal oral epithelium cells (GMSM-K). Unexpectedly, the risk allele G yielded higher promoter activity in GMSM-K. Consistent with the reporter studies, expression of IRF6 in lip tissues from NSCLO and MCL patients with the G/G phenotype was higher than in those from patients with the C/C phenotype. These results indicate that rs12403599 is tagging the risk haplotype for NSCL/P better than rs2235371 in Han Chinese and supports investigation of the mechanisms by which the allele of rs12403599 affects IRF6 expression and tests of this association in different populations. [ABSTRACT FROM AUTHOR]

Wiedemann-Steiner syndrome (WSS) is an autosomal dominant disorder caused by monoallelic variants in KMT2A and characterized by intellectual disability and hypertrichosis. We performed a retrospective, multicenter, observational study of 104 individuals with WSS from five continents to characterize the clinical and molecular spectrum of WSS in diverse populations, to identify physical features that may be more prevalent in White versus Black Indigenous People of Color individuals, to delineate genotype-phenotype correlations, to define developmental milestones, to describe the syndrome through adulthood, and to examine clinicians' differential diagnoses. Sixty-nine of the 82 variants (84%) observed in the study were not previously reported in the literature. Common clinical features identified in the cohort included: developmental delay or intellectual disability (97%), constipation (63.8%), failure to thrive (67.7%), feeding difficulties (66.3%), hypertrichosis cubiti (57%), short stature (57.8%), and vertebral anomalies (46.9%). The median ages at walking and first words were 20 months and 18 months, respectively. Hypotonia was associated with loss of function (LoF) variants, and seizures were associated with non-LoF variants. This study identifies genotype-phenotype correlations as well as race-facial feature associations in an ethnically diverse cohort, and accurately defines developmental trajectories, medical comorbidities, and long-term outcomes in individuals with WSS.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.

Prioritizing genes for translation to therapeutics for common diseases has been challenging. Here, we propose an approach to identify drug targets with high probability of success by focusing on genes with both gain of function (GoF) and loss of function (LoF) mutations associated with opposing effects on phenotype (Bidirectional Effect Selected Targets, BEST). We find 98 BEST genes for a variety of indications. Drugs targeting those genes are 3.8-fold more likely to be approved than non-BEST genes. We focus on five genes (IGF1R, NPPC, NPR2, FGFR3, and SHOX) with evidence for bidirectional effects on stature. Rare protein-altering variants in those genes result in significantly increased risk for idiopathic short stature (ISS) (OR = 2.75, p = 3.99 × 10-8). Finally, using functional experiments, we demonstrate that adding an exogenous CNP analog (encoded by NPPC) rescues the phenotype, thus validating its potential as a therapeutic treatment for ISS. Our results show the value of looking for bidirectional effects to identify and validate drug targets.

Autosomal genetic analyses of blood lipids have yielded key insights for coronary heart disease (CHD). However, X chromosome genetic variation is understudied for blood lipids in large sample sizes. We now analyze genetic and blood lipid data in a high-coverage whole X chromosome sequencing study of 65,322 multi-ancestry participants and perform replication among 456,893 European participants. Common alleles on chromosome Xq23 are strongly associated with reduced total cholesterol, LDL cholesterol, and triglycerides (min P = 8.5 × 10-72), with similar effects for males and females. Chromosome Xq23 lipid-lowering alleles are associated with reduced odds for CHD among 42,545 cases and 591,247 controls (P = 1.7 × 10-4), and reduced odds for diabetes mellitus type 2 among 54,095 cases and 573,885 controls (P = 1.4 × 10-5). Although we observe an association with increased BMI, waist-to-hip ratio adjusted for BMI is reduced, bioimpedance analyses indicate increased gluteofemoral fat, and abdominal MRI analyses indicate reduced visceral adiposity. Co-localization analyses strongly correlate increased CHRDL1 gene expression, particularly in adipose tissue, with reduced concentrations of blood lipids.

The International Mouse Phenotyping Consortium reports the generation of new mouse mutant strains for over 5,000 genes, including 2,850 novel null, 2,987 novel conditional- ready, and 4,433 novel reporter alleles.

The most recent updated classification of inborn errors of immunity/primary immunodeficiencies, compiled by the International Union of Immunological Societies Expert Committee, was published in January 2020. Within days of completing this report, it was already out of date, evidenced by the frequent publication of genetic variants proposed to cause novel inborn errors of immunity. As the next formal report from the IUIS Expert Committee will not be published until 2022, we felt it important to provide the community with a brief update of recent contributions to the field of inborn errors of immunity. Herein, we highlight studies that have identified 26 additional monogenic gene defects that reach the threshold to represent novel causes of immune defects.

Studies of acute myeloid leukemia rely on DNA sequencing and immunophenotyping by flow cytometry as primary tools for disease characterization. However, leukemia tumor heterogeneity complicates integration of DNA variants and immunophenotypes from separate measurements. Here we introduce DAb-seq, a technology for simultaneous capture of DNA genotype and cell surface phenotype from single cells at high throughput, enabling direct profiling of proteogenomic states in tens of thousands of cells. To demonstrate the approach, we analyze the disease of three patients with leukemia over multiple treatment timepoints and disease recurrences. We observe complex genotype-phenotype dynamics that illustrate the subtlety of the disease process and the degree of incongruity between blast cell genotype and phenotype in different clinical scenarios. Our results highlight the importance of combined single-cell DNA and protein measurements to fully characterize the heterogeneity of leukemia.

BackgroundThe Chr17q12-21.2 region is the strongest and most consistently associated region with asthma susceptibility. The functional genes or single nucleotide polymorphisms (SNPs) are not obvious due to linkage disequilibrium.ObjectivesWe sought to comprehensively investigate whole-genome sequence and RNA sequence from human bronchial epithelial cells to dissect functional genes/SNPs for asthma severity in the Severe Asthma Research Program.MethodsExpression quantitative trait loci analysis (n = 114), correlation analysis (n = 156) of gene expression and asthma phenotypes, and pathway analysis were performed in bronchial epithelial cells and replicated. Genetic association for asthma severity (426 severe vs 531 nonsevere asthma) and longitudinal asthma exacerbations (n = 273) was performed.ResultsMultiple SNPs in gasdermin B (GSDMB) associated with asthma severity (odds ratio, >1.25) and longitudinal asthma exacerbations (P < .05). Expression quantitative trait loci analyses identified multiple SNPs associated with expression levels of post-GPI attachment to proteins 3, GSDMB, or gasdermin A (3.1 × 10-9 -4). Higher expression levels of GSDMB correlated with asthma and greater number of exacerbations (P < .05). Expression levels of GSDMB correlated with genes involved in IFN signaling, MHC class I antigen presentation, and immune system pathways (false-discovery rate-adjusted P < .05). rs1031458 and rs3902920 in GSDMB colocalized with IFN regulatory factor binding sites and associated with GSDMB expression, asthma severity, and asthma exacerbations (P < .05).ConclusionsBy using a unique set of gene expression data from lung cells obtained using bronchoscopy from comprehensively characterized subjects with asthma, we show that SNPs in GSDMB associated with asthma severity, exacerbations, and GSDMB expression levels. Furthermore, its expression levels correlated with asthma exacerbations and antiviral pathways. Thus, GSDMB is a functional gene for both asthma susceptibility and severity.

Rationale: The 17q12-21.1 locus is one of the most highly replicated genetic associations with asthma. Individuals of African descent have lower linkage disequilibrium in this region, which could facilitate identifying causal variants.Objectives: To identify functional variants at 17q12-21.1 associated with early-onset asthma among African American individuals.Methods: We evaluated African American participants from SAPPHIRE (Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity) (n = 1,940), SAGE II (Study of African Americans, Asthma, Genes and Environment) (n = 885), and GCPD-A (Study of the Genetic Causes of Complex Pediatric Disorders-Asthma) (n = 2,805). Associations with asthma onset at ages under 5 years were meta-analyzed across cohorts. The lead signal was reevaluated considering haplotypes informed by genetic ancestry (i.e., African vs. European). Both an expression-quantitative trait locus analysis and a phenome-wide association study were performed on the lead variant.Measurements and Main Results: The meta-analyzed results from SAPPHIRE, SAGE II, and the GCPD-A identified rs11078928 as the top association for early-onset asthma. A haplotype analysis suggested that the asthma association partitioned most closely with the rs11078928 genotype. Genetic ancestry did not appear to influence the effect of this variant. In the expression-quantitative trait locus analysis, rs11078928 was related to alternative splicing of GSDMB (gasdermin-B) transcripts. The phenome-wide association study of rs11078928 suggested that this variant was predominantly associated with asthma and asthma-associated symptoms.Conclusions: A splice-acceptor polymorphism appears to be a causal variant for asthma at the 17q12-21.1 locus. This variant appears to have the same magnitude of effect in individuals of African and European descent.

Whereas large-scale statistical analyses can robustly identify disease-gene relationships, they do not accurately capture genotype-phenotype correlations or disease mechanisms. We use multiple lines of independent evidence to show that different variant types in a single gene, SATB1, cause clinically overlapping but distinct neurodevelopmental disorders. Clinical evaluation of 42 individuals carrying SATB1 variants identified overt genotype-phenotype relationships, associated with different pathophysiological mechanisms, established by functional assays. Missense variants in the CUT1 and CUT2 DNA-binding domains result in stronger chromatin binding, increased transcriptional repression, and a severe phenotype. In contrast, variants predicted to result in haploinsufficiency are associated with a milder clinical presentation. A similarly mild phenotype is observed for individuals with premature protein truncating variants that escape nonsense-mediated decay, which are transcriptionally active but mislocalized in the cell. Our results suggest that in-depth mutation-specific genotype-phenotype studies are essential to capture full disease complexity and to explain phenotypic variability.

The integration of proteomic, transcriptomic, and genetic variant annotation data will improve our understanding of genotype-phenotype associations. Due, in part, to challenges associated with accurate inter-database mapping, such multi-omic studies have not extended to chemoproteomics, a method that measures the intrinsic reactivity and potential "druggability" of nucleophilic amino acid side chains. Here, we evaluated mapping approaches to match chemoproteomic-detected cysteine and lysine residues with their genetic coordinates. Our analysis revealed that database update cycles and reliance on stable identifiers can lead to pervasive misidentification of labeled residues. Enabled by this examination of mapping strategies, we then integrated our chemoproteomics data with computational methods for predicting genetic variant pathogenicity, which revealed that codons of highly reactive cysteines are enriched for genetic variants that are predicted to be more deleterious and allowed us to identify and functionally characterize a new damaging residue in the cysteine protease caspase-8. Our study provides a roadmap for more precise inter-database mapping and points to untapped opportunities to improve the predictive power of pathogenicity scores and to advance prioritization of putative druggable sites.

Along with early-onset cancers, multiple primary cancers (MPCs) are likely resulting from increased genetic susceptibility; however, the associated predisposition genes or prevalence of the pathogenic variants genes in MPC patients are often unknown. We screened 71 patients with MPC of the stomach, colorectal, and endometrium, sequencing 65 cancer predisposition genes. A subset of 19 patients with early-onset MPC of stomach and colorectum were further evaluated for variants in cancer related genes using both normal and tumor whole exome sequencing. Among 71 patients with MPCs, variants classified to be pathogenic were observed in 15 (21.1%) patients and affected Lynch Syndrome (LS) genes: MLH1 (n = 10), MSH6 (n = 2), PMS2 (n = 2), and MSH2 (n = 1). All carriers had tumors with high microsatellite instability and 13 of them (86.7%) were early-onset, consistent with LS. In 19 patients with early-onset MPCs, loss of function (LoF) variants in RECQL5 were more prevalent in non-LS MPC than in matched sporadic cancer patients (OR = 31.6, 2.73-1700.6, p = 0.001). Additionally, there were high-confidence LoF variants at FANCG and CASP8 in two patients accompanied by somatic loss of heterozygosity in tumor, respectively. The results suggest that genetic screening should be considered for synchronous cancers and metachronous MPCs of the LS tumor spectrum, particularly in early-onset. Susceptibility variants in non-LS genes for MPC patients may exist, but evidence for their role is more elusive than for LS patients.

ObjectiveCertain copy number variants (CNVs) greatly increase the risk of autism. The authors conducted a genetics-first study to investigate whether heterogeneity in the clinical presentation of autism is underpinned by specific genotype-phenotype relationships.MethodsThis international study included 547 individuals (mean age, 12.3 years [SD=4.2], 54% male) who were ascertained on the basis of having a genetic diagnosis of a rare CNV associated with high risk of autism (82 16p11.2 deletion carriers, 50 16p11.2 duplication carriers, 370 22q11.2 deletion carriers, and 45 22q11.2 duplication carriers), as well as 2,027 individuals (mean age, 9.1 years [SD=4.9], 86% male) with autism of heterogeneous etiology. Assessments included the Autism Diagnostic Interview-Revised and IQ testing.ResultsThe four genetic variant groups differed in autism symptom severity, autism subdomain profile, and IQ profile. However, substantial variability was observed in phenotypic outcome in individual genetic variant groups (74%-97% of the variance, depending on the trait), whereas variability between groups was low (1%-21%, depending on the trait). CNV carriers who met autism criteria were compared with individuals with heterogeneous autism, and a range of profile differences were identified. When clinical cutoff scores were applied, 54% of individuals with one of the four CNVs who did not meet full autism diagnostic criteria had elevated levels of autistic traits.ConclusionsMany CNV carriers do not meet full diagnostic criteria for autism but nevertheless meet clinical cutoffs for autistic traits. Although profile differences between variants were observed, there is considerable variability in clinical symptoms in the same variant.

Arthropod-borne rickettsial pathogens cause mild and severe human disease worldwide. The tick-borne pathogen Rickettsia parkeri elicits skin lesions (eschars) and disseminated disease in humans; however, inbred mice are generally resistant to infection. We report that intradermal infection of mice lacking both interferon receptors (Ifnar1-/-;Ifngr1-/-) with as few as 10 R. parkeri elicits eschar formation and disseminated, lethal disease. Similar to human infection, eschars exhibited necrosis and inflammation, with bacteria primarily found in leukocytes. Using this model, we find that the actin-based motility factor Sca2 is required for dissemination from the skin to internal organs, and the outer membrane protein OmpB contributes to eschar formation. Immunizing Ifnar1-/-;Ifngr1-/- mice with sca2 and ompB mutant R. parkeri protects against rechallenge, revealing live-attenuated vaccine candidates. Thus, Ifnar1-/-;Ifngr1-/- mice are a tractable model to investigate rickettsiosis, virulence factors, and immunity. Our results further suggest that discrepancies between mouse and human susceptibility may be due to differences in interferon signaling.

People with schizophrenia are enriched for rare coding variants in genes associated with neurodevelopmental disorders, particularly autism spectrum disorders and intellectual disability. However, it is unclear if the same changes to gene function that increase risk to neurodevelopmental disorders also do so for schizophrenia. Using data from 3444 schizophrenia trios and 37,488 neurodevelopmental disorder trios, we show that within shared risk genes, de novo variants in schizophrenia and neurodevelopmental disorders are generally of the same functional category, and that specific de novo variants observed in neurodevelopmental disorders are enriched in schizophrenia (P = 5.0 × 10-6). The latter includes variants known to be pathogenic for syndromic disorders, suggesting that schizophrenia be included as a characteristic of those syndromes. Our findings imply that, in part, neurodevelopmental disorders and schizophrenia have shared molecular aetiology, and therefore likely overlapping pathophysiology, and support the hypothesis that at least some forms of schizophrenia lie on a continuum of neurodevelopmental disorders.

The ectodysplasin receptor (EDAR) is a tumor necrosis factor receptor (TNF) superfamily member. A substitution in an exon of EDAR at position 370 (EDARV370A) creates a gain of function mutant present at high frequencies in Asian and Indigenous American populations but absent in others. Its frequency is intermediate in populations of Mexican ancestry. EDAR regulates the development of ectodermal tissues, including mammary ducts. Obesity and type 2 diabetes mellitus are prevalent in people with Indigenous and Latino ancestry. Latino patients also have altered prevalence and presentation of breast cancer. It is unknown whether EDARV370A might connect these phenomena. The goals of this study were to determine 1) whether EDARV370A is associated with metabolic phenotypes and 2) if there is altered breast anatomy in women carrying EDARV370A. Participants were from two Latino cohorts, the Arizona Insulin Resistance (AIR) registry and Sangre por Salud (SPS) biobank. The frequency of EDARV370A was 47% in the Latino cohorts. In the AIR registry, carriers of EDARV370A (GG homozygous) had significantly (p < 0.05) higher plasma triglycerides, VLDL, ALT, 2-hour post-challenge glucose, and a higher prevalence of prediabetes/diabetes. In a subset of the AIR registry, serum levels of ectodysplasin A2 (EDA-A2) also were associated with HbA1c and prediabetes (p < 0.05). For the SPS biobank, participants that were carriers of EDARV370A had lower breast density and higher HbA1c (both p < 0.05). The significant associations with measures of glycemia remained when the cohorts were combined. We conclude that EDARV370A is associated with characteristics of the metabolic syndrome and breast density in Latinos.

BackgroundWith the increasing number of genomic sequencing studies, hundreds of genes have been implicated in neurodevelopmental disorders (NDDs). The rate of gene discovery far outpaces our understanding of genotype-phenotype correlations, with clinical characterization remaining a bottleneck for understanding NDDs. Most disease-associated Mendelian genes are members of gene families, and we hypothesize that those with related molecular function share clinical presentations.MethodsWe tested our hypothesis by considering gene families that have multiple members with an enrichment of de novo variants among NDDs, as determined by previous meta-analyses. One of these gene families is the heterogeneous nuclear ribonucleoproteins (hnRNPs), which has 33 members, five of which have been recently identified as NDD genes (HNRNPK, HNRNPU, HNRNPH1, HNRNPH2, and HNRNPR) and two of which have significant enrichment in our previous meta-analysis of probands with NDDs (HNRNPU and SYNCRIP). Utilizing protein homology, mutation analyses, gene expression analyses, and phenotypic characterization, we provide evidence for variation in 12 HNRNP genes as candidates for NDDs. Seven are potentially novel while the remaining genes in the family likely do not significantly contribute to NDD risk.ResultsWe report 119 new NDD cases (64 de novo variants) through sequencing and international collaborations and combined with published clinical case reports. We consider 235 cases with gene-disruptive single-nucleotide variants or indels and 15 cases with small copy number variants. Three hnRNP-encoding genes reach nominal or exome-wide significance for de novo variant enrichment, while nine are candidates for pathogenic mutations. Comparison of HNRNP gene expression shows a pattern consistent with a role in cerebral cortical development with enriched expression among radial glial progenitors. Clinical assessment of probands (n = 188-221) expands the phenotypes associated with HNRNP rare variants, and phenotypes associated with variation in the HNRNP genes distinguishes them as a subgroup of NDDs.ConclusionsOverall, our novel approach of exploiting gene families in NDDs identifies new HNRNP-related disorders, expands the phenotypes of known HNRNP-related disorders, strongly implicates disruption of the hnRNPs as a whole in NDDs, and supports that NDD subtypes likely have shared molecular pathogenesis. To date, this is the first study to identify novel genetic disorders based on the presence of disorders in related genes. We also perform the first phenotypic analyses focusing on related genes. Finally, we show that radial glial expression of these genes is likely critical during neurodevelopment. This is important for diagnostics, as well as developing strategies to best study these genes for the development of therapeutics.

Coat color is a trait of economic significance in horses. Variants in seven genes have been documented to cause white patterning in horses. Of the 34 variants that have been identified in KIT proto-oncogene, receptor tyrosine kinase (KIT), 27 have only been reported in a single individual or family and thus not all are routinely offered for genetic testing. Therefore, to enable proper use of marker-assisted selection, determining breed specificity for these alleles is warranted. Screening 19 unregistered all-white Shetland ponies for 16 white patterning markers identified 14 individuals whose phenotype could not be explained by testing results. In evaluating other known dominant white variants, 14 horses were heterozygous for W13. W13 was previously only reported in two quarter horses and a family of Australian miniature horses. Genotyping known white spotting variants in 30 owner-reported white animals (25 Miniature Horses and five Shetland ponies) identified two additional W13/N American Miniature Horses. The estimated allele frequency of W13 in the American Miniature Horse was 0.0063 (79 N/N, 1 W13/N) and the allele was not detected in a random sample (n = 59) of Shetland ponies. No homozygous W13 individuals were identified and W13/N ponies had a similar all-white coat with pink skin phenotype, regardless of the other white spotting variants present, demonstrating that W13 results in a Mendelian inherited dominant white phenotype and homozygosity is likely lethal. These findings document the presence of W13 in the American Miniature Horse and Shetland pony populations at a low frequency and illustrate the importance of testing for this variant in additional breeds.

Patients with inherited retinal dystrophies (IRDs) were recruited from two understudied populations: Mexico and Pakistan as well as a third well-studied population of European Americans to define the genetic architecture of IRD by performing whole-genome sequencing (WGS). Whole-genome analysis was performed on 409 individuals from 108 unrelated pedigrees with IRDs. All patients underwent an ophthalmic evaluation to establish the retinal phenotype. Although the 108 pedigrees in this study had previously been examined for mutations in known IRD genes using a wide range of methodologies including targeted gene(s) or mutation(s) screening, linkage analysis and exome sequencing, the gene mutations responsible for IRD in these 108 pedigrees were not determined. WGS was performed on these pedigrees using Illumina X10 at a minimum of 30X depth. The sequence reads were mapped against hg19 followed by variant calling using GATK. The genome variants were annotated using SnpEff, PolyPhen2, and CADD score; the structural variants (SVs) were called using GenomeSTRiP and LUMPY. We identified potential causative sequence alterations in 61 pedigrees (57%), including 39 novel and 54 reported variants in IRD genes. For 57 of these pedigrees the observed genotype was consistent with the initial clinical diagnosis, the remaining 4 had the clinical diagnosis reclassified based on our findings. In seven pedigrees (12%) we observed atypical causal variants, i.e. unexpected genotype(s), including 4 pedigrees with causal variants in more than one IRD gene within all affected family members, one pedigree with intrafamilial genetic heterogeneity (different affected family members carrying causal variants in different IRD genes), one pedigree carrying a dominant causative variant present in pseudo-recessive form due to consanguinity and one pedigree with a de-novo variant in the affected family member. Combined atypical and large structural variants contributed to about 20% of cases. Among the novel mutations, 75% were detected in Mexican and 50% found in European American pedigrees and have not been reported in any other population while only 20% were detected in Pakistani pedigrees and were not previously reported. The remaining novel IRD causative variants were listed in gnomAD but were found to be very rare and population specific. Mutations in known IRD associated genes contributed to pathology in 63% Mexican, 60% Pakistani and 45% European American pedigrees analyzed. Overall, contribution of known IRD gene variants to disease pathology in these three populations was similar to that observed in other populations worldwide. This study revealed a spectrum of mutations contributing to IRD in three populations, identified a large proportion of novel potentially causative variants that are specific to the corresponding population or not reported in gnomAD and shed light on the genetic architecture of IRD in these diverse global populations.

"Quantile-dependent expressivity" occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g., leptin) is high or low relative to its distribution. Leptin concentrations are strongly related to adiposity, whose heritability is quantile dependent. Whether inheritance of leptin concentrations is quantile dependent, and whether this explains the greater heritability in women than men in accordance with their greater adiposity, and explains other gene-environment interactions, remains to be determined. Therefore, leptin and leptin receptor concentrations from 3068 siblings in 1133 sibships from the Framingham Heart Study Third Generation Cohort were analyzed. Free leptin index (FLI) was calculated as the ratio of leptin to soluble leptin receptor concentrations. Full-sib (βFS) regression slopes were robustly estimated by quantile regression with nonparametric significance assigned from 1000 bootstrap samples. The analyses showed βFS increased significantly with increasing percentiles of the offspring's age- and sex-adjusted leptin distribution (Plinear = 0.0001), which was accelerated at the higher concentrations (Pquadratic = 0.0003). βFS at the 90th percentile (0.418 ± 0.066) was 4.7-fold greater than at the 10th percentile (0.089 ± 0.032, Pdifference = 3.6 × 10-6). Consistent with quantile-dependent expressivity, the βFS was greater in female sibs, which was attributable to their higher leptin concentrations. Reported gene-environment interactions involving adiposity and LEP, LEPR, MnSOD, PPARγ, PPARγ2, and IRS-1 polymorphisms were consistent with quantile-dependent expressivity of leptin concentrations. βFS for leptin receptor concentrations and free leptin index also increased significantly with increasing percentiles of their distributions (Plinear = 0.04 and Plinear = 8.5 × 10-6, respectively). In conclusion, inherited genetic and shared environmental effects on leptin concentrations were quantile dependent, which likely explains male-female differences in heritability and some gene-environment interactions.

Therapeutic innovation for human papilloma virus-related (HPV+) head and neck squamous cell carcinomas (HNSCCs) is impaired by inadequate preclinical models and the absence of accurate biomarkers. Our study establishes the first well-characterized panel of patient-derived xenografts (PDXs) and organoids from HPV+ HNSCCs while determining fidelity of the models to the distinguishing genetic features of this cancer type. Despite low engraftment rates, whole exome sequencing showed that PDXs retain multiple distinguishing features of HPV+ HNSCC lost in existing cell lines, including PIK3CA mutations, TRAF3 deletion and the absence of EGFR amplifications. Engrafted HPV+ tumors frequently contained NOTCH1 mutations, thus providing new models for a negatively prognostic alteration in this disease. Genotype-phenotype associations in the models were then tested for prediction of tumor progression and survival in published clinical cohorts. Observation of high tumor mutational burdens (TMBs) in the faster-growing models facilitated identification of a novel association between TMB and local progression in both HPV+ and HPV- patients that was prognostic in HPV- cases. In addition, reduced E7 and p16INK4A levels found in a PDX from an outlier case with lethal outcome led to detection of similar profiles among recurrent HPV+ HNSCCs. Transcriptional data from the Cancer Genome Atlas was used to demonstrate that the lower E2F target gene expression predicted by reduced E7 levels has potential as a biomarker of disease recurrence risk. Our findings bridge a critical gap in preclinical models for HPV+ HNSCCs and simultaneously reveal novel potential applications of quantifying mutational burden and viral oncogene functions for biomarker development.

BackgroundBeyond its structural role in the skeleton, the extracellular matrix (ECM), particularly basement membrane proteins, facilitates communication with intracellular signaling pathways and cell to cell interactions to control differentiation, proliferation, migration and survival. Alterations in extracellular proteins cause a number of skeletal disorders, yet the consequences of an abnormal ECM on cellular communication remains less well understood METHODS: Clinical and radiographic examinations defined the phenotype in this unappreciated bent bone skeletal disorder. Exome analysis identified the genetic alteration, confirmed by Sanger sequencing. Quantitative PCR, western blot analyses, immunohistochemistry, luciferase assay for WNT signaling were employed to determine RNA, proteins levels and localization, and dissect out the underlying cell signaling abnormalities.  Migration and wound healing assays examined cell migration properties.FindingsThis bent bone dysplasia resulted from biallelic mutations in LAMA5, the gene encoding the alpha-5 laminin basement membrane protein. This finding uncovered a mechanism of disease driven by ECM-cell interactions between alpha-5-containing laminins, and integrin-mediated focal adhesion signaling, particularly in cartilage. Loss of LAMA5 altered β1 integrin signaling through the non-canonical kinase PYK2 and the skeletal enriched SRC kinase, FYN. Loss of LAMA5 negatively impacted the actin cytoskeleton, vinculin localization, and WNT signaling.InterpretationThis newly described mechanism revealed a LAMA5-β1 Integrin-PYK2-FYN focal adhesion complex that regulates skeletogenesis, impacted WNT signaling and, when dysregulated, produced a distinct skeletal disorder.FundingSupported by NIH awards R01 AR066124, R01 DE019567, R01 HD070394, and U54HG006493, and Czech Republic grants INTER-ACTION LTAUSA19030, V18-08-00567 and GA19-20123S.

BackgroundTo assess the genetic architecture of hypertrophic cardiomyopathy (HCM) in patients of predominantly Chinese ancestry.MethodsWe sequenced HCM disease genes in Singaporean patients (n=224) and Singaporean controls (n=3634), compared findings with additional populations and White HCM cohorts (n=6179), and performed in vitro functional studies.ResultsSingaporean HCM patients had significantly fewer confidently interpreted HCM disease variants (pathogenic/likely pathogenic: 18%, P

Objective: Interindividual variability in the clinical progression of COVID-19 may be explained by host genetics. Emerging literature supports a potential inherited predisposition to severe forms of COVID-19. Demographic and inflammatory characteristics of COVID-19 suggest that acquired hematologic mutations leading to clonal hematopoiesis (CH) may further increase vulnerability to adverse sequelae. This review summarizes the available literature examining genetic predispositions to severe COVID-19 and describes how these findings could eventually be used to improve its clinical management. DATA SOURCES: A PubMed literature search was performed. STUDY SELECTION: Studies examining the significance of inherited genetic variation or acquired CH mutations in severe COVID-19 were selected for inclusion. DATA EXTRACTION: Relevant genetic association data and aspects of study design were qualitatively assessed and narratively synthesized. DATA SYNTHESIS: Genetic variants affecting inflammatory responses may increase susceptibility to severe COVID-19. Genome-wide association studies and candidate gene approaches have identified a list of inherited mutations, which likely alter cytokine and interferon secretion, and lung-specific mechanisms of immunity in COVID-19. The potential role of CH in COVID-19 is more uncertain at present; however, the available evidence suggests that the various types of acquired mutations and their differential influence on immune cell function must be carefully considered. CONCLUSIONS: The current literature supports the hypothesis that host genetic factors affect vulnerability to severe COVID-19. Further research is required to confirm the full scope of relevant variants and the causal mechanisms underlying these associations. Clinical approaches, which consider the genetic basis of interindividual variability in COVID-19 and potentially other causes of critical illness, could optimize hospital resource allocation, predict responsiveness to treatment, identify more efficacious drug targets, and ultimately improve outcomes. [ABSTRACT FROM AUTHOR]

Dentinogenesis imperfecta (DI) is the main orodental manifestation of osteogenesis imperfecta (OI) caused by COL1A1 or COL1A2 heterozygous pathogenic variants. Its prevalence varies according to the studied population. Here, we report the molecular analysis of 81 patients with OI followed at reference centers in Brazil and France presenting COL1A1 or COL1A2 variants. Patients were submitted to clinical and radiographic dental examinations to diagnose the presence of DI. In addition, a systematic literature search and a descriptive statistical analysis were performed to investigate OI/DI phenotype–genotype correlation in a worldwide sample. In our cohort, 50 patients had COL1A1 pathogenic variants, and 31 patients had COL1A2 variants. A total of 25 novel variants were identified. Overall, data from a total of 906 individuals with OI were assessed. Results show that DI was more frequent in severe and moderate OI cases. DI prevalence was also more often associated with COL1A2 (67.6%) than with COL1A1 variants (45.4%) because COL1A2 variants mainly lead to qualitative defects that predispose to DI more than quantitative defects. For the first time, 4 DI hotspots were identified. In addition, we showed that 1) glycine substitution by branched and charged amino acids in the α2(I) chain and 2) substitutions occurring in major ligand binding regions—MLRB2 in α1(I) and MLBR 3 in α2(I)—could significantly predict DI (P < 0.05). The accumulated variant data analysis in this study provides a further basis for increasing our comprehension to better predict the occurrence and severity of DI and appropriate OI patient management. [ABSTRACT FROM AUTHOR]

The significance of rare germline mutations in transplant-associated thrombotic microangiopathy (TA-TMA) is not well studied. We performed a genetic association study in 100 adult TA-TMA patients vs. 98 post-transplant controls after matching by race, sex, and year. We focused on 5 pathways in complement, von Willebrand factor (VWF) function and related proteins, VWF clearance, ADAMTS13 function and related proteins, and endothelial activation (3641variants in 52 genes). In the primary analysis focused on 189 functional rare variants, no differential variant enrichment was observed in any of the pathways; specifically, 29 % TA-TMA and 33 % controls had at least 1 rare complement mutation. In the secondary analysis focused on 37 rare variants predicted to be pathogenic or likely pathogenic by ClinVar, Complement Database, or REVEL in-silico prediction tool, rare variants in the VWF clearance pathway were found to be significantly associated with TA-TMA (p = 0.008). On the gene level, LRP1 was the only one with significantly increased variants in TA-TMA in both analyses (p = 0.025 and 0.015). In conclusion, we did not find a significant association between rare variants in the complement pathway and TA-TMA; however, we discovered a new signal in the VWF clearance pathway driven by the gene LRP1 among likely pathogenic variants. • It is unclear if germline mutations in complement genes are associated with TA-TMA. • We performed a genetic association case-control study with 100 cases and 98 controls. • We did not observe rare variants in complement genes in association with TA-TMA. • 29% TA-TMA and 33% controls had at least 1 rare variant in complement genes. • Rare variants in LRP1 (VWF clearance/oxidative stress) were associated with TA-TMA. [ABSTRACT FROM AUTHOR]

Common genetic association models for structured populations, including principal component analysis (PCA) and linear mixed-effects models (LMMs), model the correlation structure between individuals using population kinship matrices, also known as genetic relatedness matrices. However, the most common kinship estimators can have severe biases that were only recently determined. Here we characterize the effect of these kinship biases on genetic association. We employ a large simulated admixed family and genotypes from the 1000 Genomes Project, both with simulated traits, to evaluate key kinship estimators. Remarkably, we find practically invariant association statistics for kinship matrices of different bias types (matching all other features). We then prove using statistical theory and linear algebra that LMM association tests are invariant to these kinship biases, and PCA approximately so. Our proof shows that the intercept and relatedness effect coefficients compensate for the kinship bias, an argument that extends to generalized linear models. As a corollary, association testing is also invariant to changing the reference ancestral population of the kinship matrix. Lastly, we observed that all kinship estimators, except for popkin ratio-of-means, can give improper non-positive semidefinite matrices, which can be problematic although some LMMs handle them surprisingly well, and condition numbers can be used to choose kinship estimators. Overall, we find that existing association studies are robust to kinship estimation bias, and our calculations may help improve association methods by taking advantage of this unexpected robustness, as well as help determine the effects of kinship bias in related problems. [ABSTRACT FROM AUTHOR]

The FTO rs9939609 gene, which presents three polymorphisms (AA, AT, and TT), has been associated with the development of obesity through an increased fat accumulation; however, the associations of the gene with other physiological mechanisms, such as appetite or fat oxidation, are still unclear. Therefore, this study aims to evaluate the influence of the FTO rs9939609 gene on different obesity-related factors in young adults. The FTO rs9939609 polymorphism was genotyped in 73 participants (28 women, 22.27 ± 3.70 years). Obesity-related factors included dietary assessment, physical activity expenditure, body composition, appetite sensation, resting metabolic rate, maximal fat oxidation during exercise (MFO), and cardiorespiratory fitness. Our results showed that TT allele participants expressed higher values of hunger (p = 0.049) and appetite (p = 0.043) after exercising compared to the AT allele group. Moreover, the TT allele group showed significantly higher values of MFO (p = 0.031) compared to the AT group, regardless of sex and body mass index. Thus, our results suggest that the FTO rs9939609 gene has an influence on appetite, hunger, and fat oxidation during exercise, with TT allele participants showing significantly higher values compared to the AT allele group. These findings may have practical applications for weight loss and exercise programs. [ABSTRACT FROM AUTHOR]

Background: In studies of time-to-events, it is common to collect information about events that occurred before the inclusion in a prospective cohort. When the studied risk factors are independent of time, including both pre- and post-inclusion events in the analyses, generally referred to as relying on an ambispective design, increases the statistical power but may lead to a selection bias. In the field of venous thromboembolism (VT), ABO blood groups have been the subject of extensive research due to their substantial effect on VT risk. However, few studies have investigated their effect on the risk of VT recurrence. Motivated by the study of the association of genetically determined ABO blood groups with VT recurrence, we propose a methodology to include pre-inclusion events in the analysis of ambispective studies while avoiding the selection bias due to mortality. Methods: This work relies on two independent cohorts of VT patients, the French MARTHA study built on an ambispective design and the Dutch MEGA study built on a standard prospective design. For the analysis of the MARTHA study, a weighted Cox model was developed where weights were defined by the inverse of the survival probability at the time of data collection about the events. Thanks to the collection of information on the vital status of patients, we could estimate the survival probabilities using a delayed-entry Cox model on the death risk. Finally, results obtained in both studies were then meta-analysed. Results: In the combined sample totalling 2,752 patients including 993 recurrences, the A1 blood group has an increased risk (Hazard Ratio (HR) of 1.18, p = 4.2 × 10–3) compared with the O1 group, homogeneously in MARTHA and in MEGA. The same trend (HR = 1.19, p = 0.06) was observed for the less frequent A2 group. Conclusion: The proposed methodology increases the power of studies relying on an ambispective design which is frequent in epidemiologic studies about recurrent events. This approach allowed to clarify the association of ABO blood groups with the risk of VT recurrence. Besides, this methodology has an immediate field of application in the context of genome wide association studies. [ABSTRACT FROM AUTHOR]

We introduce a method to draw causal inferences-inferences immune to all possible confounding-from genetic data that include parents and offspring. Causal conclusions are possible with these data because the natural randomness in meiosis can be viewed as a high-dimensional randomized experiment. We make this observation actionable by developing a conditional independence test that identifies regions of the genome containing distinct causal variants. The proposed digital twin test compares an observed offspring to carefully constructed synthetic offspring from the same parents to determine statistical significance, and it can leverage any black-box multivariate model and additional nontrio genetic data to increase power. Crucially, our inferences are based only on a well-established mathematical model of recombination and make no assumptions about the relationship between the genotypes and phenotypes. We compare our method to the widely used transmission disequilibrium test and demonstrate enhanced power and localization.

The five-protein MuvB core complex is highly conserved in animals. This nuclear complex interacts with RB-family tumor suppressor proteins and E2F-DP transcription factors to form DREAM complexes that repress genes that regulate cell cycle progression and cell fate. The MuvB core complex also interacts with Myb family oncoproteins to form the Myb-MuvB complexes that activate many of the same genes. We show that animal-type Myb genes are present in Bilateria, Cnidaria and Placozoa, the latter including the simplest known animal species. However, bilaterian nematode worms lost their animal-type Myb genes hundreds of millions of years ago. Nevertheless, amino acids in the LIN9 and LIN52 proteins that directly interact with the MuvB-binding domains of human B-Myb and Drosophila Myb are conserved in Caenorhabditiselegans Here, we show that, despite greater than 500 million years since their last common ancestor, the Drosophila melanogaster Myb protein can bind to the nematode LIN9-LIN52 proteins in vitro and can cause a synthetic multivulval (synMuv) phenotype in vivo This phenotype is similar to that caused by loss-of-function mutations in C. elegans synMuvB-class genes including those that encode homologs of the MuvB core, RB, E2F and DP. Furthermore, amino acid substitutions in the MuvB-binding domain of Drosophila Myb that disrupt its functions in vitro and in vivo also disrupt these activities in C. elegans We speculate that nematodes and other animals may contain another protein that can bind to LIN9 and LIN52 in order to activate transcription of genes repressed by DREAM complexes.

BACKGROUND: African ancestry is associated with a higher prevalence and greater severity of asthma than European ancestries, yet genetic studies of the most common locus associated with childhood-onset asthma, 17q12-21, in African Americans have been inconclusive. The aim of this study was to leverage both the phenotyping of the Childrens Respiratory and Environmental Workgroup (CREW) birth cohort consortium, and the reduced linkage disequilibrium in African Americans, to fine map the 17q12-21 locus. METHODS: We first did a genetic association study and meta-analysis using 17q12-21 tag single-nucleotide polymorphisms (SNPs) for childhood-onset asthma in 1613 European American and 870 African American children from the CREW consortium. Nine tag SNPs were selected based on linkage disequilibrium patterns at 17q12-21 and their association with asthma, considering the effect allele under an additive model (0, 1, or 2 effect alleles). Results were meta-analysed with publicly available summary data from the EVE consortium (on 4303 European American and 3034 African American individuals) for seven of the nine SNPs of interest. Subsequently, we tested for expression quantitative trait loci (eQTLs) among the SNPs associated with childhood-onset asthma and the expression of 17q12-21 genes in resting peripheral blood mononuclear cells (PBMCs) from 85 African American CREW children and in upper airway epithelial cells from 246 African American CREW children; and in lower airway epithelial cells from 44 European American and 72 African American adults from a case-control study of asthma genetic risk in Chicago (IL, USA). FINDINGS: 17q12-21 SNPs were broadly associated with asthma in European Americans. Only two SNPs (rs2305480 in gasdermin-B [GSDMB] and rs8076131 in ORMDL sphingolipid biosynthesis regulator 3 [ORMDL3]) were associated with asthma in African Americans, at a Bonferroni-corrected threshold of p

Shape is a critical element of the visual appeal of strawberry fruit and is influenced by both genetic and non-genetic determinants. Current fruit phenotyping approaches for external characteristics in strawberry often rely on the human eye to make categorical assessments. However, fruit shape is an inherently multi-dimensional, continuously variable trait and not adequately described by a single categorical or quantitative feature. Morphometric approaches enable the study of complex, multi-dimensional forms but are often abstract and difficult to interpret. In this study, we developed a mathematical approach for transforming fruit shape classifications from digital images onto an ordinal scale called the Principal Progression of k Clusters (PPKC). We use these human-recognizable shape categories to select quantitative features extracted from multiple morphometric analyses that are best fit for genetic dissection and analysis. We transformed images of strawberry fruit into human-recognizable categories using unsupervised machine learning, discovered 4 principal shape categories, and inferred progression using PPKC. We extracted 68 quantitative features from digital images of strawberries using a suite of morphometric analyses and multivariate statistical approaches. These analyses defined informative feature sets that effectively captured quantitative differences between shape classes. Classification accuracy ranged from 68% to 99% for the newly created phenotypic variables for describing a shape. Our results demonstrated that strawberry fruit shapes could be robustly quantified, accurately classified, and empirically ordered using image analyses, machine learning, and PPKC. We generated a dictionary of quantitative traits for studying and predicting shape classes and identifying genetic factors underlying phenotypic variability for fruit shape in strawberry. The methods and approaches that we applied in strawberry should apply to other fruits, vegetables, and specialty crops.

Transcriptomics has become an important tool for identification of biological pathways dysregulated in Alzheimer's disease (AD). We performed a network-based gene expression analysis of blood-based microarray gene expression profiles using 2 independent cohorts, Alzheimer's Disease Neuroimaging Initiative (ADNI; N = 661) and AddNeuroMed (N = 674). Weighted gene coexpression network analysis identified 17 modules from ADNI and 13 from AddNeuroMed. Four of the modules derived in ADNI were significantly related to AD; 5 modules in AddNeuroMed were significant. Gene-set enrichment analysis of the AD-related modules identified and replicated 3 biological pathways including the Fc gamma receptor-mediated phagocytosis pathway. Module-based association analysis showed the AD-related module, which has the 3 pathways, to be associated with cognitive function and neuroimaging biomarkers. Gene-based association analysis identified PRKCD in the Fc gamma receptor-mediated phagocytosis pathway as being significantly associated with cognitive function and cerebrospinal fluid biomarkers. The identification of the Fc gamma receptor-mediated phagocytosis pathway implicates the peripheral innate immune system in the pathophysiology of AD. PRKCD is known to be related to neurodegeneration induced by amyloid-β.

Next-generation sequencing has revolutionized rare disease diagnostics, but many patients remain without a molecular diagnosis, particularly because many candidate variants usually survive despite strict filtering. Exomiser was launched in 2014 as a Java tool that performs an integrative analysis of patients' sequencing data and their phenotypes encoded with Human Phenotype Ontology (HPO) terms. It prioritizes variants by leveraging information on variant frequency, predicted pathogenicity, and gene-phenotype associations derived from human diseases, model organisms, and protein-protein interactions. Early published releases of Exomiser were able to prioritize disease-causative variants as top candidates in up to 97% of simulated whole-exomes. The size of the tested real patient datasets published so far are very limited. Here, we present the latest Exomiser version 12.0.1 with many new features. We assessed the performance using a set of 134 whole-exomes from patients with a range of rare retinal diseases and known molecular diagnosis. Using default settings, Exomiser ranked the correct diagnosed variants as the top candidate in 74% of the dataset and top 5 in 94%; not using the patients' HPO profiles (i.e., variant-only analysis) decreased the performance to 3% and 27%, respectively. In conclusion, Exomiser is an effective support tool for rare Mendelian phenotype-driven variant prioritization.

"Quantile-dependent expressivity" is a dependence of genetic effects on whether the phenotype (e.g., triglycerides) is high or low relative to its distribution in the population. Quantile-specific offspring-parent regression slopes (βOP) were estimated by quantile regression for 6227 offspring-parent pairs. Quantile-specific heritability (h2), estimated by 2βOP/(1 + rspouse), decreased 0.0047 ± 0.0007 (P = 2.9 × 10-14) for each one-percent decrement in fasting triglyceride concentrations, i.e., h2 ± SE were: 0.428 ± 0.059, 0.230 ± 0.030, 0.111 ± 0.015, 0.050 ± 0.016, and 0.033 ± 0.010 at the 90th, 75th, 50th, 25th, and 10th percentiles of the triglyceride distribution, respectively. Consistent with quantile-dependent expressivity, 11 drug studies report smaller genotype differences at lower (post-treatment) than higher (pre-treatment) triglyceride concentrations. This meant genotype-specific triglyceride changes could not move in parallel when triglycerides were decreased pharmacologically, so that subtracting pre-treatment from post-treatment triglyceride levels necessarily created a greater triglyceride decrease for the genotype with a higher pre-treatment value (purported precision-medicine genetic markers). In addition, sixty-five purported gene-environment interactions were found to be potentially attributable to triglyceride's quantile-dependent expressivity, including gene-adiposity (APOA5, APOB, APOE, GCKR, IRS-1, LPL, MTHFR, PCSK9, PNPLA3, PPARγ2), gene-exercise (APOA1, APOA2, LPL), gene-diet (APOA5, APOE, INSIG2, LPL, MYB, NXPH1, PER2, TNFA), gene-alcohol (ALDH2, APOA5, APOC3, CETP, LPL), gene-smoking (APOC3, CYBA, LPL, USF1), gene-pregnancy (LPL), and gene-insulin resistance interactions (APOE, LPL).

We describe an 11-year old boy with severe global developmental delays, failure to thrive and growth retardation, refractory seizures with recurrent status epilepticus, hypogammaglobulinemia, hypergonadotropic hypogonadism, and duodenal strictures. He had facial and skin findings compatible with trichothiodystrophy, including sparse and brittle hair, thin eyebrows, and dry skin. Exome sequencing showed a hemizygous, truncating variant in RNF113A, c.903_910delGCAGACCA, predicting p.(Gln302fs*12), that was inherited from his mother. Although his clinical features overlap closely with features described in the two previously reported male first cousins with RNF113A loss of function mutations, the duodenal strictures seen in this patient have not been reported. Interestingly, the patient's mother had short stature and 100% skewed X-inactivation as seen in other obligate female carriers. A second male with developmental delays, microcephaly, seizures, ambiguous genitalia, and facial anomalies that included sparse and brittle hair, thin eyebrows and dry skin was recently reported to have c.897_898delTG, predicting p.(Cys299*) in RNF113A and we provide additional clinical details for this patient. This report further supports deleterious variants in RNF113A as a cause of a novel trichothiodystrophy syndrome.

Knowledge of the genetic architecture of importantly agronomical traits can speed up genetic improvement in cultivated rice (Oryza sativa L.). Many recent investigations have leveraged genome-wide association studies (GWAS) to identify single nucleotide polymorphisms (SNPs), associated with agronomic traits in various rice populations. The reported trait-relevant SNPs appear to be arbitrarily distributed along the genome, including genic and nongenic regions. Whether the SNPs in different genomic regions play different roles in trait heritability and which region is more responsible for phenotypic variation remains opaque. We analyzed a natural rice population of 524 accessions with 3,616,597 SNPs to compare the genetic contributions of functionally distinct genomic regions for five agronomic traits, i.e., yield, heading date, plant height, grain length, and grain width. An analysis of heritability in the functionally partitioned rice genome showed that regulatory or intergenic regions account for the most trait heritability. A close look at the trait-associated SNPs (TASs) indicated that the majority of the TASs are located in nongenic regions, and the genetic effects of the TASs in nongenic regions are generally greater than those in genic regions. We further compared the predictabilities using the genetic variants from genic regions with those using nongenic regions. The results revealed that nongenic regions play a more important role than genic regions in trait heritability in rice, which is consistent with findings in humans and maize. This conclusion not only offers clues for basic research to disclose genetics behind these agronomic traits, but also provides a new perspective to facilitate genomic selection in rice.

MotivationHigh-throughput phenomic projects generate complex data from small treatment and large control groups that increase the power of the analyses but introduce variation over time. A method is needed to utlize a set of temporally local controls that maximizes analytic power while minimizing noise from unspecified environmental factors.ResultsHere we introduce 'soft windowing', a methodological approach that selects a window of time that includes the most appropriate controls for analysis. Using phenotype data from the International Mouse Phenotyping Consortium (IMPC), adaptive windows were applied such that control data collected proximally to mutants were assigned the maximal weight, while data collected earlier or later had less weight. We applied this method to IMPC data and compared the results with those obtained from a standard non-windowed approach. Validation was performed using a resampling approach in which we demonstrate a 10% reduction of false positives from 2.5 million analyses. We applied the method to our production analysis pipeline that establishes genotype-phenotype associations by comparing mutant versus control data. We report an increase of 30% in significant P-values, as well as linkage to 106 versus 99 disease models via phenotype overlap with the soft-windowed and non-windowed approaches, respectively, from a set of 2082 mutant mouse lines. Our method is generalizable and can benefit large-scale human phenomic projects such as the UK Biobank and the All of Us resources.Availability and implementationThe method is freely available in the R package SmoothWin, available on CRAN http://CRAN.R-project.org/package=SmoothWin.Supplementary informationSupplementary data are available at Bioinformatics online.

The bHLH transcription factor Hes1 is a key downstream effector for the Notch signaling pathway. During embryogenesis neural progenitors express low levels of Hes1 in an oscillating pattern, whereas glial brain boundary regions (e.g., isthmus) have high, sustained Hes1 levels that suppress neuronal fates. Here, we show that in the embryonic mouse retina, the optic nerve head and stalk express high Hes1, with the ONH constituting a boundary between the neural retina and glial cells that ultimately line the optic stalk. Using two Cre drivers with distinct spatiotemporal expression we conditionally inactivated Hes1, to delineate the requirements for this transcriptional repressor during retinal neurogenesis versus patterning of the optic cup and stalk. Throughout retinal neurogenesis, Hes1 maintains proliferation and blocks retinal ganglion cell formation, but surprisingly we found it also promotes cone photoreceptor genesis. In the postnatal eye, Hes1 inactivation with Rax-Cre resulted in increased bipolar neurons and a mispositioning of Müller glia. Our results indicate that Notch pathway regulation of cone genesis is more complex than previously assumed, and reveal a novel role for Hes1 in maintaining the optic cup-stalk boundary.SIGNIFICANCE STATEMENT The bHLH repressor Hes1 regulates the timing of neurogenesis, rate of progenitor cell division, gliogenesis, and maintains tissue compartment boundaries. This study expands current eye development models by showing Notch-independent roles for Hes1 in the developing optic nerve head (ONH). Defects in ONH formation result in optic nerve coloboma; our work now inserts Hes1 into the genetic hierarchy regulating optic fissure closure. Given that Hes1 acts analogously in the ONH as the brain isthmus, it prompts future investigation of the ONH as a signaling factor center, or local organizer. Embryonic development of the ONH region has been poorly studied, which is surprising given it is where the pan-ocular disease glaucoma is widely believed to inflict damage on RGC axons.

Each additional copy of the apolipoprotein E4 (APOE4) allele is associated with a higher risk of Alzheimer's dementia, while the APOE2 allele is associated with a lower risk of Alzheimer's dementia, it is not yet known whether APOE2 homozygotes have a particularly low risk. We generated Alzheimer's dementia odds ratios and other findings in more than 5,000 clinically characterized and neuropathologically characterized Alzheimer's dementia cases and controls. APOE2/2 was associated with a low Alzheimer's dementia odds ratios compared to APOE2/3 and 3/3, and an exceptionally low odds ratio compared to APOE4/4, and the impact of APOE2 and APOE4 gene dose was significantly greater in the neuropathologically confirmed group than in more than 24,000 neuropathologically unconfirmed cases and controls. Finding and targeting the factors by which APOE and its variants influence Alzheimer's disease could have a major impact on the understanding, treatment and prevention of the disease.