Sequence Conservation Analysis Research Articles

Cryo-EM structure of the CDK2-cyclin A-CDC25A complex

The cell division cycle 25 phosphatases CDC25A, B and C regulate cell cycle transitions by dephosphorylating residues in the conserved glycine-rich loop of CDKs to activate their activity. Here, we present the cryo-EM structure of CDK2-cyclin A in complex with CDC25A at 2.7 Å resolution, providing a detailed structural analysis of the overall complex architecture and key protein-protein interactions that underpin this 86 kDa complex. We further identify a CDC25A C-terminal helix that is critical for complex formation. Sequence conservation analysis suggests CDK1/2-cyclin A, CDK1-cyclin B and CDK2/3-cyclin E are suitable binding partners for CDC25A, whilst CDK4/6-cyclin D complexes appear unlikely substrates. A comparative structural analysis of CDK-containing complexes also confirms the functional importance of the conserved CDK1/2 GDSEID motif. This structure improves our understanding of the roles of CDC25 phosphatases in CDK regulation and may inform the development of CDC25-targeting anticancer strategies.

Computer‐Aided Protein Surface Modification Strategy to Improve the Thermostability of α‐Amylase

AbstractStability under high‐temperature environments is crucial for amylase to function in starch‐based industries. This study develops a method of modifying protein surfaces by combining computer‐aided tools including FoldX, PoPMuSiC, Discovery Studio, and I‐Mutant 2.0 with conserved sequence analysis. A truncated α‐amylase ∆AmyPTG from Parageobacillus thermoglucosidasius DSMZ 2542 is rationally designed through this method to improve its thermostability. Seven single‐site variants are constructed and five of them displayed enhanced thermostability. Next, three double‐site variants are constructed with one particularly successful double‐site variant N31RT213R, exhibiting a 4.3‐fold longer half‐life at 80 °C. Notably, the specific activity of N31RT213R reaches 10 567.16 U mg−1, higher than ∆AmyPTG (6645.43 U mg−1). When applied to the corn starch liquefaction reaction, the mutant N31RT213R gets a higher yield of product concentration of about 255.70 µg mL−1, compared to 190.72 µg mL−1 for ∆AmyPTG. Intramolecular forces analysis and surface electrostatic charges analysis are conducted to determine possible causes for the improvement. Also, molecular dynamics simulation is used to analyze the flexibility shifts of the entire protein. This innovative rational engineering approach has proven to be a successful strategy for the selection of hot spots for protein thermostability evolution and has the potential to be applied to other enzymes.

KDM4B mutations in human cancers

Homologous recombination (HR) is essential for repair of DNA double-strand breaks (DSBs) and restart of stalled or collapsed replication forks. Most cancers are characterized by mutations in components of the DSB repair pathways. Redundant DSB repair pathways exist in eukaryotes from yeast to humans and recent evidence has shown that complete loss of HR function appears to be lethal. Recent evidence has also shown that cancer cells with mutations in one DSB repair pathway can be killed by inhibiting one or more parallel pathways, a strategy that is currently aggressively explored as a cancer therapy. KDM4B is a histone demethylase with pleiotropic functions, which participates in preparing DSBs for repair by contributing to chromatin remodeling. In this report we carried out a pan-cancer analysis of KDM4B mutations with the goal of understanding their distribution and interaction with other DSB genes. We find that although KDM4B mutations co-occur with DSB repair genes, most KDM4B mutations are not drivers or pathogenic. A sequence conservation analysis from yeast to humans shows that highly conserved residues are resistant to mutation. Finally, all mutations occur in a heterozygous state. A single mutation, R986L, was predicted to significantly affect protein structure using computational modeling. This analysis suggests that KDM4B makes contributions to DSB repair but is not a key player.

Genotyping Hepatitis B virus by Next-Generation Sequencing: Detection of Mixed Infections and Analysis of Sequence Conservation.

