Search

Hybridization between coastal cutthroat trout (Oncorhynchus clarkii clarkii) and steelhead (O. mykiss) was assessed in the Smith River, California. Individuals were categorized as pure or as 1 of 10 hybrid classes using 30 “diagnostic” single-nucleotide polymorphisms positioned on 26 separate chromosomes. Most of the individuals examined (n = 876), were pure coastal cutthroat trout (n = 634) or pure steelhead (n = 213), and 29 individuals were identified as having hybrid ancestry. Among hybrids, first generation hybrids (n = 15) and coastal cutthroat trout backcrosses (n = 12) were the most common. No individuals were identified as backcrosses to SH, suggesting the presence of genetic or behavioral mechanisms constraining such backcrosses, or the growth and survival of their progeny. Mitochondrial DNA of 14 of 15 F1 hybrids was of steelhead origin, suggesting that hybridization was driven primarily by sneak-mating of male coastal cutthroat trout with female steelhead. Evaluation of classical phenotypic characters for coastal cutthroat trout and steelhead (i.e. jaw slash, maxillary length, and hyoid teeth) were not reliable by themselves for identification of either pure parental fish or hybrids. In contrast, analysis with geometric morphometrics revealed distinctive body shapes for pure coastal cutthroat trout and steelhead, and the combination of classical traits and geometric morphology was mostly accurate in distinguishing them. However, first generation hybrids and backcrosses overlapped completely with parental types, highlighting challenges in hybrid identification using phenotypic traits.

Pancreatic cancer still has one of the worst prognoses of all solid malignancies, despite developments in cancer knowledge and care. Research into pancreatic cancer has not fully translated into clinical improvements and as a result, fewer than 1% of patients survive 10 years post-diagnosis. This bleak outlook for patients could be improved by earlier diagnosis. The human erythrocyte phosphatidylinositol glycan class A (PIG-A) assay monitors the mutation status of the X-linked PIG-A gene by measuring glycosyl phosphatidylinositol (GPI)-anchored proteins on the extracellular surface. We have previously identified an elevated PIG-A mutant frequency in oesophageal adenocarcinoma patients and here investigate whether this could be seen in a pancreatic cancer cohort, given the urgent need for novel pancreatic cancer biomarkers. In our pilot study, an elevated PIG-A mutant frequency (5.775 × 10−6 (95% CI 4.777–10) mutants per million) was seen in pancreatic cancer patients (n = 30) when compared to the non-cancer control group (n = 14) who had an erythrocyte mutant frequency of 4.211 × 10−6 (95% CI 1.39–5.16) mutants per million (p = 0.0052). A cut-off value of 4.7 mutants per million provided an AUROC of 0.7595 with a sensitivity of 70% and specificity of 78.57%. A secondary measure of DNA damage in an alternative blood cell population also showed an increase in peripheral lymphocytes using the cytokinesis-block micronucleus assay (p = 0.0164) (AUROC = 0.77, sensitivity = 72.22%, specificity = 72.73%). The micronucleus frequency and PIG-A status show some potential as blood-based biomarkers of pancreatic cancer, but further investigations of these DNA damage tests are required to assess their utility in pancreatic cancer diagnosis.

Hibiscus liliiflorus, endemic to the Indian Ocean island of Rodrigues, is one of the rarest plant species in the world; only 2 wild individuals remain. Previously, when 4 wild individuals remained, the Mauritian Wildlife Foundation (MWF) in Rodrigues propagated cuttings of them in their nursery, then planted seedlings produced in the nursery into 3 outplanted populations on the island. Our goals were to: 1) assess whether all 4 original wild genotypes are represented in the MWF nursery; 2) determine whether ex situ living collections at international botanical gardens maintain unique genotypes of H. liliiflorus; 3) assess whether nursery individuals have crossed or self-fertilized to produce seed and quantify their relative contributions to outplanted populations; and 4) provide recommendations for future conservation actions. We used a 2b-RADseq approach to produce 2,711 genome-wide single nucleotide polymorphisms (SNPs) from 98 samples. Genotype identity analysis, principal component analysis, and model-based clustering in STRUCTURE found 4 genotypes extant in Rodrigues but no unique genotypes in ex situ botanic garden collections. Only 3 genotypes are represented in the MWF nursery; the one remaining genotype is represented by an extant wild individual. Parentage analysis showed that seeds produced in the MWF nursery resulted from both self-fertilization and crossing between genotypes, a result supported by internal relatedness and hybrid index calculations. Each outplanted population is dominated by a subset of parental genotypes, and we propose actions to balance the parental contributions to outplanted populations. Our study highlights how genetic assessments of ex situ conservation projects help conserve critically endangered species.

Orofacial cleft (OC) is a common congenital anomaly in humans, which has lifelong implications for affected individuals. This disorder can be classified as syndromic or non-syndromic depending on the presence or absence of additional physical or neurodevelopmental abnormalities, respectively. Non-syndromic cleft is often non-familial in nature and has a complex aetiology, whereas syndromic forms tend to be monogenic. Although individual OC-related syndromes have been frequently described in the medical literature, there has not been a comprehensive review across syndromes, thereby leaving a gap in our knowledge, which this paper aims to address. Six hundred and three patients with cleft-related human phenotype ontology terms were identified within the Deciphering Developmental Disorders study. Genes carrying pathogenic/likely pathogenic variants were identified and reviewed enabling a diagnostic yield of 36.5%. In total, 124 candidate genes for syndromic OC were identified, including 34 new genes that should be considered for inclusion in clinical clefting panels. Functional enrichment and gene expression analyses identified three key processes that were significantly overrepresented in syndromic OC gene lists: embryonic morphogenesis, protein stability and chromatin organization. Comparison with non-syndromic OC gene networks led us to propose that chromatin remodelling specifically contributes to the aetiology of syndromic OC. Disease-driven gene discovery is a valid approach to gene identification and curation of gene panels. Through this approach, we have started to unravel common molecular pathways contributing to syndromic orofacial clefting.

Telomere length (TL), which is maintained by human telomerase reverse transcriptase (hTERT; component of telomerase) and/or TRF1/TRF2 (core components of shelterin) via different mechanisms, is essential for chromosomal stability and cell survival. Folates comprise a group of essential B9 vitamin that involve in DNA synthesis and methylation. This study aimed to evaluate the effects of folic acid (FA) and 5-methyltetrahydrofolate (5-MeTHF) on TL, chromosome stability, and cell survival of telomerase-negative BJ and telomerase-positive A375 cells in vitro. BJ and A375 cells were cultured in modified medium with FA or 5-MeTHF (22.6 or 2260 nM) for 28 days. TL and mRNA expression were determined by RT-qPCR. Chromosome instability (CIN) and cell death were measured by CBMN-Cyt assay. Results showed that abnormal TL elongation was observed in FA- and 5-MeTHF-deficient BJ cells. The TL of A375 cells showed no obvious alterations under the FA-deficient condition but was significantly elongated under the 5-MeTHF-deficient condition. In both BJ and A375 cells, FA and 5-MeTHF deficiency caused decreased TRF1, TRF2, and hTERT expression, increased CIN and cell death; while a high concentration of 5-MeTHF induced elongated TL, elevated CIN, increased TRF1 and TRF2 expression, and decreased hTERT expression, when compared with the FA counterpart. These findings concluded that folate deficiency induced TL instability in both telomerase-negative and -positive cells, and FA was more efficient in maintaining TL and chromosome stability compared with 5-MeTHF.

The spiral gingers (Costus L.) are a pantropical genus of herbaceous perennial monocots; the Neotropical clade of Costus radiated rapidly in the past few million years into over 60 species. The Neotropical spiral gingers have a rich history of evolutionary and ecological research that can motivate and inform modern genetic investigations. Here, we present the first 2 chromosome-level genome assemblies in the genus, for C. pulverulentus and C. lasius, and briefly compare their synteny. We assembled the C. pulverulentus genome from a combination of short-read data, Chicago and Dovetail Hi-C chromatin-proximity sequencing, and alignment with a linkage map. We annotated the genome by mapping a C. pulverulentus transcriptome and querying mapped transcripts against a protein database. We assembled the C. lasius genome with Pacific Biosciences HiFi long reads and alignment to the C. pulverulentus genome. These 2 assemblies are the first published genomes for non-cultivated tropical plants. These genomes solidify the spiral gingers as a model system and will facilitate research on the poorly understood genetic basis of tropical plant diversification.

