The DNA sequence and biology of human chromosome 19

Range-limited Heaps' law for functional DNA words in the human genome

Heaps' or Herdan-Heaps' law is a linguistic law describing the relationship between the vocabulary/dictionary size (type) and word counts (token) to be a power-law function. Its existence in genomes with certain definition of DNA words is unclear partly because the dictionary size in genome could be much smaller than that in a human language. We define a DNA word as a coding region in a genome that codes for a protein domain. Using human chromosomes and chromosome arms as individual samples, we establish the existence of Heaps' law in the human genome within limited range. Our definition of words in a genomic or proteomic context is different from other definitions such as over-represented k-mers which are much shorter in length. Although an approximate power-law distribution of protein domain sizes due to gene duplication and the related Zipf's law is well known, their translation to the Heaps' law in DNA words is not automatic. Several other animal genomes are shown herein also to exhibit range-limited Heaps' law with our definition of DNA words, though with various exponents. When tokens were randomly sampled and sample sizes reach to the maximum level, a deviation from the Heaps' law was observed, but a quadratic regression in log-log type-token plot fits the data perfectly. Investigation of type-token plot and its regression coefficients could provide an alternative narrative of reusage and redundancy of protein domains as well as creation of new protein domains from a linguistic perspective.

Transcriptome analysis of cynomolgus macaques throughout their lifespan reveals age-related immune patterns

Despite the different perspectives by diverse research sectors spanning several decades, aging research remains uncharted territory for human beings. Therefore, we investigated the transcriptomic characteristics of eight male healthy cynomolgus macaques, and the annual sampling was designed with two individuals in four age groups. As a laboratory animal, the macaques were meticulously shielded from all environmental factors except aging. The results showed recent findings of certain immune response and the age-associated network of primate immunity. Three important aging patterns were identified and each gene clusters represented a different immune response. The increased expression pattern was predominantly associated with innate immune cells, such as Neutrophils and NK cells, causing chronic inflammation with aging whereas the other two decreased patterns were associated with adaptive immunity, especially "B cell activation" affecting antibody diversity of aging. Furthermore, the hub gene network of the patterns reflected transcriptomic age and correlated with human illness status, aiding in future human disease prediction. Our macaque transcriptome profiling results offer systematic insights into the age-related immunological features of primates.

High-throughput mRNA-seq atlas of human placenta shows vast transcriptome remodeling from first to third trimester†

The placenta, composed of chorionic villi, changes dramatically across gestation. Understanding differences in ongoing pregnancies are essential to identify the role of chorionic villi at specific times in gestation and develop biomarkers and prognostic indicators of maternal-fetal health. The normative mRNA profile is established using next-generation sequencing of 124 first trimester and 43 third trimester human placentas from ongoing healthy pregnancies. Stably expressed genes (SEGs) not different between trimesters and with low variability are identified. Differential expression analysis of first versus third trimester adjusted for fetal sex is performed, followed by a subanalysis with 23 matched pregnancies to control for subject variability using the same genetic and environmental background. Placenta expresses 14,979 polyadenylated genes above sequencing noise (transcripts per million > 0.66), with 10.7% SEGs across gestation. Differentially expressed genes (DEGs) account for 86.7% of genes in the full cohort [false discovery rate (FDR) < 0.05]. Fold changes highly correlate between the full cohort and subanalysis (Pearson = 0.98). At stricter thresholds (FDR < 0.001, fold change > 1.5), there remains 50.1% DEGs (3353 upregulated in first and 4155 upregulated in third trimester). This is the largest mRNA atlas of healthy human placenta across gestation, controlling for genetic and environmental factors, demonstrating substantial changes from first to third trimester in chorionic villi. Specific differences and SEGs may be used to understand the specific role of the chorionic villi throughout gestation and develop first trimester biomarkers of placental health that transpire across gestation, which can be used for future development of biomarkers for maternal-fetal health.

Clinical analysis and literature review of two paediatric cases of anti-IgLON5 antibody-related encephalitis

Introduction

Anti-IgLON5 antibody-related encephalitis is a rare autoimmune disorder of the central nervous system, predominantly occurring in middle-aged elderly individuals, with paediatric cases being exceptionally rare. This study aims to enhance the understanding of paediatric anti-IgLON5 antibody-related encephalitis by summarising its clinical and therapeutic characteristics.

Method

A retrospective analysis was conducted on two paediatric patients diagnosed with anti-IgLON5 antibody-related encephalitis at Hunan Children's Hospital from August 2022 to November 2023. This involved reviewing their medical records and follow-up data, in addition to a literature review.

Results

The study involved two patients, one male and one female, aged between 2.5 and 9.6 years, both presenting with an acute/subacute course of illness. Clinically, both exhibited movement disorders (including dystonia, involuntary movements, and ataxia), cognitive impairments, sleep disturbances, and psychiatric symptoms. Patient 1 experienced epileptic seizures, while Patient 2 exhibited brainstem symptoms and abnormal eye movements. Neither patient showed autonomic dysfunction. Patient 1 had normal cerebrospinal fluid (CSF) and Brain MRI findings, whereas Patient 2 showed moderate leukocytosis and mild protein elevation in the CSF, and Brain MRI revealed symmetrical lesions in the basal ganglia and cerebellum. Oligoclonal bands in the CSF were positive in both cases. Both patients tested negative for HLA-DQB*05:01 and HLA-DRB*10:01. They received both first-line and second-line immunotherapies, with Patient 2 showing a poor response to treatment.

Discussion

Paediatric cases of anti-IgLON5 antibody-related encephalitis similarly present sleep disturbances as a core symptom, alongside various forms of movement disorders. Immunotherapy is partially effective. Compared to adult patients, these paediatric cases tend to exhibit more pronounced psychiatric symptoms, a more rapid onset, and more evident inflammatory changes in the CSF. The condition appears to have a limited association with HLA-DQB*05:01 and HLA-DRB*10:01 polymorphisms.

FuncPEP v2.0: An Updated Database of Functional Short Peptides Translated from Non-Coding RNAs

Over the past decade, there have been reports of short novel functional peptides (less than 100 aa in length) translated from so-called non-coding RNAs (ncRNAs) that have been characterized using mass spectrometry (MS) and large-scale proteomics studies. Therefore, understanding the bivalent functions of some ncRNAs as transcripts that encode both functional RNAs and short peptides, which we named ncPEPs, will deepen our understanding of biology and disease. In 2020, we published the first database of functional peptides translated from non-coding RNAs-FuncPEP. Herein, we have performed an update including the newly published ncPEPs from the last 3 years along with the categorization of host ncRNAs. FuncPEP v2.0 contains 152 functional ncPEPs, out of which 40 are novel entries. A PubMed search from August 2020 to July 2023 incorporating specific keywords was performed and screened for publications reporting validated functional peptides derived from ncRNAs. We did not observe a significant increase in newly discovered functional ncPEPs, but a steady increase. The novel identified ncPEPs included in the database were characterized by a wide array of molecular and physiological parameters (i.e., types of host ncRNA, species distribution, chromosomal density, distribution of ncRNA length, identification methods, molecular weight, and functional distribution across humans and other species). We consider that, despite the fact that MS can now easily identify ncPEPs, there still are important limitations in proving their functionality.

Evolution of Transcript Abundance is Influenced by Indels in Protein Low Complexity Regions

Protein Protein low complexity regions (LCRs) are compositionally biased amino acid sequences, many of which have significant evolutionary impacts on the proteins which contain them. They are mutationally unstable experiencing higher rates of indels and substitutions than higher complexity regions. LCRs also impact the expression of their proteins, likely through multiple effects along the path from gene transcription, through translation, and eventual protein degradation. It has been observed that proteins which contain LCRs are associated with elevated transcript abundance (TAb), despite having lower protein abundance. We have gathered and integrated human data to investigate the co-evolution of TAb and LCRs through ancestral reconstructions and model inference using an approximate Bayesian calculation based method. We observe that on short evolutionary timescales TAb evolution is significantly impacted by changes in LCR length, with insertions driving TAb down. But in contrast, the observed data is best explained by indel rates in LCRs which are unaffected by shifts in TAb. Our work demonstrates a coupling between LCR and TAb evolution, and the utility of incorporating multiple responses into evolutionary analyses.

Genomic origin, fragmentomics, and transcriptional properties of long cell-free DNA molecules in human plasma