Our aim was to develop an accurate, highly sensitive method for HBV genotype determination and detection of genotype mixtures. We examined the preS and 5' end of the HBV X gene (5X) regions of the HBV genome using next-generation sequencing (NGS). The 1852 haplotypes obtained were subjected to genotyping via the Distance-Based discrimination method (DB Rule) using two sets of 95 reference sequences of genotypes A-H. In clinical samples from 125 patients, the main genotypes were A, D, F and H in Caucasian, B and C in Asian and A and E in Sub-Saharan patients. Genotype mixtures were identified in 28 (22.40%) cases, and potential intergenotypic recombination was observed in 29 (23.20%) cases. Furthermore, we evaluated sequence conservation among haplotypes classified into genotypes A, C, D, and E by computing the information content. The preS haplotypes exhibited limited shared conserved regions, whereas the 5X haplotypes revealed two groups of conserved regions across the genotypes assessed. In conclusion, we developed an NGS-based HBV genotyping method utilizing the DB Rule for genotype classification. We identified two regions conserved across different genotypes at 5X, offering promising targets for RNA interference-based antiviral therapies.

A novel enrofloxacin-degrading fungus, Humicola sp. KC0924g, isolated from the rhizosphere sediment of the submerged macrophyte Vallisneria spiralis L.

A novel enrofloxacin-degrading fungus was isolated from a rhizosphere sediment of the submerged macrophyte Vallisneria spiralis L.. The isolate, designated KC0924g, was identified as a member of the genus Humicola based on morphological characteristics and tandem conserved sequence analysis. The optimal temperature and pH for enrofloxacin degradation by strain KC0924g were 28°C and 9.0, respectively. Under such condition, 98.2% of enrofloxacin with an initial concentration of 1mg L-1 was degraded after 72h of incubation, with nine possible degradation products identified. Four different metabolic pathways were proposed, which were initiated by cleavage of the piperazine moiety, hydroxylation of the aromatic ring, oxidative decarboxylation, or defluorination. In addition to enrofloxacin, strain KC0924g also degraded other fluoroquinolone antibiotics (ciprofloxacin, norfloxacin, and ofloxacin), malachite green (an illegal additive in aquaculture), and leucomalachite green. Pretreatment of cells of strain KC0924g with Cu2+ accelerated ENR degradation. Furthermore, it was speculated that a flavin-dependent monooxygenase was involved in ENR degradation, based on the increased transcriptional levels of these two genes after Cu2+ induction. This work enriches strain resources for enrofloxacin remediation and, more importantly, would facilitate studies on the molecular mechanism of ENR degradation with degradation-related transcriptome available.

JAG1 Variants Confer Genetic Susceptibility to Thyroid Dysgenesis and Thyroid Dyshormonogenesis in 813 Congenital Hypothyroidism in China.

Congenital hypothyroidism (CH) is indeed a prevalent neonatal endocrine disorder, affecting approximately 1 in 2000-3000 newborns worldwide, and 1 in 2400 newborns in China. Despite its high incidence, the genetic causes of CH, particularly those related to thyroid dysgenesis (TD), are still not well understood. However, previous studies have suggested that JAG1 may be a potential susceptibility gene for congenital thyroid defects. To explore the association between JAG1 and CH, we screened JAG1 variants in a large cohort of 813 CH patients. We performed genetic analysis of JAG1 using next-generation sequencing in 813 CH cases. The pathogenicity of the variants was assessed by bioinformatics softwares, protein sequence conservation analysis, and hydrophobic analysis. Further genetic analysis was conducted targeting 20 CH-related genes in these 25 JAG1 variant carriers. We identified 10 pathogenic missense mutations (p.V45L, p.V272I, p.P552L, p.G610E, p.G852D, p.A891T, p.E1030K, p.R1060W, p.A1131T, p.P1174L) carried by 25 patients, the mutation rate of JAG1 in CH was 3.08%. Among these 25 patients, 16 with 1 variant, 6 with 2 variants, and the other 3 with 3 variants. Our findings indicated that JAG1 variants confer genetic susceptibility to both TD and DH, but with different inheritance models. JAG1 variants lead to TD mainly through monogenic model, while for DH cases, both monogenic mechanisms and oligogenic mechanisms play a pivotal role. Oligogenicity may contribute to the disease severity of DH. JAG1 is a shared genetic factor in TD and DH, with a detection rate of 3.08% in Chinese individuals with CH. A comparison between the oligogenic and monogenic groups suggests a gene dosage effect in CH. Patients with the same JAG1 mutation exhibit diverse clinical phenotypes, indicating complex mechanisms underlying phenotypic heterogeneity.

