Experiences in providing a community educational resource for the All of Us Researcher Workbench.

OBJECTIVE: Educational offerings to fill the bioinformatics knowledge gap are a key component to enhancing access and use of health data from the All of Us Research Program. We developed a Train the Trainer-based, innovative training series including project-based learning, modular on-demand demonstrations, and unstructured tutorial time as a model for educational engagement in the All of Us community. MATERIALS AND METHODS: We highlight our training modules and content, with training survey data informing cycles of development in the creation of a 6-module training series with modular demonstrations. RESULTS: We have conducted 2 public iterations of the Train the Trainer (Tx3) Series based on survey feedback while training over 300 registered researchers to access and analyze data on the All of Us Researcher Workbench. DISCUSSION AND CONCLUSION: Future directions of the Tx3 Series include enhanced focus on project-based learning and learner requests for modularity and asynchronous materials access.

Development and evaluation of a training curriculum to engage researchers on accessing and analyzing the All of Us data.

OBJECTIVE: The All of Us Evenings with Genetics (EwG) Research Program at Baylor College of Medicine (BCM), funded to engage research scholars to work with the All of Us data, developed a training curriculum for the Researcher Workbench, the platform to access and analyze All of Us data. All of Us EwG developed the curriculum so that it could teach scholars regardless of their skills and background in programming languages and cloud computing. All of Us EwG delivered this curriculum at the first annual All of Us EwG Faculty Summit in May 2022. The curriculum was evaluated both during and after the Faculty Summit so that it could be improved for future training. MATERIALS AND METHODS: Surveys were administered to assess scholars' familiarity with the programming languages and computational tools required to use the Researcher Workbench. The curriculum was developed using backward design and was informed by the survey results, a review of available resources for training users on the Researcher Workbench, and All of Us EwG members' collective experience training students. The curriculum was evaluated using feedback surveys during the Faculty Summit as well as virtual meetings and emails following the Faculty Summit. RESULTS: The evaluation results demonstrated the success of the curriculum and identified areas for improvement. DISCUSSION AND CONCLUSION: The curriculum has been adapted and improved in response to evaluations and in response to changes to the All of Us data and infrastructure to train more researchers through this program and other scholarly programs.

Diagnostic utility of DNA methylation analysis in genetically unsolved pediatric epilepsies and CHD2 episignature refinement.

Sequence-based genetic testing identifies causative variants in ~ 50% of individuals with developmental and epileptic encephalopathies (DEEs). Aberrant changes in DNA methylation are implicated in various neurodevelopmental disorders but remain unstudied in DEEs. We interrogate the diagnostic utility of genome-wide DNA methylation array analysis on peripheral blood samples from 582 individuals with genetically unsolved DEEs. We identify rare differentially methylated regions (DMRs) and explanatory episignatures to uncover causative and candidate genetic etiologies in 12 individuals. Using long-read sequencing, we identify DNA variants underlying rare DMRs, including one balanced translocation, three CG-rich repeat expansions, and four copy number variants. We also identify pathogenic variants associated with episignatures. Finally, we refine the CHD2 episignature using an 850 K methylation array and bisulfite sequencing to investigate potential insights into CHD2 pathophysiology. Our study demonstrates the diagnostic yield of genome-wide DNA methylation analysis to identify causal and candidate variants as 2% (12/582) for unsolved DEE cases.

Exploring the complexity of systemic sclerosis etiology by trio whole genome sequencing.

Systemic sclerosis (SSc) is a heterogeneous rare autoimmune fibrosing disorder affecting connective tissue. The etiology of systemic sclerosis is largely unknown and many genes have been suggested as susceptibility loci of modest impact by genome-wide association study (GWAS). Multiple factors can contribute to the pathological process of the disease, which makes it more difficult to identify possible disease-causing genetic alterations. In this study, we have applied whole genome sequencing (WGS) in 101 indexed family trios, supplemented with transcriptome sequencing on cultured fibroblast cells of four patients and five family controls where available. Single nucleotide variants (SNVs) and copy number variants (CNVs) were examined, with emphasis on de novo variants. We also performed enrichment test for rare variants in candidate genes previously proposed in association with systemic sclerosis. We identified 42 exonic and 34 ncRNA de novo SNV changes in 101 trios, from a total of over 6000 de novo variants genome wide. We observed higher than expected de novo variants in PRKXP1 gene. We also observed such phenomenon along with increased expression in patient group in NEK7 gene. Additionally, we also observed significant enrichment of rare variants in candidate genes in the patient cohort, further supporting the complexity/multi-factorial etiology of systemic sclerosis. Our findings identify new candidate genes including PRKXP1 and NEK7 for future studies in SSc. We observed rare variant enrichment in candidate genes previously proposed in association with SSc, which suggest more efforts should be pursued to further investigate possible pathogenetic mechanisms associated with those candidate genes.