Recent studies have revealed an unexplored population of long cell-free DNA (cfDNA) molecules in human plasma using long-read sequencing technologies. However, the biological properties of long cfDNA molecules (>500 bp) remain largely unknown. To this end, we have investigated the origins of long cfDNA molecules from different genomic elements. Analysis of plasma cfDNA using long-read sequencing reveals an uneven distribution of long molecules from across the genome. Long cfDNA molecules show overrepresentation in euchromatic regions of the genome, in sharp contrast to short DNA molecules. We observe a stronger relationship between the abundance of long molecules and mRNA gene expression levels, compared with short molecules (Pearson's r = 0.71 vs. -0.14). Moreover, long and short molecules show distinct fragmentation patterns surrounding CpG sites. Leveraging the cleavage preferences surrounding CpG sites, the combined cleavage ratios of long and short molecules can differentiate patients with hepatocellular carcinoma (HCC) from non-HCC subjects (AUC = 0.87). We also investigated knockout mice in which selected nuclease genes had been inactivated in comparison with wild-type mice. The proportion of long molecules originating from transcription start sites are lower in Dffb-deficient mice but higher in Dnase1l3-deficient mice compared with that of wild-type mice. This work thus provides new insights into the biological properties and potential clinical applications of long cfDNA molecules.

The Risk Genes for Neuropsychiatric Disorders negr1 and opcml Are Expressed throughout Zebrafish Brain Development

The immunoglobulin LAMP/OBCAM/NTM (IgLON) family of cell adhesion molecules comprises five members known for their involvement in establishing neural circuit connectivity, fine-tuning, and maintenance. Mutations in IgLON genes result in alterations in these processes and can lead to neuropsychiatric disorders. The two IgLON family members NEGR1 and OPCML share common links with several of them, such as schizophrenia, autism, and major depressive disorder. However, the onset and the underlying molecular mechanisms have remained largely unresolved, hampering progress in developing therapies. NEGR1 and OPCML are evolutionarily conserved in teleosts like the zebrafish (Danio rerio), which is excellently suited for disease modelling and large-scale screening for disease-ameliorating compounds. To explore the potential applicability of zebrafish for extending our knowledge on NEGR1- and OPCML-linked disorders and to develop new therapeutic strategies, we investigated the spatio-temporal expression of the two genes during early stages of development. negr1 and opcml are expressed maternally and subsequently in partially distinct domains of conserved brain regions. Other areas of expression in zebrafish have not been reported in mammals to date. Our results indicate that NEGR1 and OPCML may play roles in neural circuit development and function at stages earlier than previously anticipated. A detailed functional analysis of the two genes based on our findings could contribute to understanding the mechanistic basis of related psychiatric disorders.

Implicit bias in diagnosing mosaicism amongst preimplantation genetic testing providers: results from a multicenter study of 36 395 blastocysts

STUDY QUESTION

Does the diagnosis of mosaicism affect ploidy rates across different providers offering preimplantation genetic testing for aneuploidies (PGT-A)?

SUMMARY ANSWER

Our analysis of 36 395 blastocyst biopsies across eight genetic testing laboratories revealed that euploidy rates were significantly higher in providers reporting low rates of mosaicism.

WHAT IS KNOWN ALREADY

Diagnoses consistent with chromosomal mosaicism have emerged as a third category of possible embryo ploidy outcomes following PGT-A. However, in the era of mosaicism, embryo selection has become increasingly complex. Biological, technical, analytical, and clinical complexities in interpreting such results have led to substantial variability in mosaicism rates across PGT-A providers and clinics. Critically, it remains unknown whether these differences impact the number of euploid embryos available for transfer. Ultimately, this may significantly affect clinical outcomes, with important implications for PGT-A patients.

STUDY DESIGN, SIZE, DURATION

In this international, multicenter cohort study, we reviewed 36 395 consecutive PGT-A results, obtained from 10 035 patients across 11 867 treatment cycles, conducted between October 2015 and October 2021. A total of 17 IVF centers, across eight PGT-A providers, five countries and three continents participated in the study. All blastocysts were tested using trophectoderm biopsy and next-generation sequencing. Both autologous and donation cycles were assessed. Cycles using preimplantation genetic testing for structural rearrangements were excluded from the analysis.

PARTICIPANTS/MATERIALS, SETTING, METHODS

The PGT-A providers were randomly categorized (A to H). Providers B, C, D, E, F, G, and H all reported mosaicism, whereas Provider A reported embryos as either euploid or aneuploid. Ploidy rates were analyzed using multilevel mixed linear regression. Analyses were adjusted for maternal age, paternal age, oocyte source, number of embryos biopsied, day of biopsy, and PGT-A provider, as appropriate. We compared associations between genetic testing providers and PGT-A outcomes, including the number of chromosomally normal (euploid) embryos determined to be suitable for transfer.

MAIN RESULTS AND THE ROLE OF CHANCE

The mean maternal age (±SD) across all providers was 36.2 (±5.2). Our findings reveal a strong association between PGT-A provider and the diagnosis of euploidy and mosaicism. Amongst the seven providers that reported mosaicism, the rates varied from 3.1% to 25.0%. After adjusting for confounders, we observed a significant difference in the likelihood of diagnosing mosaicism across providers (P < 0.001), ranging from 6.5% (95% CI: 5.2-7.4%) for Provider B to 35.6% (95% CI: 32.6-38.7%) for Provider E. Notably, adjusted euploidy rates were highest for providers that reported the lowest rates of mosaicism (Provider B: euploidy, 55.7% (95% CI: 54.1-57.4%), mosaicism, 6.5% (95% CI: 5.2-7.4%); Provider H: euploidy, 44.5% (95% CI: 43.6-45.4%), mosaicism, 9.9% (95% CI: 9.2-10.6%)); and Provider D: euploidy, 43.8% (95% CI: 39.2-48.4%), mosaicism, 11.0% (95% CI: 7.5-14.5%)). Moreover, the overall chance of having at least one euploid blastocyst available for transfer was significantly higher when mosaicism was not reported, when we compared Provider A to all other providers (OR = 1.30, 95% CI: 1.13-1.50). Differences in diagnosing and interpreting mosaic results across PGT-A laboratories raise further concerns regarding the accuracy and relevance of mosaicism predictions. While we confirmed equivalent clinical outcomes following the transfer of mosaic and euploid blastocysts, we found that a significant proportion of mosaic embryos are not used for IVF treatment.

LIMITATIONS, REASONS FOR CAUTION

Due to the retrospective nature of the study, associations can be ascertained, however, causality cannot be established. Certain parameters such as blastocyst grade were not available in the dataset. Furthermore, certain platform-related and clinic-specific factors may not be readily quantifiable or explicitly captured in our dataset. As such, a full elucidation of all potential confounders accounting for variability may not be possible.

WIDER IMPLICATIONS OF THE FINDINGS

Our findings highlight the strong need for standardization and quality assurance in the industry. The decision not to transfer mosaic embryos may ultimately reduce the chance of success of a PGT-A cycle by limiting the pool of available embryos. Until we can be certain that mosaic diagnoses accurately reflect biological variability, reporting mosaicism warrants utmost caution. A prudent approach is imperative, as it may determine the difference between success or failure for some patients.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the Torres Quevedo Grant, awarded to M.P. (PTQ2019-010494) by the Spanish State Research Agency, Ministry of Science and Innovation, Spain. M.P., L.B., A.R.L., A.L.R.d.C.L., N.P.P., M.P., D.S., F.A., A.P., B.M., L.D., F.V.M., D.S., M.R., E.P.d.l.B., A.R., and R.V. have no competing interests to declare. B.L., R.M., and J.A.O. are full time employees of IB Biotech, the genetics company of the Instituto Bernabeu group, which performs preimplantation genetic testing. M.G. is a full time employee of Novagen, the genetics company of Cegyr, which performs preimplantation genetic testing.

TRIAL REGISTRATION NUMBER

N/A.

Polyamine Dysregulation and Nucleolar Disruption in Alzheimer's Disease

