Andre Luiz Martins Reis

Living systems operate far from thermodynamic equilibrium, sustaining complexity through continuous free-energy influx and heat dissipation. Yet the higher-order functional constraints imposed by non-equilibrium dynamics remain largely unexplored. To formalize how changes in cellular entropy production constrain intracellular organization across proliferative state transitions, we developed a stochastic-thermodynamic model linking reduced dissipative drive to transient increases in intracellular configurational entropy. Within this framework, we examined abundance-weighted predicted RNA folding stability as a transcriptome-level projection of the intracellular configurational remodelling. Consistent with the model, transcriptome reorganization across the cell cycle produces coordinated, phase-dependent changes in the distribution of transcript-level RNA folding free energies. Extending this analysis to stress adaptation, we find that human cancer cells exposed to molecularly targeted agents and nutrient-limited Escherichia coli converge on expressed RNA pools enriched for transcripts with less negative RNA folding free energies, despite marked heterogeneity in the transcriptional responses within each system. Opposing shifts during acquired resistance under continued drug exposure and bacterial recovery from starvation further support an association with growth state. Together, these findings suggest that transcriptome organization extends beyond gene-regulatory programs to collective biophysical properties that may reflect non-equilibrium constraints shaping cellular-state transitions.

Structural variants (SVs) are a diverse group of genetic variants defined by a minimum size of 50 base pairs. SVs account for the majority of all variant bases in a person's genome and are commonly implicated in inherited disease and cancer. However, SV analysis is complex due to their wide variation in type and size, degree of polymorphism, involvement of repetitive sequences, and the myriad ways they may elicit a functional impact, as well as technical factors like imprecise breakpoint detection, and alternative representations of the same event. Despite recent advances in the detection and characterisation of SVs, it remains difficult to assess them beyond basic annotations and comparisons. Here we introduce SVlog, a transparent and extensible meta-programming framework for SV analysis. With the logic programming language Soufflé as its engine, SVlog provides a declarative ontology describing relationships among SVs, genes and other genomic elements. Genome annotations and SV datasets — both user-provided and public reference data — are converted into relational 'facts', to which SVlog applies logical rules that define 'predicates'. Predicates are specific, transparent and deterministic, yet fully flexible and composable, enabling detailed evaluation of SVs without relying on stochastic "black box" approaches. To showcase SVlog, we have developed a ready-made predicate library for SV annotation, comparison and prioritisation in the context of rare inherited disease. Despite its compact codebase, SVlog evaluates more than 50 input predicates to generate over 70 informative output predicates. It synthesises evidence from population and clinical genomic databases, and applies a tiered filtering strategy to identify candidate pathogenic SVs in patients with inherited disease. By focusing on explainability and modularity, SVlog offers a fast, reliable library for SV analysis and is a powerful deterministic alternative to traditional bioinformatics pipelines for clinical variant curation.

The genetic variants that cause inherited myopathies vary widely in type, size and sequence context, encompassing small sequence variants, large structural variants, repeat expansions, and more complex events, such as the D4Z4 macrosatellite contraction and hypomethylation that causes facioscapulohumeral muscular dystrophy. Many of these are challenging to characterise using next-generation sequencing and other older molecular technologies. To address this, we developed a targeted long-read sequencing assay and bioinformatics analysis framework that captures the full suite of genes, variants and epigenetic signatures currently implicated in inherited myopathies. Applying this to a cohort of myopathy patients, we demonstrate the analytical validity of our approach and its improved accuracy and resolution compared to existing methods. Our assay led to new genetic diagnoses in 35.5% (11/31) of patients who remained undiagnosed after standard clinical genetic testing. This methodology constitutes a single streamlined assay for comprehensive genetic and epigenetic characterisation of inherited myopathies.

DNA methylation (5mC) is an epigenetic mark that plays a critical role in defining cell fate. Following fertilisation, DNA methylation inherited from gametes must be reprogrammed to establish totipotency and enable the parental-to-zygotic transition. To accomplish this, non-mammalian vertebrates such as zebrafish and medaka subtly reprogramme maternal 5mC profiles while maintaining high methylation levels throughout embryogenesis. In contrast, eutherian mammals such as mouse and human undergo global 5mC erasure in both embryonic and extraembryonic lineages. However, while embryonic 5mC is rapidly re-established to high levels upon implantation, the trophectoderm, which gives rise to the placenta, displays sustained and conserved DNA hypomethylation, suggesting that this drastic 5mC erasure may be functionally linked to complex placentation in mammals. To clarify whether extensive post-fertilisation 5mC erasure co-evolved with placentation, we explored embryonic methylation dynamics in marsupials, a lineage of therian mammals with a short-lived placenta. We produced a near complete telomere-to-telomere (T2T) genome and generated detailed epigenome maps of embryonic development for an Australian marsupial, the fat-tailed dunnart (Sminthopsis crassicaudata). We found the dunnart embryo exhibits genome wide DNA demethylation at the blastocyst stage, but these changes occur in the trophectoderm only, suggesting that 5mC erasure in the placenta is an ancestral state in therian mammals. Furthermore, the T2T-level dunnart genome assembly enabled identification of sex chromosomes, uncovering extensive hypomethylation of the paternally-inherited inactive X chromosome in females and revealing the previously unannotated master regulator of X chromosome inactivation, lncRNA RSX. Our data indicate that while the use of genome-wide 5mC erasure differs between eutherian and marsupial lineages, 5mC erasure in extraembryonic tissue is ancestral to therian mammals and may be necessary to support placental development.