Rational immunosilencing of a promiscuous T-cell epitope in the capsid of an adeno-associated virus.

Owing to the immunogenicity of adeno-associated viruses (AAVs), gene therapies using AAVs face considerable obstacles. Here, by leveraging ex vivo T-cell assays, the prediction of epitope binding to major histocompatibility complex class-II alleles, sequence-conservation analysis in AAV phylogeny and site-directed mutagenesis, we show that the replacement of amino acid residues in a promiscuous and most immunodominant T-cell epitope in the AAV9 capsid with AAV5 sequences abrogates the immune responses of peripheral blood mononuclear cells to the chimaeric vector while preserving its functions, potency, cellular specificity, transduction efficacy and biodistribution. This rational approach to the immunosilencing of capsid epitopes promiscuously binding to T cells may be applied to other AAV vectors and epitope regions.

An NAD+-dependent group Ⅲ alcohol dehydrogenase involved in long-chain alkane degradation in Acinetobacter venetianus RAG-1

Alcohol dehydrogenases (ADHs) are a class of key enzymes responsible for the oxidation of alkyl alcohols in the aerobic alkane metabolic pathway. Currently, the degradation mechanisms of short- and medium-chain alkanes are commonly reported, while those of long-chain alkanes have received less attention. In this work, a putative long-chain ADH was screened from Acinetobacter venetianus RAG-1 via RNA-seq with n-octacosane (C28) as the sole carbon source. Conserved sequence analysis revealed that it is a group III (Fe-containing/activated) ADH, which is widespread in the genus Acinetobacter. The deletion of adhA led to a significant reduction in the degradation of C28. AdhA exhibited optimal oxidative activity at pH 8.0 and 50 °C with NAD+ as coenzyme, while showing better tolerability to chemical reagents. Enzyme activity assay showed that AdhA owed the oxidative activity to a wide range of substrates including alkyl alcohols (C1-C32) and some isomeric alcohols, such as isopropanol, isobutanol, isoamyl alcohol, and propanetriol, and could reduce the alkyl aldehyde (C1-C12). Meanwhile, the binding of AdhA to different alkyl alcohols was mediated by different amino acids. AdhA is an ADH with an extremely broad substrate utilization range and excellent biochemical characteristics. These results provided important insights in the subsequent investigation of long-chain alkane degradation and petroleum pollution bioremediation.

Comparative analysis of the physical properties of murine and human S100A7: Insight into why zinc piracy is mediated by human but not murine S100A7

S100 proteins are a subfamily of EF-hand calcium-binding proteins found primarily in vertebrate animals. They are distinguished by binding of transition metals and functioning in both the intracellular and extracellular milieu. S100A7 functions in the protection of the skin and mucous membranes and is a biomarker in inflammatory skin disease. A recent study of Neisseria gonorrhoeae infection revealed that human but not murine S100A7 could be used to evade host nutritional immunity. To understand the molecular basis for this difference, we carried out a comparative analysis of the physical and structural properties of human and murine S100A7. The X-ray crystal structure of Ca2+-loaded mouse S100A7 (mS100A7) was determined to 1.69Å resolution, and Ca2+-induced conformational changes were assessed by NMR. Unlike human S100A7 (hS100A7), which exhibits conformational changes in response to binding of Ca2+, no significant changes in mS100A7 were detected. Dynamic light scattering, circular dichroism, and a competition chelator assay were used to compare the Zn2+ affinity and the effects of ion binding on mS100A7 versus hS100A7. Alignment of their sequences revealed a substantial difference in the C-terminal region, which is an important mediator of protein-protein interactions, suggesting a rationale for the specificity of N.gonorrhoeae for hS100A7. These data, along with more detailed analysis of S100A7 sequence conservation across different species, support the proposal that, although hS100A7 is highly conserved in many mammals, the murine protein is a distinct ortholog. Our results highlight the potential limitations of using mouse models for studying bacterial infections in humans.

