Reverse transcriptase meets DNA, again: Possible roles for transposable elements in host DNA repair.

Retrotransposons are selfish genetic elements that encode an enzyme, reverse transcriptase (RT), which converts the element-encoded RNA into DNA prior to or during genomic integration. New studies provide compelling evidence that a bacterial group II intron-like RT has adapted enzymatic activities associated with RTs to function in host DNA repair.

Genome-wide de novo L1 Retrotransposition Connects Endonuclease Activity with Replication.

L1 retrotransposon-derived sequences comprise approximately 17% of the human genome. Darwinian selective pressures alter L1 genomic distributions during evolution, confounding the ability to determine initial L1 integration preferences. Here, we generated high-confidence datasets of greater than 88,000 engineered L1 insertions in human cell lines that act as proxies for cells that accommodate retrotransposition in vivo. Comparing these insertions to a null model, in which L1 endonuclease activity is the sole determinant dictating L1 integration preferences, demonstrated that L1 insertions are not significantly enriched in genes, transcribed regions, or open chromatin. By comparison, we provide compelling evidence that the L1 endonuclease disproportionately cleaves predominant lagging strand DNA replication templates, while lagging strand 3'-hydroxyl groups may prime endonuclease-independent L1 retrotransposition in a Fanconi anemia cell line. Thus, acquisition of an endonuclease domain, in conjunction with the ability to integrate into replicating DNA, allowed L1 to become an autonomous, interspersed retrotransposon.

Poly(ADP-Ribose) Polymerase 2 Recruits Replication Protein A to Sites of LINE-1 Integration to Facilitate Retrotransposition.

Long interspersed element-1 (LINE-1 or L1) retrotransposition poses a threat to genome integrity, and cells have evolved mechanisms to restrict retrotransposition. However, how cellular proteins facilitate L1 retrotransposition requires elucidation. Here, we demonstrate that single-strand DNA breaks induced by the L1 endonuclease trigger the recruitment of poly(ADP-ribose) polymerase 2 (PARP2) to L1 integration sites and that PARP2 activation leads to the subsequent recruitment of the replication protein A (RPA) complex to facilitate retrotransposition. We further demonstrate that RPA directly binds activated PARP2 through poly(ADP-ribosyl)ation and can protect single-strand L1 integration intermediates from APOBEC3-mediated cytidine deamination in vitro. Paradoxically, we provide evidence that RPA can guide APOBEC3A, and perhaps other APOBEC3 proteins, to sites of L1 integration. Thus, the interplay of L1-encoded and evolutionarily conserved cellular proteins is required for efficient retrotransposition; however, these interactions also may be exploited to restrict L1 retrotransposition in the human genome.

Variable patterns of retrotransposition in different HeLa strains provide mechanistic insights into SINE RNA mobilization processes.

Alu elements are non-autonomous Short INterspersed Elements (SINEs) derived from the 7SL RNA gene that are present at over one million copies in human genomic DNA. Alu mobilizes by a mechanism known as retrotransposition, which requires the Long INterspersed Element-1 (LINE-1) ORF2-encoded protein (ORF2p). Here, we demonstrate that HeLa strains differ in their capacity to support Alu retrotransposition. Human Alu elements retrotranspose efficiently in HeLa-HA and HeLa-CCL2 (Alu-permissive) strains, but not in HeLa-JVM or HeLa-H1 (Alu-nonpermissive) strains. A similar pattern of retrotransposition was observed for other 7SL RNA-derived SINEs and tRNA-derived SINEs. In contrast, mammalian LINE-1s, a zebrafish LINE, a human SINE-VNTR-Alu (SVA) element, and an L1 ORF1-containing mRNA can retrotranspose in all four HeLa strains. Using an in vitro reverse transcriptase-based assay, we show that Alu RNAs associate with ORF2p and are converted into cDNAs in both Alu-permissive and Alu-nonpermissive HeLa strains, suggesting that 7SL- and tRNA-derived SINEs use strategies to 'hijack' L1 ORF2p that are distinct from those used by SVA elements and ORF1-containing mRNAs. These data further suggest ORF2p associates with the Alu RNA poly(A) tract in both Alu-permissive and Alu-nonpermissive HeLa strains, but that Alu retrotransposition is blocked after this critical step in Alu-nonpermissive HeLa strains.