DNA methylation (5mC) is an epigenetic mark that plays a critical role in defining cell fate. Following fertilisation, DNA methylation inherited from gametes must be reprogrammed to establish totipotency and enable the parental-to-zygotic transition. To accomplish this, non-mammalian vertebrates such as zebrafish and medaka subtly reprogram maternal 5mC profiles while maintaining high methylation levels throughout embryogenesis. In contrast, eutherian mammals such as mouse and human undergo global 5mC erasure in both embryonic and extraembryonic lineages. However, while embryonic 5mC is rapidly re-established to high levels upon implantation, the trophectoderm, which gives rise to the placenta, displays sustained and conserved DNA hypomethylation, suggesting that this drastic 5mC erasure may be functionally linked to complex placentation in mammals. To clarify whether extensive post-fertilisation 5mC erasure co-evolved with placentation, we explored embryonic methylation dynamics in marsupials, a lineage of therian mammals with a short-lived placenta. We produced a near complete telomere-to-telomere (T2T) genome and generated detailed epigenome maps of embryonic development for an Australian marsupial, the fat-tailed dunnart (Sminthopsis crassicaudata). We found the dunnart embryo exhibits genome wide DNA demethylation at the blastocyst stage, but these changes occur in the trophectoderm only, suggesting that 5mC erasure in the placenta is an ancestral state in therian mammals. Furthermore, the T2T-level dunnart genome assembly enabled identification of sex chromosomes, uncovering extensive hypomethylation of the paternally-inherited inactive X chromosome in females and revealing the previously unannotated master regulator of X chromosome inactivation, lncRNA Rsx. Our data indicate that while the use of genome-wide 5mC erasure differs between eutherian and marsupial lineages, 5mC erasure in extraembryonic tissue is ancestral to therian mammals and may be necessary to support placental development.

Inherited myopathies are a group of disorders with diverse and complex genetic aetiologies. The causative genetic variants vary widely in type, size and sequence context, encompassing small sequence variants and large structural variants in protein-coding genes, mitochondrial variants, repeat expansions, and more complex events, such as the chromosome 4 D4Z4 macrosatellite contraction and hypomethylation that causes facioscapulohumeral muscular dystrophy (FSHD). This poses a challenge for analysis with next generation sequencing and other molecular methods. This is further compounded by phenotypic variability between patients and phenotypic overlap of different forms of genetic myopathy. To accelerate myopathy research and improve diagnosis we have developed a targeted long-read sequencing assay and an integrated bioinformatics analysis framework that captures the full suite of genes, variants and epigenetic signatures currently implicated in all forms of inherited myopathy. Applying this to a cohort of 53 myopathy patients with and without previous genetic diagnoses, we demonstrate the analytical validity of our approach, as well as its improved accuracy and resolution compared to existing methods - especially for FSHD. Our LRS assay identified an array of new information about the genetic and epigenetic landscape of inherited myopathies and provided new diagnoses for 29% of patients who had remained undiagnosed following clinical genetic testing. In our cohort, FSHD and oculopharyngodistal myopathies were the most common new diagnoses that were missed or mis-diagnosed by standard clinical genetic testing. Our new method constitutes a single streamlined assay for comprehensive genetic and epigenetic characterisation of inherited myopathies.

Advances in long-read sequencing (LRS) and assembly algorithms have made it possible to create highly complete genome assemblies for humans, animals and plants. However, ongoing development is needed to improve accessibility, affordability, and assembly quality and completeness. ‘Cornetto’ is a new strategy in which we use programmable selective nanopore sequencing to focus LRS data production onto the unsolved regions of a nascent assembly. This improves assembly quality and streamlines the process, both for humans and non-human vertebrates. Cornetto enables us to generate highly complete diploid human genome assemblies using only nanopore LRS data, surpassing the quality of previous efforts at a fraction of the cost. Cornetto enables genome assembly from challenging sample types like human saliva. Finally, we obtain accurate assemblies for clinically-relevant repetitive loci at the extremes of the genome, demonstrating valid approaches for genetic diagnosis in facioscapulohumeral muscular dystrophy (FSHD) and MUC1-autosomal dominant tubulointerstitial kidney disease (MUC1-ADTKD).