Molecular insights into titin’s A-band

The thick filament-associated A-band region of titin is a highly repetitive component of the titin chain with important scaffolding properties that support thick filament assembly. It also has a demonstrated link to human disease. Despite its functional significance, it remains a largely uncharacterized part of the titin protein. Here, we have performed an analysis of sequence and structure conservation of A-band titin, with emphasis on poly-FnIII tandem components. Specifically, we have applied multi-dimensional sequence pairwise similarity analysis to FnIII domains and complemented this with the crystallographic elucidation of the 3D-structure of the FnIII-triplet A84-A86 from the fourth long super-repeat in the C-zone (C4). Structural models serve here as templates to map sequence conservation onto super-repeat C4, which we show is a prototypical representative of titin’s C-zone. This templating identifies positionally conserved residue clusters in C super-repeats with the potential of mediating interactions to thick-filament components. Conservation localizes to two super-repeat positions: Ig domains in position 1 and FnIII domains in position 7. The analysis also allows conclusions to be drawn on the conserved architecture of titin’s A-band, as well as revisiting and expanding the evolutionary model of titin’s A-band.

Standing in the canine precision medicine knowledge gap: Improving annotation of canine cancer genomic biomarkers through systematic comparative analysis of human cancer mutations in COSMIC.

The accrual of cancer mutation data and related functional and clinical associations have revolutionised human oncology, enabling the advancement of precision medicine and biomarker-guided clinical management. The catalogue of cancer mutations is also growing in canine cancers. However, without direct high-powered functional data in dogs, it remains challenging to interpret and utilise them in research and clinical settings. It is well-recognised that canine and human cancers share genetic, molecular and phenotypic similarities. Therefore, leveraging the massive wealth of human mutation data may help advance canine oncology. Here, we present a structured analysis of sequence conservation and conversion of human mutations to the canine genome through a 'caninisation' process. We applied this analysis to COSMIC, the Catalogue of Somatic Mutations in Cancer, the most prominent human cancer mutation database. For the project's initial phase, we focused on the subset of the COSMIC data corresponding to Cancer Gene Census (CGC) genes. A total of 670 canine orthologs were found for 721 CGC genes. In these genes, 365 K unique mutations across 160 tumour types were converted successfully to canine coordinates. We identified shared putative cancer-driving mutations, including pathogenic and hotspot mutations and mutations bearing similar biomarker associations with diagnostic, prognostic and therapeutic utility. Thus, this structured caninisation of human cancer mutations facilitates the interpretation and annotation of canine mutations and helps bridge the knowledge gap to enable canine precision medicine.

Understanding and modeling human traits and diseases: Insights from the comparative genomics resources of Zoonomia

Understanding and modeling human traits and diseases: Insights from the comparative genomics resources of Zoonomia

An In Silico Functional Analysis of Non-Synonymous Single-Nucleotide Polymorphisms of Bovine CMAH Gene and Potential Implication in Pathogenesis.

The sugar molecule N-glycolylneuraminic acid (Neu5Gc) is one of the most common sialic acids discovered in mammals. Cytidine monophospho-N-acetylneuraminic acid hydroxylase (CMAH) catalyses the conversion of N-acetylneuraminic acid (Neu5Ac) to Neu5Gc, and it is encoded by the CMAH gene. On the one hand, food metabolic incorporation of Neu5Gc has been linked to specific human diseases. On the other hand, Neu5Gc has been shown to be highly preferred by some pathogens linked to certain bovine diseases. We used various computational techniques to perform an in silico functional analysis of five non-synonymous single-nucleotide polymorphisms (nsSNPs) of the bovine CMAH (bCMAH) gene identified from the 1000 Bull Genomes sequence data. The c.1271C>T (P424L) nsSNP was predicted to be pathogenic based on the consensus result from different computational tools. The nsSNP was also predicted to be critical based on sequence conservation, stability, and post-translational modification site analysis. According to the molecular dynamic simulation and stability analysis, all variations promoted stability of the bCMAH protein, but mutation A210S significantly promoted CMAH stability. In conclusion, c.1271C>T (P424L) is expected to be the most harmful nsSNP among the five detected nsSNPs based on the overall studies. This research could pave the way for more research associating pathogenic nsSNPs in the bCMAH gene with diseases.