The clinical utility and diagnostic implementation of human subject cell transdifferentiation followed by RNA sequencing.

RNA sequencing (RNA-seq) has recently been used in translational research settings to facilitate diagnoses of Mendelian disorders. A significant obstacle for clinical laboratories in adopting RNA-seq is the low or absent expression of a significant number of disease-associated genes/transcripts in clinically accessible samples. As this is especially problematic in neurological diseases, we developed a clinical diagnostic approach that enhanced the detection and evaluation of tissue-specific genes/transcripts through fibroblast-to-neuron cell transdifferentiation. The approach is designed specifically to suit clinical implementation, emphasizing simplicity, cost effectiveness, turnaround time, and reproducibility. For clinical validation, we generated induced neurons (iNeurons) from 71 individuals with primary neurological phenotypes recruited to the Undiagnosed Diseases Network. The overall diagnostic yield was 25.4%. Over a quarter of the diagnostic findings benefited from transdifferentiation and could not be achieved by fibroblast RNA-seq alone. This iNeuron transcriptomic approach can be effectively integrated into diagnostic whole-transcriptome evaluation of individuals with genetic disorders.

A comparative study of structural variant calling in WGS from Alzheimer's disease families.

Detecting structural variants (SVs) in whole-genome sequencing poses significant challenges. We present a protocol for variant calling, merging, genotyping, sensitivity analysis, and laboratory validation for generating a high-quality SV call set in whole-genome sequencing from the Alzheimer's Disease Sequencing Project comprising 578 individuals from 111 families. Employing two complementary pipelines, Scalpel and Parliament, for SV/indel calling, we assessed sensitivity through sample replicates (N = 9) with in silico variant spike-ins. We developed a novel metric, D-score, to evaluate caller specificity for deletions. The accuracy of deletions was evaluated by Sanger sequencing. We generated a high-quality call set of 152,301 deletions of diverse sizes. Sanger sequencing validated 114 of 146 detected deletions (78.1%). Scalpel excelled in accuracy for deletions ≤100 bp, whereas Parliament was optimal for deletions >900 bp. Overall, 83.0% and 72.5% of calls by Scalpel and Parliament were validated, respectively, including all 11 deletions called by both Parliament and Scalpel between 101 and 900 bp. Our flexible protocol successfully generated a high-quality deletion call set and a truth set of Sanger sequencing-validated deletions with precise breakpoints spanning 1-17,000 bp.

De novo variants in PLCG1 are associated with hearing impairment, ocular pathology, and cardiac defects.

Phospholipase C isozymes (PLCs) hydrolyze phosphatidylinositol 4,5-bisphosphate into inositol 1,4,5-trisphosphate and diacylglycerol, important signaling molecules involved in many cellular processes. PLCG1 encodes the PLCγ1 isozyme that is broadly expressed. Hyperactive somatic mutations of PLCG1 are observed in multiple cancers, but only one germline variant has been reported. Here we describe three unrelated individuals with de novo heterozygous missense variants in PLCG1 (p.Asp1019Gly, p.His380Arg, and p.Asp1165Gly) who exhibit variable phenotypes including hearing loss, ocular pathology and cardiac septal defects. To model these variants in vivo, we generated the analogous variants in the Drosophila ortholog, small wing (sl). We created a null allele slT2A and assessed the expression pattern. sl is broadly expressed, including in wing discs, eye discs, and a subset of neurons and glia. Loss of sl causes wing size reductions, ectopic wing veins and supernumerary photoreceptors. We document that mutant flies exhibit a reduced lifespan and age-dependent locomotor defects. Expressing wild-type sl in slT2A mutant rescues the loss-of-function phenotypes whereas expressing the variants causes lethality. Ubiquitous overexpression of the variants also reduces viability, suggesting that the variants are toxic. Ectopic expression of an established hyperactive PLCG1 variant (p.Asp1165His) in the wing pouch causes severe wing phenotypes, resembling those observed with overexpression of the p.Asp1019Gly or p.Asp1165Gly variants, further arguing that these two are gain-of-function variants. However, the wing phenotypes associated with p.His380Arg overexpression are mild. Our data suggest that the PLCG1 de novo heterozygous missense variants are pathogenic and contribute to the features observed in the probands.