A hypothesis of Alzheimer's disease etiology is proposed describing how cellular stress induces excessive polyamine synthesis and recycling which can disrupt nucleoli. Polyamines are essential in nucleolar functions, such as RNA folding and ribonucleoprotein assembly. Changes in the nucleolar pool of anionic RNA and cationic polyamines acting as counterions can cause significant nucleolar dynamics. Polyamine synthesis reduces S-adenosylmethionine which, at low levels, triggers tau phosphorylation. Also, polyamine recycling reduces acetyl-CoA needed for acetylcholine, which is low in Alzheimer's disease. Extraordinary nucleolar expansion and/or contraction can disrupt epigenetic control in peri-nucleolar chromatin, such as chromosome 14 with the presenilin-1 gene; chromosome 21 with the amyloid precursor protein gene; chromosome 17 with the tau gene; chromosome 19 with the APOE4 gene; and the inactive X chromosome (Xi; aka "nucleolar satellite") with normally silent spermine synthase (polyamine synthesis) and spermidine/spermine-N1-acetyltransferase (polyamine recycling) alleles. Chromosomes 17, 19 and the Xi have high concentrations of Alu elements which can be transcribed by RNA polymerase III if positioned nucleosomes are displaced from the Alu elements. A sudden flood of Alu RNA transcripts can competitively bind nucleolin which is usually bound to Alu sequences in structural RNAs that stabilize the nucleolar heterochromatic shell. This Alu competition leads to loss of nucleolar integrity with leaking of nucleolar polyamines that cause aggregation of phosphorylated tau. The hypothesis was developed with key word searches (e.g., PubMed) using relevant terms (e.g., Alzheimer's, lupus, nucleolin) based on a systems biology approach and exploring autoimmune disease tautology, gaining synergistic insights from other diseases.

Transformation-associated recombination (TAR) cloning and its applications for gene function; genome architecture and evolution; biotechnology and biomedicine

Transformation-associated recombination (TAR) cloning represents a unique tool to selectively and efficiently recover a given chromosomal segment up to several hundred kb in length from complex genomes (such as animals and plants) and simple genomes (such as bacteria and viruses). The technique exploits a high level of homologous recombination in the yeast Sacharomyces cerevisiae. In this review, we summarize multiple applications of the pioneering TAR cloning technique, developed previously for complex genomes, for functional, evolutionary, and structural studies, and extended the modified TAR versions to isolate biosynthetic gene clusters (BGCs) from microbes, which are the major source of pharmacological agents and industrial compounds, and to engineer synthetic viruses with novel properties to design a new generation of vaccines. TAR cloning was adapted as a reliable method for the assembly of synthetic microbe genomes for fundamental research. In this review, we also discuss how the TAR cloning in combination with HAC (human artificial chromosome)- and CRISPR-based technologies may contribute to the future.

Viral and host mediators of non-suppressible HIV-1 viremia

Non-suppressible HIV-1 viremia (NSV) is defined as persistent low-level viremia on antiretroviral therapy (ART) without evidence of ART non-adherence or significant drug resistance. Unraveling the mechanisms behind NSV would broaden our understanding of HIV-1 persistence. Here we analyzed plasma virus sequences in eight ART-treated individuals with NSV (88% male) and show that they are composed of large clones without evidence of viral evolution over time in those with longitudinal samples. We defined proviruses that match plasma HIV-1 RNA sequences as 'producer proviruses', and those that did not as 'non-producer proviruses'. Non-suppressible viremia arose from expanded clones of producer proviruses that were significantly larger than the genome-intact proviral reservoir of ART-suppressed individuals. Integration sites of producer proviruses were enriched in proximity to the activating H3K36me3 epigenetic mark. CD4+ T cells from participants with NSV demonstrated upregulation of anti-apoptotic genes and downregulation of pro-apoptotic and type I/II interferon-related pathways. Furthermore, participants with NSV showed significantly lower HIV-specific CD8+ T cell responses compared with untreated viremic controllers with similar viral loads. We identified potential critical host and viral mediators of NSV that may represent targets to disrupt HIV-1 persistence.

Sequencing-based transcriptome analysis reveals diversification of immune response- and angiogenesis-related expression patterns of early-stage cervical carcinoma as compared with high-grade CIN

Background

Molecular diversity of virus-associated cervical cancer remains a relatively underexplored issue, and interrelations of immunologic and angiogenic features during the establishment of a particular landscape of the cervical cancer microenvironment are not well-characterized, especially for its earliest clinical stages, although this may provide insight into the mechanisms behind the differences in tumor aggressiveness, treatment responsiveness and prognosis. In this research, we were aimed at identifying transcriptomic landscapes of early-stage cervical carcinoma that differ substantially in their immune-related characteristics, patterns of signaling pathways and composition of the microenvironment in comparison with immediate precursor (intraepithelial) lesions.

Methods

We performed the Illumina platform-based RNA sequencing using a panel of fresh tissue samples that included human papillomavirus-positive cervical intraepithelial neoplastic lesions (CIN), invasive squamous carcinoma of the cervix of FIGO IA1-IIB stages, and morphologically normal epithelium. The derived transcriptomic profiles were bioinformatically analyzed and compared by patterns of signaling pathway activation, distribution of tumor-infiltrating cell populations, and genomic regions involved.

Result

According to hierarchical cluster analysis of the whole-transcriptome profiles, tissue samples were distributed between three groups, or gene expression patterns (the one comprising most pre-cancer cases and the other two encompassing mostly early-stage invasive cancer cases). Differentially expressed genes were retrieved in each intergroup pairwise comparison followed by Gene Ontology analysis. Gene set enrichment analysis of the two groups of tumor samples in comparison with the CIN group identified substantial differences in immunological and angiogenic properties between tumorous groups suggesting the development of different molecular phenotypes. Cell composition analysis confirmed the diverse changes in the abundancies of immune and non-immune populations and, accordingly, different impacts of the immune and stromal compartments on the tumor microenvironment in these two groups of tumors compared to CIN. Positional gene expression analysis demonstrated that the identified transcriptomic differences were linked to different chromosomal regions and co-localized with particular gene families implicated in immune regulation, inflammation, cell differentiation, and tumor invasion.

Conclusions

Overall, detection of different transcriptomic patterns of invasive cervical carcinoma at its earliest stages supports the diverse impacts of immune response- and angiogenesis-related mechanisms on the onset of tumor invasion and progression. This may provide new options for broadening the applicability and increasing the efficiency of target anti-angiogenic and immune-based therapy of virus-associated cervical carcinoma.

Ancestry-related differences in chromatin accessibility and gene expression of APOE ε4 are associated with Alzheimer's disease risk

INTRODUCTION

European local ancestry (ELA) surrounding apolipoprotein E (APOE) ε4 confers higher risk for Alzheimer's disease (AD) compared to African local ancestry (ALA). We demonstrated significantly higher APOE ε4  expression in ELA versus ALA in AD brains from APOE ε4/ε4 carriers. Chromatin accessibility differences could contribute to these expression changes.

METHODS

We performed single nuclei assays for transposase accessible chromatin sequencing from the frontal cortex of six ALA and six ELA AD brains, homozygous for local ancestry and APOE ε4.

RESULTS

Our results showed an increased chromatin accessibility at the APOE ε4  promoter area in ELA versus ALA astrocytes. This increased accessibility in ELA astrocytes extended genome wide. Genes with increased accessibility in ELA in astrocytes were enriched for synapsis, cholesterol processing, and astrocyte reactivity.

DISCUSSION

Our results suggest that increased chromatin accessibility of APOE ε4  in ELA astrocytes contributes to the observed elevated APOE ε4  expression, corresponding to the increased AD risk in ELA versus ALA APOE ε4/ε4 carriers.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

The human placenta exhibits a unique transcriptomic void

The human placenta exhibits a unique genomic architecture with an unexpectedly high mutation burden and many uniquely expressed genes. The aim of this study is to identify transcripts that are uniquely absent or depleted in the placenta. Here, we show that 40 of 46 of the other organs have no selectively depleted transcripts and that, of the remaining six, the liver has the largest number, with 26. In contrast, the term placenta has 762 depleted transcripts. Gene Ontology analysis of this depleted set highlighted multiple pathways reflecting known unique elements of placental physiology. For example, transcripts associated with neuronal function are in the depleted set-as expected given the lack of placental innervation. However, this demonstrated overrepresentation of genes involved in mitochondrial function (p = 5.8 × 10-10), including PGC-1α, the master regulator of mitochondrial biogenesis, and genes involved in polyamine metabolism (p = 2.1 × 10-4).

High-throughput mRNA-seq atlas of human placenta shows vast transcriptome remodeling from first to third trimester

Background

The placenta, composed of chorionic villi, changes dramatically across gestation. Understanding differences in ongoing pregnancies are essential to identify the role of chorionic villi at specific times in gestation and develop biomarkers and prognostic indicators of maternal- fetal health.

Methods

The normative mRNA profile is established using next-generation sequencing of 124 first trimester and 43 third trimester human placentas from ongoing healthy pregnancies. Stably expressed genes not different between trimesters and with low variability are identified. Differential expression analysis of first versus third trimester adjusted for fetal sex is performed, followed by a subanalysis with 23 matched pregnancies to control for subject variability using the same genetic and environmental background.

Results

Placenta expresses 14,979 mRNAs above sequencing noise (TPM>0.66), with 1,545 stably expressed genes across gestation. Differentially expressed genes account for 86.7% of genes in the full cohort (FDR<0.05). Fold changes highly correlate between the full cohort and subanalysis (Pearson = 0.98). At stricter thresholds (FDR<0.001, fold change>1.5), there are 6,941 differentially expressed protein coding genes (3,206 upregulated in first and 3,735 upregulated in third trimester).

Conclusion

This is the largest mRNA atlas of healthy human placenta across gestation, controlling for genetic and environmental factors, demonstrating substantial changes from first to third trimester in chorionic villi. Specific differences and stably expressed genes may be used to understand the specific role of the chorionic villi throughout gestation and develop first trimester biomarkers of placental health that transpire across gestation, which can be used for future development of biomarkers in maternal-fetal disease.

RDA coupled with deep sequencing detects somatic SVA-retrotranspositions and mosaicism in the human brain

Cells of the developing human brain are affected by the progressive acquisition of genetic and epigenetic alterations that have been reported to contribute to somatic mosaicism in the adult brain and are increasingly considered a possible cause of neurogenetic disorders. A recent work uncovered that the copy-paste transposable element (TE) LINE-1 (L1) is mobilized during brain development, and thus mobile non-autonomous TEs like AluY and SINE-VNTR-Alu (SVA) families can use L1 activity in trans, leading to de novo insertions that may influence the variability of neural cells at genetic and epigenetic levels. In contrast to SNPs and when considering substitutional sequence evolution, the presence or absence of TEs at orthologous loci represents highly informative clade markers that provide insights into the lineage relationships between neural cells and how the nervous system evolves in health and disease. SVAs, as the 'youngest' class of hominoid-specific retrotransposons preferentially found in gene- and GC-rich regions, are thought to differentially co-regulate nearby genes and exhibit a high mobility in the human germline. Therefore, we determined whether this is reflected in the somatic brain and used a subtractive and kinetic enrichment technique called representational difference analysis (RDA) coupled with deep sequencing to compare different brain regions with respect to de novo SINE-VNTR-Alu insertion patterns. As a result, we detected somatic de novo SVA integrations in all human brain regions analyzed, and the majority of de novo insertions can be attributed to lineages of telencephalon and metencephalon, since most of the examined integrations are unique to different brain regions under scrutiny. The SVA positions were used as presence/absence markers, forming informative sites that allowed us to create a maximum parsimony phylogeny of brain regions. Our results largely recapitulated the generally accepted evo-devo patterns and revealed chromosome-wide rates of de novo SVA reintegration targets and preferences for specific genomic regions, e.g., GC- and TE-rich regions as well as close proximity to genes that tend to fall into neural-specific Gene Ontology pathways. We concluded that de novo SVA insertions occur in the germline and somatic brain cells at similar target regions, suggesting that similar retrotransposition modes are effective in the germline and soma.

Polymorphic Rearrangements of Human Chromosome 9 and Male Infertility: New Evidence and Impact on Spermatogenesis

Chromosomal polymorphisms are structural variations in chromosomes that define the genomic variance of a species. These alterations are recurrent in the general population, and some of them appear to be more recurrent in the infertile population. Human chromosome 9 is highly heteromorphic, and how its rearrangement affects male fertility remains to be fully investigated. In this study, we aimed to investigate the association between the polymorphic rearrangements of chromosome 9 and male infertility via an Italian cohort of male infertile patients. Cytogenetic analysis was carried out, along with Y microdeletion screening, semen analysis, fluorescence in situ hybridization, and TUNEL assays using spermatic cells. Chromosome 9 rearrangements were observed in six patients: three of them showed a pericentric inversion, while the others showed a polymorphic heterochromatin variant 9qh. Of these, four patients exhibited oligozoospermia associated with teratozoospermia, along with a percentage of aneuploidy in the sperm of above 9%, in particular, an increase in XY disomy. Additionally, high values for sperm DNA fragmentation (≥30%) were observed in two patients. None of them had microdeletions to the AZF loci on chromosome Y. Our results suggest that polymorphic rearrangements of chromosome 9 might be associated with abnormalities in sperm quality due to incorrect spermatogenesis regulation.

Plasma cell-free DNA methylome profiling in pre- and post-surgery oral cavity squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC), a highly heterogeneous disease that involves multiple anatomic sites, is a leading cause of cancer-related mortality worldwide. Although the utility of noninvasive biomarkers based on circulating cell-free DNA (cfDNA) methylation profiling has been widely recognized, limited studies have been reported so far regarding the dynamics of cfDNA methylome in oral cavity squamous cell carcinoma (OCSCC). It is hypothesized in this study that comparison of methylation profiles in pre- and postsurgery plasma samples will reveal OCSCC-specific prognostic and diagnostic biomarkers. As a strategy to further prioritize tumor-specific targets, top differential methylated regions (DMRs) were called by reanalyzing methylation data from paired tumor and normal tissue collected in the the cancer genome atlas head-neck squamous cell carcinoma (TCGA) head and neck cancer cohort. Matched plasma samples from eight patients with OCSCC were collected at Moffitt Cancer Center before and after surgical resection. Plasma-derived cfDNA was analyzed by cfMBD-seq, which is a high-sensitive methylation profiling assay. Differential methylation analysis was then performed based on the matched samples profiled. In the top 200 HNSCC-specific DMRs detected based on the TCGA data set, a total of 23 regions reached significance in the plasma-based DMR test. The top five validated DMR regions (ranked by the significance in the plasma study) are located in the promoter regions of genes PENK, NXPH1, ZIK1, TBXT, and CDO1, respectively. The genome-wide cfDNA DMR analysis further highlighted candidate biomarkers located in genes SFRP4, SOX1, IRF4, and PCDH17. The prognostic relevance of candidate genes was confirmed by survival analysis using the TCGA data. This study supports the utility of cfDNA-based methylome profiling as a promising noninvasive biomarker source for OCSCC and HNSCC.

