Electronic health record signatures identify undiagnosed patients with common variable immunodeficiency disease

Human inborn errors of immunity include rare disorders entailing functional and quantitative antibody deficiencies due to impaired B cells called the common variable immunodeficiency (CVID) phenotype. Patients with CVID face delayed diagnoses and treatments for 5 to 15 years after symptom onset because the disorders are rare (prevalence of ~1/25,000), and there is extensive heterogeneity in CVID phenotypes, ranging from infections to autoimmunity to inflammatory conditions, overlapping with other more common disorders. The prolonged diagnostic odyssey drives excessive system-wide costs before diagnosis. Because there is no single causal mechanism, there are no genetic tests to definitively diagnose CVID. Here, we present PheNet, a machine learning algorithm that identifies patients with CVID from their electronic health records (EHRs). PheNet learns phenotypic patterns from verified CVID cases and uses this knowledge to rank patients by likelihood of having CVID. PheNet could have diagnosed more than half of our patients with CVID 1 or more years earlier than they had been diagnosed. When applied to a large EHR dataset, followed by blinded chart review of the top 100 patients ranked by PheNet, we found that 74% were highly probable to have CVID. We externally validated PheNet using >6 million records from disparate medical systems in California and Tennessee. As artificial intelligence and machine learning make their way into health care, we show that algorithms such as PheNet can offer clinical benefits by expediting the diagnosis of rare diseases.

The omics era: a nexus of untapped potential for Mendelian chromatinopathies

The OMICs cascade describes the hierarchical flow of information through biological systems. The epigenome sits at the apex of the cascade, thereby regulating the RNA and protein expression of the human genome and governs cellular identity and function. Genes that regulate the epigenome, termed epigenes, orchestrate complex biological signaling programs that drive human development. The broad expression patterns of epigenes during human development mean that pathogenic germline mutations in epigenes can lead to clinically significant multi-system malformations, developmental delay, intellectual disabilities, and stem cell dysfunction. In this review, we refer to germline developmental disorders caused by epigene mutation as "chromatinopathies". We curated the largest number of human chromatinopathies to date and our expanded approach more than doubled the number of established chromatinopathies to 179 disorders caused by 148 epigenes. Our study revealed that 20.6% (148/720) of epigenes cause at least one chromatinopathy. In this review, we highlight key examples in which OMICs approaches have been applied to chromatinopathy patient biospecimens to identify underlying disease pathogenesis. The rapidly evolving OMICs technologies that couple molecular biology with high-throughput sequencing or proteomics allow us to dissect out the causal mechanisms driving temporal-, cellular-, and tissue-specific expression. Using the full repertoire of data generated by the OMICs cascade to study chromatinopathies will provide invaluable insight into the developmental impact of these epigenes and point toward future precision targets for these rare disorders.

Increased Frequency of Clonal Hematopoiesis of Indeterminate Potential in Bloom Syndrome Probands and Carriers

Background

Bloom Syndrome (BSyn) is an autosomal recessive disorder caused by biallelic germline variants in BLM, which functions to maintain genomic stability. BSyn patients have poor growth, immune defects, insulin resistance, and a significantly increased risk of malignancies, most commonly hematologic. The malignancy risk in carriers of pathogenic variants in BLM (BLM variant carriers) remains understudied. Clonal hematopoiesis of indeterminate potential (CHIP) is defined by presence of somatic mutations in leukemia-related genes in blood of individuals without leukemia and is associated with increased risk of leukemia. We hypothesize that somatic mutations driving clonal expansion may be an underlying mechanism leading to increased cancer risk in BSyn patients and BLM variant carriers.

Methods

To determine whether de novo or somatic variation is increased in BSyn patients or carriers, we performed and analyzed exome sequencing on BSyn and control trios.

Results

We discovered that both BSyn patients and carriers had increased numbers of low-frequency, putative somatic variants in CHIP genes compared to controls. Furthermore, BLM variant carriers had increased numbers of somatic variants in DNA methylation genes compared to controls. There was no statistical difference in the numbers of de novo variants in BSyn probands compared to control probands.

Conclusion

Our findings of increased CHIP in BSyn probands and carriers suggest that one or two germline pathogenic variants in BLM could be sufficient to increase the risk of clonal hematopoiesis. These findings warrant further studies in larger cohorts to determine the significance of CHIP as a potential biomarker of aging, cancer, cardiovascular disease, morbidity and mortality.

KAT6A mutations in Arboleda-Tham syndrome drive epigenetic regulation of posterior HOXC cluster

Arboleda-Tham Syndrome (ARTHS) is a rare genetic disorder caused by heterozygous, de novo mutations in Lysine(K) acetyltransferase 6A (KAT6A). ARTHS is clinically heterogeneous and characterized by several common features, including intellectual disability, developmental and speech delay, and hypotonia, and affects multiple organ systems. KAT6A is the enzymatic core of a histone-acetylation protein complex; however, the direct histone targets and gene regulatory effects remain unknown. In this study, we use ARTHS patient (n = 8) and control (n = 14) dermal fibroblasts and perform comprehensive profiling of the epigenome and transcriptome caused by KAT6A mutations. We identified differential chromatin accessibility within the promoter or gene body of 23% (14/60) of genes that were differentially expressed between ARTHS and controls. Within fibroblasts, we show a distinct set of genes from the posterior HOXC gene cluster (HOXC10, HOXC11, HOXC-AS3, HOXC-AS2, and HOTAIR) that are overexpressed in ARTHS and are transcription factors critical for early development body segment patterning. The genomic loci harboring HOXC genes are epigenetically regulated with increased chromatin accessibility, high levels of H3K23ac, and increased gene-body DNA methylation compared to controls, all of which are consistent with transcriptomic overexpression. Finally, we used unbiased proteomic mass spectrometry and identified two new histone post-translational modifications (PTMs) that are disrupted in ARTHS: H2A and H3K56 acetylation. Our multi-omics assays have identified novel histone and gene regulatory roles of KAT6A in a large group of ARTHS patients harboring diverse pathogenic mutations. This work provides insight into the role of KAT6A on the epigenomic regulation in somatic cell types.

3D Printing as an Effective Quality Assurance Implementation in Massive-Scale SARS-CoV-2 Testing at a SwabSeq Next-Generation Sequencing Laboratory

Massive-scale SARS-CoV-2 testing using the SwabSeq diagnostic platform came with quality assurance challenges due to the novelty and scale of sequencing-based testing. The SwabSeq platform relies on accurate mapping between specimen identifiers and molecular barcodes to match a result back to a patient specimen. To identify and mitigate mapping errors, we instituted quality control using placement of negative controls within a rack of patient samples. We designed 2-dimensional paper templates to fit over a 96-position rack of specimens with holes to show the control tube placements. We designed and 3-dimensionally printed plastic templates that fit onto 4 racks of patient specimens and provide accurate indications of the correct control tube placements. The final plastic templates dramatically reduced plate mapping errors from 22.55% in January 2021 to less than 1% after implementation and training in January 2021. We show how 3D printing can be a cost-effective quality assurance tool to mitigate human error in the clinical laboratory.