In most animals, mitochondrial DNA is strictly maternally inherited and non-recombining. One exception to these assumptions is called doubly uniparental inheritance (DUI): a phenomenon involving the independent transmission of female and male mitochondrial genomes. DUI is known only from the molluscan class Bivalvia. The phylogenetic distribution of male mitochondrial DNA in bivalves is consistent with several evolutionary scenarios, including multiple independent gains, losses, and varying degrees of recombination with female mitochondrial DNA. In this study, we use phylogenetic methods to test male mitochondrial DNA origination hypotheses and infer the prevalence of mitochondrial recombination in bivalves with DUI. Phylogenetic modeling using site concordance factors supported a single origin of male mitochondrial DNA in bivalves coupled with recombination acting over long evolutionary timescales. Ongoing mitochondrial recombination is present in Mytilida and Venerida, which results in a pattern of concerted evolution of female and male mitochondrial DNA. Mitochondrial recombination could be favored to offset the deleterious effects of asexual inheritance and maintain mitonuclear compatibility across tissues. Cardiida and Unionida have gone without recent recombination, possibly due to an extension of theCOX2gene in male mitochondrial DNA. The loss of recombination may be neutral but could be connected to the role of M mtDNA in sex determination or sexual development. Our results support recombination events in DUI species may occur throughout their genomes. Future investigations may reveal more complex patterns of inheritance of recombinants, which could explain the retention of signal for a single origination of male mitochondrial DNA in protein coding genes.

The cheetah (Acinonyx jubatus, SCHREBER 1775) is a large felid and is considered the fastest land animal. Historically, it inhabited open grassland across Africa, the Arabian Peninsula, and southwestern Asia; however, only small and fragmented populations remain today. Here, we present a de novo genome assembly of the cheetah based on PacBio continuous long reads and Hi-C proximity ligation data. The final assembly (VMU_Ajub_asm_v1.0) has a total length of 2.38 Gb, of which 99.7% are anchored into the expected 19 chromosome-scale scaffolds. The contig and scaffold N50 values of 96.8 Mb and 144.4 Mb, respectively, a BUSCO completeness of 95.4% and a k-mer completeness of 98.4%, emphasize the high quality of the assembly. Furthermore, annotation of the assembly identified 23,622 genes and a repeat content of 40.4%. This new highly contiguous and chromosome-scale assembly will greatly benefit conservation and evolutionary genomic analyses and will be a valuable resource, e.g., to gain a detailed understanding of the function and diversity of immune response genes in felids.

The Aeolian wall lizard, Podarcis raffonei, is an endangered species endemic to the Aeolian archipelago, Italy, where it is present only in 3 tiny islets and a narrow promontory of a larger island. Because of the extremely limited area of occupancy, severe population fragmentation and observed decline, it has been classified as Critically Endangered by the International Union for the Conservation of Nature (IUCN). Using Pacific Biosciences (PacBio) High Fidelity (HiFi) long-read sequencing, Bionano optical mapping and Arima chromatin conformation capture sequencing (Hi-C), we produced a high-quality, chromosome-scale reference genome for the Aeolian wall lizard, including Z and W sexual chromosomes. The final assembly spans 1.51 Gb across 28 scaffolds with a contig N50 of 61.4 Mb, a scaffold N50 of 93.6 Mb, and a BUSCO completeness score of 97.3%. This genome constitutes a valuable resource for the species to guide potential conservation efforts and more generally for the squamate reptiles that are underrepresented in terms of available high-quality genomic resources.

Genome sequences can reveal the extent of inbreeding in small populations. Here, we present the first genomic characterization of type D killer whales, a distinctive eco/morphotype with a circumpolar, subantarctic distribution. Effective population size is the lowest estimated from any killer whale genome and indicates a severe population bottleneck. Consequently, type D genomes show among the highest level of inbreeding reported for any mammalian species (FROH ≥ 0.65). Detected recombination cross-over events of different haplotypes are up to an order of magnitude rarer than in other killer whale genomes studied to date. Comparison of genomic data from a museum specimen of a type D killer whale that stranded in New Zealand in 1955, with 3 modern genomes from the Cape Horn area, reveals high covariance and identity-by-state of alleles, suggesting these genomic characteristics and demographic history are shared among geographically dispersed social groups within this morphotype. Limitations to the insights gained in this study stem from the nonindependence of the 3 closely related modern genomes, the recent coalescence time of most variation within the genomes, and the nonequilibrium population history which violates the assumptions of many model-based methods. Long-range linkage disequilibrium and extensive runs of homozygosity found in type D genomes provide the potential basis for both the distinctive morphology, and the coupling of genetic barriers to gene flow with other killer whale populations.

To investigate the distribution of polymorphisms and their frequent haplotypes in the regulatory region of MGMT in residents of high background radiation area (HBRA) and their impacts on transcriptional activity, we collected DNA samples from 83 healthy Chinese residents in HBRA and searched for genetic polymorphisms in the regulatory region of MGMT. Haplotypes were characterized by Haploview analysis. Transcriptional activities of different polymorphism haplotypes were detected by using a dual-luciferase reporter assay. Six genetic polymorphisms were identified within the regulatory region (1024 bp) of MGMT. Linkage disequilibrium (LD) patterns and haplotype profiles were analyzed using the identified genetic polymorphisms. These polymorphisms we found to be in high LD, with a Dʹ of 0.928 (r2 = 0.581) for −808 T>C and −19 C>T, 0.928 (r2 = 0.581) for −797 G>A and −19 C>T in Han Chinese HBRA residents. Complete LD with a Dʹ of 1.0 (r2 = 1.0) was observed between −808 T>C and −797 G>A. Haploview analysis revealed the existence of three polymorphism haplotypes in the core region of regulatory region of MGMT. Using serially truncated regulatory region of human MGMT luciferase reporter gene constructs, we found a 1002 bp (−637 nt to +365 nt) fragment in the MGMT gene was the core region. Dual-luciferase reporter assays showed that different polymorphism haplotypes bearing different variant alleles exhibit distinct transcriptional activities, especially the polymorphism haplotype carrying −19 T has the strongest transcriptional activity. In summary, the present study obtained genetic characteristics of the six polymorphisms in the regulatory region of the MGMT gene in HBRA residents, and the results suggest that different polymorphism haplotypes have significant effects on the transcriptional activity of the MGMT and that the −19 C>T polymorphism may be a functional variant involved in the transcriptional regulation of the MGMT gene.

Variation in DNA methylation is associated with many ecological and life history traits, including niche breadth and lifespan. In vertebrates, DNA methylation occurs almost exclusively at “CpG” dinucleotides. Yet, how variation in the CpG content of the genome impacts organismal ecology has been largely overlooked. Here, we explore associations between promoter CpG content, lifespan and niche breadth among 60, amniote vertebrate species. The CpG content of 16 functionally relevant gene promoters was strongly, positively associated with lifespan in mammals and reptiles, but was not related to niche breadth. Possibly, by providing more substrate for CpG methylation to occur, high promoter CpG content extends the time taken for deleterious, age-related errors in CpG methylation patterns to accumulate, thereby extending lifespan. The association between CpG content and lifespan was driven by gene promoters with intermediate CpG enrichment—those known to be predisposed to regulation by methylation. Our findings provide novel support for the idea that high CpG content has been selected for in long-lived species to preserve the capacity for gene expression regulation by CpG methylation. Intriguingly, promoter CpG content was also dependent on gene function in our study; immune genes had on average 20% less CpG sites than metabolic- and stress-related genes.

Testing the effect of rare variants on phenotypic variation is difficult due to the need for extremely large cohorts to identify associated variants given expected effect sizes. An alternative approach is to investigate the effect of rare genetic variants on low-level genomic traits, such as gene expression or DNA methylation (DNAm), as effect sizes are expected to be larger for low-level compared to higher-order complex traits. Here, we investigate DNAm in healthy ageing populations - the Lothian Birth cohorts of 1921 and 1936 and identify both transient and stable outlying DNAm levels across the genome. We find an enrichment of rare genetic variants within 1kb of DNAm sites in individuals with stable outlying DNAm, implying genetic control of this extreme variation. Using a family-based cohort, the Brisbane Systems Genetics Study, we observed increased sharing of DNAm outliers among more closely related individuals, consistent with these outliers being driven by rare genetic variation. We demonstrated that outlying DNAm levels have a functional consequence on gene expression levels, with extreme levels of DNAm being associated with gene expression levels towards the tails of the population distribution. Overall, this study demonstrates the role of rare variants in the phenotypic variation of low-level genomic traits, and the effect of extreme levels of DNAm on gene expression.

