Alternative intronic promoters in development and disease

Associations of genome-wide structural variations with phenotypic differences in cross-bred Eurasian pigs

BACKGROUND

During approximately 10,000 years of domestication and selection, a large number of structural variations (SVs) have emerged in the genome of pig breeds, profoundly influencing their phenotypes and the ability to adapt to the local environment. SVs (≥ 50 bp) are widely distributed in the genome, mainly in the form of insertion (INS), mobile element insertion (MEI), deletion (DEL), duplication (DUP), inversion (INV), and translocation (TRA). While studies have investigated the SVs in pig genomes, genome-wide association studies (GWAS)-based on SVs have been rarely conducted.

RESULTS

Here, we obtained a high-quality SV map containing 123,151 SVs from 15 Large White and 15 Min pigs through integrating the power of several SV tools, with 53.95% of the SVs being reported for the first time. These high-quality SVs were used to recover the population genetic structure, confirming the accuracy of genotyping. Potential functional SV loci were then identified based on positional effects and breed stratification. Finally, GWAS were performed for 36 traits by genotyping the screened potential causal loci in the F2 population according to their corresponding genomic positions. We identified a large number of loci involved in 8 carcass traits and 6 skeletal traits on chromosome 7, with FKBP5 containing the most significant SV locus for almost all traits. In addition, we found several significant loci in intramuscular fat, abdominal circumference, heart weight, and liver weight, etc. CONCLUSIONS: We constructed a high-quality SV map using high-coverage sequencing data and then analyzed them by performing GWAS for 25 carcass traits, 7 skeletal traits, and 4 meat quality traits to determine that SVs may affect body size between European and Chinese pig breeds.

HOX-Gene Cluster Organization and Genome Duplications in Fishes and Mammals: Transcript Variant Distribution along the Anterior–Posterior Axis

Hox genes play a crucial role in morphogenesis, especially in anterior–posterior body axis patterning. The organization of Hox clusters in vertebrates is a result of several genome duplications: two rounds of duplication in the ancestors of all vertebrates and a third round that was specific for teleost fishes. Teleostei cluster structure has been significantly modified in the evolutionary processes by Hox gene losses and co-options, while mammals show no such tendency. In mammals, the Hox gene number in a single cluster is stable and generally large, and the numbers are similar to those in the Chondrichthyes. Hox gene alternative splicing activity slightly differs between fishes and mammals. Fishes and mammals have differences in their known alternative splicing activity for Hox gene distribution along the anterior–posterior body axis. The analyzed fish groups—the Coelacanthiformes, Chondrichthyes, and Teleostei—all have higher known alternative mRNA numbers from the anterior and posterior regions, whereas mammals have a more uniform Hox transcript distribution along this axis. In fishes, most Hox transcripts produce functioning proteins, whereas mammals have significantly more known transcripts that do not produce functioning proteins.

Concordant Androgen-Regulated Expression of Divergent Rhox5 Promoters in Sertoli Cells

Concordant transcriptional regulation can generate multiple gene products that collaborate to achieve a common goal. Here we report a case of concordant transcriptional regulation that instead drives a single protein to be produced in the same cell type from divergent promoters. This gene product-the RHOX5 homeobox transcription factor-is translated from 2 different mRNAs with different 5' untranslated regions (UTRs) transcribed from alternative promoters. Despite the fact that these 2 promoters-the proximal promoter (Pp) and the distal promoter (Pd)-exhibit different patterns of tissue-specific activity, share no obvious sequence identity, and depend on distinct transcription factors for expression, they exhibit a remarkably similar expression pattern in the testes. In particular, both depend on androgen signaling for expression in the testes, where they are specifically expressed in Sertoli cells and have a similar stage-specific expression pattern during the seminiferous epithelial cycle. We report evidence for 3 mechanisms that collaborate to drive concordant Pp/Pd expression. First, both promoters have an intrinsic ability to respond to androgen receptor and androgen. Second, the Pp acts as an enhancer to promote androgen-dependent transcription from the Pd. Third, Pd transcription is positively autoregulated by the RHOX5 protein, which is first produced developmentally from the Pp. Together, our data support a model in which the Rhox5 homeobox gene evolved multiple mechanisms to activate both of its promoters in Sertoli cells to produce Rhox5 in an androgen-dependent manner during different phases of spermatogenesis.

Intron retention and its impact on gene expression and protein diversity: A review and a practical guide

Intron retention (IR) occurs when a complete and unspliced intron remains in mature mRNA. An increasing body of literature has demonstrated a major role for IR in numerous biological functions, including several that impact human health and disease. Although experimental technologies used to study other forms of mRNA splicing can also be used to investigate IR, a specialized downstream computational analysis is optimal for IR discovery and analysis. Here we provide a review of IR and its biological implications, as well as a practical guide for how to detect and analyze it. Several methods, including long read third generation direct RNA sequencing, are described. We have developed an R package, FakIR, to facilitate the execution of the bioinformatic tasks recommended in this review and a tutorial on how to fit them to users aims. Additionally, we provide guidelines and experimental protocols to validate IR discovery and to evaluate the potential impact of IR on gene expression and protein output. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Splicing Regulation/Alternative Splicing RNA Methods > RNA Analyses in vitro and In Silico.

