Search

Type II transmembrane serine proteases (TTSPs) are a family of trypsin-like membrane-anchored serine proteases that play key roles in the regulation of some crucial processes in physiological conditions, including cardiac function, digestion, cellular iron homeostasis, epidermal differentiation, and immune responses. However, some of them, in particular TTSPs expressed in the human airways, were identified as host factors that promote the proteolytic activation and spread of respiratory viruses such as influenza virus, human metapneumovirus, and coronaviruses, including SARS-CoV-2. Given their involvement in viral priming, we hypothesized that members of the TTSP family may represent targets of positive selection, possibly as the result of virus-driven pressure. Thus, we investigated the evolutionary history of sixteen TTSP genes in mammals. Evolutionary analyses indicate that most of the TTSP genes that have a verified role in viral proteolytic activation present signals of pervasive positive selection, suggesting that viral infections represent a selective pressure driving the evolution of these proteases. We also evaluated genetic diversity in human populations and we identified targets of balancing selection in TMPRSS2 and TMPRSS4. This scenario may be the result of an ancestral and still ongoing host–pathogen arms race. Overall, our results provide evolutionary information about candidate functional sites and polymorphic positions in TTSP genes. Supplementary Information The online version contains supplementary material available at 10.1007/s00439-022-02435-y.

The molecular mechanisms that direct splice-site selection and assure orderly exon juxtaposition have not been fully clarified. The extraordinary nature of the dystrophin gene points to several hurdles in the processing of transcripts. In this study, dystrophin statistical and thermodynamic splicing parameters have been evaluated providing the first comprehensive description for a single human gene. We show that concomitant use of consensus values (CV) and DeltaDG degrees 37 values for U1snRNA annealing better discriminates between real donor sites and donor-like sequences. Evidence is also provided that, on average, out-of-frame dystrophin exons have significantly stronger CVs and more favorable DeltaDG degrees 37 values; this feature has never been reported and might reflect evolutionary-driven minimization of out-of-frame exon misplicing. Dystrophin splicing mutations have been reported to determine either Duchenne or Becker Muscular Dystrophy, but no comprehensive genotypic/phenotypic correlation has ever been investigated. We have analyzed splicing affecting single base-pair substitutions in the dystrophin gene with respect to their effect on splicing parameters; functional and clinical consequences are also reported. We have found 5'-splice-site mutation occurrence to be statistically related to mutability quotients and propose the use of DeltaDG degrees 37 values as a more effective tool than CV alone to describe donor site mutation consequences. Our analysis also indicates a nearly 100% correlation between clinical phenotype and the reading-frame rule determined at the RNA level. We consider that elucidation of the relative importance of splicing determinants might help to clarify the molecular mechanisms that direct correct splicing in complex genes and might be useful in the validation of predictive models.

Contrasting results have been reported concerning the association of a splice-site polymorphism (rs10774671) in OAS1 with multiple sclerosis (MS). We analysed two OAS1 regions encompassing alternatively spliced exons. While the region carrying the splice-site variant is neutrally evolving, a signature of long-standing balancing selection was observed across an alternative exon 7. Analysis of variants in this exon identified an insertion/deletion polymorphism (rs11352835, A/−) that originates predicted products with distinct C termini. This variant is located along the major branch of the haplotype genealogy, suggesting that it may represent the selection target. A case/control study for MS indicated that rs11352835 is associated with disease susceptibility (for an allelic model with the deleted allele predisposing to MS, OR 1.27, 95% CI 1.072–1.513, p = 0.010). No association was found between rs10774671 and MS. As the two SNPs are in linkage disequilibrium in Europeans, the previously reported association between rs10774671 and MS susceptibility might be driven by rs11352835, possibly explaining the contrasting results previously observed for the splice-site polymorphism. Thus, we describe a novel susceptibility variant for MS in OAS1 and show that population genetic analyses can be instrumental to the identification of selection targets and, consequently, of functional polymorphisms with an effect on phenotypic traits. Electronic supplementary material The online version of this article (doi:10.1007/s00439-011-1053-2) contains supplementary material, which is available to authorized users.

