Whole exome sequencing and co-expression analysis identify an SCN1A variant that modifies pathogenicity in a family with genetic epilepsy and febrile seizures plus.

OBJECTIVE: Family members carrying the same SCN1A variant often exhibit differences in the clinical severity of epilepsy. This variable expressivity suggests that other factors aside from the primary sodium channel variant influence the clinical manifestation. However, identifying such factors has proven challenging in humans. METHODS: We perform whole exome sequencing (WES) in a large family in which an SCN1A variant (p.K1372E) is segregating that is associated with a broad spectrum of phenotypes ranging from lack of epilepsy, to febrile seizures and absence seizures, to Dravet syndrome. We assessed the hypothesis that the severity of the SCN1A-related phenotype was affected by alternate alleles at a modifier locus (or loci). RESULTS: One of our top candidates identified by WES was a second variant in the SCN1A gene (p.L375S) that was shared exclusively by unaffected carriers of the K1372E allele. To test the hypothesized that L375S variant nullifies the loss-of-function effect of K1372E, we transiently expressed Nav1.1 carrying the two variants in HEK293T cells and compared their biophysical properties with the wild-type (WT) variant, and then co-expressed WT with K1372E or L375S with K1372E in equal quantity and tested the functional consequence. The data demonstrated that co-expression of the L375S and K1372E alleles reversed the loss-of-function property brought by the K1372E variant, whereas WT-K1372E co-expression remained partial loss-of-function. SIGNIFICANCE: These results support the hypothesis that L375S counteracts the loss-of-function effect of K1372E such that individuals carrying both alleles in trans do not present epilepsy-related symptoms. We demonstrate that monogenic epilepsies with wide expressivity can be modified by additional variants in the disease gene, providing a novel framework for the gene-phenotype relationship in genetic epilepsies.

Gibbon genome and the fast karyotype evolution of small apes.

Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

Great ape genetic diversity and population history.

Most great ape genetic variation remains uncharacterized; however, its study is critical for understanding population history, recombination, selection and susceptibility to disease. Here we sequence to high coverage a total of 79 wild- and captive-born individuals representing all six great ape species and seven subspecies and report 88.8 million single nucleotide polymorphisms. Our analysis provides support for genetically distinct populations within each species, signals of gene flow, and the split of common chimpanzees into two distinct groups: Nigeria-Cameroon/western and central/eastern populations. We find extensive inbreeding in almost all wild populations, with eastern gorillas being the most extreme. Inferred effective population sizes have varied radically over time in different lineages and this appears to have a profound effect on the genetic diversity at, or close to, genes in almost all species. We discover and assign 1,982 loss-of-function variants throughout the human and great ape lineages, determining that the rate of gene loss has not been different in the human branch compared to other internal branches in the great ape phylogeny. This comprehensive catalogue of great ape genome diversity provides a framework for understanding evolution and a resource for more effective management of wild and captive great ape populations.

Exome sequencing reveals new causal mutations in children with epileptic encephalopathies.

PURPOSE: The management of epilepsy in children is particularly challenging when seizures are resistant to antiepileptic medications, or undergo many changes in seizure type over time, or have comorbid cognitive, behavioral, or motor deficits. Despite efforts to classify such epilepsies based on clinical and electroencephalographic criteria, many children never receive a definitive etiologic diagnosis. Whole exome sequencing (WES) is proving to be a highly effective method for identifying de novo variants that cause neurologic disorders, especially those associated with abnormal brain development. Herein we explore the utility of WES for identifying candidate causal de novo variants in a cohort of children with heterogeneous sporadic epilepsies without etiologic diagnoses. METHODS: We performed WES (mean coverage approximately 40×) on 10 trios comprised of unaffected parents and a child with sporadic epilepsy characterized by difficult-to-control seizures and some combination of developmental delay, epileptic encephalopathy, autistic features, cognitive impairment, or motor deficits. Sequence processing and variant calling were performed using standard bioinformatics tools. A custom filtering system was used to prioritize de novo variants of possible functional significance for validation by Sanger sequencing. KEY FINDINGS: In 9 of 10 probands, we identified one or more de novo variants predicted to alter protein function, for a total of 15. Four probands had de novo mutations in genes previously shown to harbor heterozygous mutations in patients with severe, early onset epilepsies (two in SCN1A, and one each in CDKL5 and EEF1A2). In three children, the de novo variants were in genes with functional roles that are plausibly relevant to epilepsy (KCNH5, CLCN4, and ARHGEF15). The variant in KCNH5 alters one of the highly conserved arginine residues of the voltage sensor of the encoded voltage-gated potassium channel. In vitro analyses using cell-based assays revealed that the CLCN4 mutation greatly impaired ion transport by the ClC-4 2Cl(-) /H(+) -exchanger and that the mutation in ARHGEF15 reduced GEF exchange activity of the gene product, Ephexin5, by about 50%. Of interest, these seven probands all presented with seizures within the first 6 months of life, and six of these have intractable seizures. SIGNIFICANCE: The finding that 7 of 10 children carried de novo mutations in genes of known or plausible clinical significance to neuronal excitability suggests that WES will be of use for the molecular genetic diagnosis of sporadic epilepsies in children, especially when seizures are of early onset and difficult to control.

