The cytochrome oxidase subunit I gene of Tetrahymena: a 57 amino acid NH2-terminal extension and a 108 amino acid insert

Summary The gene sequence for cytochrome oxidase subunit I (COI) in the ciliate Tetrahymena mitochondrial DNA has been determined and shown to be coded by the same strand as codes the genes (in order) for 14S rRNA, tRNA^trp tRNA^glu 21S rRNA, tRNA^leu and tRNA^met. The predicted protein has 698 amino acids, including an NH₂-terminal 57 amino acid extension and a 108 amino acid insert originally found in Paramecium COI. These extension and insert segments are not highly hydrophobic but are relatively rich in lysine, arginine and serine. In analogy with the presequence of nuclear-encoded mitochondrial proteins, they might function as a transmembrane signal. The remaining poly-peptide segments show a hydrophobicity characteristic of membrane spanning proteins. TCOI shows a 64% amino acid identity with Paramecium COI but less than a 38% amino acid conservation with human COI. The Tetrahymena mitochondrial code is analogous with the mammalian mitochondrial code; but differs from the Tetrahymena nuclear genetic code; TGA is exclusively translated as tryptophan; ATA is used as an initiation codon probably for methionine, and TAA as a stop codon; the arginine codons (CGN) are not used. The use of the leucine codon TTA in TCOI is contradictory to the codon recognition pattern previously obtained from the isolated tRNA^leu isoacceptors recognizing only the CUN codons, but consistent with the tRNA^leu (anticodon UAA) gene encoded in the genome. The reason for this inconsistency has not been resolved.     

B*4440: a novel HLA-B allele identified by sequence-specific oligonucleotide hybridization and sequence-specific amplification

A novel human leukocyte antigen B (HLA-B) allele, B*4440, is described. The allele was identified in an adult stem cell donor of Caucasian origin. HLA-B*4440 most closely matches to B*4403 differing by a substitution of three nucleotides at codon 44, 45, and 50. Thus, low-resolution HLA typing using sequence-specific oligonucleotide hybridization or amplification using sequence-specific primers gave inconclusive results. DNA sequencing confirmed a variation of codons 44 and 45 (AGG AAG-->AGA GAG) and codon 50 (CCA-->CCG), resulting in an amino acid substitution Lys-->Glu at codon 45.     

Functional analysis of human MLH1 and MSH2 missense variants and hybrid human-yeast MLH1 proteins in Saccharomyces cerevisiae

Hereditary non-polyposis colorectal cancer (HNPCC) is an autosomal dominant inherited disease caused by defects in the process of DNA mismatch repair (MMR), and mutations in the hMLH1 or hMSH2 genes are responsible for the majority of HNPCC. In addition to clear loss-of-function mutations conferred by nonsense or frameshift alterations in the coding sequence or by splice variants, genetic screening has revealed a large number of missense codons with less obvious functional consequences. The ability to discriminate between a loss-of-function mutation and a silent polymorphism is important for genetic testing for inherited diseases like HNPCC where the opportunity exists for early diagnosis and preventive intervention. In this study, quantitative in vivo DNA MMR assays in the yeast Saccharomyces cerevisiae were performed to determine the functional significance of amino acid replacements observed in the human population. Missense codons previously observed in human genes were introduced at the homologous residue in the yeast MLH1 or MSH2 genes. This study also demonstrated feasibility of constructing genes that encode functional hybrid human-yeast MLH1 proteins. Three classes of missense codons were found: (i) complete loss of function, i.e. mutations; (ii) variants indistinguishable from wild-type protein, i.e. silent polymorphisms; and (iii) functional variants which support MMR at reduced efficiency, i.e. efficiency polymorphisms. There was a good correlation between the functional results in yeast and available human clinical data regarding penetrance of the missense codon. The results reported here raise the intriguing possibility that differences in the efficiency of DNA MMR exist between individuals in the human population due to common polymorphisms.     