Background: The central bearded dragon ( Pogona vitticeps) is widely distributed in central eastern Australia and adapts readily to captivity. Among other attributes, it is distinctive because it undergoes sex reversal from ZZ genotypic males to phenotypic females at high incubation temperatures. Here, we report an annotated telomere to telomere phased assembly of the genome of a female ZW central bearded dragon. Results: Genome assembly length is 1.75 Gbp with a scaffold N50 of 266.2 Mbp, N90 of 28.1 Mbp, 26 gaps and 42.2% GC content. Most (99.6%) of the reference assembly is scaffolded into 6 macrochromosomes and 10 microchromosomes, including the Z and W microchromosomes, corresponding to the karyotype. The genome assembly exceeds standard recommended by the Earth Biogenome Project (6CQ40): 0.003% collapsed sequence, 0.03% false expansions, 99.8% k-mer completeness, 97.9% complete single copy BUSCO genes and an average of 93.5% of transcriptome data mappable back to the genome assembly. The mitochondrial genome (16,731 bp) and the model rDNA repeat unit (length 9.5 Kbp) were assembled. Male vertebrate sex genes Amh and Amhr2 were discovered as copies in the small non-recombining region of the Z chromosome, absent from the W chromosome. This, coupled with the prior discovery of differential Z and W transcriptional isoform composition arising from pseudoautosomal sex gene Nr5a1 , suggests that complex interactions between these genes, their autosomal copies and their resultant transcription factors and intermediaries, determines sex in the bearded dragon. Conclusion: This high-quality assembly will serve as a resource to enable and accelerate research into the unusual reproductive attributes of this species and for comparative studies across the Agamidae and reptiles more generally. Species Taxonomy: Eukaryota; Animalia; Chordata; Reptilia; Squamata; Iguania; Agamidae; Amphibolurinae; Pogona ; Pogona vitticeps (Ahl, 1926) (NCBI:txid103695). Graphical

Background: The central bearded dragon (Pogona vitticeps) is widely distributed in central eastern Australia and adapts readily to captivity. Among other attributes, it is distinctive because it undergoes sex reversal from ZZ genotypic males to phenotypic females at high incubation temperatures. Here, we report an annotated near telomere-to-telomere phased assembly of the genome of a female ZW central bearded dragon. Results: Genome assembly length is 1.75 Gbp with a scaffold N50 of 266.2 Mbp, N90 of 28.1 Mbp, 26 gaps, and 42.2% GC content. Most (99.6%) of the reference assembly is scaffolded into 6 macrochromosomes and 10 microchromosomes, including the Z and W microchromosomes, corresponding to the karyotype. The genome assembly exceeds standard recommended by the Earth Biogenome Project (6CQ40): 0.003% collapsed sequence, 0.03% false expansions, 99.8% k-mer completeness, 97.9% complete single-copy BUSCO genes, and an average of 93.5% of transcriptome data mappable back to the genome assembly. The mitochondrial genome (16,731 bp) and the model ribosomal DNA repeat unit (length 9.5 Kbp) were assembled. Male vertebrate sex genes Amh and Amhr2 were discovered as copies in the small non-recombining region of the Z chromosome, absent from the W chromosome. This, coupled with the prior discovery of differential Z and W transcriptional isoform composition arising from pseudo-autosomal sex gene Nr5a1, suggests that complex interactions between these genes, their autosomal copies, and their resultant transcription factors and intermediaries determine sex in the bearded dragon. Conclusion: This high-quality assembly will serve as a resource to enable and accelerate research into the unusual reproductive attributes of this species and for comparative studies across the Agamidae and reptiles more generally.

Background: Loss‐of‐function variants in MME (membrane metalloendopeptidase) are a known cause of recessive Charcot–Marie–Tooth Neuropathy (CMT). A deep intronic variant, MME c.1188+428A>G (NM_000902.5), was identified through whole genome sequencing (WGS) of two Australian families with recessive inheritance of axonal CMT using the seqr platform. MME c.1188+428A>G was detected in a homozygous state in Family 1, and in a compound heterozygous state with a known pathogenic MME variant (c.467del; p.Pro156Leufs*14) in Family 2.Aims: We aimed to determine the pathogenicity of the MME c.1188+428A>G variant through segregation and splicing analysis.Methods: The splicing impact of the deep intronic MME variant c.1188+428A>G was assessed using an in vitro exon‐trapping assay.Results: The exon‐trapping assay demonstrated that the MME c.1188+428A>G variant created a novel splice donor site resulting in the inclusion of an 83 bp pseudoexon between MME exons 12 and 13. The incorporation of the pseudoexon into MME transcript is predicted to lead to a coding frameshift and premature termination codon (PTC) in MME exon 14 (p.Ala397ProfsTer47). This PTC is likely to result in nonsense mediated decay (NMD) of MME transcript leading to a pathogenic loss‐of‐function.Interpretation: To our knowledge, this is the first report of a pathogenic deep intronic MME variant causing CMT. This is of significance as deep intronic variants are missed using whole exome sequencing screening methods. Individuals with CMT should be reassessed for deep intronic variants, with splicing impacts being considered in relation to the potential pathogenicity of variants.