LINE-1 retrotransposition activity in human genomes.

Highly active (i.e., "hot") long interspersed element-1 (LINE-1 or L1) sequences comprise the bulk of retrotransposition activity in the human genome; however, the abundance of hot L1s in the human population remains largely unexplored. Here, we used a fosmid-based, paired-end DNA sequencing strategy to identify 68 full-length L1s that are differentially present among individuals but are absent from the human genome reference sequence. The majority of these L1s were highly active in a cultured cell retrotransposition assay. Genotyping 26 elements revealed that two L1s are only found in Africa and that two more are absent from the H952 subset of the Human Genome Diversity Panel. Therefore, these results suggest that hot L1s are more abundant in the human population than previously appreciated, and that ongoing L1 retrotransposition continues to be a major source of interindividual genetic variation.

Problematic meta-analyses: Bayesian and frequentist perspectives on combining randomized controlled trials and non-randomized studies.

PURPOSE: In the literature, the propriety of the meta-analytic treatment-effect produced by combining randomized controlled trials (RCT) and non-randomized studies (NRS) is questioned, given the inherent confounding in NRS that may bias the meta-analysis. The current study compared an implicitly principled pooled Bayesian meta-analytic treatment-effect with that of frequentist pooling of RCT and NRS to determine how well each approach handled the NRS bias. MATERIALS & METHODS: Binary outcome Critical-Care meta-analyses, reflecting the importance of such outcomes in Critical-Care practice, combining RCT and NRS were identified electronically. Bayesian pooled treatment-effect and 95% credible-intervals (BCrI), posterior model probabilities indicating model plausibility and Bayes-factors (BF) were estimated using an informative heavy-tailed heterogeneity prior (half-Cauchy). Preference for pooling of RCT and NRS was indicated for Bayes-factors > 3 or < 0.333 for the converse. All pooled frequentist treatment-effects and 95% confidence intervals (FCI) were re-estimated using the popular DerSimonian-Laird (DSL) random effects model. RESULTS: Fifty meta-analyses were identified (2009-2021), reporting pooled estimates in 44; 29 were pharmaceutical-therapeutic and 21 were non-pharmaceutical therapeutic. Re-computed pooled DSL FCI excluded the null (OR or RR = 1) in 86% (43/50). In 18 meta-analyses there was an agreement between FCI and BCrI in excluding the null. In 23 meta-analyses where FCI excluded the null, BCrI embraced the null. BF supported a pooled model in 27 meta-analyses and separate models in 4. The highest density of the posterior model probabilities for 0.333 < Bayes factor < 1 was 0.8. CONCLUSIONS: In the current meta-analytic cohort, an integrated and multifaceted Bayesian approach gave support to including NRS in a pooled-estimate model. Conversely, caution should attend the reporting of naïve frequentist pooled, RCT and NRS, meta-analytic treatment effects.

Post-transcriptional microRNA repression of PMP22 dose in severe Charcot-Marie-Tooth disease type 1.