A gene-to-patient approach uplifts novel disease gene discovery and identifies 18 putative novel disease genes

PurposeExome and genome sequencing have drastically accelerated novel disease gene discoveries. However, discovery is still hindered by myriad variants of uncertain significance found in genes of undetermined biological function. This necessitates intensive functional experiments on genes of equal predicted causality, leading to a major bottleneck.MethodsWe apply the loss-of-function observed/expected upper-bound fraction metric of intolerance to gene inactivation to curate a list of predicted haploinsufficient disease genes. Using data from the 100,000 Genomes Project, we adopt a gene-to-patient approach that matches de novo loss-of-function variants in constrained genes to patients with rare disease. Through large-scale aggregation of data, we reduce excess analytical noise currently hindering novel discoveries.ResultsResults from 13,949 trios revealed 643 rare, de novo predicted loss-of-function events filtered from 1044 loss-of-function observed/expected upper-bound fraction–constrained genes. A total of 168 variants occurred within 126 genes without a known disease-gene relationship. Of these, 27 genes had >1 kindred affected, and for 18 of these genes, multiple kindreds had overlapping phenotypes. Two years after initial analysis, 11 of 18 (61%) of these genes have been independently published as novel disease gene discoveries.ConclusionUsing large cohorts and adopting gene-based approaches can rapidly and objectively accelerate dominantly inherited novel gene discovery by targeting the most appropriate genes for functional validation.     

Nonrandom occurrence of multiple de novo coding variants in a proband indicates the existence of an oligogenic model in autism

PurposeElucidating the genetic architecture underlying autism spectrum disorder (ASD) will aid in the understanding of its genetic etiology and clinical diagnosis.MethodsA comprehensive set of coding de novo variants (DNVs) from 4504 trios with ASD and 3012 control/sibling trios from several large-scale sequencing studies were collected and combined. Multiple in-depth analyses including DNVs burden, clinical phenotypes, and functional networks underlying the combined data set were used to evaluate the nonrandom occurrence of multiple extreme DNVs (loss-of-function and damaging missense variants) in the same patients.ResultsWe observed a significant excess of multiple extreme DNVs among patients with ASD compared with controls. Meanwhile, patients with ASD carrying 2+ extreme DNVs had significantly lower IQs than patients carrying 0 or 1 DNV. Moreover, much closer functional connectivity than expected was observed among 2 or more genes with extreme DNVs from the same individuals. In particular, we identified 56 key genes as more confident ASD genes compared with other known ASD genes. In addition, we detected 23 new ASD candidate genes with recurrent DNVs, including VIP, ZWILCH, MSL2, LRRC4, and CAPRIN1.ConclusionsOur findings present compelling statistical evidence supporting an oligogenic model and provide new insights into the genetic architecture of ASD.     

Analysis of synonymous codon usage in the UL24 gene of duck enteritis virus

The analysis on codon usage bias of UL24 gene of duck enteritis virus (DEV) may improve our understanding of the evolution and pathogenesis of DEV and provide a basis for understanding the relevant mechanism for biased usage of synonymous codons and for selecting appropriate expression systems to improve the expression of target genes. The codon usage bias of UL24 genes of DEV and 27 reference herpesviruses were analyzed. The results showed that codon of UL24 gene of DEV was strong bias toward the synonymous codons with A and T at the third codon position. A high level of diversity in codon usage bias existed, and the effective number of codons used in a gene plot revealed that the genetic heterogeneity in UL24 gene of herpesviruses was constrained by the G + C content. The phylogentic analysis suggested that DEV was evolutionarily closer to Alphaherpesvirinae and that there was no significant deviation in codon usage in different virus strains. There were 20 codons showing distinct usage differences between DEV and Escherichia coli, 23 between DEV and Homo sapiens, but only 16 codons between DEV and yeast. Therefore the yeast expression system may be more suitable for the expression of DEV genes. Renyong Jia, Anchun Cheng, Mingshu Wang, and Hongyi Xin contributed equally to this work and should be considered as first author.     