Diagnostic Utility of Genome-wide DNA Methylation Analysis in Genetically Unsolved Developmental and Epileptic Encephalopathies and Refinement of a CHD2 Episignature.

Sequence-based genetic testing currently identifies causative genetic variants in ∼50% of individuals with developmental and epileptic encephalopathies (DEEs). Aberrant changes in DNA methylation are implicated in various neurodevelopmental disorders but remain unstudied in DEEs. Rare epigenetic variations ("epivariants") can drive disease by modulating gene expression at single loci, whereas genome-wide DNA methylation changes can result in distinct "episignature" biomarkers for monogenic disorders in a growing number of rare diseases. Here, we interrogate the diagnostic utility of genome-wide DNA methylation array analysis on peripheral blood samples from 516 individuals with genetically unsolved DEEs who had previously undergone extensive genetic testing. We identified rare differentially methylated regions (DMRs) and explanatory episignatures to discover causative and candidate genetic etiologies in 10 individuals. We then used long-read sequencing to identify DNA variants underlying rare DMRs, including one balanced translocation, three CG-rich repeat expansions, and two copy number variants. We also identify pathogenic sequence variants associated with episignatures; some had been missed by previous exome sequencing. Although most DEE genes lack known episignatures, the increase in diagnostic yield for DNA methylation analysis in DEEs is comparable to the added yield of genome sequencing. Finally, we refine an episignature for CHD2 using an 850K methylation array which was further refined at higher CpG resolution using bisulfite sequencing to investigate potential insights into CHD2 pathophysiology. Our study demonstrates the diagnostic yield of genome-wide DNA methylation analysis to identify causal and candidate genetic causes as ∼2% (10/516) for unsolved DEE cases.

Dominant negative variants in IKZF2 cause ICHAD syndrome, a new disorder characterised by immunodysregulation, craniofacial anomalies, hearing impairment, athelia and developmental delay.

BACKGROUND: Helios (encoded by IKZF2), a member of the Ikaros family of transcription factors, is a zinc finger protein involved in embryogenesis and immune function. Although predominantly recognised for its role in the development and function of T lymphocytes, particularly the CD4+ regulatory T cells (Tregs), the expression and function of Helios extends beyond the immune system. During embryogenesis, Helios is expressed in a wide range of tissues, making genetic variants that disrupt the function of Helios strong candidates for causing widespread immune-related and developmental abnormalities in humans. METHODS: We performed detailed phenotypic, genomic and functional investigations on two unrelated individuals with a phenotype of immune dysregulation combined with syndromic features including craniofacial differences, sensorineural hearing loss and congenital abnormalities. RESULTS: Genome sequencing revealed de novo heterozygous variants that alter the critical DNA-binding zinc fingers (ZFs) of Helios. Proband 1 had a tandem duplication of ZFs 2 and 3 in the DNA-binding domain of Helios (p.Gly136_Ser191dup) and Proband 2 had a missense variant impacting one of the key residues for specific base recognition and DNA interaction in ZF2 of Helios (p.Gly153Arg). Functional studies confirmed that both these variant proteins are expressed and that they interfere with the ability of the wild-type Helios protein to perform its canonical function-repressing IL2 transcription activity-in a dominant negative manner. CONCLUSION: This study is the first to describe dominant negative IKZF2 variants. These variants cause a novel genetic syndrome characterised by immunodysregulation, craniofacial anomalies, hearing impairment, athelia and developmental delay.

Clinical exome sequencing uncovers a high frequency of Mendelian disorders in infants with stroke: A retrospective analysis.

Stroke causes significant disability and is a common cause of death worldwide. Previous studies have estimated that 1%-5% of stroke is attributable to monogenic etiologies. We set out to assess the utility of clinical exome sequencing (ES) in the evaluation of stroke. We retrospectively analyzed 124 individuals who received ES at the Baylor Genetics reference lab between 2012 and 2021 who had stroke as a major part of their reported phenotype. Ages ranged from 10 days to 69 years. 8.9% of the cohort received a diagnosis, including 25% of infants less than 1 year old; an additional 10.5% of the cohort received a probable diagnosis. We identified several syndromes that predispose to stroke such as COL4A1-related brain small vessel disease, homocystinuria caused by CBS mutation, POLG-related disorders, TTC19-linked mitochondrial disease, and RNASEH2A associated Aicardi-Goutieres syndrome. We also observed pathogenic variants in NSD1, PKHD1, HRAS, and ATP13A2, which are genes rarely associated with stroke. Although stroke is a complex phenotype with varying pathologies and risk factors, these results show that use of exome sequencing can be highly relevant in stroke, especially for those presenting <1 year of age.

