Incomplete transcriptional dosage compensation of chicken and platypus sex chromosomes is balanced by post-transcriptional compensation.

Heteromorphic sex chromosomes (XY or ZW) present problems of gene dosage imbalance between sexes and with autosomes. A need for dosage compensation has long been thought to be critical in vertebrates. However, this was questioned by findings of unequal mRNA abundance measurements in monotreme mammals and birds. Here, we demonstrate unbalanced mRNA levels of X genes in platypus males and females and a correlation with differential loading of histone modifications. We also observed unbalanced transcripts of Z genes in chicken. Surprisingly, however, we found that protein abundance ratios were 1:1 between the sexes in both species, indicating a post-transcriptional layer of dosage compensation. We conclude that sex chromosome output is maintained in chicken and platypus (and perhaps many other non therian vertebrates) via a combination of transcriptional and post-transcriptional control, consistent with a critical importance of sex chromosome dosage compensation.

Origin and maintenance of large ribosomal RNA gene repeat size in mammals.

The genes encoding ribosomal RNA are highly conserved across life and in almost all eukaryotes are present in large tandem repeat arrays called the rDNA. rDNA repeat unit size is conserved across most eukaryotes but has expanded dramatically in mammals, principally through the expansion of the intergenic spacer region that separates adjacent rRNA coding regions. Here, we used long-read sequence data from representatives of the major amniote lineages to determine where in amniote evolution rDNA unit size increased. We find that amniote rDNA unit sizes fall into two narrow size classes: "normal" (∼11-20 kb) in all amniotes except monotreme, marsupial, and eutherian mammals, which have "large" (∼35-45 kb) sizes. We confirm that increases in intergenic spacer length explain much of this mammalian size increase. However, in stark contrast to the uniformity of mammalian rDNA unit size, mammalian intergenic spacers differ greatly in sequence. These results suggest a large increase in intergenic spacer size occurred in a mammalian ancestor and has been maintained despite substantial sequence changes over the course of mammalian evolution. This points to a previously unrecognized constraint on the length of the intergenic spacer, a region that was thought to be largely neutral. We finish by speculating on possible causes of this constraint.

Identification of the RSX interactome in a marsupial shows functional coherence with the Xist interactome during X inactivation.

The marsupial specific RSX lncRNA is the functional analogue of the eutherian specific XIST, which coordinates X chromosome inactivation. We characterized the RSX interactome in a marsupial representative (the opossum Monodelphis domestica), identifying 135 proteins, of which 54 had orthologues in the XIST interactome. Both interactomes were enriched for biological pathways related to RNA processing, regulation of translation, and epigenetic transcriptional silencing. This represents a remarkable example showcasing the functional coherence of independently evolved lncRNAs in distantly related mammalian lineages.

Cognitive Function Is Associated with the Genetically Determined Efficiency of DNA Repair Mechanisms.

Several modifiable risk factors for neurodegeneration and dementia have been identified, although individuals vary in their vulnerability despite a similar risk of exposure. This difference in vulnerability could be explained at least in part by the variability in DNA repair mechanisms' efficiency between individuals. Therefore, the aim of this study was to test associations between documented, prevalent genetic variation (single nucleotide polymorphism, SNP) in DNA repair genes, cognitive function, and brain structure. Community-living participants (n = 488,159; 56.54 years (8.09); 54.2% female) taking part in the UK Biobank study and for whom cognitive and genetic measures were available were included. SNPs in base excision repair (BER) genes of the bifunctional DNA glycosylases OGG1 (rs1052133, rs104893751), NEIL1 (rs7402844, rs5745906), NEIL2 (rs6601606), NEIL3 (rs10013040, rs13112390, rs13112358, rs1395479), MUTYH (rs34612342, rs200165598), NTHL1 (rs150766139, rs2516739) were considered. Cognitive measures included fluid intelligence, the symbol-digit matching task, visual matching, and trail-making. Hierarchical regression and latent class analyses were used to test the associations between SNPs and cognitive measures. Associations between SNPs and brain measures were also tested in a subset of 39,060 participants. Statistically significant associations with cognition were detected for 12 out of the 13 SNPs analyzed. The strongest effects amounted to a 1-6% difference in cognitive function detected for NEIL1 (rs7402844), NEIL2 (rs6601606), and NTHL1 (rs2516739). Associations varied by age and sex, with stronger effects detected in middle-aged women. Weaker associations with brain measures were also detected. Variability in some BER genes is associated with cognitive function and brain structure and may explain variability in the risk for neurodegeneration and dementia.