Clinical validity assessment of genes frequently tested on intellectual disability/autism sequencing panels

PurposeNeurodevelopmental disorders (NDDs), such as intellectual disability (ID) and autism spectrum disorder (ASD), exhibit genetic and phenotypic heterogeneity, making them difficult to differentiate without a molecular diagnosis. The Clinical Genome Resource Intellectual Disability/Autism Gene Curation Expert Panel (GCEP) uses systematic curation to distinguish ID/ASD genes that are appropriate for clinical testing (ie, with substantial evidence supporting their relationship to disease) from those that are not.MethodsUsing the Clinical Genome Resource gene–disease validity curation framework, the ID/Autism GCEP classified genes frequently included on clinical ID/ASD testing panels as Definitive, Strong, Moderate, Limited, Disputed, Refuted, or No Known Disease Relationship.ResultsAs of September 2021, 156 gene–disease pairs have been evaluated. Although most (75%) were determined to have definitive roles in NDDs, 22 (14%) genes evaluated had either Limited or Disputed evidence. Such genes are currently not recommended for use in clinical testing owing to the limited ability to assess the effect of identified variants.ConclusionOur understanding of gene–disease relationships evolves over time; new relationships are discovered and previously-held conclusions may be questioned. Without periodic re-examination, inaccurate gene–disease claims may be perpetuated. The ID/Autism GCEP will continue to evaluate these claims to improve diagnosis and clinical care for NDDs.     

Codon usage bias is correlated with gene expression levels in the fission yeast Schizosaccharomyces pombe

Usage of synonymous codons represents a characteristic pattern of preference in each organism. It has been inferred that such bias of codon usage has evolved as a result of adaptation for efficient synthesis of proteins. Here we examined synonymous codon usage in genes of the fission yeast Schizosaccharomyces pombe , and compared codon usage bias with expression levels of the gene. In this organism, synonymous codon usage bias was correlated with expression levels of the gene; the bias was most obvious in two-codon amino acids. A similar pattern of the codon usage bias was also observed in Saccharomyces cerevisiae , Arabidopsis thaliana and Caenorhabditis elegans , but was not obvious in Oryza sativa , Drosophila melanogaster , Takifugu rubripes and Homo sapiens . As codons of the highly expressed genes have greater influence on translational efficiency than codons of genes expressed at lower levels, it is likely that codon usage in the S. pombe genome has been optimized by translational selection through evolution.     

An extensive analysis on the global codon usage pattern of baculoviruses

Baculovirus-insect cell systems have been widely used over the past decades. However, few studies to date have addressed baculovirus codon usage. In this study, we calculated the effective number of codons (ENC) for all 5,842 ORFs from 42 completely sequenced baculoviruses. The results revealed that most of the baculoviruses lacked strong codon bias (ENC > 35). Exceptions were Lymantria dispar nucleopolyhedrovirus (LdMNPV) and Orgyia pseudotsugata nucleopolyhedrovirus (OpMNPV), which were found to have a strong codon bias (ENC < 35) in 20.9 and 11.8%, respectively, of their total genes. Comparisons of preferred codons based on taxonomic clades showed that the preferred codons were different in different clades, but nine codons (UUU, UAC, UUG, CAC, CAA, AAA, GUG, GAA, and AUU) were preferably adopted by most baculovirus genes. Correspondence analysis showed that the major trend in codon usage variation among all genes significantly correlated with the GC content of sequences. Analyses also suggested that the high condon bias of LdMNPV and OpMNPV were correlated with their high GC%.     

Dominant Leber congenital amaurosis,cone‐rod degeneration,and retinitis pigmentosa caused by mutant versions of the transcription factor CRX