KAT6A mutations in Arboleda-Tham syndrome drive epigenetic regulation of posterior HOXC cluster

Arboleda-Tham Syndrome (ARTHS) is a rare genetic disorder caused by heterozygous, de novo truncating mutations in Lysine(K) acetyltransferase 6A (KAT6A). ARTHS is clinically heterogeneous and characterized by several common features including intellectual disability, developmental and speech delay, hypotonia and affects multiple organ systems. KAT6A is highly expressed in early development and plays a key role in cell-type specific differentiation. KAT6A is the enzymatic core of a histone-acetylation protein complex, however the direct histone targets and gene regulatory effects remain unknown. In this study, we use ARTHS patient (n=8) and control (n=14) dermal fibroblasts and perform comprehensive profiling of the epigenome and transcriptome caused by KAT6A mutations. We identified differential chromatin accessibility within the promoter or gene body of 23%(14/60) of genes that were differentially expressed between ARTHS and controls. Within fibroblasts, we show a distinct set of genes from the posterior HOXC gene cluster (HOXC10, HOXC11, HOXC-AS3, HOXC-AS2, HOTAIR) that are overexpressed in ARTHS and are transcription factors critical for early development body segment patterning. The genomic loci harboring HOXC genes are epigenetically regulated with increased chromatin accessibility, high levels of H3K23ac, and increased gene-body DNA methylation compared to controls, all of which are consistent with transcriptomic overexpression. Finally, we used unbiased proteomic mass spectrometry and identified two new histone post-translational modifications (PTMs) that are disrupted in ARTHS: H2A and H3K56 acetylation. Our multi-omics assays have identified novel histone and gene regulatory roles of KAT6A in a large group of ARTHS patients harboring diverse pathogenic mutations. This work provides insight into the role of KAT6A on the epigenomic regulation in somatic cell types.

Disease risk and healthcare utilization among ancestrally diverse groups in the Los Angeles region

An individual's disease risk is affected by the populations that they belong to, due to shared genetics and environmental factors. The study of fine-scale populations in clinical care is important for identifying and reducing health disparities and for developing personalized interventions. To assess patterns of clinical diagnoses and healthcare utilization by fine-scale populations, we leveraged genetic data and electronic medical records from 35,968 patients as part of the UCLA ATLAS Community Health Initiative. We defined clusters of individuals using identity by descent, a form of genetic relatedness that utilizes shared genomic segments arising due to a common ancestor. In total, we identified 376 clusters, including clusters with patients of Afro-Caribbean, Puerto Rican, Lebanese Christian, Iranian Jewish and Gujarati ancestry. Our analysis uncovered 1,218 significant associations between disease diagnoses and clusters and 124 significant associations with specialty visits. We also examined the distribution of pathogenic alleles and found 189 significant alleles at elevated frequency in particular clusters, including many that are not regularly included in population screening efforts. Overall, this work progresses the understanding of health in understudied communities and can provide the foundation for further study into health inequities.

Whole-Exome Sequencing Identifies Homozygote Nonsense Variants in LMOD2 Gene Causing Infantile Dilated Cardiomyopathy

As an essential component of the sarcomere, actin thin filament stems from the Z-disk extend toward the middle of the sarcomere and overlaps with myosin thick filaments. Elongation of the cardiac thin filament is essential for normal sarcomere maturation and heart function. This process is regulated by the actin-binding proteins Leiomodins (LMODs), among which LMOD2 has recently been identified as a key regulator of thin filament elongation to reach a mature length. Few reports have implicated homozygous loss of function variants of LMOD2 in neonatal dilated cardiomyopathy (DCM) associated with thin filament shortening. We present the fifth case of DCM due to biallelic variants in the LMOD2 gene and the second case with the c.1193G>A (p.W398*) nonsense variant identified by whole-exome sequencing. The proband is a 4-month male infant of Hispanic descent with advanced heart failure. Consistent with previous reports, a myocardial biopsy exhibited remarkably short thin filaments. However, compared to other cases of identical or similar biallelic variants, the patient presented here has an unusually late onset of cardiomyopathy during infancy. Herein, we present the phenotypic and histological features of this variant, confirm the pathogenic impact on protein expression and sarcomere structure, and discuss the current knowledge of LMOD2-related cardiomyopathy.

Multiomics of Bohring-Opitz syndrome truncating ASXL1 mutations identify canonical and noncanonical Wnt signaling dysregulation