A multi-task SCCA method for brain imaging genetics and its application in neurodegenerative diseases

BACKGROUND AND OBJECTIVES

In brain imaging genetics, multi-task sparse canonical correlation analysis (MTSCCA) is effective to study the bi-multivariate associations between genetic variations such as single nucleotide polymorphisms (SNPs) and multi-modal imaging quantitative traits (QTs). However, most existing MTSCCA methods are neither supervised nor capable of distinguishing the shared patterns of multi-modal imaging QTs from the specific patterns.

METHODS

A new diagnosis-guided MTSCCA (DDG-MTSCCA) with parameter decomposition and graph-guided pairwise group lasso penalty was proposed. Specifically, the multi-tasking modeling paradigm enables us to comprehensively identify risk genetic loci by jointly incorporating multi-modal imaging QTs. The regression sub-task was raised to guide the selection of diagnosis-related imaging QTs. To reveal the diverse genetic mechanisms, the parameter decomposition and different constraints were utilized to facilitate the identification of modality-consistent and -specific genotypic variations. Besides, a network constraint was added to find out meaningful brain networks. The proposed method was applied to synthetic data and two real neuroimaging data sets respectively from Alzheimer's disease neuroimaging initiative (ADNI) and Parkinson's progression marker initiative (PPMI) databases.

RESULTS

Compared with the competitive methods, the proposed method exhibited higher or comparable canonical correlation coefficients (CCCs) and better feature selection results. In particular, in the simulation study, DDG-MTSCCA showed the best anti-noise ability and achieved the highest average hit rate, about 25% higher than MTSCCA. On the real data of Alzheimer's disease (AD) and Parkinson's disease (PD), our method obtained the highest average testing CCCs, about 40% ∼ 50% higher than MTSCCA. Especially, our method could select more comprehensive feature subsets, and the top five SNPs and imaging QTs were all disease-related. The ablation experimental results also demonstrated the significance of each component in the model, i.e., the diagnosis guidance, parameter decomposition, and network constraint.

CONCLUSIONS

These results on simulated data, ADNI and PPMI cohorts suggested the effectiveness and generalizability of our method in identifying meaningful disease-related markers. DDG-MTSCCA could be a powerful tool in brain imaging genetics, worthy of in-depth study.

Viral and Host Mediators of Non-Suppressible HIV-1 Viremia

Non-suppressible HIV-1 viremia (NSV) can occur in persons with HIV despite adherence to combination antiretroviral therapy (ART) and in the absence of significant drug resistance. Here, we show that plasma NSV sequences are comprised primarily of large clones without evidence of viral evolution over time. We defined proviruses that contribute to plasma viremia as "producer", and those that did not as "non-producer". Compared to ART-suppressed individuals, NSV participants had a significantly larger producer reservoir. Producer proviruses were enriched in chromosome 19 and in proximity to the activating H3K36me3 epigenetic mark. CD4+ cells from NSV participants demonstrated upregulation of anti-apoptotic genes and downregulation of pro-apoptotic and type I/II interferon-related pathways. Furthermore, NSV participants showed no elevation in HIV-specific CD8+ cell responses and producer proviruses were enriched for HLA escape mutations. We identified critical host and viral mediators of NSV that represent potential targets to disrupt HIV persistence and promote viral silencing.

Age-related methylation changes in the human sperm epigenome

Advanced paternal age is associated with increased risks for reproductive and offspring medical problems. Accumulating evidence suggests age-related changes in the sperm epigenome as one underlying mechanism. Using reduced representation bisulfite sequencing on 73 sperm samples of males attending a fertility center, we identified 1,162 (74%) regions which were significantly (FDR-adjusted) hypomethylated and 403 regions (26%) being hypermethylated with age. There were no significant correlations with paternal BMI, semen quality, or ART outcome. The majority (1,152 of 1,565; 74%) of age-related differentially methylated regions (ageDMRs) were located within genic regions, including 1,002 genes with symbols. Hypomethylated ageDMRs were closer to transcription start sites than hypermethylated DMRs, half of which reside in gene-distal regions. In this and conceptually related genome-wide studies, so far 2,355 genes have been reported with significant sperm ageDMRs, however most (90%) of them in only one study. The 241 genes which have been replicated at least once showed significant functional enrichments in 41 biological processes associated with development and the nervous system and in 10 cellular components associated with synapses and neurons. This supports the hypothesis that paternal age effects on the sperm methylome affect offspring behaviour and neurodevelopment. It is interesting to note that sperm ageDMRs were not randomly distributed throughout the human genome; chromosome 19 showed a highly significant twofold enrichment with sperm ageDMRs. Although the high gene density and CpG content have been conserved, the orthologous marmoset chromosome 22 did not appear to exhibit an increased regulatory potential by age-related DNA methylation changes.