We summarize 18 mutations in the human CRX gene that have been associated with Leber congenital amaurosis (congenital retinal blindness), cone‐rod degeneration, or retinitis pigmentosa. Except for one obviously null allele not definitely associated with a phenotype (a frameshift in codon 9), all CRX mutations appear to be completely penetrant and cause disease in heterozygotes. These dominant alleles fall into two categories. In one group are missense mutations and short, in‐frame deletions; in the second group are frameshift mutations, all of which are in the last exon. All of these dominant mutations are likely to produce stable mRNA that is translated. Mutations in the missense group preferentially affect the conserved homeobox (codons 39–98), and all frameshift mutations leave the homeodomain intact but alter the OTX motif encoded by codons 284–295 at the carboxy terminus. We could not uncover any correlation between type of disease (congenital amaurosis vs. cone‐rod degeneration or retinitis pigmentosa) and the type of mutation (missense vs. frameshift). Four of the 18 mutations (～20%) were de novo mutations, and all of these were found in isolate cases of Leber congenital amaurosis. Dominant CRX mutations have not been associated with mental retardation or developmental delay that has sometimes been found in Leber congenital amaurosis caused by other genes. Implications regarding potential future therapies are discussed. Hum Mutat 18:488–498, 2001. © 2001 Wiley‐Liss, Inc.     

A rapid solubility assay of protein domain misfolding for pathogenicity assessment of rare DNA sequence variants

PurposeDNA sequencing technology has unmasked a vast number of uncharacterized single-nucleotide variants in disease-associated genes, and efficient methods are needed to determine pathogenicity and enable clinical care.MethodsWe report an E. coli–based solubility assay for assessing the effects of variants on protein domain stability for three disease-associated proteins.ResultsFirst, we examined variants in the Kv11.1 channel PAS domain (PASD) associated with inherited long QT syndrome type 2 and found that protein solubility correlated well with reported in vitro protein stabilities. A comprehensive solubility analysis of 56 Kv11.1 PASD variants revealed that disruption of membrane trafficking, the dominant loss-of-function disease mechanism, is largely determined by domain stability. We further validated this assay by using it to identify second-site suppressor PASD variants that improve domain stability and Kv11.1 protein trafficking. Finally, we applied this assay to several cancer-linked P53 tumor suppressor DNA-binding domain and myopathy-linked Lamin A/C Ig-like domain variants, which also correlated well with reported protein stabilities and functional analyses.ConclusionThis simple solubility assay can aid in determining the likelihood of pathogenicity for sequence variants due to protein misfolding in structured domains of disease-associated genes as well as provide insights into the structural basis of disease.     

Retrospective analysis of a clinical exome sequencing cohort reveals the mutational spectrum and identifies candidate disease–associated loci for BAFopathies

PurposeBRG1/BRM-associated factor (BAF) complex is a chromatin remodeling complex that plays a critical role in gene regulation. Defects in the genes encoding BAF subunits lead to BAFopathies, a group of neurodevelopmental disorders with extensive locus and phenotypic heterogeneity.MethodsWe retrospectively analyzed data from 16,243 patients referred for clinical exome sequencing (ES) with a focus on the BAF complex. We applied a genotype-first approach, combining predicted genic constraints to propose candidate BAFopathy genes.ResultsWe identified 127 patients carrying pathogenic variants, likely pathogenic variants, or de novo variants of unknown clinical significance in 11 known BAFopathy genes. Those include 34 patients molecularly diagnosed using ES reanalysis with new gene–disease evidence (n = 21) or variant reclassifications in known BAFopathy genes (n = 13). We also identified de novo or predicted loss-of-function variants in 4 candidate BAFopathy genes, including ACTL6A, BICRA (implicated in Coffin-Siris syndrome during this study), PBRM1, and SMARCC1.ConclusionWe report the mutational spectrum of BAFopathies in an ES cohort. A genotype-driven and pathway-based reanalysis of ES data identified new evidence for candidate genes involved in BAFopathies. Further mechanistic and phenotypic characterization of additional patients are warranted to confirm their roles in human disease and to delineate their associated phenotypic spectrums.     

