DNA methylation is a major epigenetic event that affects not only cellular gene expression but that also has the potential to influence bacterial and viral DNA in their host-dependent functions. Adeno-associated virus (AAV) genome contains a high degree of CpG sequences capable of methylation in its terminal repeat sequences, which are the sole elements retained in AAV-based vectors used in gene therapy. The present study determined the influence of methylation status of the host cell on wild type (wt) AAV integration and recombinant (r) AAV transgene expression in HeLa cells. Results of the study indicated that hypo-methylation significantly enhanced both wtAAV chromosomal integration and transgene expression of rAAV. A direct influence of methylation on AAV integration was further confirmed by methylating the AAVS1 integration sites prior to viral infection with DNA trans-complementation assay. These results signify the importance of epigenetic status of target cells as one of the key factors in long-term transgene expression in AAV gene therapy.

The avian α-herpesvirus known as Marek's disease virus (MDV) linearly integrates its genomic DNA into host telomeres during infection. The resulting disease, Marek's disease (MD), is characterized by virally-induced lymphomas with high mortality. The temporal dynamics of MDV-positive (MDV+) transformed cells and expansion of MD lymphomas remain targets for further understanding. It also remains to be determined whether specific host chromosomal sites of MDV telomere integration confer an advantage to MDV-transformed cells during tumorigenesis. We applied MDV-specific fluorescence in situ hybridization (MDV FISH) to investigate virus-host cytogenomic interactions within and among a total of 37 gonad lymphomas and neoplastic splenic samples in birds infected with virulent MDV. We also determined single-cell, chromosome-specific MDV integration profiles within and among transformed tissue samples, including multiple samples from the same bird. Most mitotically-dividing cells within neoplastic samples had the cytogenomic phenotype of 'MDV telomere-integrated only', and tissue-specific, temporal changes in phenotype frequencies were detected. Transformed cell populations composing gonad lymphomas exhibited significantly lower diversity, in terms of heterogeneity of MDV integration profiles, at the latest stages of tumorigenesis (>50 days post-infection (dpi)). We further report high interindividual and lower intraindividual variation in MDV integration profiles of lymphoma cells. There was no evidence of integration hotspots into a specific host chromosome(s). Collectively, our data suggests that very few transformed MDV+ T cell populations present earlier in MDV-induced lymphomas (32-50 dpi), survive, and expand to become the dominant clonal population in more advanced MD lymphomas (51-62 dpi) and establish metastatic lymphomas.

Superior oil yield is always the top priority of the oil palm industry. Short trunk height (THT) and compactness traits have become increasingly important to improve harvesting efficiency since the industry started to suffer yield losses due to labor shortages. Breeding populations with low THT and short frond length (FL) are actually available, such as Dumpy AVROS pisifera (DAV) and Gunung Melayu dura (GM). However, multiple trait stacking still remains a challenge for oil palm breeding, which usually requires 12-20 years to complete a breeding cycle. In this study, yield and height increment in the GM × GM (GM-3341) and the GM × DAV (GM-DAV-3461) crossing programs were evaluated and palms with good yield and smaller height increment were identified. In the GM-3341 family, non-linear THT growth between THT_2008 (seven years old) and THT_2014 (13 years old) was revealed by a moderate correlation, suggesting that inter-palm competition becomes increasingly important. In total, 19 quantitative trait loci (QTLs) for THT_2008 (8), oil per palm (O/P) (7) and FL (4) were localized on the GM-3341 linkage map, with an average mapping interval of 2.01 cM. Three major QTLs for THT_2008, O/P and FL are co-located on chromosome 11 and reflect the correlation of THT_2008 with O/P and FL. Multiple trait selection for high O/P and low THT (based on the cumulative effects of positive alleles per trait) identified one palm from 100 palms, but with a large starting population of 1000-1500 seedling per cross, this low frequency could be easily compensated for during breeding selection.

