Evidence for opposing selective forces operating on human-specific duplicated TCAF genes in Neanderthals and humans

TRP channel-associated factor 1/2 (TCAF1/TCAF2) proteins antagonistically regulate the cold-sensor protein TRPM8 in multiple human tissues. Understanding their significance has been complicated given the locus spans a gap-ridden region with complex segmental duplications in GRCh38. Using long-read sequencing, we sequence-resolve the locus, annotate full-length TCAF models in primate genomes, and show substantial human-specific TCAF copy number variation. We identify two human super haplogroups, H4 and H5, and establish that TCAF duplications originated ~1.7 million years ago but diversified only in Homo sapiens by recurrent structural mutations. Conversely, in all archaic-hominin samples the fixation for a specific H4 haplotype without duplication is likely due to positive selection. Here, our results of TCAF copy number expansion, selection signals in hominins, and differential TCAF2 expression between haplogroups and high TCAF2 and TRPM8 expression in liver and prostate in modern-day humans imply TCAF diversification among hominins potentially in response to cold or dietary adaptations.

Identification and Characterization of DAMs Mutations Associated With Early Blooming in Sweet Cherry, and Validation of DNA-Based Markers for Selection

Dormancy release and bloom time of sweet cherry cultivars depend on the environment and the genotype. The knowledge of these traits is essential for cultivar adaptation to different growing areas, and to ensure fruit set in the current climate change scenario. In this work, the major sweet cherry bloom time QTL qP-BT1.1m (327 Kbs; Chromosome 1) was scanned for candidate genes in the Regina cv genome. Six MADS-box genes (PavDAMs), orthologs to peach and Japanese apricot DAMs, were identified as candidate genes for bloom time regulation. The complete curated genomic structure annotation of these genes is reported. To characterize PavDAMs intra-specific variation, genome sequences of cultivars with contrasting chilling requirements and bloom times (N = 13), were then mapped to the ‘Regina’ genome. A high protein sequence conservation (98.8–100%) was observed. A higher amino acid variability and several structural mutations were identified in the low-chilling and extra-early blooming cv Cristobalina. Specifically, a large deletion (694 bp) upstream of PavDAM1, and various INDELs and SNPs in contiguous PavDAM4 and -5 UTRs were identified. PavDAM1 upstream deletion in ‘Cristobalina’ revealed the absence of several cis-acting motifs, potentially involved in PavDAMs expression. Also, due to this deletion, a non-coding gene expressed in late-blooming ‘Regina’ seems truncated in ‘Cristobalina’. Additionally, PavDAM4 and -5 UTRs mutations revealed different splicing variants between ‘Regina’ and ‘Cristobalina’ PavDAM5. The results indicate that the regulation of PavDAMs expression and post-transcriptional regulation in ‘Cristobalina’ may be altered due to structural mutations in regulatory regions. Previous transcriptomic studies show differential expression of PavDAM genes during dormancy in this cultivar. The results indicate that ‘Cristobalina’ show significant amino acid differences, and structural mutations in PavDAMs, that correlate with low-chilling and early blooming, but the direct implication of these mutations remains to be determined. To complete the work, PCR markers designed for the detection of ‘Cristobalina’ structural mutations in PavDAMs, were validated in an F2 population and a set of cultivars. These PCR markers are useful for marker-assisted selection of early blooming seedlings, and probably low-chilling, from ‘Cristobalina’, which is a unique breeding source for these traits.

A case report of achondrogenesis type II (Langer-Saldino achondrogenesis)

Background Achondrogenesis type II is a rare autosomal dominant skeletal dysplasia with a frequency of ~0.2 per 100,000 births. It is one of the lethal short-limbed dwarfisms associated with structural mutations in type II collagen and is also known as Langer-Saldino achondrogenesis. It is characterized by severe micromelia (shortening of entire limb), narrow chest, and prominent abdomen. It shares the striking feature of partial or complete vertebral body demineralization with achondrogenesis type I. Case presentation We present a case with antenatal diagnosis of this rare entity which was confirmed by post-termination radiographs of abortus. Conclusion The imaging plays a cardinal role in the diagnosis of this condition. This case represents only the 4th case of this rare entity from India.

Defining function of wild-type and patient specific TP53 mutations in a zebrafish model of embryonal rhabdomyosarcoma