The landscape of genomic structural variation in Indigenous Australians.

Indigenous Australians harbour rich and unique genomic diversity. However, Aboriginal and Torres Strait Islander ancestries are historically under-represented in genomics research and almost completely missing from reference datasets1-3. Addressing this representation gap is critical, both to advance our understanding of global human genomic diversity and as a prerequisite for ensuring equitable outcomes in genomic medicine. Here we apply population-scale whole-genome long-read sequencing4 to profile genomic structural variation across four remote Indigenous communities. We uncover an abundance of large insertion-deletion variants (20-49 bp; n = 136,797), structural variants (50 b-50 kb; n = 159,912) and regions of variable copy number (>50 kb; n = 156). The majority of variants are composed of tandem repeat or interspersed mobile element sequences (up to 90%) and have not been previously annotated (up to 62%). A large fraction of structural variants appear to be exclusive to Indigenous Australians (12% lower-bound estimate) and most of these are found in only a single community, underscoring the need for broad and deep sampling to achieve a comprehensive catalogue of genomic structural variation across the Australian continent. Finally, we explore short tandem repeats throughout the genome to characterize allelic diversity at 50 known disease loci5, uncover hundreds of novel repeat expansion sites within protein-coding genes, and identify unique patterns of diversity and constraint among short tandem repeat sequences. Our study sheds new light on the dimensions and dynamics of genomic structural variation within and beyond Australia.

Indigenous Australian genomes show deep structure and rich novel variation.

The Indigenous peoples of Australia have a rich linguistic and cultural history. How this relates to genetic diversity remains largely unknown because of their limited engagement with genomic studies. Here we analyse the genomes of 159 individuals from four remote Indigenous communities, including people who speak a language (Tiwi) not from the most widespread family (Pama-Nyungan). This large collection of Indigenous Australian genomes was made possible by careful community engagement and consultation. We observe exceptionally strong population structure across Australia, driven by divergence times between communities of 26,000-35,000 years ago and long-term low but stable effective population sizes. This demographic history, including early divergence from Papua New Guinean (47,000 years ago) and Eurasian groups1, has generated the highest proportion of previously undescribed genetic variation seen outside Africa and the most extended homozygosity compared with global samples. A substantial proportion of this variation is not observed in global reference panels or clinical datasets, and variation with predicted functional consequence is more likely to be homozygous than in other populations, with consequent implications for medical genomics2. Our results show that Indigenous Australians are not a single homogeneous genetic group and their genetic relationship with the peoples of New Guinea is not uniform. These patterns imply that the full breadth of Indigenous Australian genetic diversity remains uncharacterized, potentially limiting genomic medicine and equitable healthcare for Indigenous Australians.

Uncovering the Genetic and Molecular Features of Huntington's Disease in Northern Colombia.

Huntington's disease (HD) is a genetic disorder caused by a CAG trinucleotide expansion in the huntingtin (HTT) gene. Juan de Acosta, Atlántico, a city located on the Caribbean coast of Colombia, is home to the world's second-largest HD pedigree. Here, we include 291 descendants of this pedigree with at least one family member with HD. Blood samples were collected, and genomic DNA was extracted. We quantified the HTT CAG expansion using an amplicon sequencing protocol. The genetic heterogeneity was measured as the ratio of the mosaicism allele's read peak and the slippage ratio of the allele's read peak from our sequence data. The statistical and bioinformatic analyses were performed with a significance threshold of p < 0.05. We found that the average HTT CAG repeat length in all participants was 21.91 (SD = 8.92). Of the 291 participants, 33 (11.3%, 18 females) had a positive molecular diagnosis for HD. Most affected individuals were adults, and the most common primary and secondary alleles were 17/7 (CAG/CCG) and 17/10 (CAG/CCG), respectively. The mosaicism increased with age in the participants with HD, while the slippage analyses revealed differences by the HD allele type only for the secondary allele. The slippage tended to increase with the HTT CAG repeat length in the participants with HD, but the increase was not statistically significant. This study analyzed the genetic and molecular features of 291 participants, including 33 with HD. We found that the mosaicism increased with age in the participants with HD, particularly for the secondary allele. The most common haplotype was 17/7_17/10. The slippage for the secondary allele varied by the HD allele type, but there was no significant difference in the slippage by sex. Our findings offer valuable insights into HD and could have implications for future research and clinical management.