Eucommia ulmoides has attracted much attention as a valuable natural rubber (Eu-rubber) production tree. As a strategic material, Eu-rubber plays a vital role in general and defence industries. However, the study of Eu-rubber biosynthesis at a molecular level is scarce, and the regulatory network between microRNAs (miRNAs) and messenger RNAs (mRNAs) in Eu-rubber biosynthesis has not been assessed. In this study, we comprehensively analyzed the transcriptomes, small RNAs (sRNAs) and degradome to reveal the regulatory network of Eu-rubber biosynthesis in E. ulmoides. A total of 82,065 unigenes and 221 miRNAs were identified using high-throughput sequencing; 20,815 targets were predicted using psRNATarget software. Of these targets, 779 miRNA-target pairs were identified via degradome sequencing. Thirty-one miRNAs were differentially expressed; 22 targets of 34 miRNAs were annotated in the terpenoid backbone biosynthesis pathway (ko00900) based on the Kyoto Encyclopedia of Genes and Genomes (KEGG). These miRNAs were putatively related to Eu-rubber biosynthesis. A regulatory network was constructed according to the expression profiles of miRNAs and their targets. These results provide a comprehensive analysis of transcriptomics, sRNAs and degradome to reveal the Eu-rubber accumulation, and provide new insights into genetic engineering techniques which may improve the content of Eu-rubber in E. ulmoides.

Newcastle Disease (ND) is a continuing global threat to domestic poultry, especially in developing countries, where severe outbreaks of velogenic ND virus (NDV) often cause major economic losses to households. Local chickens are of great importance to rural family livelihoods through provision of high-quality protein. To investigate the genetic basis of host response to NDV, three popular Tanzanian chicken ecotypes (regional populations) were challenged with a lentogenic (vaccine) strain of NDV at 28 days of age. Various host response phenotypes, including anti-NDV antibody levels (pre-infection and 10 days post-infection, dpi), and viral load (2 and 6 dpi) were measured, in addition to growth rate. We estimated genetic parameters and conducted genome-wide association study analyses by genotyping 1399 chickens using the Affymetrix 600K chicken SNP chip. Estimates of heritability of the evaluated traits were moderate (0.18-0.35). Five quantitative trait loci (QTL) associated with growth and/or response to NDV were identified by single-SNP analyses, with some regions explaining ≥1% of genetic variance based on the Bayes-B method. Immune related genes, such as ETS1, TIRAP, and KIRREL3, were located in regions associated with viral load at 6 dpi. The moderate estimates of heritability and identified QTL indicate that NDV response traits may be improved through selective breeding of chickens to enhance increased NDV resistance and vaccine efficacy in Tanzanian local ecotypes.

In the last 5 years, frequent outbreaks of infectious bronchitis virus (IBV) are observed in both broiler and layer chicken flocks in the Kingdom of Saudi Arabia (KSA) in spite of extensive usage of vaccines. The IBV is a widespread avian coronavirus affecting both vaccinated and unvaccinated chicken flocks and is attributed to significant economic losses, around the globe. In the present study, 58 (n = 58) samples were collected from four different commercial poultry flocks from 8 KSA districts during 2019. A total of nine positive isolates (9/58; 15.5%), based on real-time reverse transcriptase PCR targeting nucleocapsid (N) gene, were used for further genetic characterization and evolutionary analysis. Genetic characterization of the partial spike (S1) gene revealed the clustering of the reported isolates into three different genotypes, whereas four additional isolates were grouped within 4/91 genotype, two isolates within IS/885 genotype, one isolate was closely related to IS/1494/06, and two isolates were grouped within classic serotype (vaccine-like strains). Phylodynamic revealed clustering of four isolated viruses within GI-13 lineage, three isolates within GI-23 lineage, and two isolates within GI-1 lineage. Results indicate that there are high evolutionary distances between the newly identified IBV strains in this study and the commercially used vaccines (GI-1), suggesting that IBV strains circulating in the KSA are under constant evolutionary pressures. Selective pressure biostatistics analyses consistently demonstrate the presence of a higher positive score which highlights the role of natural selection, a mechanism of virus evolution on sites located on the protein surface, within or nearby domains involved in viral attachment or related functions. Recombination analysis revealed emergence of two isolates through recombination events resulting in new recombinant viruses. Taken together, these finding demonstrate the genetic and evolutionary insights into the currently circulating IBV genotypes in KSA, which could help to better understand the origin, spread, and evolution of infectious bronchitis viruses, and to ascertain the importance of disease monitoring as well as re-evaluation for the currently used vaccines and vaccination programs.