An improved ovine reference genome assembly to facilitate in-depth functional annotation of the sheep genome.

BACKGROUND: The domestic sheep (Ovis aries) is an important agricultural species raised for meat, wool, and milk across the world. A high-quality reference genome for this species enhances the ability to discover genetic mechanisms influencing biological traits. Furthermore, a high-quality reference genome allows for precise functional annotation of gene regulatory elements. The rapid advances in genome assembly algorithms and emergence of sequencing technologies with increasingly long reads provide the opportunity for an improved de novo assembly of the sheep reference genome. FINDINGS: Short-read Illumina (55× coverage), long-read Pacific Biosciences (75× coverage), and Hi-C data from this ewe retrieved from public databases were combined with an additional 50× coverage of Oxford Nanopore data and assembled with canu v1.9. The assembled contigs were scaffolded using Hi-C data with Salsa v2.2, gaps filled with PBsuitev15.8.24, and polished with Nanopolish v0.12.5. After duplicate contig removal with PurgeDups v1.0.1, chromosomes were oriented and polished with 2 rounds of a pipeline that consisted of freebayes v1.3.1 to call variants, Merfin to validate them, and BCFtools to generate the consensus fasta. The ARS-UI_Ramb_v2.0 assembly is 2.63 Gb in length and has improved continuity (contig NG50 of 43.18 Mb), with a 19- and 38-fold decrease in the number of scaffolds compared with Oar_rambouillet_v1.0 and Oar_v4.0. ARS-UI_Ramb_v2.0 has greater per-base accuracy and fewer insertions and deletions identified from mapped RNA sequence than previous assemblies. CONCLUSIONS: The ARS-UI_Ramb_v2.0 assembly is a substantial improvement in contiguity that will optimize the functional annotation of the sheep genome and facilitate improved mapping accuracy of genetic variant and expression data for traits in sheep.

LMOD2-related dilated cardiomyopathy presenting in late infancy.

Leiomodin-2 (LMOD2) is an important regulator of the thin filament length, known to promote elongation of actin through polymerization at pointed ends. Mice with Lmod2 deficiency die around 3 weeks of age due to severe dilated cardiomyopathy (DCM), resulting from decreased heart contractility due to shorter thin filaments. To date, there have been three infants from two families reported with biallelic variants in LMOD2, presenting with perinatal onset DCM. Here, we describe a third family with a child harboring a previously described homozygous frameshift variant, c.1243_1244delCT (p.L415Vfs*108) with DCM, presenting later in infancy at 9 months of age. Family history was relevant for a sibling who died suddenly at 1 year of age after being diagnosed with cardiomegaly. LMOD2-related cardiomyopathy is a rare form of inherited cardiomyopathy resulting from thin filament length dysregulation and should be considered in genetic evaluation of newborns and infants with suspected autosomal recessive inheritance or sporadic early onset cardiomyopathy.

The genomic basis of evolutionary differentiation among honey bees.

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Comparative genomic analysis of sifakas (Propithecus) reveals selection for folivory and high heterozygosity despite endangered status.

Sifakas (genus Propithecus) are critically endangered, large-bodied diurnal lemurs that eat leaf-based diets and show corresponding anatomical and microbial adaptations to folivory. We report on the genome assembly of Coquerel's sifaka (P. coquereli) and the resequenced genomes of Verreaux's (P. verreauxi), the golden-crowned (P. tattersalli), and the diademed (P. diadema) sifakas. We find high heterozygosity in all sifakas compared with other primates and endangered mammals. Demographic reconstructions nevertheless suggest declines in effective population size beginning before human arrival on Madagascar. Comparative genomic analyses indicate pervasive accelerated evolution in the ancestral sifaka lineage affecting genes in several complementary pathways relevant to folivory, including nutrient absorption and xenobiotic and fatty acid metabolism. Sifakas show convergent evolution at the level of the pathway, gene family, gene, and amino acid substitution with other folivores. Although sifakas have relatively generalized diets, the physiological challenges of habitual folivory likely led to strong selection.

Primate phylogenomics uncovers multiple rapid radiations and ancient interspecific introgression.