Postmating isolation is thought to be an important driver of the late stages of speciation. However, relatively little is empirically known about the process compared with other isolating mechanisms that drive the early stages of speciation, especially in non-model organisms. We characterized the genetic architecture of postmating isolation between 2 rockfishes, Sebastes schlegelii and S. trivittatus, whose reproductive isolation is complete. We examined transmission ratio distortion (TRD) patterns of genetic markers in 2 reciprocal backcross populations. Markers showing either of the 2 types of TRD was widespread across the genome, with some of the distorted markers forming extensive clusters around the recombination coldspots. These suggest that the postmating isolation effectively prevents gene flow across the genome and the recombination landscape contributes to the genetic architecture. Comparisons between 2 backcross families and 2 developmental stages showed little similarity in the distorted markers, suggesting asymmetry and stage specificity of the isolation. This may be due to hybrid incompatibility involving maternal factors or extrinsic selection. The lack of sex-ratio distortion in the mapping families suggested that Haldane’s rule in terms of hybrid inviability does not hold. Additionally, quantitative trait loci (QTL) mapping detected significant QTLs for sex and the morphological traits relevant to speciation and convergence of rockfishes, including body coloration. Genes in the melanocortin system, including agouti-signaling protein 1 (asip1) and melanocortin 1 receptor (mc1r), might underlie the horizontal and vertical color patterns on the body, respectively. These findings constitute an essential step toward a comprehensive understanding of speciation and morphological diversification of rockfishes.

Spinal muscular atrophy (SMA) is a monogenic disease that clinically manifests as severe muscle weakness owing to neurotransmission defects and motoneuron degeneration. Individuals affected by SMA experience neuromuscular weakness that impacts functional activities of daily living. We have used a mouse model of severe SMA (SMNΔ7) to test whether a calcium channel gating modifier (GV-58), alone or in combination with a potassium channel antagonist (3,4-diaminopyridine; 3,4-DAP), can improve neuromuscular function in this mouse model. Bath application of GV-58 alone or in combination with 3,4-DAP significantly restored neuromuscular transmission to control levels in both a mildly vulnerable forearm muscle and a strongly vulnerable trunk muscle in SMNΔ7 mice at postnatal days 10–12. Similarly, acute subcutaneous administration of GV-58 to postnatal day 10 SMNΔ7 mice, alone or in combination with 3,4-DAP, significantly increased a behavioral measure of muscle strength. These data suggest that GV-58 may be a promising treatment candidate that could address deficits in neuromuscular function and strength and that the addition of 3,4-DAP to GV-58 treatment could aid in restoring function in SMA.

Dysregulation of circulating lipids is a central element for the metabolic syndrome. However, it is not well established whether human subcutaneous adipose tissue is affected by or affect circulating lipids through epigenetic mechanisms. Hence, our aim was to investigate the association between circulating lipids and DNA methylation levels in human adipose tissue. DNA methylation and gene expression were analysed genome-wide in subcutaneous adipose tissue from two different cohorts, including 85 men and 93 women, respectively. Associations between DNA methylation and circulating levels of triglycerides, low-density lipoprotein, high-density lipoprotein and total cholesterol were analysed. Causal mediation analyses tested if adipose tissue DNA methylation mediates the effects of triglycerides on gene expression or insulin resistance. We found 115 novel associations between triglycerides and adipose tissue DNA methylation, e.g. in the promoter of RFS1, ARID2 and HOXA5 in the male cohort (P ≤ 1.1 × 10−7), and 63 associations, e.g. within the gene body of PTPRN2 and COL6A3 in the female cohort. We further connected these findings to altered mRNA expression levels in adipose tissue (e.g. HOXA5, IL11 and FAM45B). Interestingly, there was no overlap between methylation sites associated with triglycerides in men and the sites found in women, which points towards sex-specific effects of triglycerides on the epigenome. Finally, a causal mediation analysis provided support for adipose tissue DNA methylation as a partial mediating factor between circulating triglycerides and insulin resistance. This study identified novel epigenetic alterations in adipose tissue associated with circulating lipids. Identified epigenetic changes seem to mediate effects of triglycerides on insulin resistance.

Bi-allelic mutations in GBA1, the gene that encodes β-glucocerebrosidase (GCase), cause Gaucher disease (GD), whereas mono-allelic mutations do not cause overt pathology. Yet mono- or bi-allelic GBA1 mutations are the highest known risk factor for Parkinson’s disease (PD). GCase deficiency results in the accumulation of glucosylceramide (GluCer) and its deacylated metabolite glucosylsphingosine (GluSph). Brains from patients with neuronopathic GD have high levels of GluSph, and elevation of this lipid in GBA1-associated PD has been reported. To uncover the mechanisms involved in GBA1-associated PD, we used human induced pluripotent stem cell-derived dopaminergic (DA) neurons from patients harboring heterozygote mutations in GBA1 (GBA1/PD–DA neurons). We found that compared with gene-edited isogenic controls, GBA1/PD–DA neurons exhibit mammalian target of rapamycin complex 1 (mTORC1) hyperactivity, a block in autophagy, an increase in the levels of phosphorylated α-synuclein (129) and α-synuclein aggregation. These alterations were prevented by incubation with mTOR inhibitors. Inhibition of acid ceramidase, the lysosomal enzyme that deacylates GluCer to GluSph, prevented mTOR hyperactivity, restored autophagic flux and lowered α-synuclein levels, suggesting that GluSph was responsible for these alterations. Incubation of gene-edited wild type (WT) controls with exogenous GluSph recapitulated the mTOR/α-synuclein abnormalities of GBA1/PD neurons, and these phenotypic alterations were prevented when GluSph treatment was in the presence of mTOR inhibitors. We conclude that GluSph causes an aberrant activation of mTORC1, suppressing normal lysosomal functions, including the clearance of pathogenic α-synuclein species. Our results implicate acid ceramidase in the pathogenesis of GBA1-associated PD, suggesting that this enzyme is a potential therapeutic target for treating synucleinopathies caused by GCase deficiency.

Exposure of tobacco workers handling dried tobacco leaves has been linked to an increased risk of toxicity and respiratory illness due to the presence of nicotine and other chemicals. This study aimed to evaluate the DNA damage caused by the exposure of tobacco growers during the dry leaf classification process and the relation to cellular mechanisms. A total of 86 individuals participated in the study, divided into a group exposed to dry tobacco (n = 44) and a control group (n = 42). Genotoxicity was evaluated using the alkaline comet assay and lymphocyte micronucleus (MN) assay (CBMN-Cyt), and measurement of telomere length. The levels of oxidative and nitrosative stress were evaluated through the formation of thiobarbituric acid reactive species, and nitric oxide levels, respectively. The inorganic elements were measured in the samples using particle-induced X-ray emission method. The combination of variables was demonstrated through principal component analysis and the interactions were expanded through systems biology. Comet assay, MN, death cells, thiobarbituric acid reactive species, and nitrosative stress showed a significant increase for all exposed groups in relation to the control. Telomere length showed a significant decrease for exposed women and total exposed group in relation to men and control groups, respectively. Bromine (Br) and rubidium (Rb) in the exposed group presented higher levels than control groups. Correlations between nitrate and apoptosis; Br and MN and necrosis; and Rb and telomeres; besides age and DNA damage and death cells were observed. The systems biology analysis demonstrated that tobacco elements can increase the nuclear translocation of NFKB dimers inducing HDAC2 expression, which, associated with BRCA1 protein, can potentially repress transcription of genes that promote DNA repair. Dry tobacco workers exposed to dry leaves and their different agents showed DNA damage by different mechanisms, including redox imbalance.