In embryonal rhabdomyosarcoma (ERMS) and generally in sarcomas, the role of wild-type and loss or gain of function TP53 mutations remains largely undefined. Eliminating mutant or restoring wild-type p53 is challenging; nevertheless, understanding TP53 effects on tumorigenesis remains central to realizing better treatment outcomes. In ERMS, >70% of patients retain wild-type TP53, yet TP53 mutations when present in tumors are associated with poor prognosis. Employing a kRASG12D-driven ERMS tumor model and newly generated tp53 null (tp53-/-) zebrafish, we define both wild-type and patient-specific TP53 mutant effects on tumorigenesis. We demonstrate that tp53 is a major suppressor of tumor initiation, where tp53 loss expands tumors initiation from <35% to >97% of animals. Next, characterizing three patient-specific mutants finds that TP53C176F partially retains wild-type p53 apoptotic activity that can be exploited, while the TP53P153Δ and TP53Y220C mutants define two structural mutations that predispose to head musculature ERMS.

Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer

Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara is a broad-leaved valerian endemic to Ulleung Island, a noted hot spot of endemism in Korea. However, despite its widespread pharmacological use, this plant remains comparatively understudied. Plant cells generally contain two types of organellar genomes (the plastome and the mitogenome) that have undergone independent evolution, which accordingly can provide valuable information for elucidating the phylogenetic relationships and evolutionary histories of terrestrial plants. Moreover, the extensive mega-data available for plant genomes, particularly those of plastomes, can enable researchers to gain an in-depth understanding of the transfer of genes between different types of genomes. In this study, we analyzed two organellar genomes (the 155,179 bp plastome and the 1,187,459 bp mitogenome) of V. sambucifolia f. dageletiana and detected extensive changes throughout the plastome sequence, including rapid structural mutations associated with inverted repeat (IR) contraction and genetic variation. We also described features characterizing the first reported mitogenome sequence obtained for a plant in the order Dipsacales and confirmed frequent gene transfer in this mitogenome. We identified eight non-plastome-originated regions (NPRs) distributed within the plastome of this endemic plant, for six of which there were no corresponding sequences in the current nucleotide sequence databases. Indeed, one of these unidentified NPRs unexpectedly showed certain similarities to sequences from bony fish. Although this is ostensibly difficult to explain, we suggest that this surprising association may conceivably reflect the occurrence of gene transfer from a bony fish to the plastome of an ancestor of V. sambucifolia f. dageletiana mediated by either fungi or bacteria.

Paradoxical androgen receptor regulation by small molecule enantiomers

Small molecules that target the androgen receptor (AR) are the mainstay of therapy for lethal castration-resistant prostate cancer (CRPC), yet existing drugs lose their efficacy during continued treatment. This evolution of resistance is due to heterogenous mechanisms which include AR mutations causing the identical drug to activate instead of inhibit the receptor. Understanding in molecular detail the paradoxical phenomenon wherein an AR antagonist is transformed into an agonist by structural mutations in the target receptor is thus of paramount importance. Herein, we describe a reciprocal paradox: opposing antagonist and agonist AR regulation determined uniquely by enantiomeric forms of the same drug structure. The antiandrogen BMS-641988, which has (R)-chirality at C-5 encompasses a previously uncharacterized (S)-stereoisomer that is, surprisingly, a potent agonist of AR, as demonstrated by transcriptional assays supported by cell imaging studies. This duality was reproduced in a series of novel compounds derived from the BMS-641988 scaffold. Coupled with in silico modeling studies, the results inform an AR model that explains the switch from potent antagonist to high-affinity agonist in terms of C-5 substituent steric interactions with helix 12 of the ligand binding site. They imply strategies to overcome AR drug resistance and demonstrate that insufficient enantiopurity in this class of AR antagonist can confound efforts to correlate structure with function.

Gene disruption by structural mutations drives selection in US rice breeding over the last century

The genetic basis of general plant vigor is of major interest to food producers, yet the trait is recalcitrant to genetic mapping because of the number of loci involved, their small effects, and linkage. Observations of heterosis in many crops suggests that recessive, malfunctioning versions of genes are a major cause of poor performance, yet we have little information on the mutational spectrum underlying these disruptions. To address this question, we generated a long-read assembly of a tropical japonica rice (Oryza sativa) variety, Carolina Gold, which allowed us to identify structural mutations (>50 bp) and orient them with respect to their ancestral state using the outgroup, Oryza glaberrima. Supporting prior work, we find substantial genome expansion in the sativa branch. While transposable elements (TEs) account for the largest share of size variation, the majority of events are not directly TE-mediated. Tandem duplications are the most common source of insertions and are highly enriched among 50-200bp mutations. To explore the relative impact of various mutational classes on crop fitness, we then track these structural events over the last century of US rice improvement using 101 resequenced varieties. Within this material, a pattern of temporary hybridization between medium and long-grain varieties was followed by recent divergence. During this long-term selection, structural mutations that impact gene exons have been removed at a greater rate than intronic indels and single-nucleotide mutations. These results support the use of ab initio estimates of mutational burden, based on structural data, as an orthogonal predictor in genomic selection.