Copy number variation (CNV) may lead to pathological traits, and Charcot-Marie-Tooth disease type 1A (CMT1A), the commonest inherited peripheral neuropathy, is due to a genomic duplication encompassing the dosage-sensitive PMP22 gene. MicroRNAs act as repressors on post-transcriptional regulation of gene expression and in rodent models of CMT1A, overexpression of one such microRNA (miR-29a) has been shown to reduce the PMP22 transcript and protein level. Here we present genomic and functional evidence, for the first time in a human CNV-associated phenotype, of the 3' untranslated region (3'-UTR)-mediated role of microRNA repression on gene expression. The proband of the family presented with an early-onset, severe sensorimotor demyelinating neuropathy and harboured a novel de novo deletion in the PMP22 3'-UTR. The deletion is predicted to include the miR-29a seed binding site and transcript analysis of dermal myelinated nerve fibres using a novel platform, revealed a marked increase in PMP22 transcript levels. Functional evidence from Schwann cell lines harbouring the wild-type and mutant 3'-UTR showed significantly increased reporter assay activity in the latter, which was not ameliorated by overexpression of a miR-29a mimic. This shows the importance of miR-29a in regulating PMP22 expression and opens an avenue for therapeutic drug development.

Condensin I and condensin II proteins form a LINE-1 dependent super condensin complex and cooperate to repress LINE-1.

Condensin I and condensin II are multi-subunit complexes that are known for their individual roles in genome organization and preventing genomic instability. However, interactions between condensin I and condensin II subunits and cooperative roles for condensin I and condensin II, outside of their genome organizing functions, have not been reported. We previously discovered that condensin II cooperates with Gamma Interferon Activated Inhibitor of Translation (GAIT) proteins to associate with Long INterspersed Element-1 (LINE-1 or L1) RNA and repress L1 protein expression and the retrotransposition of engineered L1 retrotransposition in cultured human cells. Here, we report that the L1 3'UTR is required for condensin II and GAIT association with L1 RNA, and deletion of the L1 RNA 3'UTR results in increased L1 protein expression and retrotransposition. Interestingly, like condensin II, we report that condensin I also binds GAIT proteins, associates with the L1 RNA 3'UTR, and represses L1 retrotransposition. We provide evidence that the condensin I protein, NCAPD2, is required for condensin II and GAIT protein association with L1 RNA. Furthermore, condensin I and condensin II subunits interact to form a L1-dependent super condensin complex (SCC) which is located primarily within the cytoplasm of both transformed and primary epithelial cells. These data suggest that increases in L1 expression in epithelial cells promote cytoplasmic condensin protein associations that facilitate a feedback loop in which condensins may cooperate to mediate L1 repression.

Long-read assembly of a Great Dane genome highlights the contribution of GC-rich sequence and mobile elements to canine genomes.

Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3' end of LINE-1_Cfs (i.e., LINE-1_Cf 3'-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

p53 genes function to restrain mobile elements.

Throughout the animal kingdom, p53 genes govern stress response networks by specifying adaptive transcriptional responses. The human member of this gene family is mutated in most cancers, but precisely how p53 functions to mediate tumor suppression is not well understood. Using Drosophila and zebrafish models, we show that p53 restricts retrotransposon activity and genetically interacts with components of the piRNA (piwi-interacting RNA) pathway. Furthermore, transposon eruptions occurring in the p53(-) germline were incited by meiotic recombination, and transcripts produced from these mobile elements accumulated in the germ plasm. In gene complementation studies, normal human p53 alleles suppressed transposons, but mutant p53 alleles from cancer patients could not. Consistent with these observations, we also found patterns of unrestrained retrotransposons in p53-driven mouse and human cancers. Furthermore, p53 status correlated with repressive chromatin marks in the 5' sequence of a synthetic LINE-1 element. Together, these observations indicate that ancestral functions of p53 operate through conserved mechanisms to contain retrotransposons. Since human p53 mutants are disabled for this activity, our findings raise the possibility that p53 mitigates oncogenic disease in part by restricting transposon mobility.

Modelling of intensive care unit (ICU) length of stay as a quality measure: a problematic exercise.