Heart failure (HF) is a complex clinical condition characterized by functional and structural defects of the myocardium, but genetic and environmental factors are considered to play an important role in the development of the disease. In the present study, we investigated the genome instability (DNA and chromosomal damage) in patients with heart failure with reduced ejection fraction (HFrEF) ≤40% and its association with risk factors. The studied population included 48 individuals, of which 29 HFrEF patients (mean age 57.41 ± 5.74 years) and 19 healthy controls (mean age 57.63 ± 6.09 years). The genetic damage index in peripheral blood lymphocytes was analyzed using the comet assay, while micronuclei frequency and nuclear division index were analyzed using the cytokinesis-block micronucleus assay. Our results showed that HFrEF patients had a significantly higher genetic damage index compared with the healthy controls (P < .001). Cytokinesis-block micronucleus assay showed that the average micronucleus frequency in peripheral blood lymphocytes of patients was significantly higher, while the nuclear division index values were significantly lower than in controls (P < .01). Using multiple linear regression analysis, pathological state, ejection fraction, creatinine, glucose, associated disease, residence, proBNP, troponin, urea, ACE-inhibitors, and length of the drug therapy were identified as predictors of DNA and/or chromosomal damage in HF patients. We can conclude that DNA and chromosomal damage was increased in patients with HF, which may be a consequence of disease and/or drug therapy.

Breast cancer (bc) is the second most common type of human malignancies with highest morbidity and mortality in the female population. Therefore, it is essential to develop novel and effective therapies for bc treatment. The main aim of the current study is to investigate the functions of CEBPB and THBS2 in bc and the underlying mechanism. Reverse transcription–quantitative real-time polymerase chain reaction and western blot were performed for the measurement of ribonucleic acids and proteins. Function and mechanism assays were, respectively, conducted for the evaluation of bc biological behaviors and exploration of the potential correlation of genes. According to bioinformatics analyses and experimental results, THBS2, up-regulated in bc tissues and cell lines, could facilitate cell migration, invasion and EMT in bc. CEBPB was validated to facilitate miR-29a-3p transcription, thus negatively modulating THBS2 expression. The results of rescue experiments reflected that CEBPB could regulate the malignant behaviors of bc cells via THBS2. Furthermore, CEBPB was ascertained to inhibit the transcription of B3GALTL to affect THBS2 protein O-fucosylation and secretion. The interaction between THBS2 and ITGB1 was confirmed, and THBS2 was found to activate the PI3K/AKT signaling pathway. To conclude, CEBPB could restrain bc cell migration, invasion and EMT via inhibition on THBS2 expression and O-fucosylation.

Human DUX4 and its mouse ortholog Dux are normally expressed in the early embryo—the 4-cell or 2-cell cleavage stage embryo, respectively—and activate a portion of the first wave of zygotic gene expression. DUX4 is epigenetically suppressed in nearly all somatic tissue, whereas facioscapulohumeral dystrophy (FSHD)-causing mutations result in its aberrant expression in skeletal muscle, transcriptional activation of the early embryonic program and subsequent muscle pathology. Although DUX4 and Dux both activate an early totipotent transcriptional program, divergence of their DNA binding domains limits the use of DUX4 expressed in mice as a preclinical model for FSHD. In this study, we identify the porcine DUXC messenger ribonucleic acid expressed in early development and show that both pig DUXC and human DUX4 robustly activate a highly similar early embryonic program in pig muscle cells. These results support further investigation of pig preclinical models for FSHD.

Pathogenic variations in the sodium voltage-gated channel alpha subunit 1 (SCN1A) gene are responsible for multiple epilepsy phenotypes, including Dravet syndrome, febrile seizures (FS) and genetic epilepsy with FS plus. Phenotypic heterogeneity is a hallmark of SCN1A-related epilepsies, the causes of which are yet to be clarified. Genetic variation in the non-coding regulatory regions of SCN1A could be one potential causal factor. However, a comprehensive understanding of the SCN1A regulatory landscape is currently lacking. Here, we summarized the current state of knowledge of SCN1A regulation, providing details on its promoter and enhancer regions. We then integrated currently available data on SCN1A promoters by extracting information related to the SCN1A locus from genome-wide repositories and clearly defined the promoter and enhancer regions of SCN1A. Further, we explored the cellular specificity of differential SCN1A promoter usage. We also reviewed and integrated the available human brain-derived enhancer databases and mouse-derived data to provide a comprehensive computationally developed summary of SCN1A brain-active enhancers. By querying genome-wide data repositories, extracting SCN1A-specific data and integrating the different types of independent evidence, we created a comprehensive catalogue that better defines the regulatory landscape of SCN1A, which could be used to explore the role of SCN1A regulatory regions in disease.

The testis-specific adenosine deaminase domain-containing (ADAD) protein family, including ADAD1 and ADAD2, has been confirmed to be essential in mouse male fertility. However, the roles of ADAD1 and ADAD2 in human reproductive biology are unclear. Herein, whole-exome sequencing was conducted for 337 infertile patients to detect pathogenic variants in ADAD1 and ADAD2. Importantly, a novel deleterious biallelic variant of NM_001159285.2:c.1408G > T (p.V470F) and NM_001159285.2:c.1418A > G (p.E473G) in ADAD1 and a pathogenic homozygous missense variant of NM_001145400.2:c.1381C > T (p.R461W) in ADAD2 were identified in this infertile cohort with frequencies of 0.29 (1/337) and 0.59% (2/337), respectively. Electron microscopy revealed an abnormal morphology and severely disorganized ultrastructure of sperm from the patients. Immunofluorescence and western blotting showed a sharp decrease in ADAD1 and ADAD2 expression in sperm from the patients. Mechanistically, bioinformatics analysis suggested that ADAD2 interacts with DNAH17. Furthermore, we demonstrated that the expression of DNAH17 was markedly downregulated in the sperm of patients harboring ADAD2 variants. In addition, the expression of several autophagy regulators was significantly disrupted in the sperm of patients harboring ADAD2 variants. In conclusion, we identified novel ADAD1 and ADAD2 variants in three infertile patients from a large infertile cohort, first providing evidence that ADAD1 and ADAD2 variants might be a candidate genetic cause of human male infertility. Moreover, an important new dimension to our understanding of the genotype–phenotype correlations between the ADAD gene family and male infertility in humans has been uncovered, providing valuable information for the genetic diagnosis of male infertility.

Oligoasthenoteratozoospermia (OAT) can result in male infertility owing to reduced sperm motility and abnormal spermatozoan morphology. The Tektins are a family of highly conserved filamentous proteins expressed in the axoneme and associated structures in many different metazoan species. Earlier studies on mice identified Tektin3 (Tekt3) as a testis-enriched gene, and knockout of Tekt3 resulted in asthenozoospermia in the mice. Here, whole-exome sequencing of 100 males with asthenozoospermia from unrelated families was performed, followed by Sanger sequencing, leading to the identification of TEKT3 as a candidate gene in two of these patients and their associated family members. In total, three mutations in the TEKT3 gene were identified in both these patients, including one homozygous deletion–insertion mutation (c.543_547delinsTTGAT: p.Glu182*) and one compound heterozygous mutation (c.[548G > A]; [752A > C], p.[Arg183Gln]; [Gln251Pro]). Both of these mutations resulted in the complete loss of TEKT3 expression. The patients were both found to produce sperm that, although those showed no apparent defects in the flagellar structure, had reduced progressive motility. In contrast to mice, most sperm from these two patients exhibited acrosomal hypoplasia, although this did not prevent the use of the sperm for in vitro fertilization through an ICSI approach. TEKT3 was found to bind to other TEKT proteins, suggesting that these proteins form a complex within human spermatozoa. Overall, these results suggest that a loss of TEKT3 function can contribute to OAT incidence in humans. TEKT3 deficiencies can reduce sperm motility and contribute to severe acrosomal hypoplasia in spermatozoa, compromising their normal function.

Robert (Bob) Keith Wayne lost his battle with cancer in his home on 26 December 2022 with his wife, Dr. Blaire Van Valkenburgh, by his side. This essay, written by his former graduate students, highlights the foundation that survives Bob and our vision of continuing his efforts in building a future for endangered species.