ATP5D Is a Potential Biomarker for Male Fertility

Background

Infertility is a global medical and social problem that affects human health and social development. At present, about 15% of couples of the right age in the world are infertile. As all we know, genetic defects are the most likely underlying cause of the pathology. ATP5D is also known as the delta subunit of mitochondrial ATP synthase. Mitochondria maintain sperm vitality, capacitation, acrosome reaction, and DNA integrity through ATP. Mitochondrial damage can trigger energy synthesis disorders, resulting in decreased sperm quality and function or even disappearance. The specific role of ATP5D in regulation of the male reproductive system remains elusive.

Methods

In this study, semen from normal and infertile males were collected and their indicators were examined by analysis of routine sperm parameters; ATP5D protein content in semen was examined by ELISA. Singer sequencing was used to detect whether there was a mutated of ATP5D in semen. Meanwhile, ATP5D knockout (KO) and knockin (KI) male mice were selected at 8-12 weeks of age and mated with adult wild-type (WT) female mice for more than two months to assess their fertility and reproductive ability. Morphological changes in tissues such as testes and epididymis were observed by HE staining; spermatozoa were taken from the epididymis of the mice; sperm counts were performed and morphological changes were observed by Diff-Quik staining.

Results

The results showed that the expression of ATP5D in infertile males was significantly lower than that in normal males (P < 0.001) and the normal morphology rate of spermatozoa was much lower than that of normal males, and the sequencing results showed no mutations. The animal reproductive experiments showed no significant changes in the number of fertility in KO/KI mice compared with WT mice, but the duration of fertility was significantly longer (P = 0.02). The testicular cells in KO mice were loosely arranged and disorganized, the lumen was larger, the interstitial cells were atrophied, and the number of spermatozoa was reduced and the malformation rate was higher in WT males. This suggests that ATP5D is an essential protein for sperm formation and fertility in male mice and may be used as a biomarker of male fertility.

Conclusion

This study found ATP5D correlated with male infertility and the expression levels were significantly reduced in the seminal plasma of all male infertile patients without gene mutations. KO male significantly prolonged fertility time and impaired testicular histomorphology. This suggests that ATP5D may be associated with spermatogenic function and fertility in male mice and may be used as a biomarker for male fertility. Future studies are required to elucidate the potential mechanisms. The trial registration number is KLL-2021-266.

Clinical and genomic delineation of the new proximal 19p13.3 microduplication syndrome

A small but growing body of scientific literature is emerging about clinical findings in patients with 19p13.3 microdeletion or duplication. Recently, a proximal 19p13.3 microduplication syndrome was described, associated with growth delay, microcephaly, psychomotor delay and dysmorphic features. The aim of our study was to better characterize the syndrome associated with duplications in the proximal 19p13.3 region (prox 19p13.3 dup), and to propose a comprehensive analysis of the underlying genomic mechanism. We report the largest cohort of patients with prox 19p13.3 dup through a collaborative study. We collected 24 new patients with terminal or interstitial 19p13.3 duplication characterized by array-based Comparative Genomic Hybridization (aCGH). We performed mapping, phenotype-genotype correlations analysis, critical region delineation and explored three-dimensional chromatin interactions by analyzing Topologically Associating Domains (TADs). We define a new 377 kb critical region (CR 1) in chr19: 3,116,922-3,494,377, GRCh37, different from the previously described critical region (CR 2). The new 377 kb CR 1 includes a TAD boundary and two enhancers whose common target is PIAS4. We hypothesize that duplications of CR 1 are responsible for tridimensional structural abnormalities by TAD disruption and misregulation of genes essentials for the control of head circumference during development, by breaking down the interactions between enhancers and the corresponding targeted gene.

Microdeletions at 19p13.11p12 in five individuals with neurodevelopmental delay

Only few copy number variants at chromosome 19p13.11 have been reported, thus associated clinical information is scarce. Proximal to these copy number losses, we now identified deletions in five unrelated individuals with neurodevelopmental disorders. They presented with psychomotor delay as well as behavioral and sleeping disorders, while complex cardiovascular, skeletal, and various other malformations were more variable. Dysmorphic features were rather unspecific and not considered as a recognizable gestalt. Neither of the analyzed parents carried their offsprings' deletions, indicating de novo occurrence. The deletion sizes ranged between 0.7 and 5.2 Mb, were located between 18 and 24 megabases from the telomere, and contained a variable number of protein-coding genes (n = 25-68). Although not all microdeletions shared a common region, the smallest common overlap of some of the deletions provided interesting insights in the chromosomal region 19p13.11p12. Diligent literature review using OMIM and Pubmed did not identify a satisfying candidate gene for neurodevelopmental disorders. In the literature, a de novo in-frame deletion in MAU2 was considered pathogenic in an individual with Cornelia de Lange syndrome. Therefore, the clinical differential diagnosis of this latter syndrome in one individual and the encompassment of MAU2 in three individuals' deletions suggest clinical and genetic overlap with this specific syndrome. Three of the four here reported individuals with deletion encompassing GDF1 had different congenital heart defects, suggesting that this gene's haploinsufficiency might contribute to the cardiovascular phenotype, however, with reduced penetrance. Our findings indicate an association of microdeletions at 19p13.11/ 19p13.11p12 with neurodevelopmental disorders, variable symptoms, and malformations, and delineate the phenotypic spectrum of deletions within this genomic region.

A case with cardiac, skeletal, speech, and motor traits narrows the subtelomeric 19p13.3 microdeletion region to 46 kb

Subtelomeric 19p13.3 deletions have been associated with diverse anatomical and developmental phenotypes. A recent study of eight patients with subtelomeric interstitial 19p13.3 microdeletions at 0.3-1.4 Mb (hg 19) showed associations with growth restrictions, skeletal deformities, craniofacial anomalies, congenital heart defects, renal malformations, hernias, immune system deficits, fine and gross motor delays, speech delays, and developmental and learning delays. The authors defined two small regions of overlap containing four and 11 genes, respectively, with potential haploinsufficiency. Here, we present a new case with a de novo 184 kb deletion containing eight genes, three of which fall into the second previously identified small region of overlap, reducing the shared region to 46 kb. Phenotypic traits include most of the core findings in the previously reported cases but not growth restrictions, craniofacial anomalies, renal malformation, and learning disability. A closer look at the speech and motor delays reveals apraxic speech and discoordination in the fine and gross motor domain, consistent with cerebellar involvement across motor systems. Findings are consistent with a role of AZU1 in the observed immune deficiencies and PTBP1 in the observed skeletal, abdominal, speech, language, motor, and sensory traits. This case thus contributes to a more nuanced understanding of the subtelomeric 19p13.3 deletion region.

Maternal Copy Number Imbalances in Non-Invasive Prenatal Testing: Do They Matter?

Non-invasive prenatal testing (NIPT) has become a routine practice in screening for common aneuploidies of chromosomes 21, 18, and 13 and gonosomes X and Y in fetuses worldwide since 2015 and has even expanded to include smaller subchromosomal events. In fact, the fetal fraction represents only a small proportion of cell-free DNA on a predominant background of maternal DNA. Unlike fetal findings that have to be confirmed using invasive testing, it has been well documented that NIPT provides information on maternal mosaicism, occult malignancies, and hidden health conditions due to copy number variations (CNVs) with diagnostic resolution. Although large duplications or deletions associated with certain medical conditions or syndromes are usually well recognized and easy to interpret, very little is known about small, relatively common copy number variations on the order of a few hundred kilobases and their potential impact on human health. We analyzed data from 6422 NIPT patient samples with a CNV detection resolution of 200 kb for the maternal genome and identified 942 distinct CNVs; 328 occurred repeatedly. We defined them as multiple occurring variants (MOVs). We scrutinized the most common ones, compared them with frequencies in the gnomAD SVs v2.1, dbVar, and DGV population databases, and analyzed them with an emphasis on genomic content and potential association with specific phenotypes.

Paralogous Genes Involved in Embryonic Development: Lessons from the Eye and Other Tissues

During evolution, gene duplications lead to a naturally increased gene dosage. Duplicated genes can be further retained or eliminated over time by purifying selection pressure. The retention probability is increased by functional diversification and by the acquisition of novel functions. Interestingly, functionally diverged paralogous genes can maintain a certain level of functional redundancy and at least a partial ability to replace each other. In such cases, diversification probably occurred at the level of transcriptional regulation. Nevertheless, some duplicated genes can maintain functional redundancy after duplication and the ability to functionally compensate for the loss of each other. Many of them are involved in proper embryonic development. The development of particular tissues/organs and developmental processes can be more or less sensitive to the overall gene dosage. Alterations in the gene dosage or a decrease below a threshold level may have dramatic phenotypic consequences or even lead to embryonic lethality. The number of functional alleles of particular paralogous genes and their mutual cooperation and interactions influence the gene dosage, and therefore, these factors play a crucial role in development. This review will discuss individual interactions between paralogous genes and gene dosage sensitivity during development. The eye was used as a model system, but other tissues are also included.