BACKGROUND: Intensive care unit (ICU) length of stay (LOS) and the risk adjusted equivalent (RALOS) have been used as quality metrics. The latter measures entail either ratio or difference formulations or ICU random effects (RE), which have not been previously compared. METHODS: From calendar year 2016 data of an adult ICU registry-database (Australia & New Zealand Intensive Care Society (ANZICS) CORE), LOS predictive models were established using linear (LMM) and generalised linear (GLMM) mixed models. Model fixed effects quality-metric formulations were estimated as RALOSR for LMM (geometric mean derived from log(ICU LOS)) and GLMM (day) and observed minus expected ICU LOS (OMELOS from GLMM). Metric confidence intervals (95%CI) were estimated by bootstrapping; random effects (RE) were predicted for LMM and GLMM. Forest-plot displays of ranked quality-metric point-estimates (95%CI) were generated for ICU hospital classifications (metropolitan, private, rural/regional, and tertiary). Robust rank confidence sets (point estimate and 95%CI), both marginal (pertaining to a singular ICU) and simultaneous (pertaining to all ICU differences), were established. RESULTS: The ICU cohort was of 94,361 patients from 125 ICUs (metropolitan 16.9%, private 32.8%, rural/regional 6.4%, tertiary 43.8%). Age (mean, SD) was 61.7 (17.5) years; 58.3% were male; APACHE III severity-of-illness score 54.6 (25.7); ICU annual patient volume 1192 (702) and ICU LOS 3.2 (4.9). There was no concordance of ICU ranked model predictions, GLMM versus LMM, nor for the quality metrics used, RALOSR, OMELOS and site-specific RE for each of the ICU hospital classifications. Furthermore, there was no concordance between ICU ranking confidence sets, marginal and simultaneous for models or quality metrics. CONCLUSIONS: Inference regarding adjusted ICU LOS was dependent upon the statistical estimator and the quality index used to quantify any LOS differences across ICUs. That is, there was no "one best model"; thus, ICU "performance" is determined by model choice and any rankings thereupon should be circumspect.

A 3' Poly(A) Tract Is Required for LINE-1 Retrotransposition.

L1 retrotransposons express proteins (ORF1p and ORF2p) that preferentially mobilize their encoding RNA in cis, but they also can mobilize Alu RNA and, more rarely, cellular mRNAs in trans. Although these RNAs differ in sequence, each ends in a 3' polyadenosine (poly(A)) tract. Here, we replace the L1 polyadenylation signal with sequences derived from a non-polyadenylated long non-coding RNA (MALAT1), which can form a stabilizing triple helix at the 3' end of an RNA. L1/MALAT RNAs accumulate in cells, lack poly(A) tails, and are translated; however, they cannot retrotranspose in cis. Remarkably, the addition of a 16 or 40 base poly(A) tract downstream of the L1/MALAT triple helix restores retrotransposition in cis. The presence of a poly(A) tract also allows ORF2p to bind and mobilize RNAs in trans. Thus, a 3' poly(A) tract is critical for the retrotransposition of sequences that comprise approximately one billion base pairs of human DNA.

H3K27 demethylases are dispensable for activation of Polycomb-regulated injury response genes in peripheral nerve.

The induction of nerve injury response genes in Schwann cells depends on both transcriptional and epigenomic reprogramming. The nerve injury response program is regulated by the repressive histone mark H3K27 trimethylation (H3K27me3), deposited by Polycomb repressive complex 2 (PRC2). Loss of PRC2 function leads to early and augmented induction of the injury response gene network in peripheral nerves, suggesting H3K27 demethylases are required for derepression of Polycomb-regulated nerve injury genes. To determine the function of H3K27 demethylases in nerve injury, we generated Schwann cell-specific knockouts of H3K27 demethylase Kdm6b and double knockouts of Kdm6b/Kdm6a (encoding JMJD3 and UTX). We found that H3K27 demethylases are largely dispensable for Schwann cell development and myelination. In testing the function of H3K27 demethylases after injury, we found early induction of some nerve injury genes was diminished compared with control, but most injury genes were largely unaffected at 1 and 7 days post injury. Although it was proposed that H3K27 demethylases are required to activate expression of the cyclin-dependent kinase inhibitor Cdkn2a in response to injury, Schwann cell-specific deletion of H3K27 demethylases affected neither the expression of this gene nor Schwann cell proliferation after nerve injury. To further characterize the regulation of nerve injury response genes, we found that injury genes are associated with repressive histone H2AK119 ubiquitination catalyzed by PRC1, which declines after injury. Overall, our results indicate H3K27 demethylation is not required for induction of injury response genes and that other mechanisms likely are involved in activating Polycomb-repressed injury genes in peripheral nerve.