Gene expression of male pathway genes sox9 and amh during early sex differentiation in a reptile departs from the classical amniote model.

BACKGROUND: Sex determination is the process whereby the bipotential embryonic gonads become committed to differentiate into testes or ovaries. In genetic sex determination (GSD), the sex determining trigger is encoded by a gene on the sex chromosomes, which activates a network of downstream genes; in mammals these include SOX9, AMH and DMRT1 in the male pathway, and FOXL2 in the female pathway. Although mammalian and avian GSD systems have been well studied, few data are available for reptilian GSD systems. RESULTS: We conducted an unbiased transcriptome-wide analysis of gonad development throughout differentiation in central bearded dragon (Pogona vitticeps) embryos with GSD. We found that sex differentiation of transcriptomic profiles occurs at a very early stage, before the gonad consolidates as a body distinct from the gonad-kidney complex. The male pathway genes dmrt1 and amh and the female pathway gene foxl2 play a key role in early sex differentiation in P. vitticeps, but the central player of the mammalian male trajectory, sox9, is not differentially expressed in P. vitticeps at the bipotential stage. The most striking difference from GSD systems of other amniotes is the high expression of the male pathway genes amh and sox9 in female gonads during development. We propose that a default male trajectory progresses if not repressed by a W-linked dominant gene that tips the balance of gene expression towards the female trajectory. Further, weighted gene expression correlation network analysis revealed novel candidates for male and female sex differentiation. CONCLUSION: Our data reveal that interpretation of putative mechanisms of GSD in reptiles cannot solely depend on lessons drawn from mammals.

Arctic introgression and chromatin regulation facilitated rapid Qinghai-Tibet Plateau colonization by an avian predator.

The Qinghai-Tibet Plateau (QTP), possesses a climate as cold as that of the Arctic, and also presents uniquely low oxygen concentrations and intense ultraviolet (UV) radiation. QTP animals have adapted to these extreme conditions, but whether they obtained genetic variations from the Arctic during cold adaptation, and how genomic mutations in non-coding regions regulate gene expression under hypoxia and intense UV environment, remain largely unknown. Here, we assemble a high-quality saker falcon genome and resequence populations across Eurasia. We identify female-biased hybridization with Arctic gyrfalcons in the last glacial maximum, that endowed eastern sakers with alleles conveying larger body size and changes in fat metabolism, predisposing their QTP cold adaptation. We discover that QTP hypoxia and UV adaptations mainly involve independent changes in non-coding genomic variants. Our study highlights key roles of gene flow from Arctic relatives during QTP hypothermia adaptation, and cis-regulatory elements during hypoxic response and UV protection.

Kmer2SNP: Reference-Free Heterozygous SNP Calling Using k-mer Frequency Distributions.

DNA sequencing technologies enable the generation of genetic profiles from many individuals at a rapid rate. Identifying single-nucleotide polymorphism (SNP) between biological samples is fundamental in genetics with various applications, such as disease diagnosis and associations and ancestry and relationship inference. Most methods use a species-specific reference genome for aligning raw sequenced reads for accurate SNP calling. However, high-quality reference genomes may not be available for all species. Therefore, we developed a reference-free algorithm, Kmer2SNP, to identify heterozygous SNPs from raw sequenced reads to facilitate genetic studies in species without the reference genome. Kmer2SNP first calculates the k-mer frequency distribution from reads to determine k-mers containing heterozygous SNPs. Next, these k-mers are rapidly matched with each other to identify pairs of exact heterozygous k-mers that belong to one of the two possible haplotypes in a diploid organism. Finally, using overlapping neighboring k-mers, weights are assigned for SNP assignments; higher weights increase SNP discovery confidence.