Current Advances in Adeno-Associated Virus-Mediated Gene Therapy to Prevent Acquired Hearing Loss

Adeno-associated viruses (AAVs) are viral vectors that offer an excellent platform for gene therapy due to their safety profile, persistent gene expression in non-dividing cells, target cell specificity, lack of pathogenicity, and low immunogenicity. Recently, gene therapy for genetic hearing loss with AAV transduction has shown promise in animal models. However, AAV transduction for gene silencing or expression to prevent or manage acquired hearing loss is limited. This review provides an overview of AAV as a leading gene delivery vector for treating genetic hearing loss in animal models. We highlight the advantages and shortcomings of AAV for investigating the mechanisms and preventing acquired hearing loss. We predict that AAV-mediated gene manipulation will be able to prevent acquired hearing loss.

Maternal obesity alters methylation level of cytosine in CpG island for epigenetic inheritance in fetal umbilical cord blood

BACKGROUND

Over the past few decades, global maternal obesity prevalence has rapidly increased. This condition may induce long-lasting pathophysiological effects on either fetal or infant health that could be attributable to unknown unique changes in the umbilical blood composition.

METHODS

A total of 34 overweight/obese and 32 normal-weight pregnant women were recruited. Fifteen umbilical blood samples including 8 overweight/obese subjects and 7 normal weight women were sequenced using Targeted Bisulfite Sequencing technology to detect the average methylation level of cytosine and identify the differentially methylated region (DMR). GO and KEGG analyses were then employed to perform pathway enrichment analysis of DMR-related genes and promoters. Moreover, the mRNA levels of methylation-related genes histone deacetylases (HDACs) and DNA methyltransferases (DNMTs) were characterized in the samples obtained from these two groups.

RESULTS

Average methylated cytosine levels in both the CpG islands (CGI) and promoter significantly decreased in overweight/obese groups. A total of 1669 DMRs exhibited differences in their DNA methylation status between the overweight/obese and control groups. GO and KEGG analyses revealed that DMR-related genes and promoters were enriched in the metabolism, cancer and cardiomyopathy signaling pathways. Furthermore, the HDACs and DNMTs mRNA levels trended to decline in overweight/obese groups.

CONCLUSIONS

Decreased methylated cytosine levels in overweight/obese women induce the gene expression activity at a higher level than in the control group. DMRs between these two groups in the fetal blood may contribute to the changes in gene transcription that underlie the increased risk of metabolic disorders, cancers and cardiomyopathy in their offspring.

Functional validation of variants of unknown significance using CRISPR gene editing and transcriptomics: A Kleefstra syndrome case study

There are an estimated > 400 million people living with a rare disease globally, with genetic variants the cause of approximately 80% of cases. Next Generation Sequencing (NGS) rapidly identifies genetic variants however they are often of unknown significance. Low throughput functional validation in specialist laboratories is the current ad hoc approach for functional validation of genetic variants, which creating major bottlenecks in patient diagnosis. This study investigates the application of CRISPR gene editing followed by genome wide transcriptomic profiling to facilitate patient diagnosis. As proof-of-concept, we introduced a variant in the Euchromatin histone methyl transferase (EHMT1) gene into HEK293T cells. We identified changes in the regulation of the cell cycle, neural gene expression and suppression of gene expression changes on chromosome 19 and chromosome X, that are in keeping with Kleefstra syndrome clinical phenotype and/or provide insight into disease mechanism. This study demonstrates the utility of genome editing followed by functional readouts to rapidly and systematically validating the function of variants of unknown significance in patients suffering from rare diseases.

Comprehensive analysis of the prognostic and role in immune cell infiltration of MSR1 expression in lower-grade gliomas

BACKGROUND

The therapeutic effects of conventional treatment on gliomas are not promising. The tumor microenvironment (TME) has a close association with the invasiveness of multiple types of tumors, including low-grade gliomas (LGG). This study aims to validate the prognostic and immune-related role of macrophage scavenger receptor 1 (MSR1) in LGG patients.

METHODS

Data in this study were obtained from public databases. The differential expression of MSR1 was analyzed in LGG patients with different clinicopathological characteristics. Kaplan-Meier survival analysis, a time-dependent receiver operating characteristic (ROC) curve, and Cox regression analysis were used to assess the prognostic value of MSR1. Differentially expressed genes (DEGs) were screened between the high and low expression groups of MSR1. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to annotate the function of these DEGs. Hallmark gene sets were identified based on MSR1 by Gene Set Enrichment Analysis (GSEA). Difference analysis and correlation analysis were used to study the relationship between MSR1 and TME-related scores, tumor-infiltrating immune cells (TIICs), immune-related gene sets, and immune checkpoints (ICPs). The single-cell sequencing data were processed to identify the cell types expressing MSR1. The quantification of TIICs in TME was calculated by single-sample gene set enrichment analysis (ssGSEA). The differential expression of MSR1 in LGG and control brain tissues was verified by experiments.

RESULTS

There were significant differences in the expression level of MSR1 in different types of tissues and cells. MSR1 has a high prognostic value in LGG patients and can be used as an independent prognostic factor. MSR1 is closely related to TME and may play an important role in the immunotherapy of LGG patients.

CONCLUSIONS

The result of our study demonstrated that MSR1 is an independent prognostic biomarker in LGG patients and may play an important role in the TME of LGGs.

Recombination, selection and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.

Copy Number Variation Analysis of Euploid Pregnancy Loss

Objectives: Copy number variant (CNV) is believed to be the potential genetic cause of pregnancy loss. However, CNVs less than 3 Mb in euploid products of conceptions (POCs) remain largely unexplored. The aim of this study was to investigate the features of CNVs less than 3 Mb in POCs and their potential clinical significance in pregnancy loss/fetal death.

Methods: CNV data were extracted from a cohort in our institution and 19 peer-reviewed publications, and only those CNVs less than 3 Mb detected in euploid pregnancy loss/fetal death were included. We conducted a CNV map to analyze the distribution of CNVs in chromosomes using R packages karyoploteR_1.10.5. Gene names and annotated gene types covered by those CNVs were mined from the human Release 19 reference genome file and GENECODE database. We assessed the expression patterns and the consequences of murine knock-out of those genes using TiGER and Mouse Genome Informatics (MGI) databases. Functional enrichment and pathway analysis for genes in CNVs were performed using clusterProfiler V3.12.0.

Result: Breakpoints of 564 CNVs less than 3 Mb were obtained from 442 euploid POCs, with 349 gains and 185 losses. The CNV map showed that CNVs were distributed in all chromosomes, with the highest frequency detected in chromosome 22 and the lowest frequency in chromosome Y, and CNVs showed a higher density in the pericentromeric and sub-telomeric regions. A total of 5,414 genes mined from the CNV regions (CNVRs), Gene Ontology (GO), and pathway analysis showed that the genes were significantly enriched in multiple terms, especially in sensory perception, membrane region, and tight junction. A total of 995 protein-coding genes have been reported to present mammalian phenotypes in MGI, and 276 of them lead to embryonic lethality or abnormal embryo/placenta in knock-out mouse models. CNV located at 19p13.3 was the most common CNV of all POCs.

Conclusion: CNVs less than 3 Mb in euploid POCs distribute unevenly in all chromosomes, and a higher density was seen in the pericentromeric and sub-telomeric regions. The genes in those CNVRs are significantly enriched in biological processes and pathways that are important to embryonic/fetal development. CNV in 19p13.3 and the variations of ARID3A and FSTL3 might contribute to pregnancy loss.

Preparative capillary electrophoresis (CE) fractionation of protein digests improves protein and peptide identification in bottom-up proteomics

Reversed-phase liquid chromatography (RPLC) is widely used to reduce sample complexity prior to mass spectrometry (MS) analysis in bottom-up proteomics. Improving peptide separation in complex samples enables lower-abundance proteins to be identified. Multidimensional separations that combine orthogonal separation modes improve protein and peptide identifications over RPLC alone. Here we report a preparative capillary electrophoresis (CE) fractionation method that combines CE and RPLC separations. Using this method, we demonstrate improved protein and peptide identification in a tryptic digest of E. coli cell lysate, with 132 ± 33% more protein identifications and 185 ± 65% more peptide identifications over non-fractionated samples. Fractionation enables detection of lower-abundance proteins in this complex sample. We demonstrate improved coverage of ovarian cancer biomarker MUC16 isolated from conditioned cell media, with 6.73% sequence coverage using CE fractionation compared to 2.74% coverage without preparative fractionation. This new method will allow researchers performing bottom-up proteomics to harness the advantages of CE separations while using widely available LC-MS/MS instrumentation.