Pmp22 super-enhancer deletion causes tomacula formation and conduction block in peripheral nerves.

Copy number variation of the peripheral nerve myelin gene Peripheral Myelin Protein 22 (PMP22) causes multiple forms of inherited peripheral neuropathy. The duplication of a 1.4 Mb segment surrounding this gene in chromosome 17p12 (c17p12) causes the most common form of Charcot-Marie-Tooth disease type 1A, whereas the reciprocal deletion of this gene causes a separate neuropathy termed hereditary neuropathy with liability to pressure palsies (HNPP). PMP22 is robustly induced in Schwann cells in early postnatal development, and several transcription factors and their cognate regulatory elements have been implicated in coordinating the gene's proper expression. We previously found that a distal super-enhancer domain was important for Pmp22 expression in vitro, with particular impact on a Schwann cell-specific alternative promoter. Here, we investigate the consequences of deleting this super-enhancer in vivo. We find that loss of the super-enhancer in mice reduces Pmp22 expression throughout development and into adulthood, with greater impact on the Schwann cell-specific promoter. Additionally, these mice display tomacula formed by excessive myelin folding, a pathological hallmark of HNPP, as have been previously observed in heterozygous Pmp22 mice as well as sural biopsies from patients with HNPP. Our findings demonstrate a mechanism by which smaller copy number variations, not including the Pmp22 gene, are sufficient to reduce gene expression and phenocopy a peripheral neuropathy caused by the HNPP-associated deletion encompassing PMP22.

Hospital-acquired complications: the relative importance of hospital- and patient-related factors.

OBJECTIVE: To quantify the prevalence of hospital-acquired complications; to determine the relative influence of patient- and hospital-related factors on complication rates. DESIGN, PARTICIPANTS: Retrospective analysis of administrative data (Integrated South Australian Activity Collection; Victorian Admitted Episodes Dataset) for multiple-day acute care episodes for adults in public hospitals. SETTING: Thirty-eight major public hospitals in South Australia and Victoria, 2015-2018. MAIN OUTCOME MEASURES: Hospital-acquired complication rates, overall and by complication class, by hospital and hospital type (tertiary referral, major metropolitan service, major regional service); variance in rates (intra-class correlation coefficient, ICC) at the patient, hospital, and hospital type levels as surrogate measures of their influence on rates. RESULTS: Of 1 558 978 public hospital episodes (10 029 918 bed-days), 151 486 included a total of 214 286 hospital-acquired complications (9.72 [95% CI, 9.67-9.77] events per 100 episodes; 2.14 [95% CI, 2.13-2.15] events per 100 bed-days). Complication rates were highest in tertiary referral hospitals (12.7 [95% CI, 12.6-12.8] events per 100 episodes) and for episodes including intensive care components (37.1 [95% CI, 36.7-37.4] events per 100 episodes). For all complication classes, inter-hospital variation was determined more by patient factors (overall ICC, 0.55; 95% CI, 0.53-0.57) than by hospital factors (ICC, 0.04; 95% CI, 0.02-0.07) or hospital type (ICC, 0.01; 95% CI, 0.001-0.03). CONCLUSIONS: Hospital-acquired complications were recorded for 9.7% of hospital episodes, but patient-related factors played a greater role in determining their prevalence than the treating hospital.

Identification and characterization of occult human-specific LINE-1 insertions using long-read sequencing technology.