The genomic landscape of blood groups in Indigenous Australians in remote communities.

BACKGROUND: Red blood cell (RBC) membrane-associated blood group systems are clinically significant. Alloimmunisation is a persistent risk associated with blood transfusion owing to the antigen polymorphisms among these RBC-associated blood groups. Next-generation sequencing (NGS) offers an opportunity to characterize the blood group variant profile of a given individual. Australia comprises a large multiethnic population where most blood donors are Caucasian and blood group variants remain poorly studied among Indigenous Australians. In this study, we focused on the Tiwi Islanders, who have lived in relative isolation for thousands of years. METHODS AND MATERIALS: We predicted the blood group phenotype profiles in the Tiwi (457) and 1000 Genomes Phase 3 (1KGP3-2504) cohort individuals using RBCeq (https://www.rbceq.org/). The predicted phenotype prevalence was compared with the previous literature report. RESULTS: We report, for the first time, comprehensive blood group profiles corresponding to the 35 known blood group systems among the Indigenous Tiwi islander population and identify possible novel antigen variants therein. Our results demonstrate that the genetic makeup of the Tiwi participants is distinct from that of other populations, with a low prevalence of LU (Au[a-b+]) and ABO (A2) and D+C+c+E+e- phenotype, an absence of Diego blood group variants, and a unique RHD (DIII type4) variant. CONCLUSION: Our results may contribute to the development of a database of predicted phenotype donors among the Tiwi population and aid in improving transfusion safety for the ~2800 Tiwi people and the ~800,000 other Indigenous Australians throughout the nation.

Sex-specific splicing of Z- and W-borne nr5a1 alleles suggests sex determination is controlled by chromosome conformation.

Pogona vitticeps has female heterogamety (ZZ/ZW), but the master sex-determining gene is unknown, as it is for all reptiles. We show that nr5a1 (Nuclear Receptor Subfamily 5 Group A Member 1), a gene that is essential in mammalian sex determination, has alleles on the Z and W chromosomes (Z-nr5a1 and W-nr5a1), which are both expressed and can recombine. Three transcript isoforms of Z-nr5a1 were detected in gonads of adult ZZ males, two of which encode a functional protein. However, ZW females produced 16 isoforms, most of which contained premature stop codons. The array of transcripts produced by the W-borne allele (W-nr5a1) is likely to produce truncated polypeptides that contain a structurally normal DNA-binding domain and could act as a competitive inhibitor to the full-length intact protein. We hypothesize that an altered configuration of the W chromosome affects the conformation of the primary transcript generating inhibitory W-borne isoforms that suppress testis determination. Under this hypothesis, the genetic sex determination (GSD) system of P. vitticeps is a W-borne dominant female-determining gene that may be controlled epigenetically.

Balancing openness with Indigenous data sovereignty: An opportunity to leave no one behind in the journey to sequence all of life.

The field of genomics has benefited greatly from its "openness" approach to data sharing. However, with the increasing volume of sequence information being created and stored and the growing number of international genomics efforts, the equity of openness is under question. The United Nations Convention of Biodiversity aims to develop and adopt a standard policy on access and benefit-sharing for sequence information across signatory parties. This standardization will have profound implications on genomics research, requiring a new definition of open data sharing. The redefinition of openness is not unwarranted, as its limitations have unintentionally introduced barriers of engagement to some, including Indigenous Peoples. This commentary provides an insight into the key challenges of openness faced by the researchers who aspire to protect and conserve global biodiversity, including Indigenous flora and fauna, and presents immediate, practical solutions that, if implemented, will equip the genomics community with both the diversity and inclusivity required to respectfully protect global biodiversity.

Pharmacogenomic analysis of a genetically distinct Indigenous population.

Indigenous Australians face a disproportionately severe burden of chronic disease relative to other Australians, with elevated rates of morbidity and mortality. While genomics technologies are slowly gaining momentum in personalised treatments for many, a lack of pharmacogenomic research in Indigenous peoples could delay adoption. Appropriately implementing pharmacogenomics in clinical care necessitates an understanding of the frequencies of pharmacologically relevant genetic variants within Indigenous populations. We analysed whole-genome sequence data from 187 individuals from the Tiwi Islands and characterised the pharmacogenomic landscape of this population. Specifically, we compared variant profiles and allelic distributions of previously described pharmacologically significant genes and variants with other population groups. We identified 22 translationally relevant pharmacogenomic variants and 18 clinically actionable guidelines with implications for drug dosing and treatment of conditions including heart disease, diabetes and cancer. We specifically observed increased poor and intermediate metabolizer phenotypes in the CYP2C9 (PM:19%, IM:44%) and CYP2C19 (PM:18%, IM:44%) genes.

Genomic attributes differ between Vibrio parahaemolyticus environmental and clinical isolates including pathotypes.

Vibrio parahaemolyticus is a marine bacterium and causes opportunistic gastroenteritis in humans. Clinical strains of V. parahaemolyticus contain haemolysin and type III secretion systems (T3SS) that define their pathotype. A growing number of strains isolated recently from the environment have acquired these virulence genes constituting a pool of potential pathogens. This study used comparative genomics to identify genetic factors that delineate environmental and clinical V. parahaemolyticus population and understand the similarities and differences between the T3SS2 phylotypes. The comparative analysis revealed the presence of a cluster of genes belonging to bacterial cellulose synthesis (bcs) in isolates of environmental origin. This cluster, previously unreported in V. parahaemolyticus, exhibit significant similarity to that of Aliivibrio fischeri, and might dictate a potentially new mechanism of its environmental adaptation and persistence. The study also identified many genes predicted in silico to be T3SS effectors that are unique to T3SS2ß of tdh- trh+ and tdh+ trh+ pathotype and having no identifiable homologue in tdh+ trh- T3SS2α. Overall, these findings highlight the importance of understanding the genes and strategies V. parahaemolyticus utilize for the myriad interactions with its hosts, either marine invertebrates or humans.

Microchromosomes are building blocks of bird, reptile, and mammal chromosomes.

Microchromosomes, once considered unimportant shreds of the chicken genome, are gene-rich elements with a high GC content and few transposable elements. Their origin has been debated for decades. We used cytological and whole-genome sequence comparisons, and chromosome conformation capture, to trace their origin and fate in genomes of reptiles, birds, and mammals. We find that microchromosomes as well as macrochromosomes are highly conserved across birds and share synteny with single small chromosomes of the chordate amphioxus, attesting to their origin as elements of an ancient animal genome. Turtles and squamates (snakes and lizards) share different subsets of ancestral microchromosomes, having independently lost microchromosomes by fusion with other microchromosomes or macrochromosomes. Patterns of fusions were quite different in different lineages. Cytological observations show that microchromosomes in all lineages are spatially separated into a central compartment at interphase and during mitosis and meiosis. This reflects higher interaction between microchromosomes than with macrochromosomes, as observed by chromosome conformation capture, and suggests some functional coherence. In highly rearranged genomes fused microchromosomes retain most ancestral characteristics, but these may erode over evolutionary time; surprisingly, de novo microchromosomes have rapidly adopted high interaction. Some chromosomes of early-branching monotreme mammals align to several bird microchromosomes, suggesting multiple microchromosome fusions in a mammalian ancestor. Subsequently, multiple rearrangements fueled the extraordinary karyotypic diversity of therian mammals. Thus, microchromosomes, far from being aberrant genetic elements, represent fundamental building blocks of amniote chromosomes, and it is mammals, rather than reptiles and birds, that are atypical.

A globally diverse reference alignment and panel for imputation of mitochondrial DNA variants.

BACKGROUND: Variation in mitochondrial DNA (mtDNA) identified by genotyping microarrays or by sequencing only the hypervariable regions of the genome may be insufficient to reliably assign mitochondrial genomes to phylogenetic lineages or haplogroups. This lack of resolution can limit functional and clinical interpretation of a substantial body of existing mtDNA data. To address this limitation, we developed and evaluated a large, curated reference alignment of complete mtDNA sequences as part of a pipeline for imputing missing mtDNA single nucleotide variants (mtSNVs). We call our reference alignment and pipeline MitoImpute. RESULTS: We aligned the sequences of 36,960 complete human mitochondrial genomes downloaded from GenBank, filtered and controlled for quality. These sequences were reformatted for use in imputation software, IMPUTE2. We assessed the imputation accuracy of MitoImpute by measuring haplogroup and genotype concordance in data from the 1000 Genomes Project and the Alzheimer's Disease Neuroimaging Initiative (ADNI). The mean improvement of haplogroup assignment in the 1000 Genomes samples was 42.7% (Matthew's correlation coefficient = 0.64). In the ADNI cohort, we imputed missing single nucleotide variants. CONCLUSION: These results show that our reference alignment and panel can be used to impute missing mtSNVs in existing data obtained from using microarrays, thereby broadening the scope of functional and clinical investigation of mtDNA. This improvement may be particularly useful in studies where participants have been recruited over time and mtDNA data obtained using different methods, enabling better integration of early data collected using less accurate methods with more recent sequence data.