Trypanothione Reductase Gene Mutations in Meglumine Anti-moniate Resistant Isolates from Cutaneous Leishmaniasis Pa-tients Using Molecular Dynamics Method

Background: In this study, mutations in the tripanothion reductase of Leishmania tropica isolated from Iran was investigated using sequencing and simulation of the enzyme by the molecular dynamic method. Methods: Fifteen susceptible and 15 clinical resistant L. tropica specimens were collected from skin lesions from different regions of Iran in 2017. After DNA extraction, trypanothione reductase (TRYR or TPR), gene fragment was amplified using PCR and sequencing methods. In the case of structural mutations, the components were simulated by molecular dynamics using the GROMACS software. Results: Some structural mutations were observed in 9 amino acids surrounding the active site of the TRYR gene of L. tropica with three-dimensional trypanothione reductase alteration. Conclusion: Change in the active site of TRYR of L. tropica, could probably contribute to the development of resistant L. tropica to glucantime. Because of the likely occurrence of mutations in glucantime as well as the ease of development of L. tropica resistant populations, more samples are needed to demonstrate the relationship between mutations in this enzyme and clinical resistance to glucantime. On the other hand, it is recommended that enzymatic studies be performed to confirm the role of mutation in the function and expression of trypanothione reductase in glucantime resistant and susceptible populations.

KinaseMD: kinase mutations and drug response database

Mutations in kinases are abundant and critical to study signaling pathways and regulatory roles in human disease, especially in cancer. Somatic mutations in kinase genes can affect drug treatment, both sensitivity and resistance, to clinically used kinase inhibitors. Here, we present a newly constructed database, KinaseMD (kinase mutations and drug response), to structurally and functionally annotate kinase mutations. KinaseMD integrates 679 374 somatic mutations, 251 522 network-rewiring events, and 390 460 drug response records curated from various sources for 547 kinases. We uniquely annotate the mutations and kinase inhibitor response in four types of protein substructures (gatekeeper, A-loop, G-loop and αC-helix) that are linked to kinase inhibitor resistance in literature. In addition, we annotate functional mutations that may rewire kinase regulatory network and report four phosphorylation signals (gain, loss, up-regulation and down-regulation). Overall, KinaseMD provides the most updated information on mutations, unique annotations of drug response especially drug resistance and functional sites of kinases. KinaseMD is accessible at https://bioinfo.uth.edu/kmd/, having functions for searching, browsing and downloading data. To our knowledge, there has been no systematic annotation of these structural mutations linking to kinase inhibitor response. In summary, KinaseMD is a centralized database for kinase mutations and drug response.

De novo structural mutation rates and gamete-of-origin biases revealed through genome sequencing of 2,396 families

Each human genome includes de novo mutations that arose during gametogenesis. While these germline mutations represent a fundamental source of new genetic diversity, they can also create deleterious alleles that impact fitness. The germline mutation rate for single nucleotide variants and factors that significantly influence this rate, such as parental age, are now well established. However, far less is known about the frequency, distribution, and features that impact de novo structural mutations. We report a large, family-based study of germline mutations, excluding aneuploidy, that affect genome structure among 572 genomes from 33 families in a multigenerational CEPH-Utah cohort and 2,363 cases of non-familial autism spectrum disorder (ASD), 1,938 unaffected siblings, and both parents (9,599 genomes in total). We find that de novo structural mutations detected by alignment-based, short-read WGS occurred at an overall rate of at least 0.160 events per genome in unaffected individuals and was significantly higher (0.206 per genome) in ASD cases. In both probands and unaffected samples, nearly 73% of de novo structural mutations arose in paternal gametes, and predict most de novo structural mutations to be caused by mutational mechanisms that do not require sequence homology. After multiple testing correction we did not observe a statistically significant correlation between parental age and the rate of de novo structural variation in offspring. These results highlight that a spectrum of mutational mechanisms contribute to germline structural mutations, and that these mechanisms likely have markedly different rates and selective pressures than those leading to point mutations.