Isolated hypogonadotropic hypogonadism (IHH) is a rare disease with hypogonadism and infertility caused by the defects in embryonic migration of hypothalamic gonadotropin-releasing hormone (GnRH) neurons, hypothalamic GnRH secretion or GnRH signal transduction. PROKR2 gene, encoding a G-protein coupled receptor PROKR2, is one of the most frequently mutated genes identified in IHH patients. However, the functional consequences of several PROKR2 mutants remain elusive. In this study, we systematically analyzed the Gαq, Gαs and ERK1/2 signaling of 23 IHH-associated PROKR2 mutations which are yet to be functionally characterized. We demonstrate that blockage of Gαq, instead of MAPK/ERK pathway, inhibited PROK2-induced migration of PROKR2-expressing cells, implying that PROKR2-related IHH results primarily due to Gαq signaling pathway disruption. Combined with previous reports, we categorized a total of 63 IHH-associated PROKR2 mutations into four distinct groups according Gαq pathway functionality: (i) neutral (N, >80% activity); (ii) low pathogenicity (L, 50–80% activity); (iii) medium pathogenicity (M, 20–50% activity) and (iv) high pathogenicity (H

Nemaline myopathy (NM) is a rare neuromuscular disorder associated with congenital or childhood-onset of skeletal muscle weakness and hypotonia, which results in limited motor function. NM is a genetic disorder and mutations in 12 genes are known to contribute to autosomal dominant or recessive forms of the disease. Recessive mutations in nebulin (NEB) are the most common cause of NM affecting about 50% of patients. Because of the large size of the NEB gene and lack of mutational hot spots, developing therapies that can benefit a wide group of patients is challenging. Although there are several promising therapies under investigation, there is no cure for NM. Therefore, targeting disease modifiers that can stabilize or improve skeletal muscle function may represent alternative therapeutic strategies. Our studies have identified Nrap upregulation in nebulin deficiency that contributes to structural and functional deficits in NM. We show that genetic ablation of nrap in nebulin deficiency restored sarcomeric disorganization, reduced protein aggregates and improved skeletal muscle function in zebrafish. Our findings suggest that Nrap is a disease modifier that affects skeletal muscle structure and function in NM; thus, therapeutic targeting of Nrap in nebulin-related NM and related diseases may be beneficial for patients.

Despite increasing sequencing efforts, numerous fish families still lack a reference genome, which complicates genetic research. One such understudied family is the sand lances (Ammodytidae, literally: “sand burrower”), a globally distributed clade of over 30 fish species that tend to avoid tidal currents by burrowing into the sand. Here, we present the first annotated chromosome-level genome assembly of the great sand eel (Hyperoplus lanceolatus). The genome assembly was generated using Oxford Nanopore Technologies long sequencing reads and Illumina short reads for polishing. The final assembly has a total length of 808.5 Mbp, of which 97.1% were anchored into 24 chromosome-scale scaffolds using proximity-ligation scaffolding. It is highly contiguous with a scaffold and contig N50 of 33.7 and 31.3 Mbp, respectively, and has a BUSCO completeness score of 96.9%. The presented genome assembly is a valuable resource for future studies of sand lances, as this family is of great ecological and commercial importance and may also contribute to studies aiming to resolve the suprafamiliar taxonomy of bony fishes.

Neuroinflammation and immune dysregulation play a key role in Alzheimer’s disease (AD) and are also associated with severe Covid-19 and neurological symptoms. Also, genome-wide association studies found many risk single nucleotide polymorphisms (SNPs) for AD and Covid-19. However, our understanding of underlying gene regulatory mechanisms from risk SNPs to AD, Covid-19 and phenotypes is still limited. To this end, we performed an integrative multi-omics analysis to predict gene regulatory networks for major brain regions from population data in AD. Our networks linked transcription factors (TFs) to TF binding sites (TFBSs) on regulatory elements to target genes. Comparative network analyses revealed cross-region-conserved and region-specific regulatory networks, in which many immunological genes are present. Furthermore, we identified a list of AD–Covid genes using our networks involving known and Covid-19 genes. Our machine learning analysis prioritized 36 AD–Covid candidate genes for predicting Covid severity. Our independent validation analyses found that these genes outperform known genes for classifying Covid-19 severity and AD. Finally, we mapped genome-wide association study SNPs of AD and severe Covid that interrupt TFBSs on our regulatory networks, revealing potential mechanistic insights of those disease risk variants. Our analyses and results are open-source available, providing an AD–Covid functional genomic resource at the brain region level.

Age-related macular degeneration (AMD) is the most prevalent cause of blindness in the developed world. Vision loss in the advanced stages of the disease is caused by atrophy of retinal photoreceptors, overlying retinal pigment epithelium (RPE) and choroidal endothelial cells. The molecular events that underline the development of these cell types from in utero to adult as well as the progression to intermediate and advanced stages AMD are not yet fully understood. We performed single-cell RNA-sequencing (RNA-Seq) of human fetal and adult RPE–choroidal tissues, profiling in detail all the cell types and elucidating cell type-specific proliferation, differentiation and immunomodulation events that occur up to midgestation. Our data demonstrate that progression from the fetal to adult state is characterized by an increase in expression of genes involved in the oxidative stress response and detoxification from heavy metals, suggesting a better defence against oxidative stress in the adult RPE–choroid tissue. Single-cell comparative transcriptional analysis between a patient with intermediate AMD and an unaffected subject revealed a reduction in the number of RPE cells and melanocytes in the macular region of the AMD patient. Together these findings may suggest a macular loss of RPE cells and melanocytes in the AMD patients, but given the complex processing of tissues required for single-cell RNA-Seq that is prone to technical artefacts, these findings need to be validated by additional techniques in a larger number of AMD patients and controls.

Membrane fusion is mediated by soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins. During neurotransmitter exocytosis, SNARE proteins on a synaptic vesicle and the target membrane form a complex, resulting in neurotransmitter release. N-ethylmaleimide-sensitive factor (NSF), a homohexameric ATPase, disassembles the complex, allowing individual SNARE proteins to be recycled. Recently, the association between pathogenic NSF variants and developmental and epileptic encephalopathy (DEE) was reported; however, the molecular pathomechanism of NSF-related DEE remains unclear. Here, three patients with de novo heterozygous NSF variants were presented, of which two were associated with DEE and one with a very mild phenotype. One of the DEE patients also had hypocalcemia from parathyroid hormone deficiency and neuromuscular junction impairment. Using PC12 cells, a neurosecretion model, we show that NSF with DEE-associated variants impaired the recycling of vesicular membrane proteins and vesicle enlargement in response to exocytotic stimulation. In addition, DEE-associated variants caused neurodegenerative change and defective autophagy through overactivation of the mammalian/mechanistic target of rapamycin (mTOR) pathway. Treatment with rapamycin, an mTOR inhibitor or overexpression of wild-type NSF ameliorated these phenotypes. Furthermore, neurons differentiated from patient-derived induced pluripotent stem cells showed neurite degeneration, which was also alleviated by rapamycin treatment or gene correction using genome editing. Protein structure analysis of NSF revealed that DEE-associated variants might disrupt the transmission of the conformational change of NSF monomers and consequently halt the rotation of ATP hydrolysis, indicating a dominant negative mechanism. In conclusion, this study elucidates the pathomechanism underlying NSF-related DEE and identifies a potential therapeutic approach.

The abnormal expansion of GGGGCC hexanucleotide repeats within the C9orf72 gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The accumulation of GGGGCC repeat-containing RNAs as RNA foci, and the deposition of dipeptide repeat proteins (DPR) produced from these repeat RNAs by unconventional translation are major pathological hallmarks of C9orf72-linked ALS/FTD (C9-ALS/FTD), and are both thought to play a crucial role in the pathogenesis of these diseases. Because GGGGCC repeat RNA is likely to be the most upstream therapeutic target in the pathogenic cascade of C9-ALS/FTD, lowering the cellular level of GGGGCC repeat RNA is expected to mitigate repeat RNA toxicity, and will therefore be a disease-modifying therapeutic strategy for the treatment of C9-ALS/FTD. In this study, we demonstrated using a Drosophila model of C9-ALS/FTD that elevated expression of a subset of human RNA-binding proteins that bind to GGGGCC repeat RNA, including hnRNPA3, IGF2BP1, hnRNPA2B1, hnRNPR and SF3B3, reduces the level of GGGGCC repeat RNA, resulting in the suppression of neurodegeneration. We further showed that hnRNPA3-mediated reduction of GGGGCC repeat RNA suppresses disease pathology, such as RNA foci and DPR accumulation. These results demonstrate that hnRNPA3 and other RNA-binding proteins negatively regulate the level of GGGGCC repeat RNA, and mitigate repeat RNA toxicity in vivo, indicating the therapeutic potential of the repeat RNA-lowering approach mediated by endogenous RNA-binding proteins for the treatment of C9-ALS/FTD.

Atopic dermatitis is a chronically recurrent dermatologic disease affected by complex pathophysiology with limited therapeutic options. To identify promising biomarkers for atopic dermatitis, we conducted a Mendelian randomization (MR) study to systematically screen blood metabolome for potential causal mediators of atopic dermatitis and further predict target-mediated side effects. We selected 128 unique blood metabolites from three European-descent metabolome genome-wide association studies (GWASs) with a total of 147 827 participants. Atopic dermatitis dataset originated from a large-scale GWAS including 10 788 cases and 30 047 controls of European ancestry. MR analyses were performed to estimate the associations of blood metabolites with atopic dermatitis. We then applied a phenome-wide MR analysis to ascertain potential on-target side effects of metabolite intervention. Three metabolites were identified as potential causal mediators for atopic dermatitis, including docosahexaenoic acid (odds ratio [OR], 0.87; 95% confidence interval [CI], 0.81–0.94; P = 3.45 × 10−4), arachidonate (OR, 0.30; 95% CI, 0.17–0.53; P = 4.09 × 10−5) and 1-arachidonoylglycerophosphoethanolamine (1-arachidonoyl-GPE) (OR, 0.25; 95% CI, 0.12–0.53; P = 2.58 × 10−4). In the phenome-wide MR analysis, docosahexaenoic acid and arachidonate were also identified to have beneficial or detrimental effects on multiple diseases beyond atopic dermatitis, respectively. No adverse side effects were found for 1-arachidonoyl-GPE. In this systematic MR study, docosahexaenoic acid, arachidonate and 1-arachidonoyl-GPE were identified as potential causal and beneficial mediators in the development of atopic dermatitis. Side-effect profiles were characterized to help inform drug target prioritization, and 1-arachidonoyl-GPE was a promising target for prevention and treatment of atopic dermatitis with no predicted adverse side effects.

As the auditory and balance receptor cells in the inner ear, hair cells are responsible for converting mechanical stimuli into electrical signals, a process referred to as mechano-electrical transduction. Hair cell development and function are tightly regulated, and hair cell deficits are the main reasons for hearing loss and balance disorders. TMCC2 is an endoplasmic reticulum (ER)-residing transmembrane protein whose physiological function largely remains unknown. In the present work, we show that Tmcc2 is specifically expressed in the auditory hair cells of mouse inner ear. Tmcc2 knockout mice were then established to investigate its physiological role in hearing. Auditory brainstem responses measurements show that Tmcc2 knockout mice suffer from congenital hearing loss. Further investigations reveal progressive auditory hair cell loss in the Tmcc2 knockout mice. The general morphology and function of ER are unaffected in Tmcc2 knockout hair cells. However, increased ER stress was observed in Tmcc2 knockout mice and knockdown cells, suggesting that loss of TMCC2 leads to auditory hair cell death through elevated ER stress.

Despite the increasing feasibility of sequencing whole genomes from diverse taxa, a persistent problem in phylogenomics is the selection of appropriate genetic markers or loci for a given taxonomic group or research question. In this review, we aim to streamline the decision-making process when selecting specific markers to use in phylogenomic studies by introducing commonly used types of genomic markers, their evolutionary characteristics, and their associated uses in phylogenomics. Specifically, we review the utilities of ultraconserved elements (including flanking regions), anchored hybrid enrichment loci, conserved nonexonic elements, untranslated regions, introns, exons, mitochondrial DNA, single nucleotide polymorphisms, and anonymous regions (nonspecific regions that are evenly or randomly distributed across the genome). These various genomic elements and regions differ in their substitution rates, likelihood of neutrality or of being strongly linked to loci under selection, and mode of inheritance, each of which are important considerations in phylogenomic reconstruction. These features may give each type of marker important advantages and disadvantages depending on the biological question, number of taxa sampled, evolutionary timescale, cost effectiveness, and analytical methods used. We provide a concise outline as a resource to efficiently consider key aspects of each type of genetic marker. There are many factors to consider when designing phylogenomic studies, and this review may serve as a primer when weighing options between multiple potential phylogenomic markers.

Monoallelic or biallelic RAD51C germline mutations results in chromosome instability disorders such as fanconi anemia (FA) and cancers. The bona fide function of RAD51C is to assist RAD51 nucleoprotein filament onto single strand DNA to complete Homologous Recombination (HR) repair. In addition to HR repair, the role of RAD51C in DNA replication is emerging when replication forks are transiently or irreversibly stalled. We identified novel RAD51C variants of uncertain significance (VUS) from breast, ovarian, pancreatic, and gastric cancer patients and functionally characterized the effect of these variants in replication fork protection and double-strand breaks (DSB's) repair. In RAD51C deficient chinese hamster CL-V4B cells, expression of RAD51C F164S, A87E, L134S and E49K variants heightened sensitivity to mitomycin C (MMC), etoposide and PARP inhibition. Differently, expression of a subset RAD51C variants R24L, R24W, and R212H displayed mild sensitivity to MMC, etoposide and PARP inhibition. Further functional characterization of a subset of variants revealed that Rad51C F164S, A87E, L134S and E49K variants displayed reduced RAD51 foci formation and increased overall nuclear single strand DNA levels in the presence of replication stress. Additionally, DNA fiber assay revealed that RAD51C F164S, A87E, L134S and E49K variants displayed defective replication fork protection upon prolonged fork stalling. Investigations using patient-derived lymphoblastoid cell line carrying heterozygous RAD51C L134S variant showed an impairment in RAD51 chromatin association and replication fork protection, suggestive of deleteriousness of this VUS variant. Overall, our findings provide more insights into molecular roles of RAD51C in replication fork integrity maintenance and in DSB repair. Disclosure Summary The authors report no relevant conflicts of interest in this work. Ministry of Education (MOE) National Medical Research Council (NMRC) Singapore Ministry of Education (MOE) under its Academic Research Fund ( AcRF) Tier 1 [2019-T1-001-018]; National Cancer Centre Research Fund Terry Fox Grant [NCCRF-YR2018-NOV-1]; National Medical Research Council (NMRC) Clinician Scientist Award [MOH-000654].

The Tibetan partridge (Perdix hodgsoniae) is a widely distributed endemic species in high-altitude areas across the Tibetan Plateau where the hypoxia, lower temperature and high ultraviolet radiation are pivotal factors influencing survival. However, the underlying genetic adaptation of the Tibetan partridge to extreme environments remains uncertain due to limited genomic resources. Similarly, the phylogenetic position of Perdix within Phasianidae remains controversial due to lacking information. Consequently, we de novo assembled and annotated the whole genome of the Tibetan partridge. The genome size was 1.15 Gb with contig N50 of 3.70 Mb. A total of 202.30 Mb (17.61%) repetitive elements and 445,876 perfect microsatellites were identified. A total of 16,845 functionally annotated protein-coding genes were identified in the Tibetan partridge. Genomic phylogenetic analysis across 30 Galliformes species indicated a close relationship between Perdix and typical pheasants composed of Chrysolophus, Symaticus, Phasianus, Crossopilon, and Lophura. However, the phylogenetic relationship of (Perdix + (Chrysolophus + (Syrmaticus + other pheasants))) was different from those of (Perdix + (Syrmaticus + (Chrysolophus + other pheasants))) in previous studies. Comparative genomic results identified NFKB1 and CREBBP positively selected genes related to hypoxia with 3 and 2 Tibetan partridge-specific missense mutations, respectively. Expanded gene families were mainly associated with energy metabolism and steroid hydroxylase activity, meanwhile, contracted gene families were mainly related to immunity and olfactory perception. Our genomic data considerably contribute to the phylogeny of Perdix and the underlying adaptation strategies of the Tibetan partridge to a high-altitude environment.

A functional nerve growth factor NGF–Tropomyosin Receptor kinase A (TrkA) system is an essential requisite for the generation and maintenance of long-lasting thermal and mechanical hyperalgesia in adult mammals. Indeed, mutations in the gene encoding for TrkA are responsible for a rare condition, named Hereditary Sensory and Autonomic Neuropathy type IV (HSAN IV), characterized by the loss of response to noxious stimuli, anhidrosis and cognitive impairment. However, to date, there is no available mouse model to properly understand how the NGF–TrkA system can lead to pathological phenotypes that are distinctive of HSAN IV. Here, we report the generation of a knock-in mouse line carrying the HSAN IV TrkAR649W mutation. First, by in vitro biochemical and biophysical analyses, we show that the pathological R649W mutation leads to kinase-inactive TrkA also affecting its membrane dynamics and trafficking. In agreement with the HSAN IV human phenotype, TrkAR649W/m mice display a lower response to thermal and chemical noxious stimuli, correlating with reduced skin innervation, in addition to decreased sweating in comparison to TrkAh/m controls. Moreover, the R649W mutation decreases anxiety-like behavior and compromises cognitive abilities, by impairing spatial-working and social memory. Our results further uncover unexplored roles of TrkA in thermoregulation and sociability. In addition to accurately recapitulating the clinical manifestations of HSAN IV patients, our findings contribute to clarifying the involvement of the NGF–TrkA system in pain sensation.

Nuclear mitochondrial pseudogenes (numts) may hinder the reconstruction of mtDNA genomes and affect the reliability of mtDNA datasets for phylogenetic and population genetic comparisons. Here, we present the program Numt Parser, which allows for the identification of DNA sequences that likely originate from numt pseudogene DNA. Sequencing reads are classified as originating from either numt or true cytoplasmic mitochondrial (cymt) DNA by direct comparison against cymt and numt reference sequences. Classified reads can then be parsed into cymt or numt datasets. We tested this program using whole genome shotgun-sequenced data from two ancient Cape lions (Panthera leo), because mtDNA is often the marker of choice for ancient DNA studies and the genus Panthera is known to have numt pseudogenes. Numt Parser decreased sequence disagreements that were likely due to numt pseudogene contamination and equalized read coverage across the mitogenome by removing reads that likely originated from numts. We compared the efficacy of Numt Parser to two other bioinformatic approaches that can be used to account for numt contamination. We found that Numt Parser outperformed approaches that rely only on read alignment or Basic Local Alignment Search Tool (BLAST) properties, and was effective at identifying sequences that likely originated from numts while having minimal impacts on the recovery of cymt reads. Numt Parser therefore improves the reconstruction of true mitogenomes, allowing for more accurate and robust biological inferences.

Non-syndromic sensorineural hearing loss (NSHL) is a group of genetically heterogeneous conditions with broad phenotypic heterogeneity. There is, at present, no curative treatment for genetic hearing loss (HL). Early molecular diagnosis of progressive disorders and elucidation of the causes and pathomechanisms are essential for developing therapeutic strategies. Here, we identified a novel rare frameshift variant of LMX1A (c.915dup), which resulted in the C-terminal-altered and -truncated LMX1A (p.Val306Cysfs*32). This C-terminal frameshift mutation co-segregated with autosomal dominant (AD) NSHL in a four-generation Chinese family, suggesting that the LMX1A non-missense mutation is also contributed to ADNSHL. In this family, the affected individuals exhibited the variable auditory phenotypes ranging from profound congenital deafness at birth or to mild/moderate HL in adulthood. We also found that the embryonic cells carrying with the heterozygous variant significantly expressed several upregulated HL-associated genes at transcriptional level. In vitro splicing assay suggested that the LMX1A mRNA with c.915dup did not cause nonsense-mediated decay and was translated into a truncated LMX1A. In addition, electrophoresis mobility shift assay and luciferase assays have shown that the highly conserved C-terminal domain (amino acid 306–382) of the LMX1A was required for regulating the protein–DNA interaction and transactivation in vitro. Furthermore, apoptosis assays suggested that the C-terminal domain of the LMX1A was important for mediating apoptosis in the cochlear hair cells. Our work provided the multiline of the evidence to support that non-missense mutation of LMX1A leads to ADNSHL and the C-terminal domain of LMX1A is important for mediating transcriptional activity and associated with promoting apoptosis in the cells.

Inferences made from molecular data support regional stock assessment goals by providing insights into the genetic population dynamics of enigmatic species. Population genomics metrics, such as genetic diversity and population connectivity, serve as useful proxies for species health and stability. Sleeper sharks (genus Somniosus) are ecologically important deep-sea predators, estimated to reach ages of 250 to 300 yr and taking decades to reach sexual maturity. The subgenus Somniosus (Somniosus) is comprised of 3 species: S. pacificus, S. microcephalus, and S. antarcticus. Given the life history strategy of somniosids, they are vulnerable to overfishing and population declines. Further, data to assess the stocks of these species are limited. To address this deficiency, we used the reduced representation library method Restriction-site Associated DNA sequencing (RADseq) to conduct phylogenomic and population genomics analyses, providing novel information for use in stock assessments. Our results strongly support the species status of S. microcephalus (N = 79), but recover S. antarcticus (N = 2) intermixed within the S. pacificus (N = 170) clade. Population genomics analyses reveal genetic homogeneity within S. pacificus and S. microcephalus, and estimates of effective population size were in the hundreds for both species. Kinship analysis identified 2 first-degree relative pairs within our dataset (1 within each species). Our results contribute new information for stock assessments of these uniquely long-lived species by providing the strongest molecular evidence to date for the synonymization of S. antarcticus and S. pacificus, as well as estimating population genomic metrics for each supported species within the Somniosus (Somniosus) subgenus.

Red abalone, Haliotis rufescens, are herbivorous marine gastropods that primarily feed on kelp. They are the largest and longest-lived of abalone species with a range distribution in North America from central Oregon, United States, to Baja California, MEX. Recently, red abalone have been in decline as a consequence of overharvesting, disease, and climate change, resulting in the closure of the commercial fishery in the 1990s and the recreational fishery in 2018. Protecting this ecologically and economically important species requires an understanding of their current population dynamics and connectivity. Here, we present a new red abalone reference genome as part of the California Conservation Genomics Project (CCGP). Following the CCGP genome strategy, we used Pacific Biosciences HiFi long reads and Dovetail Omni-C data to generate a scaffold-level assembly. The assembly comprises 616 scaffolds for a total size of 1.3 Gb, a scaffold N50 of 45.7 Mb, and a BUSCO complete score of 97.3%. This genome represents a significant improvement over a previous assembly and will serve as a powerful tool for investigating seascape genomic diversity, local adaptation to temperature and ocean acidification, and informing management strategies.

Incorporating measures of taxonomic diversity into research and management plans has long been a tenet of conservation science. Increasingly, active conservation programs are turning toward multispecies landscape and regional conservation actions, and away from single species approaches. This is both a reflection of changing trends in conservation science and advances in foundational technologies, including genomics and geospatial science. Multispecies approaches may provide more fundamental insights into evolutionary processes and equip managers with a more holistic understanding of the landscapes under their jurisdiction. Central to this approach are data generation and analyses which embrace and reflect a broad range of taxonomic diversity. Here, we examine the family-level phylogenetic breadth of the California Conservation Genomics Project (CCGP) based on family-level phylogenetic diversity (PD), family-level phylogenetic distinctness, and family richness. We place this in the context of the diversity present in California and compare it to the 35-plus years of genetic research compiled in the CaliPopGen Database. We found that the family-level PD in the CCGP reflected that of California very well, slightly overrepresenting chordates and underrepresenting arthropods, and that 42% of CCGP PD represented new contributions to genetic data for the state. In one focused effort, the CCGP was able to achieve roughly half the family-level PD studied over the last several decades. To maximize studied PD, future work should focus on arthropods, a conclusion that likely reflects the overall lack of attention to this hyperdiverse clade.

The California ribbed mussel, Mytilus californianus, is an ecosystem engineer crucial for the survival of many marine species inhabiting the intertidal zone of California. Here, we describe the first reference genome for M. californianus and compare it to previously published genomes from three other Mytilus species: M. edulis, M. coruscus, and M. galloprovincialis. The M. californianus reference genome is 1.65 Gb in length, with N50 sequence length of 118 Mb, and an estimated 86.0% complete single copy genes. Compared with the other three Mytilus species, the M. californianus genome assembly is the longest, has the highest N50 value, and the highest percentage complete single copy genes. This high-quality genome assembly provides a foundation for population genetic analyses that will give insight into future conservation work along the coast of California.

The rubber boa, Charina bottae is a semi-fossorial, cold-temperature adapted snake that ranges across the wetter and cooler ecoregions of the California Floristic Province. The rubber boa is 1 of 2 species in the family Boidae native to California and currently has 2 recognized subspecies, the Northern rubber boa C. bottae bottae and the Southern rubber boa C. bottae umbratica. Recent genomic work on C. bottae indicates that these 2 subspecies are collectively composed of 4 divergent lineages that separated during the late Miocene. Analysis of habitat suitability indicates that C. bottae umbratica montane sky-island populations from southern California will lose the majority of their habit over the next 70 yr, and is listed as Threatened under the California Endangered Species Act. Here, we report a new, chromosome-level assembly of C. bottae bottae as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 289 scaffolds covering 1,804,944,895 bp, has a contig N50 of 37.3 Mb, a scaffold N50 of 97 Mb, and BUSCO completeness score of 96.3%, and represents the first reference genome for the Boidae snake family. This genome will enable studies of genetic differentiation and connectivity among C. bottae bottae and C. bottae umbratica populations across California and help manage locally endemic lineages as they confront challenges from human-induced climate warming, droughts, and wildfires across California.

Habitat loss is threatening natural communities worldwide. Small and isolated populations suffer from inbreeding and genetic drift, which jeopardize their long-term survival and adaptive capacities. However, the consequences of habitat loss for reciprocal coevolutionary interactions remain poorly studied. In this study, we investigated the effects of decreasing habitat patch size and connectivity associated with habitat loss on molecular signatures of coevolution in the Alcon blue butterfly (Phengaris alcon) and its most limited host, the marsh gentian (Gentiana pneumonanthe). Because reciprocal coevolution is characterized by negative frequency-dependent selection as a particular type of balancing selection, we investigated how signatures of balancing selection vary along a gradient of patch size and connectivity, using single nucleotide polymorphisms (SNPs). We found that signatures of coevolution were unaffected by patch characteristics in the host plants. On the other hand, more pronounced signatures of coevolution were observed in both spatially isolated and in large Alcon populations, together with pronounced spatial variation in SNPs that are putatively involved in coevolution. These findings suggest that habitat loss can facilitate coevolution in large butterfly populations through limiting swamping of locally beneficial alleles by maladaptive ones. We also found that allelic richness (Ar) of the coevolutionary SNPs is decoupled from neutral Ar in the butterfly, indicating that habitat loss has different effects on coevolutionary as compared with neutral processes. We conclude that this specialized coevolutionary system requires particular conservation interventions aiming at generating a spatial mosaic of both connected and of isolated habitat to maintain coevolutionary dynamics.

Abnormal lipid homeostasis has been observed in the brain of Parkinson's disease (PD) patients and experimental models, although the mechanism underlying this phenomenon is unclear. Notably, previous studies have reported that the PD-linked protein Parkin functionally interacts with important lipid regulators, including Sterol Regulatory Element-Binding Proteins (SREBPs) and cluster of differentiation 36 (CD36). Here, we demonstrate a functional relationship between Parkin and lipoprotein lipase (LPL), a triglyceride lipase that is widely expressed in the brain. Using a human neuroblastoma cell line and a Parkin knockout mouse model, we demonstrate that Parkin expression level positively correlates with neuronal LPL protein level and activity. Importantly, our study identified SREBP2, a major regulator of sterol and fatty acid synthesis, as a potential mediator between Parkin and LPL. Supporting this, SREBP2 genetic ablation abolished Parkin effect on LPL expression. We further demonstrate that Parkin-LPL pathway regulates the formation of intracellular lipid droplets, and that this pathway is upregulated upon exposure to PD-linked oxidative stress induced by rotenone. Finally, we show that inhibition of either LPL or SREBP2 exacerbates rotenone-induced cell death. Taken together, our findings reveal a novel pathway linking Parkin, SREBP2 and LPL in neuronal lipid homeostasis that may be relevant to the pathogenesis of PD. Ministry of Education (MOE) Ministry of Health (MOH) National Medical Research Council (NMRC) Published version Singapore Ministry of Health Open Fund-Large Collaborative Grant (MOH-OFLCG18May-0002), Lee Kong Chian School of Medicine Start Up Grant from Singapore Ministry of Education (#002642-00001), and National Medical Research CouncilTranslational and Clinical Research (NMRC/TCR/013-NNI/2014) awarded to K.L.L. and Khoo Pilot Award (Duke/Duke-NUS/RECA (Pilot)/2016/0023) awarded to K.L.L. and T.P.Y.

The mitochondrial DNA mutation m.9032T>C was previously identified in patients presenting with NARP (Neuropathy Ataxia Retinitis Pigmentosa). Their clinical features had a maternal transmission and patient’s cells showed a reduced oxidative phosphorylation capacity, elevated reactive oxygen species (ROS) production and hyperpolarization of the mitochondrial inner membrane, providing evidence that m.9032T>C is truly pathogenic. This mutation leads to replacement of a highly conserved leucine residue with proline at position 169 of ATP synthase subunit a (L169P). This protein and a ring of identical c-subunits (c-ring) move protons through the mitochondrial inner membrane coupled to ATP synthesis. We herein investigated the consequences of m.9032T>C on ATP synthase in a strain of Saccharomyces cerevisiae with an equivalent mutation (L186P). The mutant enzyme assembled correctly but was mostly inactive as evidenced by a > 95% drop in the rate of mitochondrial ATP synthesis and absence of significant ATP-driven proton pumping across the mitochondrial membrane. Intragenic suppressors selected from L186P yeast restoring ATP synthase function to varying degrees (30–70%) were identified at the original mutation site (L186S) or in another position of the subunit a (H114Q, I118T). In light of atomic structures of yeast ATP synthase recently described, we conclude from these results that m.9032T>C disrupts proton conduction between the external side of the membrane and the c-ring, and that H114Q and I118T enable protons to access the c-ring through a modified pathway.

Beta amyloid cleaving enzyme 1 (BACE1) is largely expressed by neurons and is the sole β-secretase for initiating the production of neuronal β-amyloid peptides (Aβ). To fully understand the physiological functions of neuronal BACE1, we used mouse genetic approach coupled with unbiased single nucleus RNA sequencing (snRNAseq) to investigate how targeted deletion of Bace1 in neurons, driven by Thy-1-Cre recombinase, would affect functions in the nervous system. Our transcriptome results revealed that BACE1 is essential for maturation of neural precursor cells and oligodendrocytes in mice. RNA velocity analysis confirmed deficit in the trajectory of neuroblasts in reaching the immature granule neuron state in young Bace1fl/fl; Thy1-cre mice. Further analysis of differential gene expression indicated changes in genes important for SNARE signaling, tight junction signaling, synaptogenesis and insulin secretion pathways. Morphological studies revealed a hypomyelination in Bace1fl/fl;Thy1-cre sciatic nerves, but no detectable myelination changes in the corpus callosum, despite clear reduction in myelination proteins in the brain. Functional studies showed reduction in long-term potential, defects in synaptogenesis and learning behavioral. Altogether, our results show that neuronal BACE1 is critical for optimal development of central and peripheral nervous system, and inhibition of neuronal BACE1 will result in deficits in synaptic functions and cognitive behaviors.

ADP-ribosylation factor 1 (ARF1) is a small GTPase that regulates membrane traffic at the Golgi apparatus and endosomes through recruitment of several coat proteins and lipid-modifying enzymes. Here, we report a pediatric patient with an ARF1-related disorder because of a monoallelic de novo missense variant (c.296 G > A; p.R99H) in the ARF1 gene, associated with developmental delay, hypotonia, intellectual disability and motor stereotypies. Neuroimaging revealed a hypoplastic corpus callosum and subcortical white matter abnormalities. Notably, this patient did not exhibit periventricular heterotopias previously observed in other patients with ARF1 variants (including p.R99H). Functional analysis of the R99H-ARF1 variant protein revealed that it was expressed at normal levels and properly localized to the Golgi apparatus; however, the expression of this variant caused swelling of the Golgi apparatus, increased the recruitment of coat proteins such as coat protein complex I, adaptor protein complex 1 and GGA3 and altered the morphology of recycling endosomes. In addition, we observed that the expression of R99H-ARF1 prevented dispersal of the Golgi apparatus by the ARF1-inhibitor brefeldin A. Finally, protein interaction analyses showed that R99H-ARF1 bound more tightly to the ARF1-effector GGA3 relative to wild-type ARF1. These properties were similar to those of the well-characterized constitutively active Q71L-ARF1 mutant, indicating that the pathogenetic mechanism of the R99H-ARF1 variant involves constitutive activation with resultant Golgi and endosomal alterations. The absence of periventricular nodular heterotopias in this R99H-ARF1 subject also indicates that this finding may not be a consistent phenotypic expression of all ARF1-related disorders.