ASXL1 (additional sex combs-like 1) plays key roles in epigenetic regulation of early developmental gene expression. De novo protein-truncating mutations in ASXL1 cause Bohring-Opitz syndrome (BOS; OMIM #605039), a rare neurodevelopmental condition characterized by severe intellectual disabilities, distinctive facial features, hypertrichosis, increased risk of Wilms tumor, and variable congenital anomalies, including heart defects and severe skeletal defects giving rise to a typical BOS posture. These BOS-causing ASXL1 variants are also high-prevalence somatic driver mutations in acute myeloid leukemia. We used primary cells from individuals with BOS (n = 18) and controls (n = 49) to dissect gene regulatory changes caused by ASXL1 mutations using comprehensive multiomics assays for chromatin accessibility (ATAC-seq), DNA methylation, histone methylation binding, and transcriptome in peripheral blood and skin fibroblasts. Our data show that regardless of cell type, ASXL1 mutations drive strong cross-tissue effects that disrupt multiple layers of the epigenome. The data showed a broad activation of canonical Wnt signaling at the transcriptional and protein levels and upregulation of VANGL2, which encodes a planar cell polarity pathway protein that acts through noncanonical Wnt signaling to direct tissue patterning and cell migration. This multiomics approach identifies the core impact of ASXL1 mutations and therapeutic targets for BOS and myeloid leukemias.

Leveraging genomic diversity for discovery in an electronic health record linked biobank: the UCLA ATLAS Community Health Initiative

BACKGROUND

Large medical centers in urban areas, like Los Angeles, care for a diverse patient population and offer the potential to study the interplay between genetic ancestry and social determinants of health. Here, we explore the implications of genetic ancestry within the University of California, Los Angeles (UCLA) ATLAS Community Health Initiative-an ancestrally diverse biobank of genomic data linked with de-identified electronic health records (EHRs) of UCLA Health patients (N=36,736).

METHODS

We quantify the extensive continental and subcontinental genetic diversity within the ATLAS data through principal component analysis, identity-by-descent, and genetic admixture. We assess the relationship between genetically inferred ancestry (GIA) and >1500 EHR-derived phenotypes (phecodes). Finally, we demonstrate the utility of genetic data linked with EHR to perform ancestry-specific and multi-ancestry genome and phenome-wide scans across a broad set of disease phenotypes.

RESULTS

We identify 5 continental-scale GIA clusters including European American (EA), African American (AA), Hispanic Latino American (HL), South Asian American (SAA) and East Asian American (EAA) individuals and 7 subcontinental GIA clusters within the EAA GIA corresponding to Chinese American, Vietnamese American, and Japanese American individuals. Although we broadly find that self-identified race/ethnicity (SIRE) is highly correlated with GIA, we still observe marked differences between the two, emphasizing that the populations defined by these two criteria are not analogous. We find a total of 259 significant associations between continental GIA and phecodes even after accounting for individuals' SIRE, demonstrating that for some phenotypes, GIA provides information not already captured by SIRE. GWAS identifies significant associations for liver disease in the 22q13.31 locus across the HL and EAA GIA groups (HL p-value=2.32×10-16, EAA p-value=6.73×10-11). A subsequent PheWAS at the top SNP reveals significant associations with neurologic and neoplastic phenotypes specifically within the HL GIA group.

CONCLUSIONS

Overall, our results explore the interplay between SIRE and GIA within a disease context and underscore the utility of studying the genomes of diverse individuals through biobank-scale genotyping linked with EHR-based phenotyping.

DMRscaler: a scale-aware method to identify regions of differential DNA methylation spanning basepair to multi-megabase features

BACKGROUND

Pathogenic mutations in genes that control chromatin function have been implicated in rare genetic syndromes. These chromatin modifiers exhibit extraordinary diversity in the scale of the epigenetic changes they affect, from single basepair modifications by DNMT1 to whole genome structural changes by PRM1/2. Patterns of DNA methylation are related to a diverse set of epigenetic features across this full range of epigenetic scale, making DNA methylation valuable for mapping regions of general epigenetic dysregulation. However, existing methods are unable to accurately identify regions of differential methylation across this full range of epigenetic scale directly from DNA methylation data.

RESULTS

To address this, we developed DMRscaler, a novel method that uses an iterative windowing procedure to capture regions of differential DNA methylation (DMRs) ranging in size from single basepairs to whole chromosomes. We benchmarked DMRscaler against several DMR callers in simulated and natural data comparing XX and XY peripheral blood samples. DMRscaler was the only method that accurately called DMRs ranging in size from 100 bp to 1 Mb (pearson's r = 0.94) and up to 152 Mb on the X-chromosome. We then analyzed methylation data from rare-disease cohorts that harbor chromatin modifier gene mutations in NSD1, EZH2, and KAT6A where DMRscaler identified novel DMRs spanning gene clusters involved in development.

CONCLUSION

Taken together, our results show DMRscaler is uniquely able to capture the size of DMR features across the full range of epigenetic scale and identify novel, co-regulated regions that drive epigenetic dysregulation in human disease.

Polygenic risk scores of endo-phenotypes identify the effect of genetic background in congenital heart disease

Congenital heart disease (CHD) is a rare structural defect that occurs in ∼1% of live births. Studies on CHD genetic architecture have identified pathogenic single-gene mutations in less than 30% of cases. Single-gene mutations often show incomplete penetrance and variable expressivity. Therefore, we hypothesize that genetic background may play a role in modulating disease expression. Polygenic risk scores (PRSs) aggregate effects of common genetic variants to investigate whether, cumulatively, these variants are associated with disease penetrance or severity. However, the major limitations in this field have been in generating sufficient sample sizes for these studies. Here we used CHD-phenotype matched genome-wide association study (GWAS) summary statistics from the UK Biobank (UKBB) as our base study and whole-genome sequencing data from the CHD cohort (n₁ = 711 trios, n₂ = 362 European trios) of the Gabriella Miller Kids First dataset as our target study to develop PRSs for CHD. PRSs estimated using a GWAS for heart valve problems and heart murmur explain 2.5% of the variance in case-control status of CHD (all SNVs, p = 7.90 × 10^-3; fetal cardiac SNVs, p = 8.00 × 10^-3) and 1.8% of the variance in severity of CHD (fetal cardiac SNVs, p = 6.20 × 10^-3; all SNVs, p = 0.015). These results show that common variants captured in CHD phenotype-matched GWASs have a modest but significant contribution to phenotypic expression of CHD. Further exploration of the cumulative effect of common variants is necessary for understanding the complex genetic etiology of CHD and other rare diseases.

DNA methylation signature associated with Bohring-Opitz syndrome: a new tool for functional classification of variants in ASXL genes

The additional sex combs-like (ASXL) gene family-encoded by ASXL1, ASXL2, and ASXL3-is crucial for mammalian development. Pathogenic variants in the ASXL gene family are associated with three phenotypically distinct neurodevelopmental syndromes. Our previous work has shown that syndromic conditions caused by pathogenic variants in epigenetic regulatory genes show consistent patterns of genome-wide DNA methylation (DNAm) alterations, i.e., DNAm signatures in peripheral blood. Given the role of ASXL1 in chromatin modification, we hypothesized that pathogenic ASXL1 variants underlying Bohring-Opitz syndrome (BOS) have a unique DNAm signature. We profiled whole-blood DNAm for 17 ASXL1 variants, and 35 sex- and age-matched typically developing individuals, using Illumina's Infinium EPIC array. We identified 763 differentially methylated CpG sites in individuals with BOS. Differentially methylated sites overlapped 323 unique genes, including HOXA5 and HOXB4, supporting the functional relevance of DNAm signatures. We used a machine-learning classification model based on the BOS DNAm signature to classify variants of uncertain significance in ASXL1, as well as pathogenic ASXL2 and ASXL3 variants. The DNAm profile of one individual with the ASXL2 variant was BOS-like, whereas the DNAm profiles of three individuals with ASXL3 variants were control-like. We also used Horvath's epigenetic clock, which showed acceleration in DNAm age in individuals with pathogenic ASXL1 variants, and the individual with the pathogenic ASXL2 variant, but not in individuals with ASXL3 variants. These studies enhance our understanding of the epigenetic dysregulation underpinning ASXL gene family-associated syndromes.

Balancing the transcriptome: leveraging sample similarity to improve measures of gene specificity

The spatial and temporal domain of a gene's expression can range from ubiquitous to highly specific. Quantifying the degree to which this expression is unique to a specific tissue or developmental timepoint can provide insight into the etiology of genetic diseases. However, quantifying specificity remains challenging as measures of specificity are sensitive to similarity between samples in the sample set. For example, in the Gene-Tissue Expression project (GTEx), brain subregions are overrepresented at 13 of 54 (24%) unique tissues sampled. In this dataset, existing specificity measures have a decreased ability to identify genes specific to the brain relative to other organs. To solve this problem, we leverage sample similarity information to weight samples such that overrepresented tissues do not have an outsized effect on specificity estimates. We test this reweighting procedure on 4 measures of specificity, Z-score, Tau, Tsi and Gini, in the GTEx data and in single cell datasets for zebrafish and mouse. For all of these measures, incorporating sample similarity information to weight samples results in greater stability of sets of genes called as specific and decreases the overall variance in the change of specificity estimates as sample sets become more unbalanced. Furthermore, the genes with the largest improvement in their specificity estimate's stability are those with functions related to the overrepresented sample types. Our results demonstrate that incorporating similarity information improves specificity estimates' stability to the choice of the sample set used to define the transcriptome, providing more robust and reproducible measures of specificity for downstream analyses.

Genomic epidemiology of the Los Angeles COVID-19 outbreak and the early history of the B.1.43 strain in the USA

Background

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused global disruption of human health and activity. Being able to trace the early outbreak of SARS-CoV-2 within a locality can inform public health measures and provide insights to contain or prevent viral transmission. Investigation of the transmission history requires efficient sequencing methods and analytic strategies, which can be generally useful in the study of viral outbreaks.

Methods

The County of Los Angeles (hereafter, LA County) sustained a large outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To learn about the transmission history, we carried out surveillance viral genome sequencing to determine 142 viral genomes from unique patients seeking care at the University of California, Los Angeles (UCLA) Health System. 86 of these genomes were from samples collected before April 19, 2020.

Results

We found that the early outbreak in LA County, as in other international air travel hubs, was seeded by multiple introductions of strains from Asia and Europe. We identified a USA-specific strain, B.1.43, which was found predominantly in California and Washington State. While samples from LA County carried the ancestral B.1.43 genome, viral genomes from neighboring counties in California and from counties in Washington State carried additional mutations, suggesting a potential origin of B.1.43 in Southern California. We quantified the transmission rate of SARS-CoV-2 over time, and found evidence that the public health measures put in place in LA County to control the virus were effective at preventing transmission, but might have been undermined by the many introductions of SARS-CoV-2 into the region.

Conclusion

Our work demonstrates that genome sequencing can be a powerful tool for investigating outbreaks and informing the public health response. Our results reinforce the critical need for the USA to have coordinated inter-state responses to the pandemic.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08488-7.

Lean Principles to Improve Quality in High-Throughput COVID-19 Testing Using SwabSeq: A Barcoded Sequencing-Based Testing Platform

OBJECTIVE

To describe and quantify the effect of quality control (QC) metrics to increase testing efficiency in a high-complexity, Clinical Laboratory Improvement Amendments-certified laboratory that uses amplicon-based, next generation sequencing for the clinical detection of SARS-CoV-2. To enable rapid scalability to several thousands of specimens per day without fully automated platforms, we developed internal QC methods to ensure high-accuracy testing.

METHODS

We implemented procedures to increase efficiency by applying the Lean Six Sigma model into our sequencing-based COVID-19 detection.

RESULTS

The application of the Lean Six Sigma model increased laboratory efficiency by reducing errors, allowing for a higher testing volume to be met with minimal staffing. Furthermore, these improvements resulted in an improved turnaround time.

CONCLUSION

Lean Six Sigma model execution has increased laboratory efficiency by decreasing critical testing errors and has prepared the laboratory for future scaling up to 50,000 tests per day.

Lower Severe Acute Respiratory Syndrome Coronavirus 2 Viral Shedding Following Coronavirus Disease 2019 Vaccination Among Healthcare Workers in Los Angeles, California

Among 880 healthcare workers with a positive severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) test, 264 (30.0%) infections were identified following receipt of at least 1 vaccine dose. Median SARS-CoV-2 cycle threshold values were highest among individuals receiving 2 vaccine doses, corresponding to lower viral shedding. Vaccination might lead to lower transmissibility of SARS-CoV-2.

Novel variants in KAT6B spectrum of disorders expand our knowledge of clinical manifestations and molecular mechanisms

The phenotypic variability associated with pathogenic variants in Lysine Acetyltransferase 6B (KAT6B, a.k.a. MORF, MYST4) results in several interrelated syndromes including Say-Barber-Biesecker-Young-Simpson Syndrome and Genitopatellar Syndrome. Here we present 20 new cases representing 10 novel KAT6B variants. These patients exhibit a range of clinical phenotypes including intellectual disability, mobility and language difficulties, craniofacial dysmorphology, and skeletal anomalies. Given the range of features previously described for KAT6B-related syndromes, we have identified additional phenotypes including concern for keratoconus, sensitivity to light or noise, recurring infections, and fractures in greater numbers than previously reported. We surveyed clinicians to qualitatively assess the ways families engage with genetic counselors upon diagnosis. We found that 56% (10/18) of individuals receive diagnoses before the age of 2 years (median age = 1.96 years), making it challenging to address future complications with limited accessible information and vast phenotypic severity. We used CRISPR to introduce truncating variants into the KAT6B gene in model cell lines and performed chromatin accessibility and transcriptome sequencing to identify key dysregulated pathways. This study expands the clinical spectrum and addresses the challenges to management and genetic counseling for patients with KAT6B-related disorders.

Massively scaled-up testing for SARS-CoV-2 RNA via next-generation sequencing of pooled and barcoded nasal and saliva samples

Frequent and widespread testing of members of the population who are asymptomatic for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is essential for the mitigation of the transmission of the virus. Despite the recent increases in testing capacity, tests based on quantitative polymerase chain reaction (qPCR) assays cannot be easily deployed at the scale required for population-wide screening. Here, we show that next-generation sequencing of pooled samples tagged with sample-specific molecular barcodes enables the testing of thousands of nasal or saliva samples for SARS-CoV-2 RNA in a single run without the need for RNA extraction. The assay, which we named SwabSeq, incorporates a synthetic RNA standard that facilitates end-point quantification and the calling of true negatives, and that reduces the requirements for automation, purification and sample-to-sample normalization. We used SwabSeq to perform 80,000 tests, with an analytical sensitivity and specificity comparable to or better than traditional qPCR tests, in less than two months with turnaround times of less than 24 h. SwabSeq could be rapidly adapted for the detection of other pathogens.

Retrospective Detection of SARS-CoV-2 in Symptomatic Patients prior to Widespread Diagnostic Testing in Southern California

Background

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused one of the worst pandemics in recent history. Few reports have revealed that SARS-CoV-2 was spreading in the United States as early as the end of January. In this study, we aimed to determine if SARS-CoV-2 had been circulating in the Los Angeles (LA) area at a time when access to diagnostic testing for coronavirus disease 2019 (COVID-19) was severely limited.

Methods

We used a pooling strategy to look for SARS-CoV-2 in remnant respiratory samples submitted for regular respiratory pathogen testing from symptomatic patients from November 2019 to early March 2020. We then performed sequencing on the positive samples.

Results

We detected SARS-CoV-2 in 7 specimens from 6 patients, dating back to mid-January. The earliest positive patient, with a sample collected on January 13, 2020 had no relevant travel history but did have a sibling with similar symptoms. Sequencing of these SARS-CoV-2 genomes revealed that the virus was introduced into the LA area from both domestic and international sources as early as January.

Conclusions

We present strong evidence of community spread of SARS-CoV-2 in the LA area well before widespread diagnostic testing was being performed in early 2020. These genomic data demonstrate that SARS-CoV-2 was being introduced into Los Angeles County from both international and domestic sources in January 2020.

Swab-Seq: A high-throughput platform for massively scaled up SARS-CoV-2 testing

The rapid spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is due to the high rates of transmission by individuals who are asymptomatic at the time of transmission^1,2. Frequent, widespread testing of the asymptomatic population for SARS-CoV-2 is essential to suppress viral transmission. Despite increases in testing capacity, multiple challenges remain in deploying traditional reverse transcription and quantitative PCR (RT-qPCR) tests at the scale required for population screening of asymptomatic individuals. We have developed SwabSeq, a high-throughput testing platform for SARS-CoV-2 that uses next-generation sequencing as a readout. SwabSeq employs sample-specific molecular barcodes to enable thousands of samples to be combined and simultaneously analyzed for the presence or absence of SARS-CoV-2 in a single run. Importantly, SwabSeq incorporates an in vitro RNA standard that mimics the viral amplicon, but can be distinguished by sequencing. This standard allows for end-point rather than quantitative PCR, improves quantitation, reduces requirements for automation and sample-to-sample normalization, enables purification-free detection, and gives better ability to call true negatives. After setting up SwabSeq in a high-complexity CLIA laboratory, we performed more than 80,000 tests for COVID-19 in less than two months, confirming in a real world setting that SwabSeq inexpensively delivers highly sensitive and specific results at scale, with a turn-around of less than 24 hours. Our clinical laboratory uses SwabSeq to test both nasal and saliva samples without RNA extraction, while maintaining analytical sensitivity comparable to or better than traditional RT-qPCR tests. Moving forward, SwabSeq can rapidly scale up testing to mitigate devastating spread of novel pathogens.

From chemoproteomic-detected amino acids to genomic coordinates: insights into precise multi-omic data integration

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.

Metabolic reprogramming and epigenetic changes of vital organs in SARS-CoV-2-induced systemic toxicity

Extrapulmonary manifestations of COVID-19 are associated with a much higher mortality rate than pulmonary manifestations. However, little is known about the pathogenesis of systemic complications of COVID-19. Here, we create a murine model of SARS-CoV-2-induced severe systemic toxicity and multiorgan involvement by expressing the human ACE2 transgene in multiple tissues via viral delivery, followed by systemic administration of SARS-CoV-2. The animals develop a profound phenotype within 7 days with severe weight loss, morbidity, and failure to thrive. We demonstrate that there is metabolic suppression of oxidative phosphorylation and the tricarboxylic acid (TCA) cycle in multiple organs with neutrophilia, lymphopenia, and splenic atrophy, mirroring human COVID-19 phenotypes. Animals had a significantly lower heart rate, and electron microscopy demonstrated myofibrillar disarray and myocardial edema, a common pathogenic cardiac phenotype in human COVID-19. We performed metabolomic profiling of peripheral blood and identified a panel of TCA cycle metabolites that served as biomarkers of depressed oxidative phosphorylation. Finally, we observed that SARS-CoV-2 induces epigenetic changes of DNA methylation, which affects expression of immune response genes and could, in part, contribute to COVID-19 pathogenesis. Our model suggests that SARS-CoV-2-induced metabolic reprogramming and epigenetic changes in internal organs could contribute to systemic toxicity and lethality in COVID-19.

Variants in SCAF4 Cause a Neurodevelopmental Disorder and Are Associated with Impaired mRNA Processing

RNA polymerase II interacts with various other complexes and factors to ensure correct initiation, elongation, and termination of mRNA transcription. One of these proteins is SR-related CTD-associated factor 4 (SCAF4), which is important for correct usage of polyA sites for mRNA termination. Using exome sequencing and international matchmaking, we identified nine likely pathogenic germline variants in SCAF4 including two splice-site and seven truncating variants, all residing in the N-terminal two thirds of the protein. Eight of these variants occurred de novo, and one was inherited. Affected individuals demonstrated a variable neurodevelopmental disorder characterized by mild intellectual disability, seizures, behavioral abnormalities, and various skeletal and structural anomalies. Paired-end RNA sequencing on blood lymphocytes of SCAF4-deficient individuals revealed a broad deregulation of more than 9,000 genes and significant differential splicing of more than 2,900 genes, indicating an important role of SCAF4 in mRNA processing. Knockdown of the SCAF4 ortholog CG4266 in the model organism Drosophila melanogaster resulted in impaired locomotor function, learning, and short-term memory. Furthermore, we observed an increased number of active zones in larval neuromuscular junctions, representing large glutamatergic synapses. These observations indicate a role of CG4266 in nervous system development and function and support the implication of SCAF4 in neurodevelopmental phenotypes. In summary, our data show that heterozygous, likely gene-disrupting variants in SCAF4 are causative for a variable neurodevelopmental disorder associated with impaired mRNA processing.

KAT6A Syndrome: genotype-phenotype correlation in 76 patients with pathogenic KAT6A variants

PURPOSE

Pathogenic variants in KAT6A have recently been identified as a cause of syndromic developmental delay. Within 2 years, the number of patients identified with pathogenic KAT6A variants has rapidly expanded and the full extent and variability of the clinical phenotype has not been reported.

METHODS

We obtained data for patients with KAT6A pathogenic variants through three sources: treating clinicians, an online family survey distributed through social media, and a literature review.

RESULTS

We identified 52 unreported cases, bringing the total number of published cases to 76. Our results expand the genotypic spectrum of pathogenic variants to include missense and splicing mutations. We functionally validated a pathogenic splice-site variant and identified a likely hotspot location for de novo missense variants. The majority of clinical features in KAT6A syndrome have highly variable penetrance. For core features such as intellectual disability, speech delay, microcephaly, cardiac anomalies, and gastrointestinal complications, genotype- phenotype correlations show that late-truncating pathogenic variants (exons 16-17) are significantly more prevalent. We highlight novel associations, including an increased risk of gastrointestinal obstruction.

CONCLUSION

Our data expand the genotypic and phenotypic spectrum for individuals with genetic pathogenic variants in KAT6A and we outline appropriate clinical management.

An Overview of DNA Analytical Methods

The development of rapid parallel sequencing in the last 20 years has begun a revolution in the field of genetics that is changing nearly all disciplines within biology and medicine. Genomic sequencing has become crucial to the diagnosis and clinical management of patients with constitutional diseases and cancer and has quickly become an integral part of the new era of personalized and precision medicine. The precision medicine initiative, released by the NIH in 2015, has catapulted genomic technologies to the forefront of the practice of medicine and biomedical research.This chapter focuses on the core technologies driving the genomic revolution from first generation (Sanger) sequencing to microarray-based technologies, to second, commonly referred to as next-generation sequencing (NGS) methods, and finally to the emerging third generation technologies capable of performing single-molecule and long-read sequencing. The goal of the chapter is to provide a broad overview of these methods of DNA analysis and highlight their strengths and weaknesses. Furthermore, with a knowledge of the different mutation types, we seek to provide the basis for understanding how these technologies work, and can be adopted, to explore other type of nucleic acids and epigenetic changes.

Identification of novel PIEZO1 variants using prenatal exome sequencing and correlation to ultrasound and autopsy findings of recurrent hydrops fetalis

Nonimmune hydrops fetalis (NIHF) is a rare disorder with a high perinatal mortality of at least 50%. One cause of NIHF is generalized lymphatic dysplasia (GLD), a rare form of primary lymphedema of the extremities and systemic involvement including chylothoraces and pericardial effusions. An autosomal recessive form of GLD has been described, caused by variants in the PIEZO1 gene. It has been reported clinically to cause NIHF and childhood onset of facial and limb lymphedema, most of which were diagnosed postnatally. We present a case of a woman with recurrent pregnancies affected by NIHF because of novel compound heterozygous variants in the PIEZO1 gene diagnosed prenatally using exome sequencing (ES). Two variants in PIEZO1 (c.3206G>A and c.6208A>C) were identified that were inherited from the father and mother, and are predicted to cause a nonsense and missense change, respectively, in the PIEZO1 subunits. Ultrasound demonstrated severe bilateral pleural effusions, whole body edema and polyhydramnios. Histopathology revealed an increased number of lymphatic channels, many of which showed failure of luminal canalization. Sanger sequencing confirmed the same variants in a prior fetal demise. We provide phenotypic correlation with ultrasound and autopsy finding, review PIEZO1 variants as a cause of GLD and discuss the uses of prenatal ES to date.

Identification of novel candidate genes for 46,XY disorders of sex development (DSD) using a C57BL/6J-Y POS mouse model

BACKGROUND

Disorders of sex development (DSD) have an estimated frequency of 0.5% of live births encompassing a variety of urogenital anomalies ranging from mild hypospadias to a discrepancy between sex chromosomes and external genitalia. In order to identify the underlying genetic etiology, we had performed exome sequencing in a subset of DSD cases with 46,XY karyotype and were able to identify the causative genetic variant in 35% of cases. While the genetic etiology was not ascertained in more than half of the cases, a large number of variants of unknown clinical significance (VUS) were identified in those exomes.

METHODS

To investigate the relevance of these VUS in regards to the patient's phenotype, we utilized a mouse model in which the presence of a Y chromosome from the poschiavinus strain (Y POS ) on a C57BL/6J (B6) background results in XY undervirilization and sex reversal, a phenotype characteristic to a large subset of human 46,XY DSD cases. We assessed gene expression differences between B6-Y B6 and undervirilized B6-Y POS gonads at E11.5 and identified 515 differentially expressed genes (308 underexpressed and 207 overexpressed in B6-Y POS males).

RESULTS

We identified 15 novel candidate genes potentially involved in 46,XY DSD pathogenesis by filtering the list of human VUS-carrying genes provided by exome sequencing with the list of differentially expressed genes from B6-Y POS mouse model. Additionally, we identified that 7 of the 15 candidate genes were significantly underexpressed in the XY gonads of mice with suppressed Sox9 expression in Sertoli cells suggesting that some of the candidate genes may be downstream of a well-known sex determining gene, Sox9.

CONCLUSION

The use of a DSD-specific animal model improves variant interpretation by correlating human sequence variants with transcriptome variation.

Ensuring the Quality of Point-of-Care Testing in a Large and Decentralized Ambulatory Care Setting

OBJECTIVES

In this project, we assessed the breadth, quality, trends, and outcomes of point-of-care (POC) testing and regulatory compliance in 200 University of California, Los Angeles (UCLA) Health system outpatient clinics.

METHODS

We retrospectively extracted clinic POC test orders, results, and e-prescription data from the UCLA electronic health record over a 1-year period.

RESULTS

Over 100,000 total tests were performed, encompassing 10 POC tests. Initially, 40% of clinics did not have complete licensure, but after implementation of the POC team, this metric improved to 100% licensure within 6 months. Most clinics used two or fewer POC tests, resulted fewer than 200 tests per year, and performed little to no external quality control measures. Our data analytics approach showed that peak POC testing occurred in January 2015, driven by influenza and urinalysis testing, and that both the testing and resulting clinical decision making do not routinely follow society guidelines.

CONCLUSIONS

This decentralization of laboratory testing presents challenges to ensuring quality POC testing. Optimization and analysis of informatics data allow for the identification of POC test utilization trends, areas of improvement for clinical workflows, and increased education on national guidelines.

Loss of Nardilysin, a Mitochondrial Co-chaperone for α-Ketoglutarate Dehydrogenase, Promotes mTORC1 Activation and Neurodegeneration

We previously identified mutations in Nardilysin (dNrd1) in a forward genetic screen designed to isolate genes whose loss causes neurodegeneration in Drosophila photoreceptor neurons. Here we show that NRD1 is localized to mitochondria, where it recruits mitochondrial chaperones and assists in the folding of α-ketoglutarate dehydrogenase (OGDH), a rate-limiting enzyme in the Krebs cycle. Loss of Nrd1 or Ogdh leads to an increase in α-ketoglutarate, a substrate for OGDH, which in turn leads to mTORC1 activation and a subsequent reduction in autophagy. Inhibition of mTOR activity by rapamycin or partially restoring autophagy delays neurodegeneration in dNrd1 mutant flies. In summary, this study reveals a novel role for NRD1 as a mitochondrial co-chaperone for OGDH and provides a mechanistic link between mitochondrial metabolic dysfunction, mTORC1 signaling, and impaired autophagy in neurodegeneration.

MicroRNA-146a modulates B-cell oncogenesis by regulating Egr1

miR-146a is a NF-κB induced microRNA that serves as a feedback regulator of this critical pathway. In mice, deficiency of miR-146a results in hematolymphoid cancer at advanced ages as a consequence of constitutive NF-κB activity. In this study, we queried whether the deficiency of miR-146a contributes to B-cell oncogenesis. Combining miR-146a deficiency with transgenic expression of c-Myc led to the development of highly aggressive B-cell malignancies. Mice transgenic for c-Myc and deficient for miR-146a were characterized by significantly shortened survival, increased lymph node involvement, differential involvement of the spleen and a mature B-cell phenotype. High-throughput sequencing of the tumors revealed significant dysregulation of approximately 250 genes. Amongst these, the transcription factor Egr1 was consistently upregulated in mice deficient for miR-146a. Interestingly, transcriptional targets of Egr1 were enriched in both the high-throughput dataset and in a larger set of miR-146a-deficient tumors. miR-146a overexpression led to downregulation of Egr1 and downstream targets with concomitant decrease in cell growth. Direct targeting of the human EGR1 by miR-146a was seen by luciferase assay. Together our findings illuminate a bona fide role for miR-146a in the modulation of B-cell oncogenesis and reveal the importance of understanding microRNA function in a cell- and disease-specific context.

A novel follicle-stimulating hormone receptor mutation causing primary ovarian failure: a fertility application of whole exome sequencing

STUDY QUESTION

Can whole exome sequencing (WES) and in vitro validation studies be used to find the causative genetic etiology in a patient with primary ovarian failure and infertility?

SUMMARY ANSWER

A novel follicle-stimulating hormone receptor (FSHR) mutation was found by WES and shown, via in vitro flow cytometry studies, to affect membrane trafficking.

WHAT IS KNOWN ALREADY

WES may diagnose up to 25-35% of patients with suspected disorders of sex development (DSD). FSHR mutations are an extremely rare cause of 46, XX gonadal dysgenesis with primary amenorrhea due to hypergonadotropic ovarian failure.

STUDY DESIGN, SIZE, DURATION

A WES study was followed by flow cytometry studies of mutant protein function.

PARTICIPANTS/MATERIALS, SETTING, METHODS

The study subjects were two Turkish sisters with hypergonadotropic primary amenorrhea, their parents and two unaffected sisters. The affected siblings and both parents were sequenced (trio-WES). Transient transfection of HEK 293T cells was performed with a vector containing wild-type FSHR as well as the novel FSHR variant that was discovered by WES. Cellular localization of FSHR protein as well as FSH-stimulated cyclic AMP (cAMP) production was evaluated using flow cytometry.

MAIN RESULTS AND THE ROLE OF CHANCE

Both affected sisters were homozygous for a previously unreported missense mutation (c.1222G>T, p.Asp408Tyr) in the second transmembrane domain of FSHR. Modeling predicted disrupted secondary structure. Flow cytometry demonstrated an average of 48% reduction in cell-surface signal detection (P < 0.01). The mean fluorescent signal for cAMP (second messenger of FSHR), stimulated by FSH, was reduced by 50% in the mutant-transfected cells (P < 0.01).

LIMITATIONS, REASONS FOR CAUTION

This is an in vitro validation. All novel purported genetic variants can be clinically reported only as 'variants of uncertain significance' until more patients with a similar phenotype are discovered with the same variant.

WIDER IMPLICATIONS OF THE FINDINGS

We report the first WES-discovered FSHR mutation, validated by quantitative flow cytometry. WES is a valuable tool for diagnosis of rare genetic diseases, and flow cytometry allows for quantitative characterization of purported variants. WES-assisted diagnosis allows for treatments aimed at the underlying molecular etiology of disease. Future studies should focus on pharmacological and assisted reproductive treatments aimed at the disrupted FSHR, so that patients with FSH resistance can be treated by personalized medicine.

STUDY FUNDING/COMPETING INTERESTS

E.V. is partially funded by the DSD Translational Research Network (NICHD 1R01HD068138). M.S.B. is funded by the Neuroendocrinology, Sex Differences and Reproduction training grant (NICHD 5T32HD007228). The authors have no competing interests to disclose.

Mutations in the PCNA-binding site of CDKN1C inhibit cell proliferation by impairing the entry into S phase

CDKN1C (also known as P57 (kip2) ) is a cyclin-dependent kinase inhibitor that functions as a negative regulator of cell proliferation through G1 phase cell cycle arrest. Recently, our group described gain-of-function mutations in the PCNA-binding site of CDKN1C that result in an undergrowth syndrome called IMAGe Syndrome (Intrauterine Growth Restriction, Metaphyseal dysplasia, Adrenal hypoplasia, and Genital anomalies), with life-threatening consequences. Loss-of-function mutations in CDKN1C have been identified in 5-10% of individuals with Beckwith-Wiedemann syndrome (BWS), an overgrowth disorder with features that are the opposite of IMAGe syndrome. Here, we investigate the effects of IMAGe-associated mutations on protein stability, cell cycle progression and cell proliferation. Mutations in the PCNA-binding site of CDKN1C significantly increase CDKN1C protein stability and prevent cell cycle progression into the S phase. Overexpression of either wild-type or BWS-mutant CDKN1C inhibited cell proliferation. However, the IMAGe-mutant CDKN1C protein decreased cell growth significantly more than both the wild-type or BWS protein. These findings bring new insights into the molecular events underlying IMAGe syndrome.

De novo nonsense mutations in KAT6A, a lysine acetyl-transferase gene, cause a syndrome including microcephaly and global developmental delay

Chromatin remodeling through histone acetyltransferase (HAT) and histone deactylase (HDAC) enzymes affects fundamental cellular processes including the cell-cycle, cell differentiation, metabolism, and apoptosis. Nonsense mutations in genes that are involved in histone acetylation and deacetylation result in multiple congenital anomalies with most individuals displaying significant developmental delay, microcephaly and dysmorphism. Here, we report a syndrome caused by de novo heterozygous nonsense mutations in KAT6A (a.k.a., MOZ, MYST3) identified by clinical exome sequencing (CES) in four independent families. The same de novo nonsense mutation (c.3385C>T [p.Arg1129∗]) was observed in three individuals, and the fourth individual had a nearby de novo nonsense mutation (c.3070C>T [p.Arg1024∗]). Neither of these variants was present in 1,815 in-house exomes or in public databases. Common features among all four probands include primary microcephaly, global developmental delay including profound speech delay, and craniofacial dysmorphism, as well as more varied features such as feeding difficulties, cardiac defects, and ocular anomalies. We further demonstrate that KAT6A mutations result in dysregulation of H3K9 and H3K18 acetylation and altered P53 signaling. Through histone and non-histone acetylation, KAT6A affects multiple cellular processes and illustrates the complex role of acetylation in regulating development and disease.

Exome sequencing for the diagnosis of 46,XY disorders of sex development

CONTEXT

Disorders of sex development (DSD) are clinical conditions where there is a discrepancy between the chromosomal sex and the phenotypic (gonadal or genital) sex of an individual. Such conditions can be stressful for patients and their families and have historically been difficult to diagnose, especially at the genetic level. In particular, for cases of 46,XY gonadal dysgenesis, once variants in SRY and NR5A1 have been ruled out, there are few other single gene tests available.

OBJECTIVE

We used exome sequencing followed by analysis with a list of all known human DSD-associated genes to investigate the underlying genetic etiology of 46,XY DSD patients who had not previously received a genetic diagnosis.

DESIGN

Samples were either submitted to the research laboratory or submitted as clinical samples to the UCLA Clinical Genomic Center. Sequencing data were filtered using a list of genes known to be involved in DSD.

RESULTS

We were able to identify a likely genetic diagnosis in more than a third of cases, including 22.5% with a pathogenic finding, an additional 12.5% with likely pathogenic findings, and 15% with variants of unknown clinical significance.

CONCLUSIONS

Early identification of the genetic cause of a DSD will in many cases streamline and direct the clinical management of the patient, with more focused endocrine and imaging studies and better-informed surgical decisions. Exome sequencing proved an efficient method toward such a goal in 46,XY DSD patients.

DSDs: genetics, underlying pathologies and psychosexual differentiation

Mammalian sex determination is the unique process whereby a single organ, the bipotential gonad, undergoes a developmental switch that promotes its differentiation into either a testis or an ovary. Disruptions of this complex genetic process during human development can manifest as disorders of sex development (DSDs). Sex development can be divided into two distinct processes: sex determination, in which the bipotential gonads form either testes or ovaries, and sex differentiation, in which the fully formed testes or ovaries secrete local and hormonal factors to drive differentiation of internal and external genitals, as well as extragonadal tissues such as the brain. DSDs can arise from a number of genetic lesions, which manifest as a spectrum of gonadal (gonadal dysgenesis to ovotestis) and genital (mild hypospadias or clitoromegaly to ambiguous genitalia) phenotypes. The physical attributes and medical implications associated with DSDs confront families of affected newborns with decisions, such as gender of rearing or genital surgery, and additional concerns, such as uncertainty over the child's psychosexual development and personal wishes later in life. In this Review, we discuss the underlying genetics of human sex determination and focus on emerging data, genetic classification of DSDs and other considerations that surround gender development and identity in individuals with DSDs.

Mutant cohesin in premature ovarian failure

Premature ovarian failure is a major cause of female infertility. The genetic causes of this disorder remain unknown in most patients. Using whole-exome sequence analysis of a large consanguineous family with inherited premature ovarian failure, we identified a homozygous 1-bp deletion inducing a frameshift mutation in STAG3 on chromosome 7. STAG3 encodes a meiosis-specific subunit of the cohesin ring, which ensures correct sister chromatid cohesion. Female mice devoid of Stag3 are sterile, and their fetal oocytes are arrested at early prophase I, leading to oocyte depletion at 1 week of age.

Targeted massively parallel sequencing provides comprehensive genetic diagnosis for patients with disorders of sex development

Disorders of sex development (DSD) are rare disorders in which there is discordance between chromosomal, gonadal, and phenotypic sex. Only a minority of patients clinically diagnosed with DSD obtains a molecular diagnosis, leaving a large gap in our understanding of the prevalence, management, and outcomes in affected patients. We created a novel DSD-genetic diagnostic tool, in which sex development genes are captured using RNA probes and undergo massively parallel sequencing. In the pilot group of 14 patients, we determined sex chromosome dosage, copy number variation, and gene mutations. In the patients with a known genetic diagnosis (obtained either on a clinical or research basis), this test identified the molecular cause in 100% (7/7) of patients. In patients in whom no molecular diagnosis had been made, this tool identified a genetic diagnosis in two of seven patients. Targeted sequencing of genes representing a specific spectrum of disorders can result in a higher rate of genetic diagnoses than current diagnostic approaches. Our DSD diagnostic tool provides for first time, in a single blood test, a comprehensive genetic diagnosis in patients presenting with a wide range of urogenital anomalies.

The evolution of the search for novel genes in mammalian sex determination: from mice to men

Disorders of sex determination are a genetically heterogeneous group of rare disorders, presenting with sex-specific phenotypes and variable expressivity. Prior to the advent of the Human Genome Project, the identification of novel mammalian sex determination genes was hindered by the rarity of disorders of sex determination and small family sizes that made traditional linkage approaches difficult, if not impossible. This article reviews the revolutionary role of the Human Genome Project in the history of sex determination research and highlights the important role of inbred mouse models in elucidating the role of identified sex determination genes in mammalian sex determination. Next generation sequencing technologies has made it possible to sequence complete human genomes or exomes for the purpose of providing a genetic diagnosis to more patients with unexplained disorders of sex determination and identifying novel sex determination genes. However, beyond novel gene discovery, these tools have the power to inform us on more intricate and complex regulation-taking place within the heterogeneous cells that make up the testis and ovary.

Highly efficient small interfering RNA delivery to primary mammalian neurons induces MicroRNA-like effects before mRNA degradation

The study of protein function in neurons has been hindered by the lack of highly efficient, nontoxic methods of inducing RNA interference in such cells. Here we show that application of synthetic small interfering RNA (siRNA) linked to the vector peptide Penetratin1 results in rapid, highly efficient uptake of siRNA by entire populations of cultured primary mammalian hippocampal and sympathetic neurons. This treatment leads to specific knock-down of targeted proteins within hours without the toxicity associated with transfection. In contrast to current methods, our technique permits study of protein function across entire populations with minimal disturbance of complex cellular networks. Using this technique, we found that protein knock-down (evident after 6 hr) precedes any decrease in targeted message (evident after 24 hr), suggesting an early, translational repression by perfectly targeted siRNAs.

Caspase Function in Neuronal Death: Delineation of the Role of Caspases in Ischemia

Cerebral ischemia is one of the major causes of morbidity and mortality in the Western world. Despite extensive research, adequate therapies are still elusive. Neuronal degeneration and death are hallmarks of stroke/ischemia. Understanding how the death machinery executes neuronal death in ischemia will provide therapeutic targets. Key to the death machinery are caspases: the family of cell death proteases. While much data has been published regarding caspase involvement in models of ischemia, the pathways have not been thoroughly defined. The specification of the caspases critical for death has been hampered by the use of non-specific reagents. Thus many conclusions about specificity are unwarranted. In this review we discuss how caspases can be measured and review the existing knowledge of the roles of specific caspases in ischemia. We also discuss approaches to determining the molecules that execute ischemic death.