Long Interspersed Element-1 (LINE-1) retrotransposition contributes to inter- and intra-individual genetic variation and occasionally can lead to human genetic disorders. Various strategies have been developed to identify human-specific LINE-1 (L1Hs) insertions from short-read whole genome sequencing (WGS) data; however, they have limitations in detecting insertions in complex repetitive genomic regions. Here, we developed a computational tool (PALMER) and used it to identify 203 non-reference L1Hs insertions in the NA12878 benchmark genome. Using PacBio long-read sequencing data, we identified L1Hs insertions that were absent in previous short-read studies (90/203). Approximately 81% (73/90) of the L1Hs insertions reside within endogenous LINE-1 sequences in the reference assembly and the analysis of unique breakpoint junction sequences revealed 63% (57/90) of these L1Hs insertions could be genotyped in 1000 Genomes Project sequences. Moreover, we observed that amplification biases encountered in single-cell WGS experiments led to a wide variation in L1Hs insertion detection rates between four individual NA12878 cells; under-amplification limited detection to 32% (65/203) of insertions, whereas over-amplification increased false positive calls. In sum, these data indicate that L1Hs insertions are often missed using standard short-read sequencing approaches and long-read sequencing approaches can significantly improve the detection of L1Hs insertions present in individual genomes.

RNA ligation precedes the retrotransposition of U6/LINE-1 chimeric RNA.

Long interspersed element-1 (LINE-1 or L1) amplifies via retrotransposition. Active L1s encode 2 proteins (ORF1p and ORF2p) that bind their encoding transcript to promote retrotransposition in cis The L1-encoded proteins also promote the retrotransposition of small-interspersed element RNAs, noncoding RNAs, and messenger RNAs in trans Some L1-mediated retrotransposition events consist of a copy of U6 RNA conjoined to a variably 5'-truncated L1, but how U6/L1 chimeras are formed requires elucidation. Here, we report the following: The RNA ligase RtcB can join U6 RNAs ending in a 2',3'-cyclic phosphate to L1 RNAs containing a 5'-OH in vitro; depletion of endogenous RtcB in HeLa cell extracts reduces U6/L1 RNA ligation efficiency; retrotransposition of U6/L1 RNAs leads to U6/L1 pseudogene formation; and a unique cohort of U6/L1 chimeric RNAs are present in multiple human cell lines. Thus, these data suggest that U6 small nuclear RNA (snRNA) and RtcB participate in the formation of chimeric RNAs and that retrotransposition of chimeric RNA contributes to interindividual genetic variation.

Heatwave-induced synchrony within forage fish portfolio disrupts energy flow to top pelagic predators.

During the Pacific marine heatwave of 2014-2016, abundance and quality of several key forage fish species in the Gulf of Alaska were simultaneously reduced throughout the system. Capelin (Mallotus catervarius), sand lance (Ammodytes personatus), and herring (Clupea pallasii) populations were at historically low levels, and within this community abrupt declines in portfolio effects identify trophic instability at the onset of the heatwave. Although compensatory changes in age structure, size, growth or energy content of forage fish were observed to varying degrees among all these forage fish, none were able to fully mitigate adverse impacts of the heatwave, which likely included both top-down and bottom-up forcing. Notably, changes to the demographic structure of forage fish suggested size-selective removals typical of top-down regulation. At the same time, changes in zooplankton communities may have driven bottom-up regulation as copepod community structure shifted toward smaller, warm water species, and euphausiid biomass was reduced owing to the loss of cold-water species. Mediated by these impacts on the forage fish community, an unprecedented disruption of the normal pelagic food web was signaled by higher trophic level disruptions during 2015-2016, when seabirds, marine mammals, and groundfish experienced shifts in distribution, mass mortalities, and reproductive failures. Unlike decadal-scale variability underlying ecosystem regime shifts, the heatwave appeared to temporarily overwhelm the ability of the forage fish community to buffer against changes imposed by warm water anomalies, thereby eliminating any ecological advantages that may have accrued from having a suite of coexisting forage species with differing life-history compensations.

Spliced integrated retrotransposed element (SpIRE) formation in the human genome.