Insect salivary glands play an important role for host feeding, specifically by secreting salivary proteins for digestion and potentially modulating host defenses. Compared to other hemipterans, the significance of salivary glands is less studied in the black-faced leafhopper, Graminella nigrifrons, a crop pest that vectors several agronomically important plant viruses. To identify functionally important genes in the salivary glands of the black-faced leafhopper, we compared transcriptomes between adult salivary glands (SG) and the remaining carcasses. We identified 14,297 salivary gland-enriched transcripts and 195 predicted secretory peptides (i.e., with a signal peptide and extracellular localization characteristics). Overall, the SG transcriptome included functions such as 'oxidoreduction', 'membrane transport', and 'ATP-binding', which might be important for the fundamental physiology of this tissue. We further evaluated transcripts with potential contributions in host feeding using RT-qPCR. Two SG-enriched transcripts (log2 fold change > 5), GnP19 and GnE63 (a putative calcium binding protein), were significantly upregulated in maize-fed adults relative to starved adults, validating their importance in feeding. The SG-enriched transcripts of the black-faced leafhopper could play a potential role for interacting with maize and could be targets of interest for further functional studies and improve pest control and disease transmission.

Newcastle disease virus (NDV) is a highly contagious avian pathogen that poses a tremendous threat to poultry producers in endemic zones due to its epidemic potential. To investigate host genetic resistance to NDV while under the effects of heat stress, a genome-wide association study (GWAS) was performed on Hy-Line Brown layer chickens that were challenged with NDV while under high ambient temperature to identify regions associated with host viral titer, circulating anti-NDV antibody titer, and body weight change. A single nucleotide polymorphism (SNP) on chromosome 1 was associated with viral titer at two days post-infection (dpi), while 30 SNPs spanning a quantitative trait loci (QTL) on chromosome 24 were associated with viral titer at 6 dpi. Immune related genes, such as CAMK1d and CCDC3 on chromosome 1, associated with viral titer at 2 dpi, and TIRAP, ETS1, and KIRREL3, associated with viral titer at 6 dpi, were located in two QTL regions for viral titer that were identified in this study. This study identified genomic regions and candidate genes that are associated with response to NDV during heat stress in Hy-Line Brown layer chickens. Regions identified for viral titer on chromosome 1 and 24, at 2 and 6 dpi, respectively, included several genes that have key roles in regulating the immune response.

Highly pathogenic avian influenza viruses (HPAIVs) in gallinaceous poultry are associated with viral infection of the endothelium, the induction of a 'cytokine storm, and severe disease. In contrast, in Pekin ducks, HPAIVs are rarely endothelial tropic, and a cytokine storm is not observed. To date, understanding these species-dependent differences in pathogenesis has been hampered by the absence of a pure culture of duck and chicken endothelial cells. Here, we use our recently established in vitro cultures of duck and chicken aortic endothelial cells to investigate species-dependent differences in the response of endothelial cells to HPAIV H5N1 infection. We demonstrate that chicken and duck endothelial cells display a different transcriptional response to HPAI H5N1 infection in vitro-with chickens displaying a more pro-inflammatory response to infection. As similar observations were recorded following in vitro stimulation with the viral mimetic polyI:C, these findings were not specific to an HPAIV H5N1 infection. However, similar species-dependent differences in the transcriptional response to polyI:C were not observed in avian fibroblasts. Taken together, these data demonstrate that chicken and duck endothelial cells display a different response to HPAIV H5N1 infection, and this may help account for the species-dependent differences observed in inflammation in vivo.

Cardiomyopathy syndrome is a viral disease of Atlantic salmon, mostly affecting fish during the late stages of production, resulting in significant losses to the industry. It has been shown that resistance to this disease has a strong genetic component, with quantitative trait loci (QTL) on chromosomes 27 (Ssa27) and Ssa12 to explain most of the additive genetic variance. Here, by analysing animals from a different year-class and a different population, we further aimed to confirm and narrow down the locations of these QTL. The data support the existence of the two QTL and suggest that the causative mutation on Ssa27 is most likely within the 10-10.5 Mbp segment of this chromosome. This region contains a cluster of major histocompatibility complex class I (MHC I) genes with the most strongly associated marker mapped to one of these loci. On Ssa12, the data confirmed the previous finding that the location of the causative mutation is within the 61.3 to 61.7 Mbp region. This segment contains several immune-related genes, but of particular interest are genes related to MHC II. Together, these findings highlight the likely key role of MHC genes in Atlantic salmon following infection with Piscine myocarditis virus (PMCV) and their potential impact on influencing the trajectory of this disease.