Systematic Investigation of TCP Gene Family: Genome-Wide Identification and Light-Regulated Gene Expression Analysis in Pepino (Solanum Muricatum)

Plant-specific transcription factors such as the TCP family play crucial roles in light responses and lateral branching. The commercial development of S. muricatum has been influenced by the ease with which its lateral branches can be germinated, especially under greenhouse cultivation during the winter with supplemented LED light. The present study examined the TCP family genes in S. muricatum using bioinformatics analysis (whole-genome sequencing and RNA-seq) to explore the response of this family to different light treatments. Forty-one TCP genes were identified through a genome-wide search; phylogenetic analysis revealed that the CYC/TB1, CIN and Class I subclusters contained 16 SmTCP, 11 SmTCP and 14 SmTCP proteins, respectively. Structural and conserved sequence analysis of SmTCPs indicated that the motifs in the same subcluster were highly similar in structure and the gene structure of SmTCPs was simpler than that in Arabidopsis thaliana; 40 of the 41 SmTCPs were localized to 12 chromosomes. In S. muricatum, 17 tandem repeat sequences and 17 pairs of SmTCP genes were found. We identified eight TCPs that were significantly differentially expressed (DETCPs) under blue light (B) and red light (R), using RNA-seq. The regulatory network of eight DETCPs was preliminarily constructed. All three subclusters responded to red and blue light treatment. To explore the implications of regulatory TCPs in different light treatments for each species, the TCP regulatory gene networks and GO annotations for A. thaliana and S. muricatum were compared. The regulatory mechanisms suggest that the signaling pathways downstream of the TCPs may be partially conserved between the two species. In addition to the response to light, functional regulation was mostly enriched with auxin response, hypocotyl elongation, and lateral branch genesis. In summary, our findings provide a basis for further analysis of the TCP gene family in other crops and broaden the functional insights into TCP genes regarding light responses.

Systematic Identification and Comparison of the Expressed Profiles of Exosomal MiRNAs in Pigs Infected with NADC30-like PRRSV Strain.

Exosomes are biological vesicles secreted and released by cells that act as mediators of intercellular communication and play a unique role in virus infection, antigen presentation, and suppression/promotion of body immunity. Porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most damaging pathogens in the pig industry and can cause reproductive disorders in sows, respiratory diseases in pigs, reduced growth performance, and other diseases leading to pig mortality. In this study, we used the PRRSV NADC30-like CHsx1401 strain to artificially infect 42-day-old pigs and isolate serum exosomes. Based on high-throughput sequencing technology, 305 miRNAs were identified in serum exosomes before and after infection, among which 33 miRNAs were significantly differentially expressed between groups (13 relatively upregulated and 20 relatively downregulated). Sequence conservation analysis of the CHsx1401 genome identified 8 conserved regions, of which a total of 16 differentially expressed (DE) miRNAs were predicted to bind to the conserved region closest to the 3' UTR of the CHsx1401 genome, including 5 DE miRNAs capable of binding to the CHsx1401 3' UTR (ssc-miR-34c, ssc-miR-375, ssc-miR-378, ssc-miR-486, ssc-miR-6529). Further analysis revealed that the target genes of differentially expressed miRNAs were widely involved in exosomal function-related and innate immunity-related signaling pathways, and 18 DE miRNAs (ssc-miR-4331-3p, ssc-miR-744, ssc-miR-320, ssc-miR-10b, ssc-miR-124a, ssc-miR-128, etc.) associated with PRRSV infection and immunity were screened as potential functional molecules involved in the regulation of PRRSV virus infection by exosomes.