Human Long interspersed element-1 (L1) retrotransposons contain an internal RNA polymerase II promoter within their 5' untranslated region (UTR) and encode two proteins, (ORF1p and ORF2p) required for their mobilization (i.e., retrotransposition). The evolutionary success of L1 relies on the continuous retrotransposition of full-length L1 mRNAs. Previous studies identified functional splice donor (SD), splice acceptor (SA), and polyadenylation sequences in L1 mRNA and provided evidence that a small number of spliced L1 mRNAs retrotransposed in the human genome. Here, we demonstrate that the retrotransposition of intra-5'UTR or 5'UTR/ORF1 spliced L1 mRNAs leads to the generation of spliced integrated retrotransposed elements (SpIREs). We identified a new intra-5'UTR SpIRE that is ten times more abundant than previously identified SpIREs. Functional analyses demonstrated that both intra-5'UTR and 5'UTR/ORF1 SpIREs lack Cis-acting transcription factor binding sites and exhibit reduced promoter activity. The 5'UTR/ORF1 SpIREs also produce nonfunctional ORF1p variants. Finally, we demonstrate that sequence changes within the L1 5'UTR over evolutionary time, which permitted L1 to evade the repressive effects of a host protein, can lead to the generation of new L1 splicing events, which, upon retrotransposition, generates a new SpIRE subfamily. We conclude that splicing inhibits L1 retrotransposition, SpIREs generally represent evolutionary "dead-ends" in the L1 retrotransposition process, mutations within the L1 5'UTR alter L1 splicing dynamics, and that retrotransposition of the resultant spliced transcripts can generate interindividual genomic variation.

Multivariate meta-analysis of critical care meta-analyses: a meta-epidemiological study.

BACKGROUND: Meta-analyses typically consider multiple outcomes and report univariate effect sizes considered as independent. Multivariate meta-analysis (MVMA) incorporates outcome correlation and synthesises direct evidence and related outcome estimates within a single analysis. In a series of meta-analyses from the critically ill literature, the current study contrasts multiple univariate effect estimates and their precision with those derived from MVMA. METHODS: A previous meta-epidemiological study was used to identify meta-analyses with either one or two secondary outcomes providing sufficient detail to structure bivariate or tri-variate MVMA, with mortality as primary outcome. Analysis was performed using a random effects model for both odds ratio (OR) and risk ratio (RR); borrowing of strength (BoS) between multivariate outcome estimates was reported. Estimate comparisons, ß coefficients, standard errors (SE) and confidence interval (CI) width, univariate versus multivariate, were performed using Lin's concordance correlation coefficient (CCC). RESULTS: In bivariate meta-analyses, for OR (n = 49) and RR (n = 48), there was substantial concordance (≥ 0.69) between estimates; but this was less so for tri-variate meta-analyses for both OR (n = 25; ≥ 0.38) and RR (≥ -0.10; n = 22). A variable change in the multivariate precision of primary mortality outcome estimates compared with univariate was present for both bivariate and tri-variate meta-analyses and for metrics. For second outcomes, precision tended to decrease and CI width increase for bivariate meta-analyses, but was variable in the tri-variate. For third outcomes, precision increased and CI width decreased. In bivariate meta-analyses, OR coefficient significance reversal, univariate versus MVMA, occurred once for mortality and 6 cases for second outcomes. RR coefficient significance reversal occurred in 4 cases; 2 were discordant with OR. For tri-variate OR meta-analyses reversal of coefficient estimate significance occurred in two cases for mortality, nine cases for second and 7 cases for third outcomes. In RR meta-analyses significance reversals occurred for mortality in 2 cases, 6 cases for second and 3 cases for third; there were 7 discordances with OR. BoS was greater in trivariate MVMAs compared with bivariate and for OR versus RR. CONCLUSIONS: MVMA would appear to be the preferred solution to multiple univariate analyses; parameter significance changes may occur. Analytic metric appears to be a determinant.