Hepatitis B virus (HBV) infection is a leading cause of hepatocellular carcinoma (HCC) worldwide. The integration of HBV genomic DNA into the host genome occurs randomly, early after infection, and is associated with hepatocarcinogenesis in HBV-infected patients. Therefore, it is important to analyze HBV genome integration. We analyzed HBV genome integration in human hepatoma PLC/PRF/5 cells by HBV sequence capture-based next-generation sequencing (NGS) methods. We confirmed the results by using Sanger sequencing methods. We observed that HBV genotype A is integrated into the genome of PLC/PRF/5 cells. HBV sequence capture-based NGS is useful for the analysis of HBV genome integrants and their locations in the human genome. Among the HBV genome integrants, we performed functional analysis and demonstrated the automatic expression of some HBV proteins encoded by HBV integrants from chromosomes 3 and 11 in Huh7 cells transfected with these DNA sequences. HBV sequence capture-based NGS may be a useful tool for the assessment of HBV genome integration into the human genome in clinical samples and suggests new strategies for hepatocarcinogenesis in HBV infection.

Single-cell RNA sequencing of the bronchoalveolar lavage fluid (BALF) samples from COVID-19 patients has enabled us to examine gene expression changes of human tissue in response to the SARS-CoV-2 virus infection. However, the underlying mechanisms of COVID-19 pathogenesis at single-cell resolution, its transcriptional drivers, and dynamics require further investigation. In this study, we applied machine learning algorithms to infer the trajectories of cellular changes and identify their transcriptional programs. Our study generated cellular trajectories that show the COVID-19 pathogenesis of healthy-to-moderate and healthy-to-severe on macrophages and T cells, and we observed more diverse trajectories in macrophages compared to T cells. Furthermore, our deep-learning algorithm DrivAER identified several pathways (e.g., xenobiotic pathway and complement pathway) and transcription factors (e.g., MITF and GATA3) that could be potential drivers of the transcriptomic changes for COVID-19 pathogenesis and the markers of the COVID-19 severity. Moreover, macrophages-related functions corresponded more to the disease severity compared to T cells-related functions. Our findings more proficiently dissected the transcriptomic changes leading to the severity of a COVID-19 infection.

Plants have a lot of diversity in epigenetic modifications such as DNA methylation in their natural populations or cultivars. Although many studies observing the epigenetic diversity within and among species have been reported, the mechanisms how these variations are generated are still not clear. In addition to the de novo spontaneous epi-mutation, the intra- and inter-specific crossing can also cause a change of epigenetic modifications in their progenies. Here we report an example of diversification of DNA methylation by crossing and succeeding selfing. We traced the inheritance pattern of epigenetic modification during the crossing experiment between two natural strains Columbia (Col), and Landsberg electa (Ler) in model plant Arabidopsis thaliana to observe the inheritance of DNA methylation in two organellar DNA-like sequence regions in the nuclear genome. Because organellar DNA integration to the nuclear genome is common in flowering plants and these sequences are occasionally methylated, such DNA could be the novel source of plant genome evolution. The amplicon sequencing, using bisulfite-converted DNA and a next-generation auto-sequencer, was able to efficiently track the heredity of DNA methylation in F₁ and F₂ populations. One region showed hypomethylation in the F₁ population and succeeding elevation of DNA methylation with large variance in the F₂ population. The methylation level of Col and Ler alleles in F₂ heterozygotes showed a significant positive correlation, implying the trans-chromosomal effect on DNA methylation. The results may suggest the possible mechanism causing the natural epigenetic diversity within plant populations.