Intronic variation of the SOHLH2 gene confers risk to male reproductive impairment

OBJECTIVE

To evaluate whether SOHLH2 intronic variation contributes to the genetic predisposition to male infertility traits, including severe oligospermia (SO) and different nonobstructive azoospermia (NOA) clinical phenotypes.

DESIGN

Genetic association study.

SETTING

Not applicable.

PATIENT(S)

Five hundred five cases (455 infertile patients diagnosed with NOA and 50 with SO) and 1,050 healthy controls from Spain and Portugal.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

Genomic DNA extraction from peripheral blood mononuclear cells, genotyping of the SOHLH2 polymorphisms rs1328626 and rs6563386 using the TaqMan allelic discrimination technology, case-control association analyses using logistic regression models, and exploration of functional annotations in publicly available databases.

RESULT(S)

Evidence of association was observed for both rs6563386 with SO and rs1328626 with unsuccessful sperm retrieval after testicular sperm extraction (TESE-) in the context of NOA. A dominant effect of the minor alleles was suggested in both associations, either when the subset of patients with the manifestation were compared against the control group (rs6563386/SO: P=.021, odds ratio [OR] = 0.51; rs1328626/TESE-: P=.066, OR = 1.46) or against the group of patients without the manifestation (rs6563386/SO: P=.014, OR = 0.46; rs1328626/TESE-: P=.012, OR = 2.43). The haplotype tests suggested a combined effect of both polymorphisms. In silico analyses evidenced that this effect could be due to alteration of the isoform population.

CONCLUSION(S)

Our data suggest that intronic variation of SOHLH2 is associated with spermatogenic failure. The genetic effect is likely caused by different haplotypes of rs6563386 and rs1328626, which may predispose to SO or TESE- depending on the specific allelic combination.

Alternative UNC13D Promoter Encodes a Functional Munc13-4 Isoform Predominantly Expressed in Lymphocytes and Platelets

Autosomal recessive mutations in genes required for cytotoxicity are causative of a life-threatening, early-onset hyperinflammatory syndrome termed familial hemophagocytic lymphohistiocytosis (FHL). Mutations in UNC13D cause FHL type 3. UNC13D encodes Munc13-4, a member of the Unc13 protein family which control SNARE complex formation and vesicle fusion. We have previously identified FHL3-associated mutations in the first intron of UNC13D which control transcription from an alternative transcriptional start site. Using isoform specific antibodies, we demonstrate that this alternative Munc13-4 isoform with a unique N-terminus is preferentially expressed in human lymphocytes and platelets, as compared to the conventional isoform that was mostly expressed in monocytes and neutrophils. The distinct N-terminal of the two isoforms did not impact on Munc13-4 localization or trafficking to the immunological synapse of cytotoxic T cells. Moreover, ectopic expression of both isoforms efficiently restored exocytosis by FHL3 patient-derived Munc13-4 deficient T cells. Thus, we demonstrate that the conventional and alternative Munc13-4 isoforms have different expression pattern in hematopoietic cell subsets, but display similar localization and contribution to T cell exocytosis. The use of an alternative transcriptional starting site (TSS) in lymphocytes and platelets could be selected for increasing the overall levels of Munc13-4 expression for efficient secretory granule release.

CRISPR/Cas9-based editing of a sensitive transcriptional regulatory element to achieve cell type-specific knockdown of the NEMO scaffold protein

The use of alternative promoters for the cell type-specific expression of a given mRNA/protein is a common cell strategy. NEMO is a scaffold protein required for canonical NF-κB signaling. Transcription of the NEMO gene is primarily controlled by two promoters: one (promoter B) drives NEMO transcription in most cell types and the second (promoter D) is largely responsible for NEMO transcription in liver cells. Herein, we have used a CRISPR/Cas9-based approach to disrupt a core sequence element of promoter B, and this genetic editing essentially eliminates expression of NEMO mRNA and protein in 293T human kidney cells. By cell subcloning, we have isolated targeted 293T cell lines that express no detectable NEMO protein, have defined genomic alterations at promoter B, and do not support activation of canonical NF-κB signaling in response to treatment with tumor necrosis factor. Nevertheless, non-canonical NF-κB signaling is intact in these NEMO-deficient cells. Expression of ectopic wild-type NEMO, but not certain human NEMO disease mutants, in the edited cells restores downstream NF-κB signaling in response to tumor necrosis factor. Targeting of the promoter B element does not substantially reduce NEMO expression (from promoter D) in the human SNU-423 liver cancer cell line. Thus, we have created a strategy for selectively eliminating cell type-specific expression from an alternative promoter and have generated 293T cell lines with a functional knockout of NEMO. The implications of these findings for further studies and for therapeutic approaches to target canonical NF-κB signaling are discussed.