Our understanding of the evolutionary history of primates is undergoing continual revision due to ongoing genome sequencing efforts. Bolstered by growing fossil evidence, these data have led to increased acceptance of once controversial hypotheses regarding phylogenetic relationships, hybridization and introgression, and the biogeographical history of primate groups. Among these findings is a pattern of recent introgression between species within all major primate groups examined to date, though little is known about introgression deeper in time. To address this and other phylogenetic questions, here, we present new reference genome assemblies for 3 Old World monkey (OWM) species: Colobus angolensis ssp. palliatus (the black and white colobus), Macaca nemestrina (southern pig-tailed macaque), and Mandrillus leucophaeus (the drill). We combine these data with 23 additional primate genomes to estimate both the species tree and individual gene trees using thousands of loci. While our species tree is largely consistent with previous phylogenetic hypotheses, the gene trees reveal high levels of genealogical discordance associated with multiple primate radiations. We use strongly asymmetric patterns of gene tree discordance around specific branches to identify multiple instances of introgression between ancestral primate lineages. In addition, we exploit recent fossil evidence to perform fossil-calibrated molecular dating analyses across the tree. Taken together, our genome-wide data help to resolve multiple contentious sets of relationships among primates, while also providing insight into the biological processes and technical artifacts that led to the disagreements in the first place.

Global Analysis of Transcription Start Sites in the New Ovine Reference Genome (Oar rambouillet v1.0).

The overall aim of the Ovine FAANG project is to provide a comprehensive annotation of the new highly contiguous sheep reference genome sequence (Oar rambouillet v1.0). Mapping of transcription start sites (TSS) is a key first step in understanding transcript regulation and diversity. Using 56 tissue samples collected from the reference ewe Benz2616, we have performed a global analysis of TSS and TSS-Enhancer clusters using Cap Analysis Gene Expression (CAGE) sequencing. CAGE measures RNA expression by 5' cap-trapping and has been specifically designed to allow the characterization of TSS within promoters to single-nucleotide resolution. We have adapted an analysis pipeline that uses TagDust2 for clean-up and trimming, Bowtie2 for mapping, CAGEfightR for clustering, and the Integrative Genomics Viewer (IGV) for visualization. Mapping of CAGE tags indicated that the expression levels of CAGE tag clusters varied across tissues. Expression profiles across tissues were validated using corresponding polyA+ mRNA-Seq data from the same samples. After removal of CAGE tags with <10 read counts, 39.3% of TSS overlapped with 5' ends of 31,113 transcripts that had been previously annotated by NCBI (out of a total of 56,308 from the NCBI annotation). For 25,195 of the transcripts, previously annotated by NCBI, no TSS meeting stringent criteria were identified. A further 14.7% of TSS mapped to within 50 bp of annotated promoter regions. Intersecting these predicted TSS regions with annotated promoter regions (±50 bp) revealed 46% of the predicted TSS were "novel" and previously un-annotated. Using whole-genome bisulfite sequencing data from the same tissues, we were able to determine that a proportion of these "novel" TSS were hypo-methylated (32.2%) indicating that they are likely to be reproducible rather than "noise". This global analysis of TSS in sheep will significantly enhance the annotation of gene models in the new ovine reference assembly. Our analyses provide one of the highest resolution annotations of transcript regulation and diversity in a livestock species to date.

Sequence analysis in Bos taurus reveals pervasiveness of X-Y arms races in mammalian lineages.

Studies of Y Chromosome evolution have focused primarily on gene decay, a consequence of suppression of crossing-over with the X Chromosome. Here, we provide evidence that suppression of X-Y crossing-over unleashed a second dynamic: selfish X-Y arms races that reshaped the sex chromosomes in mammals as different as cattle, mice, and men. Using super-resolution sequencing, we explore the Y Chromosome of Bos taurus (bull) and find it to be dominated by massive, lineage-specific amplification of testis-expressed gene families, making it the most gene-dense Y Chromosome sequenced to date. As in mice, an X-linked homolog of a bull Y-amplified gene has become testis-specific and amplified. This evolutionary convergence implies that lineage-specific X-Y coevolution through gene amplification, and the selfish forces underlying this phenomenon, were dominatingly powerful among diverse mammalian lineages. Together with Y gene decay, X-Y arms races molded mammalian sex chromosomes and influenced the course of mammalian evolution.

Correction to: Genome-enabled insights into the biology of thrips as crop pests.

An amendment to this paper has been published and can be accessed via the original article.

Genome-enabled insights into the biology of thrips as crop pests.

BACKGROUND: The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set. RESULTS: We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta. CONCLUSIONS: Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.