Escherichia coli are among the most common foodborne pathogens associated with infections reported from meat sources. This study investigated the virulome, pathogenicity, stress response factors, clonal lineages, and the phylogenomic relationship of E. coli isolated from different meat sources in Ghana using whole-genome sequencing. Isolates were screened from five meat sources (beef, chevon, guinea fowl, local chicken, and mutton) and five areas (Aboabo, Central market, Nyorni, Victory cinema, and Tishegu) based in the Tamale Metropolis, Ghana. Following microbial identification, the E. coli strains were subjected to whole-genome sequencing. Comparative visualisation analyses showed different DNA synteny of the strains. The isolates consisted of diverse sequence types (STs) with the most common being ST155 (n = 3/14). Based Upon Related Sequence Types (eBURST) analyses of the study sequence types identified four similar clones, five single-locus variants, and two satellite clones (more distantly) with global curated E. coli STs. All the isolates possessed at least one restriction-modification (R-M) and CRISPR defence system. Further analysis revealed conserved stress response mechanisms (detoxification, osmotic, oxidative, and periplasmic stress) in the strains. Estimation of pathogenicity predicted a higher average probability score (Pscore ≈ 0.937), supporting their pathogenic potential to humans. Diverse virulence genes that were clonal-specific were identified. Phylogenomic tree analyses coupled with metadata insights depicted the high genetic diversity of the E. coli isolates with no correlation with their meat sources and areas. The findings of this bioinformatic analyses further our understanding of E. coli in meat sources and are broadly relevant to the design of contamination control strategies in meat retail settings in Ghana.

Heavy metals (HMs) are compounds that can be hazardous and impair growth of living organisms. Bacteria have evolved the capability not only to cope with heavy metals but also to detoxify polluted environments. Three heavy metal-resistant strains of Mucilaginibacer rubeus and one of Mucilaginibacter kameinonensis were isolated from the gold/copper Zijin mining site, Longyan, Fujian, China. These strains were shown to exhibit high resistance to heavy metals with minimal inhibitory concentration reaching up to 3.5 mM Cu(II), 21 mM Zn(II), 1.2 mM Cd(II), and 10.0 mM As(III). Genomes of the four strains were sequenced by Illumina. Sequence analyses revealed the presence of a high abundance of heavy metal resistance (HMR) determinants. One of the strain, M. rubeus P2, carried genes encoding 6 putative PIB-1-ATPase, 5 putative PIB-3-ATPase, 4 putative Zn(II)/Cd(II) PIB-4 type ATPase, and 16 putative resistance-nodulation-division (RND)-type metal transporter systems. Moreover, the four genomes contained a high abundance of genes coding for putative metal binding chaperones. Analysis of the close vicinity of these HMR determinants uncovered the presence of clusters of genes potentially associated with mobile genetic elements. These loci included genes coding for tyrosine recombinases (integrases) and subunits of mating pore (type 4 secretion system), respectively allowing integration/excision and conjugative transfer of numerous genomic islands. Further in silico analyses revealed that their genetic organization and gene products resemble the Bacteroides integrative and conjugative element CTnDOT. These results highlight the pivotal role of genomic islands in the acquisition and dissemination of adaptive traits, allowing for rapid adaption of bacteria and colonization of hostile environments.

High-throughput structure profiling (SP) experiments that provide information at nucleotide resolution are revolutionizing our ability to study RNA structures. Of particular interest are RNA elements whose underlying structures are necessary for their biological functions. We previously introduced patteRNA, an algorithm for rapidly mining SP data for patterns characteristic of such motifs. This work provided a proof-of-concept for the detection of motifs and the capability of distinguishing structures displaying pronounced conformational changes. Here, we describe several improvements and automation routines to patteRNA. We then consider more elaborate biological situations starting with the comparison or integration of results from searches for distinct motifs and across datasets. To facilitate such analyses, we characterize patteRNA’s outputs and describe a normalization framework that regularizes results. We then demonstrate that our algorithm successfully discerns between highly similar structural variants of the human immunodeficiency virus type 1 (HIV-1) Rev response element (RRE) and readily identifies its exact location in whole-genome structure profiles of HIV-1. This work highlights the breadth of information that can be gleaned from SP data and broadens the utility of data-driven methods as tools for the detection of novel RNA elements.

Lentiviral-mediated integration of a CFTR transgene cassette into airway basal cells is a strategy being considered for cystic fibrosis (CF) cell-based therapies. However, CFTR expression is highly regulated in differentiated airway cell types and a subset of intermediate basal cells destined to differentiate. Since basal stem cells typically do not express CFTR, suppressing the CFTR expression from the lentiviral vector in airway basal cells may be beneficial for maintaining their proliferative capacity and multipotency. We identified miR-106b as highly expressed in proliferating airway basal cells and extinguished in differentiated columnar cells. Herein, we developed lentiviral vectors with the miR-106b-target sequence (miRT) to both study miR-106b regulation during basal cell differentiation and detarget CFTR expression in basal cells. Given that miR-106b is expressed in the 293T cells used for viral production, obstacles of viral genome integrity and titers were overcome by creating a 293T-B2 cell line that inducibly expresses the RNAi suppressor B2 protein from flock house virus. While miR-106b vectors effectively detargeted reporter gene expression in proliferating basal cells and following differentiation in the air-liquid interface and organoid cultures, the CFTR-miRT vector produced significantly less CFTR-mediated current than the non-miR-targeted CFTR vector following transduction and differentiation of CF basal cells. These findings suggest that miR-106b is expressed in certain airway cell types that contribute to the majority of CFTR anion transport in airway epithelium.