Codon usage bias of the phosphoprotein gene of spring viraemia of carp virus and high codon adaptation to the host

In this study, we calculated the relative synonymous codon usage (RSCU) value and the effective number of codons (ENC) value to carry out principal component analysis (PCA) and correlation analysis of the codon usage pattern of the phosphoprotein gene (P gene) of spring viraemia of carp virus (SVCV). The synonymous codon usage pattern in P genes is geography-specific, based on PCA analysis. The high correlation between (G + C)_1,2 % and (G + C)₃ % suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage and base components in P genes. At least 40 out of 59 synonymous codons are similarly selected in all functional genes within five complete SVCV genomes, and the hosts based on the RSCU data. These results not only provide insight into variations in the codon usage pattern of SVCV but also may help in understanding the processes governing the evolution of SVCV.     

Caution in interpretation of disease causality for heterozygous loss-of-function variants in the MYH8 gene associated with autosomal dominant disorder

To date, the NM_002472.2(MYH8):c.2021G>A (p.Arg674Gln) missense variant in the MYH8 gene is the only known genetic change in individuals with autosomal dominant trismus-pseudocamptodactyly syndrome with unknown molecular mechanism. Next-generation sequencing (NGS), including targeted gene panels and whole-exome sequencing, is routinely performed in many clinical diagnostic laboratories as standard-of-care testing aimed at identifying disease-causing genomic variants. Whole-exome sequencing has revealed loss-of-function variants in the MYH8 gene. To properly classify the MYH8 loss-of-function variants, we either retrieved them from public databases or retrospectively collected them from individuals genetically tested by custom NGS panels or by whole-exome sequencing and confirmed using Sanger sequencing. We further evaluated the respective clinical presentations of these individuals with the MYH8 loss-of-function variants. Heterozygous loss-of-function variants in the MYH8 gene were detected in 16 individuals without trismus-pseudocamptodactyly syndrome. Four of these 16 individuals had a pathogenic or likely pathogenic variant detected in another gene that could explain their clinical presentation. Moreover, there are ∼100 MYH8 heterozygous protein-truncating and splice site variants in the ExAC database in different populations. Our results, combined with the population data, indicate that loss-of-function variants in the MYH8 gene do not cause autosomal dominant trismus-pseudocamptodactyly syndrome, and the clinical significance of these variants remains unknown at present. This result highlights the importance of considering the molecular mechanism of disease, variants published in the medical literature, and population genomic data for the correct interpretation of loss-of-function variants in genes associated with autosomal dominant diseases.     

Cladistic Analysis of Human Apolipoprotein A4 Polymorphisms in Relation to Quantitative Plasma Lipid Risk Factors of Coronary Heart Disease

Genetic variation in several genes involved in lipid metabolism is known to affect population variation in quantitative lipid risk factor profiles for coronary heart disease (CHD). The apolipoprotein A‐IV gene (APOA4) is one such candidate gene. We genotyped five polymorphisms in the APOA4 gene (codon 127, codon 130, codon347, codon 360 and 3’ VNTR) and investigated their impact on plasma lipid trait levels in three populations comprising 604 U.S. non‐Hispanic Whites (NHWs), 408 U.S. Hispanics and 708 Nigerian Blacks. Cladistic analysis was carried out to identify 5‐site haplotypes that were associated with significant phenotypic differences in each population. The distribution of APOA4 genotypes was significantly different between ethnic groups. The Africans were monomorphic for two of the five sites (codons 130 and 360), but possess a unique 12 bp insertion that was not observed in NHWs and Hispanics. Due to linkage disequilibrium between the sites, only 6 haplotypes were observed in NHWs and Hispanics, and 4 in Africans. Several gender‐and ethnic‐specific associations between genotypes and plasma lipid traits were observed when single sites were used. Several haplotypes were identified by cladistic analysis that may carry functional mutations that affect plasma lipid trait levels.     