Crosstalk between Hepatitis B Virus and the 3D Genome Structure

Viruses that transcribe their DNA within the nucleus have to adapt to the existing cellular mechanisms that govern transcriptional regulation. Recent technological breakthroughs have highlighted the highly hierarchical organization of the cellular genome and its role in the regulation of gene expression. This review provides an updated overview on the current knowledge on how the hepatitis B virus interacts with the cellular 3D genome and its consequences on viral and cellular gene expression. We also briefly discuss the strategies developed by other DNA viruses to co-opt and sometimes subvert cellular genome spatial organization.

Circular DNA in the human germline and its association with recombination

Extrachromosomal circular DNA (eccDNA) is common in somatic tissue, but its existence and effects in the human germline are unexplored. We used microscopy, long-read DNA sequencing, and new analytic methods to document thousands of eccDNAs from human sperm. EccDNAs derived from all genomic regions and mostly contained a single DNA fragment, although some consisted of multiple fragments. The generation of eccDNA inversely correlates with the meiotic recombination rate, and chromosomes with high coding-gene density and Alu element abundance form the least eccDNA. Analysis of insertions in human genomes further indicates that eccDNA can persist in the human germline when the circular molecules reinsert themselves into the chromosomes. Our results suggest that eccDNA has transient and permanent effects on the germline. They explain how differences in the physical and genetic map might arise and offer an explanation of how Alu elements coevolved with genes to protect genome integrity against deleterious mutations producing eccDNA.

Mitotic Antipairing of Homologous Chromosomes

Chromosome organization is highly dynamic and plays an essential role during cell function. It was recently found that pairs of the homologous chromosomes are continuously separated at mitosis and display a haploid (1n) chromosome set, or "antipairing," organization in human cells. Here, we provide an introduction to the current knowledge of homologous antipairing in humans and its implications in human disease.

Pathophysiological role of growth differentiation factor 15 (GDF15) in obesity, cancer, and cachexia

Growth differentiation factor 15 or macrophage inhibitory cytokine-1 (GDF15/MIC-1) is a divergent member of the transforming growth factor β superfamily and has a diverse pathophysiological roles in cancers, cardiometabolic disorders, and other diseases. GDF15 controls hematopoietic growth, energy homeostasis, adipose tissue metabolism, body growth, bone remodeling, and response to stress signals. The role of GDF15 in cancer development and progression is complicated and depends on the specific cancer type, stage, and tumor microenvironment. Recently, research on GDF15 and GDF15-associated signaling has accelerated due to the identification of the GDF15 receptor: glial cell line-derived neurotrophic factor (GDNF) family receptor α-like (GFRAL). Therapeutic interventions to target GDF15 and/or GFRAL revealed the mechanisms that drive its activity and might improve overall outcomes of patients with metabolic disorders and cancer. This review highlights the structure and functions of GDF15 and its receptor, emphasizing the pleiotropic role of GDF15 in obesity, tumorigenesis, metastasis, immunomodulation, and cachexia.

Untangling the clinicopathological significance of MRE11-RAD50-NBS1 complex in sporadic breast cancers

The MRE11-RAD50-NBS1 (MRN) complex is critical for genomic stability. Although germline mutations in MRN may increase breast cancer susceptibility, such mutations are extremely rare. Here, we have conducted a comprehensive clinicopathological study of MRN in sporadic breast cancers. We have protein expression profiled for MRN and a panel of DNA repair factors involved in double-strand break repair (BRCA1, BRCA2, ATM, CHK2, ATR, Chk1, pChk1, RAD51, γH2AX, RPA1, RPA2, DNA-PKcs), RECQ DNA helicases (BLM, WRN, RECQ1, RECQL4, RECQ5), nucleotide excision repair (ERCC1) and base excision repair (SMUG1, APE1, FEN1, PARP1, XRCC1, Pol β) in 1650 clinical breast cancers. The prognostic significance of MRE11, RAD50 and NBS1 transcripts and their microRNA regulators (hsa-miR-494 and hsa-miR-99b) were evaluated in large clinical datasets. Expression of MRN components was analysed in The Cancer Genome Atlas breast cancer cohort. We show that low nuclear MRN is linked to aggressive histopathological phenotypes such as high tumour grade, high mitotic index, oestrogen receptor- and high-risk Nottingham Prognostic Index. In univariate analysis, low nuclear MRE11 and low nuclear RAD50 were associated with poor survival. In multivariate analysis, low nuclear RAD50 remained independently linked with adverse clinical outcomes. Low RAD50 transcripts were also linked with reduced survival. In contrast, overexpression of hsa-miR-494 and hsa-miR-99b microRNAs was associated with poor survival. We observed large-scale genome-wide alterations in MRN-deficient tumours contributing to aggressive behaviour. We conclude that MRN status may be a useful tool to stratify tumours for precision medicine strategies.

A catalog of curated breast cancer genes

Purpose

Decades of research have identified multiple genetic variants associated with breast cancer etiology. However, there is no database that archives breast cancer genes and variants responsible for predisposition. We set out to build a dynamic repository of curated breast cancer genes.

Methods

A comprehensive literature search was performed in PubMed and Google Scholar, followed by data extraction and harmonization for downstream analysis.

Results

Using a subset of 345 studies, we cataloged 652 breast cancer-associated loci across the genome. A majority of these were present in the non-coding region (i.e., intergenic (101) and intronic (345)), whereas only 158 were located within an exon. Using the odds ratio, we identified 429 loci to increase the disease risk and 198 to confer protection against breast cancer, whereas 25 were identified to both increase disease risk and confer protection against breast cancer. Chromosomal ideogram analysis indicated that chromosomes 17 and 19 have the highest density of breast cancer loci. We manually annotated and collated breast cancer genes in which a previous association between rare-monogenic variant and breast cancer has been documented. Finally, network and functional enrichment analysis revealed that steroid metabolism and DNA repair pathways were predominant among breast cancer genes and variants.

Conclusions