Notch is a key factor of a signaling cascade which regulates cell differentiation in all multicellular organisms. Numerous investigations have been directed mainly at studying the mechanism of Notch protein action; however, very little is known about the regulation of activity of the gene itself. Here, we provide the results of targeted 5'-end editing of the Drosophila Notch gene in its native environment and genetic and cytological effects of these changes. Using the Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR associated protein 9 (CRISPR/Cas9) system in combination with homologous recombination, we obtained a founder fly stock in which a 4-kb fragment, including the 5' nontranscribed region, the first exon, and a part of the first intron of Notch, was replaced by an attachment Phage (attP) site. Then, fly lines carrying a set of six deletions within the 5'untranscribed region of the gene were obtained by ΦC31-mediated integration of transgenic constructs. Part of these deletions does not affect gene activity, but their combinations with transgenic construct in the first intron of the gene cause defects in the Notch target tissues. At the polytene chromosome level we defined a DNA segment (~250 bp) in the Notch5'-nontranscribed region which when deleted leads to disappearance of the 3C6/C7 interband and elimination of CTC-Factor (CTCF) and Chromator (CHRIZ) insulator proteins in this region.

Gene networks have arisen as a promising tool in the comprehensive modeling and analysis of complex diseases. Particularly in viral infections, the understanding of the host-pathogen mechanisms, and the immune response to these, is considered a major goal for the rational design of appropriate therapies. For this reason, the use of gene networks may well encourage therapy-associated research in the context of the coronavirus pandemic, orchestrating experimental scrutiny and reducing costs. In this work, gene co-expression networks were reconstructed from RNA-Seq expression data with the aim of analyzing the time-resolved effects of gene Ly6E in the immune response against the coronavirus responsible for murine hepatitis (MHV). Through the integration of differential expression analyses and reconstructed networks exploration, significant differences in the immune response to virus were observed in Ly6E Δ H S C compared to wild type animals. Results show that Ly6E ablation at hematopoietic stem cells (HSCs) leads to a progressive impaired immune response in both liver and spleen. Specifically, depletion of the normal leukocyte mediated immunity and chemokine signaling is observed in the liver of Ly6E Δ H S C mice. On the other hand, the immune response in the spleen, which seemed to be mediated by an intense chromatin activity in the normal situation, is replaced by ECM remodeling in Ly6E Δ H S C mice. These findings, which require further experimental characterization, could be extrapolated to other coronaviruses and motivate the efforts towards novel antiviral approaches.

Telomeres are considered to be an internal biological clock, and their progressive shortening has been associated with the risk of age-related diseases and reproductive alterations. Over recent years, an increasing number of studies have focused on the association between telomere length and fertility, identifying sperm telomere length (STL) as a novel biomarker of male fertility. Although typically considered to be repeated DNA sequences, telomeres have recently been shown to also include a long non-coding RNA (lncRNA) known as TERRA (telomeric repeat-containing RNAs). Interestingly, males with idiopathic infertility show reduced testicular TERRA expression, suggesting a link between TERRA and male fertility. The aim of this study was to investigate the role of seminal TERRA expression in embryo quality. To this end, STL and TERRA expression were quantified by Real Time qPCR in the semen of 35 men who underwent assisted reproductive technologies (ART) and 30 fertile men. We found that TERRA expression in semen and STL was reduced in patients that underwent ART (both p < 0.001). Interestingly, TERRA and STL expressions were positively correlated (p = 0.010), and TERRA expression was positively associated with embryo quality (p < 0.001). These preliminary findings suggest a role for TERRA in the maintenance of sperm telomere integrity during gametogenesis, and for the first time, TERRA expression was found as a predictive factor for embryo quality in the setting of assisted reproduction.