Functional microRNA targetome undergoes degeneration-induced shift in the retina.

BACKGROUND: MicroRNA (miRNA) play a significant role in the pathogenesis of complex neurodegenerative diseases including age-related macular degeneration (AMD), acting as post-transcriptional gene suppressors through their association with argonaute 2 (AGO2) - a key member of the RNA Induced Silencing Complex (RISC). Identifying the retinal miRNA/mRNA interactions in health and disease will provide important insight into the key pathways miRNA regulate in disease pathogenesis and may lead to potential therapeutic targets to mediate retinal degeneration. METHODS: To identify the active miRnome targetome interactions in the healthy and degenerating retina, AGO2 HITS-CLIP was performed using a rodent model of photoreceptor degeneration. Analysis of publicly available single-cell RNA sequencing (scRNAseq) data was performed to identify the cellular location of AGO2 and key members of the microRNA targetome in the retina. AGO2 findings were verified by in situ hybridization (RNA) and immunohistochemistry (protein). RESULTS: Analysis revealed a similar miRnome between healthy and damaged retinas, however, a shift in the active targetome was observed with an enrichment of miRNA involvement in inflammatory pathways. This shift was further demonstrated by a change in the seed binding regions of miR-124-3p, the most abundant retinal AGO2-bound miRNA, and has known roles in regulating retinal inflammation. Additionally, photoreceptor cluster miR-183/96/182 were all among the most highly abundant miRNA bound to AGO2. Following damage, AGO2 expression was localized to the inner retinal layers and more in the OLM than in healthy retinas, indicating a locational miRNA response to retinal damage. CONCLUSIONS: This study provides important insight into the alteration of miRNA regulatory activity that occurs as a response to retinal degeneration and explores the miRNA-mRNA targetome as a consequence of retinal degenerations. Further characterisation of these miRNA/mRNA interactions in the context of the degenerating retina may provide an important insight into the active role these miRNA may play in diseases such as AMD.

PACIFIC: a lightweight deep-learning classifier of SARS-CoV-2 and co-infecting RNA viruses.

Viral co-infections occur in COVID-19 patients, potentially impacting disease progression and severity. However, there is currently no dedicated method to identify viral co-infections in patient RNA-seq data. We developed PACIFIC, a deep-learning algorithm that accurately detects SARS-CoV-2 and other common RNA respiratory viruses from RNA-seq data. Using in silico data, PACIFIC recovers the presence and relative concentrations of viruses with > 99% precision and recall. PACIFIC accurately detects SARS-CoV-2 and other viral infections in 63 independent in vitro cell culture and patient datasets. PACIFIC is an end-to-end tool that enables the systematic monitoring of viral infections in the current global pandemic.