Interpreting the molecular mechanisms of disease variants in human transmembrane proteins

Next-generation sequencing of human genomes reveals millions of missense variants, some of which may lead to loss of protein function and ultimately disease. Here, we investigate missense variants in membrane proteins—key drivers in cell signaling and recognition. We find enrichment of pathogenic variants in the transmembrane region across 19,000 functionally classified variants in human membrane proteins. To accurately predict variant consequences, one fundamentally needs to understand the underlying molecular processes. A key mechanism underlying pathogenicity in missense variants of soluble proteins has been shown to be loss of stability. Membrane proteins, however, are widely understudied. Here, we interpret variant effects on a larger scale by performing structure-based estimations of changes in thermodynamic stability using a membrane-specific energy function and analyses of sequence conservation during evolution of 15 transmembrane proteins. We find evidence for loss of stability being the cause of pathogenicity in more than half of the pathogenic variants, indicating that this is a driving factor also in membrane-protein-associated diseases. Our findings show how computational tools aid in gaining mechanistic insights into variant consequences for membrane proteins. To enable broader analyses of disease-related and population variants, we include variant mappings for the entire human proteome.

Sequence Conservation Analysis and Gene Relationships of Nucleocapsid (N) Gene in Orthocoronavirinae Subfamily

Coronavirus (CoV) is a virus that causes respiratory and gastrointestinal diseases in animals and humans. It belongs to the Orthocoronavirina. The nucleocapsid protein (N) plays multiple roles in virus assembly, RNA transcription, and interaction with host cells. This study aimed to analyse the N protein by identifying conserved residues and exploring the gene and protein relationships within the Orthocoronavirinae. Therefore the results of this study are expected to help identify conserved regions of N protein in SARS-CoV-2 which can be used as probes for the virus identification process and can be used as target areas in vaccine development. We used 159 N gene and protein sequences, including 64 from Alpha, 51 from Beta-, 11 from Delta-, and 20 from Gammacoronavirus genera of the Orthocoronavirinae. Three sequences from Tobaniviridae were used as outgroups. Multiple sequence alignment (MSA) and phylogenetic tree analysis were performed using the neighbour-joining and Maximum Likelihood. The MSA results revealed several conserved residues, ranging from 18 to 41, were located in the N-terminal and Cterminal domains, the linker region, Nuclear Localization Signal (NLS), Nuclear Export Signal (NES) motifs, and Packing Signal (PS) binding sites. The phylogenetic tree analysis indicated that Gammacoronavirus and Deltacoronavirus were closely related to Betacoronavirus, while Alfacoronavirus showed the most distant relationship. Furthermore, the study identified 23 conserved residues involved in RNA binding, including amino acids such as Ser89, Val111, Pro112, Gly124, Tyr125, Phe150, Tyr151, Gly154, Thr155, Gly156, Trp180, Val181, Gly409, Arg411, Asn419, Gly421, and Pro443. These residues interacted with phosphate groups, nitrogenous bases, and pentose sugars and exhibited non-specific interactions with RNA. In summary, this study investigated the N protein in the Orthocoronavirinae subfamily, providing insights into its function, structure, and evolutionary relationships.

Molecular diversity and phenotypic pleiotropy of ancient genomic regulatory loci derived from human endogenous retrovirus type H (HERVH) promoter LTR7 and HERVK promoter LTR5_Hs and their contemporary impacts on pathophysiology of Modern Humans