We have built an online interactive catalog of curated breast cancer genes (https://cbcg.dk). This will expedite clinical diagnostics and support the ongoing efforts in managing breast cancer etiology. Moreover, the database will serve as an essential repository when designing new breast cancer multigene panels.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10549-021-06441-y.

The role of distinct BRD4 isoforms and their contribution to high-grade serous ovarian carcinoma pathogenesis

High-grade serous ovarian carcinoma (HGSOC) is the most aggressive type of ovarian cancer, often diagnosed at advanced stages. Molecularly, HGSOC shows high degree of genomic instability associated with large number of genetic alterations. BRD4 is the 4th most amplified gene in HGSOC, which correlates with poor patients' prognosis. BRD4 is constitutively expressed and generates two proteins, BRD4 long (BRD4-L) and BRD4 short (BRD4-S). Both isoforms contain bromodomains that bind to lysine-acetylated histones. Amongst other functions, BRD4 participates in chromatin organization, acetylation of histones, transcriptional control and DNA damage repair. In cancer patients with amplified BRD4, the increased activity of BRD4 is associated with higher expression of oncogenes, such as MYC, NOTCH3 and NRG1. BRD4-driven oncogenes promote increased tumor cells proliferation, genetic instability, epithelial-mesenchymal transition, metastasis and chemoresistance. Ablation of BRD4 activity can be successfully achieved with bromodomain inhibitors (BETi) and degraders, and it has been applied in pre-clinical and clinical settings. Inhibition of BRD4 function has an effective anti-cancer effect, reducing tumor growth whether ablated by single agents or in combination with other drugs. When combined with standard chemotherapy, BETi are capable of sensitizing highly resistant ovarian cancer cell lines to platinum drugs. Despite the evidence that BRD4 amplification in ovarian cancer contributes to poor patient prognosis, little is known about the specific mechanisms by which BRD4 drives tumor progression. In addition, newly emerging data revealed that BRD4 isoforms exhibit contradicting functions in cancer. Therefore, it is paramount to expand studies elucidating distinct roles of BRD4-L and BRD4-S in HGSOC, which has important implications on development of therapeutic approaches targeting BRD4.

During evolution from the earliest tetrapoda, newly-recruited genes are increasingly paralogues of existing genes and distribute non-randomly among the chromosomes

BACKGROUND

The present availability of full genome sequences of a broad range of animal species across the whole range of evolutionary history enables one to ask questions as to the distribution of genes across the chromosomes. Do newly recruited genes, as new clades emerge, distribute at random or at non-random locations?

RESULTS

We extracted values for the ages of the human genes and for their current chromosome locations, from published sources. A quantitative analysis showed that the distribution of newly-added genes among and within the chromosomes appears to be increasingly non-random if one observes animals along the evolutionary series from the precursors of the tetrapoda through to the great apes, whereas the oldest genes are randomly distributed.

CONCLUSIONS

Randomization will result from chromosome evolution, but less and less time is available for this process as evolution proceeds. Much of the bunching of recently-added genes arises from new gene formation as paralogues in gene families, near the location of genes that were recruited in the preceding phylostratum. As examples we cite the KRTAP, ZNF, OR and some minor gene families. We show that bunching can also result from the evolution of the chromosomes themselves when, as for the KRTAP genes, blocks of genes that had previously been on disparate chromosomes become linked together.

Detection of Neoplasms by Metagenomic Next-Generation Sequencing of Cerebrospinal Fluid

Importance

Cerebrospinal fluid (CSF) cytologic testing and flow cytometry are insensitive for diagnosing neoplasms of the central nervous system (CNS). Such clinical phenotypes can mimic infectious and autoimmune causes of meningoencephalitis.

Objective

To ascertain whether CSF metagenomic next-generation sequencing (mNGS) can identify aneuploidy, a hallmark of malignant neoplasms, in difficult-to-diagnose cases of CNS malignant neoplasm.

Design, Setting, and Participants

Two case-control studies were performed at the University of California, San Francisco (UCSF). The first study used CSF specimens collected at the UCSF Clinical Laboratories between July 1, 2017, and December 31, 2019, and evaluated test performance in specimens from patients with a CNS malignant neoplasm (positive controls) or without (negative controls). The results were compared with those from CSF cytologic testing and/or flow cytometry. The second study evaluated patients who were enrolled in an ongoing prospective study between April 1, 2014, and July 31, 2019, with presentations that were suggestive of neuroinflammatory disease but who were ultimately diagnosed with a CNS malignant neoplasm. Cases of individuals whose tumors could have been detected earlier without additional invasive testing are discussed.

Main Outcomes and Measures

The primary outcome measures were the sensitivity and specificity of aneuploidy detection by CSF mNGS. Secondary subset analyses included a comparison of CSF and tumor tissue chromosomal abnormalities and the identification of neuroimaging characteristics that were associated with test performance.

Results

Across both studies, 130 participants were included (median [interquartile range] age, 57.5 [43.3-68.0] years; 72 men [55.4%]). The test performance study used 125 residual laboratory CSF specimens from 47 patients with a CNS malignant neoplasm and 56 patients with other neurological diseases. The neuroinflammatory disease study enrolled 12 patients and 17 matched control participants. The sensitivity of the CSF mNGS assay was 75% (95% CI, 63%-85%), and the specificity was 100% (95% CI, 96%-100%). Aneuploidy was detected in 64% (95% CI, 41%-83%) of the patients in the test performance study with nondiagnostic cytologic testing and/or flow cytometry, and in 55% (95% CI, 23%-83%) of patients in the neuroinflammatory disease study who were ultimately diagnosed with a CNS malignant neoplasm. Of the patients in whom aneuploidy was detected, 38 (90.5%) had multiple copy number variations with tumor fractions ranging from 31% to 49%.

Conclusions and Relevance

This case-control study showed that CSF mNGS, which has low specimen volume requirements, does not require the preservation of cell integrity, and was orginally developed to diagnose neurologic infections, can also detect genetic evidence of a CNS malignant neoplasm in patients in whom CSF cytologic testing and/or flow cytometry yielded negative results with a low risk of false-positive results.

Developing ZNF Gene Signatures Predicting Radiosensitivity of Patients with Breast Cancer

Adjuvant radiotherapy is one of the main treatment methods for breast cancer, but its clinical benefit depends largely on the characteristics of the patient. This study aimed to explore the relationship between the expression of zinc finger (ZNF) gene family proteins and the radiosensitivity of breast cancer patients. Clinical and gene expression data on a total of 976 breast cancer samples were obtained from The Cancer Genome Atlas (TCGA) database. ZNF gene expression was dichotomized into groups with a higher or lower level than the median level of expression. Univariate and multivariate Cox regression analyses were used to evaluate the relationship between ZNF gene expression levels and radiosensitivity. The Molecular Taxonomy Data of the International Federation of Breast Cancer (METABRIC) database was used for validation. The results revealed that 4 ZNF genes were possible radiosensitivity markers. High expression of ZNF644 and low expression levels of the other 3 genes (ZNF341, ZNF541, and ZNF653) were related to the radiosensitivity of breast cancer. Hierarchical cluster, Cox, and CoxBoost analysis based on these 4 ZNF genes indicated that patients with a favorable 4-gene signature had better overall survival on radiotherapy. Thus, this 4-gene signature may have value for selecting those patients most likely to benefit from radiotherapy. ZNF gene clusters could act as radiosensitivity signatures for breast cancer patients and may be involved in determining the radiosensitivity of cancer.

Altered DNA Methylation Profiles in SF3B1 Mutated CLL Patients

Mutations in splicing factor genes have a severe impact on the survival of cancer patients. Splicing factor 3b subunit 1 (SF3B1) is one of the most frequently mutated genes in chronic lymphocytic leukemia (CLL); patients carrying these mutations have a poor prognosis. Since the splicing machinery and the epigenome are closely interconnected, we investigated whether these alterations may affect the epigenomes of CLL patients. While an overall hypomethylation during CLL carcinogenesis has been observed, the interplay between the epigenetic stage of the originating B cells and SF3B1 mutations, and the subsequent effect of the mutations on methylation alterations in CLL, have not been investigated. We profiled the genome-wide DNA methylation patterns of 27 CLL patients with and without SF3B1 mutations and identified local decreases in methylation levels in SF3B1^mut CLL patients at 67 genomic regions, mostly in proximity to telomeric regions. These differentially methylated regions (DMRs) were enriched in gene bodies of cancer-related signaling genes, e.g., NOTCH1, HTRA3, and BCL9L. In our study, SF3B1 mutations exclusively emerged in two out of three epigenetic stages of the originating B cells. However, not all the DMRs could be associated with the methylation programming of B cells during development, suggesting that mutations in SF3B1 cause additional epigenetic aberrations during carcinogenesis.

Highly Selective, CRISPR/Cas9-Mediated Isolation of Genes and Genomic Loci from Complex Genomes by TAR Cloning in Yeast

Here we describe an updated TAR cloning protocol for the selective and efficient isolation of any genomic fragment or gene of interest up to 280 kb in size from genomic DNA. The method exploits the special recombination machinery of the yeast Saccharomyces cerevisiae. TAR cloning is based on the high level of in vivo recombination that occurs between a specific genomic DNA fragment of interest and targeting sequences (hooks) in a TAR vector that are homologous to the 5' and 3' ends of the targeted region. Upon co-transformation into yeast, this results in the isolation of the chromosomal region of interest as a circular YAC molecule, which then propagates and segregates in yeast cells and can be selected for. In the updated TAR cloning protocol described here, the fraction of region-positive clones typically obtained is increased from 1% up to 35% by pre-treatment of the genomic DNA with specifically designed CRISPR/Cas9 endonucleases that create double-strand breaks (DSBs) bracketing the target genomic DNA sequence, thereby making the ends of the chromosomal region of interest highly recombinogenic. In addition, a new TAR vector was constructed that contains YAC and BAC cassettes, permitting direct transfer of a TAR-cloned DNA from yeast to bacterial cells. Once the TAR vector with the hooks is constructed and genomic DNA is prepared, the entire procedure takes 3 weeks to complete. The updated TAR protocol does not require significant yeast experience or extensively time-consuming yeast work because screening only about a dozen yeast transformants is typically enough to find a clone with the region of interest. TAR cloning of chromosomal fragments, individual genes, or gene families can be used for functional, structural, and population studies, for comparative genomics, and for long-range haplotyping, and has potential for gene therapy. Published 2021. This article is a U.S. Government work and is in the public domain in the USA. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Preparation of CRISPR/Cas9-treated genomic DNA for TAR cloning Basic Protocol 2: Isolation of a gene or genomic locus by TAR cloning Basic Protocol 3: Transfer of TAR/YAC/BAC isolates from yeast to E. coli.

Anti-IgLON5 disease in a pediatric patient with Langerhans cell histiocytosis

INTRODUCTION

Anti-IgLON5 disease is a rare autoimmune disease of the central nervous system. Different from the previous autoimmune encephalitis, the disease is a chronic progressive disease characterized by abnormal sleep, sleep apnea and motor disorders, which is prone to misdiagnosis and missed diagnosis.

METHODS

We report a unique case of anti-IgLON5 disease in a pediatric patient with Langerhans cell histiocytosis (LCH). He gradually developed increased muscle tone and nystagmus during chemotherapy and showed signs of meningeal enhancement on cranial imaging. Due to insufficient evidence of LCH invasion of the central nervous system, the presence of autoimmune encephalitis-related antibodies was investigated by using cell-based assay (CBA) experiment in indirect immuno-fluorescence assay (IFA).

RESULTS

Clinical manifestations of sleep disorders and motor disorders, plus the presence of IgLON5 IgG antibodies (1:30) in the serum leading to a confirmed diagnosis of anti-IgLON5 disease.

CONCLUSION

Anti-IgLON5 disease is rare and almost no cases of children have been reported. In view of the difficult to recognize symptoms in pediatric patients, especially those with other comorbidities. Clinicians should raise their awareness of this disease and pay attention to the detection of autoimmune antibodies.