Mutation frequencies of EXT1 and EXT2 in 43 Japanese families with hereditary multiple exostoses

Hereditary multiple exostoses (EXT) is an autosomal dominant bone disease characterized by the formation of cartilage‐capped prominences. EXT is genetically heterogeneous with at least four chromosomal loci. Among the four loci, the exostosis type 1 gene (EXT1) and type 2 gene (EXT2) have been cloned. Previous studies have shown that disease‐type‐specific frequency of mutations is different among various ethnic populations. To determine those frequencies in the Japanese, we conducted a large‐scale mutation screening on both genes. In 23 of 43 Japanese families examined, we found 21 different mutations, of which 18 are novel. Seventeen (40%) of the 23 families had a mutation in EXT1 and six (14%) had a mutation in EXT2, suggesting that the former mutations are more frequent than the latter in Japanese EXT families. Of the 17 families with EXT1 mutations, 13 had those causing premature termination of the EXT1 protein and four showed missense mutations, whereas five of the six families with EXT2 mutations had those causing premature termination and one showed missense mutation. Interestingly, all four EXT1 missense mutations occurred in an arginine residue at codon 340 (R340) that is known as a critical site for expression of heparan sulfate glycosaminoglycans, suggesting that the region encompassing the arginine residue may play an important role in the function of the EXT1 protein. These results expand our knowledge of the ethnic difference of EXT and the structure‐function relationship of the EXT genes. © Wiley‐Liss. Inc.     

Prion protein gene (PRNP) polymorphisms in Xinjiang local sheep breeds in China

Summary. Amino acid polymorphisms of the prion protein (PrP) gene (PRNP), particularly those occurring at codons 136, 154, and 171 have a significant influence on scrapie pathogenesis in many sheep breeds. We isolated blood samples from 222 sheep representing the eight main local sheep breeds in the Xinjiang Autonomous Region, the territory with the most abundant local sheep breeds in China, to identify the PRNP polymorphisms and to determine whether these breeds were at risk for developing scrapie. A new PRNP polymorphism encoding either glycine (G) or arginine (R) at codon 85 as well as eight previously reported polymorphisms at codons 101, 112, 127, 141, 146, 154, 171, and 189 in other sheep breeds were detected. Interestingly, the alanine (A)/V polymorphism at codon 136 was not observed in this study, all sheep being homozygous for A at this position. While the previously identified polymorphism of argine (R) or histidine (H) at codon 154 was detected, the H polymorphism was rare (2.25%). Four polymorphisms at codon 171 encoding glutamine (Q), R, H, or lysine (K) were detected. The predominant ARQ allele occurred at a high frequency of 77.48%, suggesting an increased risk of scrapie in Xinjiang Autonomous Region.     

Detection of mutation in isoniazid-resistant Mycobacterium tuberculosis isolates from tuberculosis patients in Belarus

The aim of this study was to investigate the frequency, location and type of katG mutations in Mycobacterium tuberculosis strains isolated from patients in Belarus. Forty two isoniazid-resistant isolates were identified from sputum of 163 patients with active pulmonary tuberculosis. Drug susceptibility testing was determined by using CDC standard conventional proportional method and BACTEC system. Standard PCR method for detection of isoniazid resistance associated mutations was performed by katG gene amplification and DNA sequencing. Most mutations were found in katG gene codons 315, 316 and 309. Four types of mutations were identified in codon 315: AGC-->ACC (n=36) 85%, AGC-->AGG (n=1) 2.3%, AGC-->AAC (n=2) 4.7%, AGC-->GGC (n=1) 2.3%. One type of mutation was found in codon 316: GGC-->AGC (n=18) 41.4%, four types of mutations were detected in codon 309: GGT-->GGT (n=7) 16.1%, GGT-->GCT (n=4) 9.2%, GGT-->GTC (n=3)6.9%, GGT-->GGG (n=1) 2.7%. The highest frequency of mutations sharing between primary and secondary infections was found in codon 315.