Identification of an allele-specific transcription factor binding interaction that regulates PLA2G2A gene expression.

The secreted phospholipase A 2 (sPLA 2 ) isoform, sPLA 2 -IIA, has been implicated in a variety of diseases and conditions, including bacteremia, cardiovascular disease, COVID-19, sepsis, adult respiratory distress syndrome, and certain cancers. Given its significant role in these conditions, understanding the regulatory mechanisms impacting its levels is crucial. Genome-wide association studies (GWAS) have identified several single nucleotide polymorphisms (SNPs), including rs11573156, that are associated with circulating levels of sPLA 2 -IIA. Through Genotype-Tissue Expression (GTEx), 234 expression quantitative trait loci (eQTLs) were identified for the gene that encodes for sPLA 2 -IIA, PLA2G2A . SNP2TFBS ( https://ccg.epfl.ch/snp2tfbs/ ) was utilized to ascertain the binding affinities between transcription factors (TFs) to both the reference and alternative alleles of identified SNPs. Subsequently, ChIP-seq peaks highlighted the TF combinations that specifically bind to the SNP, rs11573156. SP1 emerged as a significant TF/SNP pair in liver cells, with rs11573156/SP1 interaction being most prominent in liver, prostate, ovary, and adipose tissues. Further analysis revealed that the upregulation of PLA2G2A transcript levels through the rs11573156 variant was affected by tissue SP1 protein levels. By leveraging an ordinary differential equation, structured upon Michaelis-Menten enzyme kinetics assumptions, we modeled the PLA2G2A transcription's dependence on SP1 protein levels, incorporating the SNP's influence. Collectively, these data strongly suggest that the binding affinity differences of SP1 for the different rs11573156 alleles can influence PLA2G2A expression. This, in turn, can modulate sPLA2-IIA levels, impacting a wide range of human diseases.

Genome sequence of the facultative anaerobic arsenite-oxidizing and nitrate-reducing bacterium Acidovorax sp. strain NO1.

Acidovorax sp. strain NO1, isolated from gold mine soil, was shown to be a facultative anaerobic arsenite-oxidizing and nitrate-reducing bacterium. The reported draft genome predicts the presence of genes involved in arsenic metabolism, nitrate reduction, phosphate transport, and multiple metal resistances and indicates putative horizontal gene transfer events.

Genome sequence of the highly efficient arsenite-oxidizing bacterium Achromobacter arsenitoxydans SY8.

We report the draft genome sequence of Achromobacter arsenitoxydans SY8, the first reported arsenite-oxidizing bacterium belonging to the genus Achromobacter and containing a genomic arsenic island, an intact type III secretion system, and multiple metal(loid) transporters. The genome may be helpful to explore the mechanisms intertwining metal(loid) resistance and pathogenicity.

Genome sequence of deep-sea manganese-oxidizing bacterium Marinobacter manganoxydans MnI7-9.

Here we report the draft genome of Marinobacter manganoxydans MnI7-9, isolated from a deep-sea hydrothermal vent in the Indian Ocean and capable of oxidizing manganese even when there is a very high concentration of Mn(2+). The strain also displayed high resistance and adsorption ability toward many metal(loid)s.

Genome sequence of the arsenite-oxidizing strain Agrobacterium tumefaciens 5A.

Microbial transformations of arsenic influence its mobility and toxicity. We report the draft genome sequence of the arsenite-oxidizing strain Agrobacterium tumefaciens 5A isolated from an As-contaminated soil in the Madison River Valley, MT. A large number of metal (or metalloid) resistance genes, especially contributing to arsenite oxidation, were identified.

Draft genome sequence of plant growth-promoting rhizobium Mesorhizobium amorphae, isolated from zinc-lead mine tailings.

Here, we describe the draft genome sequence of Mesorhizobium amorphae strain CCNWGS0123, isolated from nodules of Robinia pseudoacacia growing on zinc-lead mine tailings. A large number of metal(loid) resistance genes, as well as genes reported to promote plant growth, were identified, presenting a great future potential for aiding phytoremediation in metal(loid)-contaminated soil.

Genome sequence of the moderately halotolerant, arsenite-oxidizing bacterium Pseudomonas stutzeri TS44.

We present the draft genome sequence of Pseudomonas stutzeri TS44, a moderately halotolerant, arsenite-oxidizing bacterium isolated from arsenic-contaminated soil. The genome contains genes for arsenite oxidation, arsenic resistance, and ectoine/hydroxyectoine biosynthesis. The genome information will be useful for exploring adaptation of P. stutzeri TS44 to an arsenic-contaminated environment.

Draft genome sequence of Halomonas sp. strain HAL1, a moderately halophilic arsenite-oxidizing bacterium isolated from gold-mine soil.

We report the draft genome sequence of arsenite-oxidizing Halomonas sp. strain HAL1, isolated from the soil of a gold mine. Genes encoding proteins involved in arsenic resistance and transformation, phosphate utilization and uptake, and betaine biosynthesis were identified. Their identification might help in understanding how arsenic and phosphate metabolism are intertwined.

Genome sequence and mutational analysis of plant-growth-promoting bacterium Agrobacterium tumefaciens CCNWGS0286 Isolated from a zinc-lead mine tailing.

The plant-growth-promoting bacterium Agrobacterium tumefaciens CCNWGS0286, isolated from the nodules of Robinia pseudoacacia growing in zinc-lead mine tailings, both displayed high metal resistance and enhanced the growth of Robinia plants in a metal-contaminated environment. Our goal was to determine whether bacterial metal resistance or the capacity to produce phytohormones had a larger impact on the growth of host plants under zinc stress. Eight zinc-sensitive mutants and one zinc-sensitive mutant with reduced indole-3-acetic acid (IAA) production were obtained by transposon mutagenesis. Analysis of the genome sequence and of transcription via reverse transcriptase PCR (RT-PCR) combined with transposon gene disruptions revealed that ZntA-4200 and the transcriptional regulator ZntR1 played important roles in the zinc homeostasis of A. tumefaciens CCNWGS0286. In addition, interruption of a putative oligoketide cyclase/lipid transport protein reduced IAA synthesis and also showed reduced zinc and cadmium resistance but had no influence on copper resistance. In greenhouse studies, R. pseudoacacia inoculated with A. tumefaciens CCNWGS0286 displayed a significant increase in biomass production over that without inoculation, even in a zinc-contaminated environment. Interestingly, the differences in plant biomass improvement among A. tumefaciens CCNWGS0286, A. tumefaciens C58, and zinc-sensitive mutants 12-2 (zntA::Tn5) and 15-6 (low IAA production) revealed that phytohormones, rather than genes encoding zinc resistance determinants, were the dominant factor in enhancing plant growth in contaminated soil.

Draft genome of Streptomyces zinciresistens K42, a novel metal-resistant species isolated from copper-zinc mine tailings.

A draft genome sequence of Streptomyces zinciresistens K42, a novel Streptomyces species displaying a high level of resistance to zinc and cadmium, is presented here. The genome contains a large number of genes encoding proteins predicted to be involved in conferring metal resistance. Many of these genes appear to have been acquired through horizontal gene transfer.

Interpreting Clinical Trials With Omega-3 Supplements in the Context of Ancestry and FADS Genetic Variation.

Human diets in developed countries such as the US have changed dramatically over the past 75 years, leading to increased obesity, inflammation, and cardiometabolic dysfunction. Evidence over the past decade indicates that the interaction of genetic variation with changes in the intake of 18-carbon essential dietary omega-6 (n-6) and omega-3 (n-3) polyunsaturated fatty acids (PUFA), linoleic acid (LA) and α-linolenic acid (ALA), respectively, has impacted numerous molecular and clinical phenotypes. Interactions are particularly relevant with the FADS1 and FADS2 genes, which encode key fatty acid desaturases in the pathway that converts LA and ALA to their long chain (≥20 carbons), highly unsaturated fatty acid (HUFA) counterparts. These gene by nutrient interactions affect the levels and balance of n-6 and n-3 HUFA that in turn are converted to a wide array of lipids with signaling roles, including eicosanoids, docosanoids, other oxylipins and endocannabinoids. With few exceptions, n-6 HUFA are precursors of pro-inflammatory/pro-thrombotic signaling lipids, and n-3 HUFA are generally anti-inflammatory/anti-thrombotic. We and others have demonstrated that African ancestry populations have much higher frequencies (vs. European-, Asian- or indigenous Americas-ancestry populations) of a FADS "derived" haplotype that is associated with the efficient conversion of high levels of dietary n-6 PUFA to pro-inflammatory n-6 HUFA. By contrast, an "ancestral" haplotype, carrying alleles associated with a limited capacity to synthesize HUFA, which can lead to n-3 HUFA deficiency, is found at high frequency in certain Hispanic populations and is nearly fixed in several indigenous populations from the Americas. Based on these observations, a focused secondary subgroup analysis of the VITAL n-3 HUFA supplementation trial stratifying the data based on self-reported ancestry revealed that African Americans may benefit from n-3 HUFA supplementation, and both ancestry and FADS variability should be factored into future clinical trials design.

Group IIA secreted phospholipase A2 is associated with the pathobiology leading to COVID-19 mortality.

There is an urgent need to identify the cellular and molecular mechanisms responsible for severe COVID-19 that results in death. We initially performed both untargeted and targeted lipidomics as well as focused biochemical analyses of 127 plasma samples and found elevated metabolites associated with secreted phospholipase A2 (sPLA2) activity and mitochondrial dysfunction in patients with severe COVID-19. Deceased COVID-19 patients had higher levels of circulating, catalytically active sPLA2 group IIA (sPLA2-IIA), with a median value that was 9.6-fold higher than that for patients with mild disease and 5.0-fold higher than the median value for survivors of severe COVID-19. Elevated sPLA2-IIA levels paralleled several indices of COVID-19 disease severity (e.g., kidney dysfunction, hypoxia, multiple organ dysfunction). A decision tree generated by machine learning identified sPLA2-IIA levels as a central node in the stratification of patients who died from COVID-19. Random forest analysis and least absolute shrinkage and selection operator-based (LASSO-based) regression analysis additionally identified sPLA2-IIA and blood urea nitrogen (BUN) as the key variables among 80 clinical indices in predicting COVID-19 mortality. The combined PLA-BUN index performed significantly better than did either one alone. An independent cohort (n = 154) confirmed higher plasma sPLA2-IIA levels in deceased patients compared with levels in plasma from patients with severe or mild COVID-19, with the PLA-BUN index-based decision tree satisfactorily stratifying patients with mild, severe, or fatal COVID-19. With clinically tested inhibitors available, this study identifies sPLA2-IIA as a therapeutic target to reduce COVID-19 mortality.

Group IIA Secreted Phospholipase A 2 Plays a Central Role in the Pathobiology of COVID-19.

There is an urgent need to identify cellular and molecular mechanisms responsible for severe COVID-19 disease accompanied by multiple organ failure and high mortality rates. Here, we performed untargeted/targeted lipidomics and focused biochemistry on 127 patient plasma samples, and showed high levels of circulating, enzymatically active secreted phospholipase A 2 Group IIA (sPLA 2 -IIA) in severe and fatal COVID-19 disease compared with uninfected patients or mild illness. Machine learning demonstrated that sPLA 2 -IIA effectively stratifies severe from fatal COVID-19 disease. We further introduce a PLA-BUN index that combines sPLA 2 -IIA and blood urea nitrogen (BUN) threshold levels as a critical risk factor for mitochondrial dysfunction, sustained inflammatory injury and lethal COVID-19. With the availability of clinically tested inhibitors of sPLA 2 -IIA, our study opens the door to a precision intervention using indices discovered here to reduce COVID-19 mortality.

Impact of Amerind ancestry and FADS genetic variation on omega-3 deficiency and cardiometabolic traits in Hispanic populations.

Long chain polyunsaturated fatty acids (LC-PUFAs) have critical signaling roles that regulate dyslipidemia and inflammation. Genetic variation in the FADS gene cluster accounts for a large portion of interindividual differences in circulating and tissue levels of LC-PUFAs, with the genotypes most strongly predictive of low LC-PUFA levels at strikingly higher frequencies in Amerind ancestry populations. In this study, we examined relationships between genetic ancestry and FADS variation in 1102 Hispanic American participants from the Multi-Ethnic Study of Atherosclerosis. We demonstrate strong negative associations between Amerind genetic ancestry and LC-PUFA levels. The FADS rs174537 single nucleotide polymorphism (SNP) accounted for much of the AI ancestry effect on LC-PUFAs, especially for low levels of n-3 LC-PUFAs. Rs174537 was also strongly associated with several metabolic, inflammatory and anthropomorphic traits including circulating triglycerides (TGs) and E-selectin in MESA Hispanics. Our study demonstrates that Amerind ancestry provides a useful and readily available tool to identify individuals most likely to have FADS-related n-3 LC-PUFA deficiencies and associated cardiovascular risk.

Genomics in Personalized Nutrition: Can You "Eat for Your Genes"?

Genome-wide single nucleotide polymorphism (SNP) data are now quickly and inexpensively acquired, raising the prospect of creating personalized dietary recommendations based on an individual's genetic variability at multiple SNPs. However, relatively little is known about most specific gene-diet interactions, and many molecular and clinical phenotypes of interest (e.g., body mass index [BMI]) involve multiple genes. In this review, we discuss direct to consumer genetic testing (DTC-GT) and the current potential for precision nutrition based on an individual's genetic data. We review important issues such as dietary exposure and genetic architecture addressing the concepts of penetrance, pleiotropy, epistasis, polygenicity, and epigenetics. More specifically, we discuss how they complicate using genotypic data to predict phenotypes as well as response to dietary interventions. Then, several examples (including caffeine sensitivity, alcohol dependence, non-alcoholic fatty liver disease, obesity/appetite, cardiovascular, Alzheimer's disease, folate metabolism, long-chain fatty acid biosynthesis, and vitamin D metabolism) are provided illustrating how genotypic information could be used to inform nutritional recommendations. We conclude by examining ethical considerations and practical applications for using genetic information to inform dietary choices and the future role genetics may play in adopting changes beyond population-wide healthy eating guidelines.

Rare Tumor-Normal Matched Whole Exome Sequencing Identifies Novel Genomic Pathogenic Germline and Somatic Aberrations.

Whole exome sequencing (WES) of matched tumor-normal pairs in rare tumors has the potential to identify genome-wide mutations and copy number alterations (CNAs). We evaluated 27 rare cancer patients with tumor-normal matching by WES and tumor-only next generation sequencing (NGS) as a comparator. Our goal was to: 1) identify known and novel variants and CNAs in rare cancers with comparison to common cancers; 2) examine differences between germline and somatic variants and how that functionally impacts rare tumors; 3) detect and characterize alleles in biologically relevant genes-pathways that may be of clinical importance but not represented in classical cancer genes. We identified 3343 germline single nucleotide variants (SNVs) and small indel variants-1670 in oncogenes and 1673 in tumor suppressor genes-generating an average of 124 germline variants/case. The number of somatic SNVs and small indels detected in all cases was 523:306 in oncogenes and 217 in tumor suppressor genes. Of the germline variants, six were identified to be pathogenic or likely pathogenic. In the 27 analyzed rare cancer cases, CNAs are variable depending on tumor type, germline pathogenic variants are more common. Cell fate pathway mutations (e.g., Hippo, Notch, Wnt) dominate pathogenesis and double hit (mutation + CNV) represent ~18% cases.