Identification of a Novel Role for Foxo3 Isoform2 in Osteoclastic Inhibition

Foxo3 acts as an important central regulator that integrates signaling pathways and coordinates cellular responses to environmental changes. Recent studies show the involvement of Foxo3 in osteoclastogenesis and rheumatoid arthritis, which prompted us to further investigate the FOXO3 locus. Several databases document FOXO3 isoform2, an N-terminal truncated mutation of the full-length FOXO3 However, the biological function of FOXO3 isoform2 is unclear. In this study, we established a conditional allele of Foxo3 in mice that deletes the full-length Foxo3 except isoform2, a close ortholog of the human FOXO3 isoform2. Expression of Foxo3 isoform2 specifically in macrophage/osteoclast lineage suppresses osteoclastogenesis and leads to the osteopetrotic phenotype in mice. Mechanistically, Foxo3 isoform2 enhances the expression of type I IFN response genes to RANKL stimulation and thus inhibits osteoclastogenesis via endogenous IFN-β-mediated feedback inhibition. Our findings identify, to our knowledge, the first known biological function of Foxo3 isoform2 that acts as a novel osteoclastic inhibitor in bone remodeling.

KDM2A/B lysine demethylases and their alternative isoforms in development and disease

Aberrant levels of histone modifications lead to chromatin malfunctioning and consequently to various developmental defects and human diseases. Therefore, the proteins bearing the ability to modify histones have been extensively studied and the molecular mechanisms of their action are now fairly well understood. However, little attention has been paid to naturally occurring alternative isoforms of chromatin modifying proteins and to their biological roles. In this review, we focus on mammalian KDM2A and KDM2B, the only two lysine demethylases whose genes have been described to produce also an alternative isoform lacking the N-terminal demethylase domain. These short KDM2A/B-SF isoforms arise through alternative promoter usage and seem to play important roles in development and disease. We hypothesise about the biological significance of these alternative isoforms, which might represent a more common evolutionarily conserved regulatory mechanism.

Revealing the alternative promoter usage of SAF/MAZ gene by bichromatic fluorescent reporter construct

The large-scale identification of putative alternative promoters study shows more than 52% of human genes are regulated by alternative promoters. The human myc-associated zinc finger protein (SAF/MAZ) gene have SAF-1 and SAF-3 variants transcripted from two transcription start sites (TSSs). By using SAF/MAZ promoter as a model, we set up an approach to probe how the alternative promoters are regulated in real time. We have constructed the bichromatic fluorescent reporter driven by SAF/MAZ 5'-proximal promoter plasmids from which transactivation status of SAF-1 and SAF-3 alternative promoter could be monitored by EGFP and DsRed expression respectively. The results showed that the SAF-3 expression is regulated by alternative promoters. When the bichromatic fluorescent reporter was driven by -1692/+277 or -1401/+277 SAF/MAZ promoter the dominant expression of SAF-3 would be observed in comparison with SAF-1 expression. We also identified that Elk-1 is an inhibitory transcription factor for SAF-3 expression. The temporal diversity of SAF-1 and SAF-3 expressions can be observed via bichromatic fluorescent reporters. These imply that the bichromatic fluorescent reporter driven by alternative promoter construct might be a useful tool for decoding the temporal regulatory repertoire of alternative promoter in human genes.

A demethylation deficient isoform of the lysine demethylase KDM2A interacts with pericentromeric heterochromatin in an HP1a-dependent manner

Histone modifications have a profound impact on the chromatin structure and gene expression and their correct establishment and recognition is essential for correct cell functioning. Malfunction of histone modifying proteins is associated with developmental defects and diseases and detailed characterization of these proteins is therefore very important. The lysine specific demethylase KDM2A is a CpG island binding protein that has been studied predominantly for its ability to regulate CpG island-associated gene promoters by demethylating their H3K36me2. However, very little attention has been paid to the alternative KDM2A isoform that lacks the N-terminal demethylation domain, KDM2A-SF. Here we characterized KDM2A-SF more in detail and we found that, unlike the canonical full length KDM2A-LF isoform, KDM2A-SF forms distinct nuclear heterochromatic bodies in an HP1a dependent manner. Our chromatin immunoprecipitation experiments further showed that KDM2A binds to transcriptionally silent pericentromeric regions that exhibit high levels of H3K36me2. H3K36me2 is the substrate of the KDM2A demethylation activity and the high levels of this histone modification in the KDM2A-bound pericentromeric regions imply that these regions are occupied by the demethylation deficient KDM2A-SF isoform.