The human TRIM5 genes encodes a retroviral restriction factor (TRIM5α). Evolutionary analyses of this gene in mammals have revealed a complex and multifaceted scenario, suggesting that TRIM5 has been the target of exceptionally strong selective pressures, possibly exerted by recurrent waves of retroviral infections. TRIM5 displays inter-individual expression variability in humans and high levels of TRIM5 mRNA have been associated with a reduced risk of HIV-1 infection. We resequenced TRIM5 in chimpanzees and identified two polymorphisms in intron 1 that are shared with humans. Analysis of the gene region encompassing the two trans-specific variants in human populations identified exceptional nucleotide diversity levels and an excess of polymorphism compared to fixed divergence. Most tests rejected the null hypothesis of neutral evolution for this region and haplotype analysis revealed the presence of two deeply separated clades. Calculation of the time to the most recent common ancestor (TMRCA) for TRIM5 haplotypes yielded estimates ranging between 4 and 7 million years. Overall, these data indicate that long-term balancing selection, an extremely rare process outside MHC genes, has maintained trans-specific polymorphisms in the first intron of TRIM5. Bioinformatic analyses indicated that variants in intron 1 may affect transcription factor-binding sites and, therefore, TRIM5 transcriptional activity. Data herein confirm an extremely complex evolutionary history of TRIM5 genes in primates and open the possibility that regulatory variants in the gene modulate the susceptibility to HIV-1.

Duchenne and Becker Muscular Dystrophy (DMD and BMD) are caused, in the majority of cases, by deletions in the dystrophin gene ( DMD). Here we describe the unprecedented case of a BMD patient carrying a large out-of-frame intragenic deletion, together with an inversion in the DMD gene, resulting in the inclusion of a novel exon in the transcript. Multiplex PCR amplification revealed the presence of a 48-52 exon deletion, but transcript analysis identified two unexpected products, neither of them including exon 53. The shorter mRNA derived from the juxtaposition of exons 47-54 (in-frame), while the longer one resulted from the inclusion of a novel 73-bp exon between exons 47 and 54. Sequence analysis revealed that the inserted sequence derived from an inverted portion of intron 53; its inclusion is predicted to determine protein truncation. The presence of a genomic inversion involving exon 53 and flanking regions was confirmed, and inversion/deletion breakpoints were sequenced. The inverted 73-bp sequence displays splicing signals at both ends and thus it is probably recognized as a novel exon when the partially inverted hnRNA is processed. These findings highlight the importance of mRNA analysis on patients that, based on routine DNA screenings, do not follow the reading-frame rule. This is the first reported patient carrying both an intragenic deletion and inversion in the DMD locus. This case might provide further insight into both the mechanisms that determine genomic rearrangements in the DMD locus and the molecular signals that drive exon inclusion.

Haptoglobin, coded by the HP gene, is a plasma protein that acts as a scavenger for free heme, and haptoglobin-related protein (coded by the HPR gene) forms part of the trypanolytic factor TLF-1, together with apolipoprotein L1 (ApoL1). We analyse the polymorphic small intragenic duplication of the HP gene, with alleles Hp1 and Hp2, in 52 populations, and find no evidence for natural selection either from extended haplotype analysis or from correlation with pathogen richness matrices. Using fiber-FISH, the paralog ratio test, and array-CGH data, we also confirm that the HPR gene is copy number variable, with duplication of the whole HPR gene at polymorphic frequencies in west and central Africa, up to an allele frequency of 15 %. The geographical distribution of the HPR duplication allele overlaps the region where the pathogen causing chronic human African trypanosomiasis, Trypanosoma brucei gambiense, is endemic. The HPR duplication has occurred on one SNP haplotype, but there is no strong evidence of extended homozygosity, a characteristic of recent natural selection. The HPR duplication shows a slight, non-significant undertransmission to human African trypanosomiasis-affected children of unaffected parents in the Democratic Republic of Congo. However, taken together with alleles of APOL1, there is an overall significant undertransmission of putative protective alleles to human African trypanosomiasis-affected children. Electronic supplementary material The online version of this article (doi:10.1007/s00439-013-1352-x) contains supplementary material, which is available to authorized users.