Timelines of population-level effects of viruses on humans varied from the evolutionary scale of million years to contemporary spread of viral infections. Correspondingly, these events are exemplified by: (i) emergence of human endogenous retroviruses (HERVs) from ancient germline infections leading to stable integration of viral genomes into human chromosomes; and (ii) wide-spread viral infections reaching a global pandemic state such as the COVID-19 pandemic. Despite significant efforts, understanding of HERV’s roles in governance of genomic regulatory networks, their impacts on primate evolution and development of human-specific physiological and pathological phenotypic traits remains limited. Remarkably, present analyses revealed that expression of a dominant majority of genes (1696 of 1944 genes; 87%) constituting high-confidence down-steam regulatory targets of defined HERV loci was significantly altered in cells infected with the SARS-CoV-2 coronavirus, a pathogen causing the global COVID-19 pandemic. This study focused on defined sub-sets of DNA sequences derived from HERVs that are expressed at specific stages of human preimplantation embryogenesis and exert regulatory actions essential for self-renewal and pluripotency. Evolutionary histories of LTR7/HERVH and LTR5_Hs/HERVK were charted based on evidence of the earliest presence and expansion of highly conserved (HC) LTR sequences. Sequence conservation analyses of most recent releases 17 primate species’ genomes revealed that LTR7/HERVH have entered germlines of primates in Africa after the separation of the New World Monkey lineage, while LTR5_Hs/HERVK successfully colonized primates’ germlines after the segregation of Gibbons’ species. Subsequently, both LTR7 and LTR5_Hs undergo a marked ~ fourfold–fivefold expansion in genomes of Great Apes. Timelines of quantitative expansion of both LTR7 and LTR5_Hs loci during evolution of Great Apes appear to replicate the consensus evolutionary sequence of increasing cognitive and behavioral complexities of non-human primates, which seems particularly striking for LTR7 loci and 11 distinct LTR7 subfamilies. Consistent with previous reports, identified in this study, 351 human-specific (HS) insertions of LTR7 (175 loci) and LTR5_Hs (176 loci) regulatory sequences have been linked to genes implicated in establishment and maintenance of naïve and primed pluripotent states and preimplantation embryogenesis phenotypes. Unexpectedly, HS-LTRs manifest regulatory connectivity to genes encoding markers of 12 distinct cells’ populations of fetal gonads, as well as genes implicated in physiology and pathology of human spermatogenesis, including Y-linked spermatogenic failure, oligo- and azoospermia. Granular interrogations of genes linked with 11 distinct LTR7 subfamilies revealed that mammalian offspring survival (MOS) genes seem to remain one of consistent regulatory targets throughout ~ 30 MYA of the divergent evolution of LTR7 loci. Differential GSEA of MOS versus non-MOS genes identified clearly discernable dominant enrichment patterns of phenotypic traits affected by MOS genes linked with LTR7 (562 MOS genes) and LTR5_Hs (126 MOS genes) regulatory loci across the large panel of genomics and proteomics databases reflecting a broad spectrum of human physiological and pathological traits. GSEA of LTR7-linked MOS genes identified more than 2200 significantly enriched records of human common and rare diseases and gene signatures of 466 significantly enriched records of Human Phenotype Ontology traits, including Autosomal Dominant (92 genes) and Autosomal Recessive (93 genes) Inheritance. LTR7 regulatory elements appear linked with genes implicated in functional and morphological features of central nervous system, including synaptic transmission and protein–protein interactions at synapses, as well as gene signatures differentially regulated in cells of distinct neurodevelopmental stages and morphologically diverse cell types residing and functioning in human brain. These include Neural Stem/Precursor cells, Radial Glia cells, Bergman Glia cells, Pyramidal cells, Tanycytes, Immature neurons, Interneurons, Trigeminal neurons, GABAergic neurons, and Glutamatergic neurons. GSEA of LTR7-linked genes identified significantly enriched gene sets encoding markers of more than 80 specialized types of neurons and markers of 521 human brain regions, most prominently, subiculum and dentate gyrus. Identification and characterization of 1944 genes comprising high-confidence down-steam regulatory targets of LTR7 and/or LTR5_Hs loci validated and extended these observations by documenting marked enrichments for genes implicated in neoplasm metastasis, intellectual disability, autism, multiple cancer types, Alzheimer’s, schizophrenia, and other brain disorders. Overall, genes representing down-stream regulatory targets of ancient retroviral LTRs exert the apparently cooperative and exceedingly broad phenotypic impacts on human physiology and pathology. This is exemplified by altered expression of 93% high-confidence LTR targets in cells infected by contemporary viruses, revealing a convergence of virus-inflicted aberrations on genomic regulatory circuitry governed by ancient retroviral LTR elements and interference with human cells’ differentiation programs.

T cell receptor repertoire analysis in HTLV-1-associated diseases.

Human T lymphotropic virus 1 (HTLV-1) is a human retrovirus identified as the causative agent in adult T-cell leukemia/lymphoma (ATL) and chronic-progressive neuroinflammatory disorder HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP). HTLV-1 is estimated to infect between 5-20 million people worldwide, although most infected individuals remain asymptomatic. HTLV-1 infected persons carry an estimated lifetime risk of approximately 5% of developing ATL, and between 0.25% and 1.8% of developing HAM/TSP. Most HTLV-1 infection is detected in CD4+ T cells in vivo which causes the aggressive malignancy in ATL. In HAM/TSP, the increase of HTLV-1 provirus induces immune dysregulation to alter inflammatory milieu, such as expansion of HTLV-1-specific CD8+ T cells, in the central nervous system of the infected subjects, which have been suggested to underlie the pathogenesis of HAM/TSP. Factors contributing to the conversion from asymptomatic carrier to disease state remain poorly understood. As such, the identification and tracking of HTLV-1-specific T cell biomarkers that may be used to monitor the progression from primary infection to immune dysfunction and disease are of great interest. T cell receptor (TCR) repertoires have been extensively investigated as a mechanism of monitoring adaptive T cell immune response to viruses and tumors. Breakthrough technologies such as single-cell RNA sequencing have increased the specificity with which T cell clones may be characterized and continue to improve our understanding of TCR signatures in viral infection, cancer, and associated treatments. In HTLV-1-associated disease, sequencing of TCR repertoires has been used to reveal repertoire patterns, diversity, and clonal expansions of HTLV-1-specific T cells capable of immune evasion and dysregulation in ATL as well as in HAM/TSP. Conserved sequence analysis has further been used to identify CDR3 motif sequences and exploit disease- or patient-specificity and commonality in HTLV-1-associated disease. In this article we review current research on TCR repertoires and HTLV-1-specific clonotypes in HTLV-1-associated diseases ATL and HAM/TSP and discuss the implications of TCR clonal expansions on HTLV-1-associated disease course and treatments.

Aspartyl Protease Inhibitors as Anti-Filarial Drugs.

The current treatments for lymphatic filariasis and onchocerciasis do not effectively kill the adult parasitic nematodes, allowing these chronic and debilitating diseases to persist in millions of people. Thus, the discovery of new drugs with macrofilaricidal potential to treat these filarial diseases is critical. To facilitate this need, we first investigated the effects of three aspartyl protease inhibitors (APIs) that are FDA-approved as HIV antiretroviral drugs on the adult filarial nematode, Brugia malayi and the endosymbiotic bacteria, Wolbachia. From the three hits, nelfinavir had the best potency with an IC50 value of 7.78 µM, followed by ritonavir and lopinavir with IC50 values of 14.3 µM and 16.9 µM, respectively. The three APIs have a direct effect on killing adult B. malayi after 6 days of exposure in vitro and did not affect the Wolbachia titers. Sequence conservation and stage-specific gene expression analysis identified Bm8660 as the most likely primary aspartic protease target for these drug(s). Immunolocalization using antibodies raised against the Bm8660 ortholog of Onchocerca volvulus showed it is strongly expressed in female B. malayi, especially in metabolically active tissues such as lateral and dorsal/ventral chords, hypodermis, and uterus tissue. Global transcriptional response analysis using adult female B. pahangi treated with APIs identified four additional aspartic proteases differentially regulated by the three effective drugs, as well as significant enrichment of various pathways including ubiquitin mediated proteolysis, protein kinases, and MAPK/AMPK/FoxO signaling. In vitro testing against the adult gastro-intestinal nematode Trichuris muris suggested broad-spectrum potential for these APIs. This study suggests that APIs may serve as new leads to be further explored for drug discovery to treat parasitic nematode infections.