Antisense oligonucleotide (ASO)-mediated exon skipping has become a valuable tool for investigating gene function and developing gene therapy. Machine learning-based computational methods such as eSkip-Finder have been developed to predict the efficacy of ASOs via exon skipping. However, these methods are computationally demanding, and the accuracy of predictions remains suboptimal. In this study, we propose a new approach to reduce computational burden and improve prediction performance by using feature selection within machine learning algorithms and ensemble learning techniques. We evaluated our approach using a dataset of experimentally validated exon skipping events, dividing it into training and testing sets. Our results demonstrate that using a 3-way voting approach with random forest, gradient boosting, and XGBoost can significantly reduce computation time to under ten seconds while improving prediction performance, as measured by R2 for both 2’-O-methyl nucleotides (2OMe) and phosphorodiamidate morpholino oligomers (PMOs). Additionally, the feature importance ranking derived from our approach is in good agreement with previously published results. Our findings suggest that our approach has the potential to enhance the accuracy and efficiency of predicting ASO efficacy via exon skipping. It could also facilitate the development of novel therapeutic strategies. This study could contribute to the ongoing efforts to improve ASO design and optimize gene therapy approaches.

Differential RNA editing by adenosine deaminases that act on RNA (ADARs) has been implicat-ed in several neurological disorders, including Parkinson’s disease (PD). Here, we report results of an RNAi screen of genes differentially regulated in adr-2 mutants, normally encoding the only catalytically active ADAR in Caenorhabditis elegans, ADR-2. Subsequent analysis of candidate genes that alter the misfolding of human α-synuclein (α-syn) and dopaminergic neurodegenera-tion, two PD pathologies, reveal that reduced expression of xdh-1, the ortholog of human xanthine dehydrogenase (XDH), is protective against α-synuclein-induced dopaminergic neurodegenera-tion. Further RNAi experiments show that WHT-2, the worm ortholog of the human ABCG2 transporter and a predicted interactor of XDH-1, is the rate-limiting factor in the ADR-2, XDH-1, WHT-2 system for dopaminergic neuroprotection. Proteomic analysis indicates that the editing of one nucleotide in wht-2 RNA leads to the substitution of threonine with alanine at residue 124 in the WHT-2 protein, changing its structure. Thus, we propose a model where wht-2, is edited by ADR-2 which promotes optimal export of uric acid, a known substrate of WHT-2 and a product of XDH-1 activity. In the absence of editing, uric acid export is limited provoking a reduction in xdh-1 transcription to limit uric acid production and maintain cellular homeostasis. In turn, ele-vation of uric acid is protective against dopaminergic neuronal cell death. These data indicate that modifying specific targets of RNA editing may represent a promising therapeutic strategy for PD.

Selectins are responsible for the early stages of cell migration as they control cell adhesion. They make the microenvironment permissive for metastatic events by promoting the activation of other cell adhesion molecules (CAMs). We employed several robust bioinformatics tools to evaluate gene sequence; single nucleotide polymorphisms (SNPs) in intronic and UTR regions; and missense SNPs with amino acid change in L-selectin (SELL), E-selectin (SELE), P-selectin (SELP) and PSGL-1 (SELPLG). We demonstrated that gene polymorphisms rs2229569, rs1131498, rs4987360, rs4987301 and rs2205849; polymorphisms rs3917777, rs2205894 and rs2205893 of SELP gene; rs7138370, rs7300972 and rs2228315 variants of SELPLG gene; and rs1534904 and rs5368 polymorphisms of SELE gene may produce important alterations in the DNA structure and consequent alterations in the morphology and function of the corresponding proteins. These Selectin polymorphic variants deserve further investigation in cancer patients, as they may become useful clinical markers for risk determination, diagnostic and prognostic biomarkers or even targets for targeted therapies.

Spinal muscular atrophy 5q (SMA) is the most common neuromuscular inherited disease in neonates which is associated with homozygous deletion of exon 7 in the SMN1 gene. Recently established drugs can improve the motor functions of SMA infants when treated in the pre-symptomatic stage. With aim of providing an early diagnosis, newborn screening (NBS) for SMA using real-time PCR assay with dried blood spots (DBS) was performed in Saint-Petersburg from January 2022 through November 2022. Here, 36,140 newborns were screened by GenomeX real-time PCR-based screening test, and three genotypes were determined. Homozygous deletion carriers (4 newborns), heterozygous carriers (772 newborns) and wild-type individuals (35,364 newborns) were identified. All four newborns screened positive for homozygous SMN1 deletion were confirmed by alternate methods. Two of the newborns had two copies of SMN2, and two of the newborns had three copies. We determined Saint-Petersburg spinal muscular atrophy incidence of 1 in 9,035 and SMA carrier frequency of 1 in 47. In conclusion, it can be summarized that providing both timely SMN1 information and confirmation along with SMN2 copy number as part of SMA newborn screening algorithm can significantly improve clinical follow-up, family members testing, and SMA patients' treatment.

Deep learning has already revolutionised the way we process a wide range of data, in many areas of our daily life. The ability to learn abstractions and relationships from heterogeneous data, has provided impressively accurate prediction and classification tools to handle increasingly big datasets. This has a significant impact on the growing wealth of omics datasets, with the unprecedented opportunity for a better understanding of the complexity of living organisms. While this revolution is transforming the way we analyse these data, explainable deep learning is emerging as an additional tool with the potential to change the way we interpret biological data. Explainability addresses critical issues such as transparency, so important when computational tools are introduced especially in clinical environments. Moreover, it empowers artificial intelligence with the capability to provide new insights in the input data, thus adding an element of discovery to these already powerful resources. In this review we provide an overview of the transformative effects explainable deep learning is having on multiple sectors, ranging from genome engineering and genomics, from radiomics to drug design and clinical trials. We offer a perspective to life scientists, to better understand the potential of these tools, and a motivation to implement them in their research, by suggesting learning resources they can use to move their first steps in this field.

Chicory, a horticultural crop cultivated worldwide, presents many botanical varieties and local biotypes. Among these, Italian radicchio group cultivars of pure Cichorium intybus L. and interspecific hybrids with Cichorium endivia L., as the “Red of Chioggia” biotype which in turn includes several phenotypes. This study deals with a pipeline for the marker-assisted breeding of F1 hybrids: it presents the genotyping-by-sequencing results of four elite inbred lines using a RADseq approach and an original molecular assay based on CAPS markers for screening mutants with nuclear male-sterility in radicchio of Chioggia. Two thousand nine hundred fifty-three SNP-carrying RADtags were identified and used to compute the actual estimates of homozygosity, and overall genetic similarity and uniformity of the populations, as well as to determine their genetic distinctiveness and differentiation. Molecular data were further used to investigate the genomic distribution of the RADtags among the two Cichorium species, allowing their mapping in 1,131 and 1,071 coding sequences in chicory and endive, respectively. In parallel, an assay to screen the genotype at the male-sterility locus Cims-1 was developed to discriminate wild-type and mutant alleles of the causative gene myb80-like. Moreover, a RADtag mapped close to this genomic region proved the potential application of this method for future marker-assisted selection tools. Finally, combining the genotype information of the core collection, the best ten individuals from each inbred line were selected to compute the observed genetic similarity as a measure of uniformity as well as the expected homozygosity and heterozygosity estimates scorable by the putative progenies derived from selfing (pollen parent) and full-sibling (seed parent) or pair-wise crossing (F1 hybrids). This predictive approach was conducted as a pilot study for understanding the potential application of RADseq in improving molecular breeding strategies aimed at the development of inbred lines and F1 hybrids in leaf chicory.

Lupinus mutabilis Sweet (Fabaceae), “tarwi” or “chocho”, is an important grain legume in the Andean region. In Peru, studies on tarwi have been mainly focused on morphological features, however, the have not been molecularly characterized. Currently, it is possible to explore genetic parameters of plants with reliable and modern methods like genotyping-by-sequencing (GBS). We here for the first time used single nucleotide polymorphisms (SNPs) markers to infer the genetic diversity and population structure of 89 accessions of tarwi from nine Andean regions of Peru. A total of 5922 SNPs distributed along all chromosomes of tarwi were identified. STRUCTURE analysis revealed that this crop is grouped into two clusters. A dendrogram was generated using the UPGMA clustering algorithm and, similar to the principal coordinate analysis (PCoA), it showed two groups that correspond to the geographic origin of the tarwi samples. AMOVA showed a reduced variation between clusters (7.59 %) and indicated that variability within populations is 92.41 %. Population divergence (Fst) between clusters 1 and 2 revealed low genetic difference (0.019). We also detected a negative Fis for both populations, demonstrating that, similar to other Lupinus species, tarwi also depends on cross-pollination. SNPs markers were powerful and effective for the genotyping process in this germplasm. We hope that this information is the beginning of the path towards a modern genetic improvement and conservation strategies of this important Andean legume.

Prostate cancer (PC) is the most frequently diagnosed non-skin cancer in the world. Previous studies showed that genomic alterations represent the most common mechanism for molecular alterations that cause the development and progression of PC. Great efforts have been done to identify common protein-coding genetic variations; however, the impact of non-coding variations including regulatory genetic variants is not still well understood. To gain an understanding of the functional impact of genetic variants, particularly, regulatory variants in PC, we developed an integrative pipeline (AGV) that used whole genome/exome sequences, GWAS SNPs, chromosome conformation capture data, and ChIP-Seq signals to investigate the potential impact of genomic variants on the underlying target genes in PC. We identified 646 putative regulatory variants, of which 30 of them significantly altered the expression of at least one protein-coding gene. Our analysis of chromatin interactions data (Hi-C) revealed that the 30 putative regulatory variants may affect 131 coding and non-coding genes. Interestingly, our study showed the 131 protein-coding genes are involved in disease-related pathways including Reactome and MSigDB in which for most of them targeted treatment options are currently available. Together, our results provide a comprehensive map of genomic variants in PC and revealed their potential contribution to prostate cancer progression and development.

Milk is a complex liquid, and the concentrations of many of its components are under genetic control. Many genes and pathways are known to regulate milk compositon, and the purpose of this review is to highlight how the discoveries of quantitative trait loci (QTL) for milk phenotypes can elucidate these pathways. The main body of this review focusses primarily on QTL discovered in cattle (Bos taurus) as a model species for lactation biology, with occasional references to other ruminant species, and some comparisons with human milk composition are also presented. The following section describes a range of techniques that can be used to help identify the causative genes underlying QTL when the underlying mechanism involves the regulation of gene expression. As genotype and phenotype data bases continue to grow and diversify, new QTL will continue to be discovered, and, although proving the causality of underlying genes and variants remains difficult, these new data sets will further enhance our understanding of lactation biology.

AIM: The study aims to identify if mutations in the SLC4A11 gene are present in Filipino families affected with Congenital Hereditary Endothelial Dystrophy (CHED). METHODS: This is a family cohort study that investigated the genetic profile of a selected family in Northern Luzon, Philippines, whose members were diagnosed with Congenital Hereditary Endothelial Dystrophy (CHED). A patient who was diagnosed with CHED prior to this study, served as the proband for this family. A detailed family history and a complete ophthalmologic examination were done to each of the family members. A total of six affected members and three unaffected members were included in this study. DNA was isolated from peripheral blood samples of family members, and polymerase chain reaction (PCR) was used to amplify the gene’s entire coding region, (19 exons and 2 putative promoter regions), and finally, the amplified regions were analyzed using DNA sequencing. RESULTS: Consanguinity was not present in the family. Corneal haze was reported to be present since birth or shortly thereafter in all affected patients. Slit-lamp examination showed edematous corneas. Molecular studies of the SLC4A11 gene revealed four novel homozygous point mutations variably present in six affected members as well as three unaffected members. One unaffected family member (I-1) had a novel sense mutation absent in other family members. All affected siblings showed little phenotypic variability. CONCLUSION: This is the first report that gives us a genetic profile of a Northern Luzon family with members affected by CHED. This study supports earlier findings that the mutations in the SLC4A11 gene are not consistently the same among different ethnic groups worldwide, probably due to the disease’s genetic heterogeneity. Our study documented five novel mutations adding to the growing list of mutations probably responsible for acquiring the CHED phenotype. It is possible that there are more novel mutations waiting to be discovered in this hereditary disease. Screening for these specific mutations in other families may prove useful for genetic counseling, prenatal diagnosis and the future development of gene therapy.

Fanconi anemia (FA) patients display an exacerbated risk of oral squamous cell carcinoma (OSCC) and precursor lesions at young ages, mainly at the oral cavity. As patients have defects in DNA repair mechanisms, standard-of-care treatments to OSCC such as radiotherapy and chemotherapy give rise to severe toxicities. New methods for early diagnosis are urgently necessary to allow treatments in early disease stages and achieve better clinical outcomes. We have conducted a prospective, longitudinal study whereby liquid biopsies from sixteen lesion/tumor-free patients were analyzed for the presence of mutations in cancer genes. DNA from saliva and plasma were sequentially collected and deep-sequenced, and the clinical evolution followed during a median time of around 2 years. In 9/16 FA patients we detected mutations in cancer genes (mainly TP53) with molecular allele frequencies (MAF) down to 0.07 %. Importantly, all patients having mutations and clinical follow-up data after mutation detection (n=6) developed oral precursor lesions or OSCC. Lead-time between mutation detection and tumor diagnosis ranged from 23 to 630 days. Strikingly, FA patients without mutations display significantly lower risk of developing precursor lesions or OSCC. Therefore, our diagnostic approach could help to stratify FA patients into risk groups, which would allow closer surveillance for OSCC or precursor lesions.

The ZRS enhancer, a regulatory sequence found in numerous organisms, plays an important role in early embryonic limb development. ZRS controls the expression of the Sonic Hedgehog gene (Shh) and, therefore, early limb development in an organism. Shh has been shown to control the width of the limb bud by stimulating mesenchyme cell proliferation due to its ability to regulate the anterior-posterior length of the apical ectodermal ridge. Several transcription factors, acting as repressors or activators of the Shh gene, coordinate this limb development process in tandem with the ZRS enhancer. While the significance of normal ZRS activity is evident, this study looks deeper into the effects of pathogenic changes to the ZRS enhancer and the development of associated limb disorders such as preaxial polydactyly (PPD) by focusing on several aspects of ZRS regulation and its relation to Shh expression. This was accomplished by characterizing the expression of Shh and mCherry, an introduced luminescence gene regulated by ZRS, through single-cell RNA-sequenced cells from a developing limb bud of a mouse embryo. Additionally, this study characterized specific transcription factors as potential repressors or activators of ZRS by determining TF enrichment or depletion in highly expressive Shh and mCherry cells through novel deep-learning algorithms. Classifying such TFs is vital in identifying the regulatory elements that control the formation of limb malformation disorders such as preaxial polydactyly which can aid in the development of future therapeutic interventions.

Prostate cancer (PCa) is one of the most prevalent cancers among men in India. Although studies on PCa have dealt with the genetics, genomics, and the environmental influence in causality of PCa, not many studies employing the next generation sequencing (NGS) approaches of PCa have been carried out. In our previous study, we have identified some causal genes and mutations specific to Indian PCa using Whole-Exome Sequencing (WES). In the recent past, with the help of different cancer consortiums such as The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), along with differentially expressed genes (DEGs), many cancer-associated novel non-coding RNAs have been identified as biomarkers. In this work, we attempt to identify DEGs as well as long non-coding RNAs (lncRNAs) associated with signature pathways from an Indian PCa cohort using RNA-sequencing (RNA-seq) approach. From a cohort of 60, we screened 6 patients who underwent prostatectomy; we performed a whole transcriptome shotgun sequencing (WTSS)/RNA-sequencing to decipher the DEGs. We further normalized the read counts using fragments per kilobase of transcript per million mapped reads (FPKM) and analyzed the DEGs using a cohort of downstream regulatory tools, viz. GeneMANIA Stringdb, Cytoscape-Cytohubba, cbioportal to map the inherent signatures associated with PCa. By comparing the RNA-seq data obtained from the pairs of normal and PCa tissue samples using our benchmarked in-house cuffdiff pipeline, we observed some important genes specific to PCa such as STEAP2, APP, PMEPA1, PABPC1, NFE2L2, HN1L and some other important genes known to be involved in different cancer pathways such as, COL6A1, DOK5, STX6, BCAS1, BACE1, BACE2, LMOD1, SNX9, CTNND1 etc. We also identified a few novel lncRNAs such as LINC01440, SOX2OT, ENSG00000232855, and ENST00000647843.1 that need to be characterized further. Deregulation of SOX2OT is observed in various tumors, including lung cancer, gastric cancer, esophageal cancer, breast cancer, hepatocellular carcinoma, ovarian cancer, pancreatic, laryngeal squamous cell carcinoma, osteosarcoma, nasopharyngeal carcinoma, and glioblastoma. It would be interesting to characterize its function in PCa as well. In comparison with publicly available datasets, we have identified characteristic DEGs and novel lncRNAs implicated in signature PCa pathways in an Indian PCa cohort which have perhaps not been reported. As a pilot study, this has set a precedent for us to validate further experimentally, and we firmly believe this will pave a way towards discovery of biomarkers and development of novel therapies.

Background: High blood pressure (BP) has been implicated as a major risk factor for cardiovascular diseases in several global populations, including in individuals of African ancestry. Despite the elevated burden of high BP-induced cardiovascular diseases in Africa and other global populations with African ancestry, limited genetic studies have been carried out to explore the genetic machinery driving this phenomenon. Methods: We performed univariate and multivariate analyses using Genome-wide association studies (GWAS) and summary statistics data of 77,850 individuals of African ancestry for systolic (SBP) and diastolic blood pressure (DBP) traits. The six independent cohorts used included individuals derived from the African Partnership for Chronic Disease Research (APCDR), the UK Biobank, and the Million Veteran Program (MVP). Subsequently, we annotated, prioritized, visualized, and interpreted our meta-analyses results using FUMA, to gain further insight into the molecular mechanism(s) that contribute to the genetics of BP traits. Finally, loci attaining genome-wide significance, GWS (p<5x10-8) were also followed up with Bayesian fine-mapping to identify potential causal variants. Results: Our meta-analyses altogether identified 350 GWAS SNPs for SBP (166 SNPs) and DBP (184 SNPs, including two novel loci) whilst our multivariate GWAS method identified 166 SNPs (including three novel loci). Interestingly, in FUMA there was significant tissue enrichment of up-regulated differentially expressed genes (DEGs) in the sigmoid and transverse colon for SBP, as well as 10 significant gene sets from MAGMA gene set analyses, However, for DBP, no significant DEGs nor gene sets in MAGMA were found; instead, in DBP for gene property analysis for tissue specificity nine candidates were found to be significant and all nine were in different brain regions. Finally, Bayesian fine-mapping revealed that only 11 variants from the lead SNPs had >50% posterior probability (PP) of being causal and they included novel variant rs562545 (MOBP, PP = 77%) and 10 other previously published variants. Conclusion: Our results demonstrate the importance of performing GWAS in large sample sizes of global populations of African ancestry, including continental Africans; which yield novel insights, from novel loci to novel pathways/tissue expression candidates. Large-scale genomic datasets are required to enhance further discovery and fine-mapping of high-risk loci/variants in highly susceptible groups for cardiovascular disease and other related traits. Our study highlights the need for diversity in genetic research and the importance of expanding large GWASs to include ancestrally diverse populations.

(1) Background: Obesity and its comorbidities can cause burdens and limitations. Bariatric surgery (BS) is indicated as a safe procedure to reduce body mass and improve present comorbidities. Yet, several complications were reported, like vitamin D [25(OH)D] deficiency. We evaluated if 25(OH)D serum levels relate to clinical characteristics, symptoms, or habits in women after their BS and whether the vitamin D receptor (VDR) gene's TaqI and FokI polymorphisms affected 25(OH)D levels and the total body bone mineral density (TBBMD). (2) Methods: This cohort cross-sectional comparative analytical prospective study consisted of 27 women, 61.6 ± 5.0 years, submitted to BS one year prior at a public reference hospital, DF-Brazil. All participants were asked to follow the physical and dietary activity recommendations and received vitamin D3 supplements. Their anthropometric, biochemical, and immunological measurements and blood samples were obtained. (3) Results: 73.3% of participants had low 25(OH)D levels, and their levels correlated positively with TBBMD and negatively with systolic pressure. VDR TaqI did not affect 25(OH)D levels, whereas VDR FokI's allele f presence correlated to a median rise in 25(OH)D levels. Neither polymorphism correlated to TBBMD. (4) Conclusions: 25(OH)D levels were positively correlated with TBBMD, negatively with systolic blood pressure, and were higher in those with the VDR FokI allele f.

Inteins, often referred to as protein introns, are highly mobile genetic elements that invade conserved genes throughout the tree of life. Inteins have been found to invade a wide variety of key genes within actinophages. While in the process of conducting a survey of these inteins in actinophages we discovered that one protein family of methylases contained a putative intein, and two other unique insertion elements. These methylases are known to occur commonly in phages as orphan methylases (possibly as a form of resistance to restriction-modification systems). We found that the methylase family is not conserved within phage clusters and has a disparate distribution across divergent phage groups. We determined that two of the three insertion elements have a patchy distribution within the methylase protein family. We also found that the third insertion element is likely a second homing endonuclease, and that all three elements (the intein, the homing endonuclease, and what we refer to as the ShiLan domain) all have different insertion sites that are conserved in the methylase gene family. Furthermore, we find strong evidence that both the intein and ShiLan domain are partaking in long distance horizontal gene transfer events between divergent methylases in disparate phage hosts within the already dispersed methylase distribution. The reticulate evolutionary history of methylases and their insertion elements reveals high rates of gene transfer and within-gene recombination in actinophages.

Genotype imputation is widely used to enrich genetic datasets. The operation relies of panels of known reference haplotypes with typically whole-genome sequencing data. How to choose a reference panel has been widely studied and it is essential to have a panel that is well matched to the individuals who require imputation of missing genotypes. However, it is broadly accepted that such an imputation panel will have an enhanced performance with the inclusion of diversity; haplotypes from many different populations. We investigate this observation in this work by examining in fine detail exactly which reference haplotypes are contributing at different regions of the genome. This is achieved using a novel method of inserting synthetic genetic variation into the reference panel in order to track the performance of leading imputation algorithms. We show that while diversity may globally improve imputation accuracy, there can be occasions where incorrect genotypes are imputed following the inclusion of more diverse haplotypes in the reference panel. We however demonstrate a technique for retaining and benefitting from the diversity in the reference panel whilst avoiding the occasional adverse effects on imputation accuracy. What is more, our results elucidate more clearly the role of the diversity in a reference panel than has been shown in previous studies.

This article is devoted to results of in-depth analysis of the system of binary-oppositional structures in DNA n-plet alphabets and their algebraic-matrix representations. These results show that the molecular complementary replication of DNA strands is accompanied by the presence of an algebraic version of the principle "like begets like" in matrix representations of DNA alphabets having internal structures. This algebraic version is based on binary-oppositional structures in the genetic molecular system, which can be represented by bynary numbers and corresponding matrices of DNA alphabets. The received results allow thinking that the phenomenon "like begets like" (or a complementary replication in a wide sense) is a systemic in the genetic organization and is connected with algebraic features of biological organization. Correspondingly, the biological principle "like begets like" can be additionally modeled by algebraic-matrix methods and approaches. Such algebraic-matrix modeling of the genetic coding system gives new ways for studying and understanding a key role of the named principle in genetic and other inherited physiological complexes. The author believes that further study of the algebraic relationships of the genetic system and inherited physiological complexes will be increasingly revealing the key biological role of the ancient principle "like begets like" at different levels of biological organization.

Meiotic crossovers/chiasmata are not randomly distributed and strictly controlled. The mechanisms behind COs patterning remain largely unknown. In many species COs predominantly occur in the distal 2/3 of the chromosome arm and suppression of COs is observed in the proximal regions. The exceptions are some species, among which Allium fistulosum has strictly localized COs in the proximal region. The ability to manipulate the COs localization can be useful for onion breeding. We investigated the factors that may contribute to the pattern of COs in two closely related onion species A. cepa and A. fistulosum, which differ significantly in the localization of chiasmata, and their F1 diploid and triploid hybrids in pollen mother cells. We demonstrate a significant shift in the COs localization to the distal and interstitial region in F1 triploid hybrid, which has a complete diploid set of A. fistulosum chromosomes and haploid set of A. cepa chromosomes. This observation points to a possible genetic control of COs distribution. We did not find the differences in the assembly and disassembly of ASY1 and ZYP1 between A. cepa and A. fistulosum while the difference between parental species and their hybrids was observed. In diploid F1 hybrids at pachytene the chromosome pairing delay marked by ASY1 was revealed. Immunolocalization of MLH1, a marker for class I of the interference-dependent COs, and MUS81, a marker for class II of the interference-free COs, and mlh1/mus81 gene expression profiling in different stages of meiosis showed a spatiotemporal asymmetry among A. cepa, A. fistulosum and their diploid and triploid F1 hybrids. Our results support the hypothesis of genetic control of CO distribution as one of the players that affect meiotic recombination and exchange of genomic material.

A recent paradigm shift in diagnostics of medulloblastoma allows the distinction of four major groups defined by genetic data rather than histology. This new molecular classification correlates better with prognosis and will allow better clinical management for therapies targeting druggable mutations, but also offers a new combination of monitoring tumor development in real-time and treatment response by sequential liquid biopsy. This review highlights recent developments after a century of milestones in neurosurgery, radio- and chemotherapy, but also controversial theories on the cell of origin, animal models and the use of liquid biopsy.

High resolution melting (HRM) analysis is a powerful detection method for fast, high-throughput post-PCR analysis. A two-step HRM marker system was developed for identification of the N-, S-, R- and T-cytoplasms of onion. In the first step for identification of N-, S-, and R-cytoplasms, one forward primer was designed to the identical sequences of both cox1 and orf725 genes and two reverse primers specific to the polymorphic sequences of cox1 and orf725 genes were used. For the second step breeding lines with N-cytoplasm were evaluated with primers developed from the orfA501 sequence to distinguish between N- and T-cytoplasms. An amplicon with primers to the mitocondrial atp9 gene was used as an internal control. The two-step HRM marker system was tested using 246 onion plants. HRM analysis showed that the most common source of CMS, often used by Russian breeders, is S-cytoplasm, the rarest type of CMS is R-cytoplasm, and the proportion of T-cytoplasm among the analyzed breeding lines was 20.5%.

The hexaploid bread wheat Triticum aestivum L. (2n = 6x = 42) and the tetraploid durum wheat Triticum durum (2n = 4x = 28) are two commercially the most important wheat species worldwide. These two species are components of the same gene pool, which evolved under different conditions, which made it possible to acquire different adaptations to a wider range of environmental conditions. Identification and characterization of introgressed fragments of one species in the genome of another can help in the search for new loci responsible for economically valuable traits. SNP microarrays provide information on the genotypes of thousands of positions throughout the genome. This information can be used to obtain the distribution of alleles of one species in the genome of another and to determine the location of introgressed fragments. In our work, we used genotyping data from the SNP microarray to identify durum wheat chromosome fragments containing bread wheat alleles in the offspring from crossing bread and durum wheat. We have studied the distribution of the genetic material of bread wheat in 15 combinations from crosses of bread and durum wheat. Introgressed fragments of bread wheat into the durum wheat genome were identified on chromosomes 1A, 1B, 4A, 5B, 6B, and 7A. The introgressed fragments showed polymorphism in different breeding lines. A functional enrichment analysis of genes in the introgression fragments was performed, which showed the presence of statistically significant enrichment of genes responsible for various molecular functions, biological processes and cellular components. The obtained data can be used to identify new QTL of economically valuable traits for obtaining durum wheat varieties with improved agronomic traits.

Background: Tuberculosis (TB) manifests itself primarily in the lungs as pulmonary disease (PTB) and sometimes disseminates to other organs to cause extra-pulmonary TB, such as lymph node TB (LNTB). This study aimed to investigate the role of host genetic polymorphism in immunity related genes to find a genetic basis for such differences. Methods: Sixty-three, Single nucleotide polymorphisms (SNPs) in twenty-three, TB-immunity related genes including eleven innate immunity (SLCA11, VDR, TLR2, TLR4, TLR8, IRGM, P2RX7, LTA4H, SP110, DCSIGN and NOS2A) and twelve cytokine (TNFA, IFNG, IL2, Il12, IL18, IL1B, IL10, IL6, IL4, IL1RA, IL8 and TNFB) genes were investigated to find genetic associations in both PTB and LNTB as compared to healthy community controls. The serum cytokine levels were correlated for association with the genotypes. Results: PTB and LNTB showed differential genetic associations. The genetic variants in the cytokine genes (IFNG, IL12, IL4, TNFB and IL1RA and TLR2,4 associated with PTB susceptibility and cytokine levels but not LNTB (p &lt; 0.05). Similarly, genetic variants in LTA4H, P2RX7, DCSIGN and SP110 showed susceptibility to LNTB and not PTB. Pathway analysis showed abundance of cytokine related variants for PTB and apoptosis related variants for LNTB. Conclusions: PTB and LNTB outcomes of TB infection have a genetic component and should be considered for any future susceptibility and functional studies.

We attempted to understand the cervical cancer patient samples through Whole Exome Sequencing. We derived the variants from raw reads via our in-house benchmarked pipeline and validated the variants by IGV. This is the first cervical cancer exome data from the Indian cohort. Background: Cervical cancer (CC) is caused mainly by persistent infections of high-risk HPV, reduced parity, and factors like a decrease in average socioeconomic levels. We keep this because no cervical cancer exome data is available from the Indian cohort. Methods: The CC patient clinical samples were initially subjected to preparation using Qiagen DNA extraction, quantified, and the library preparation using the Agilent target enrichment system. Further, exome capture by Illumina platform (100X), quality check, alignment, and variant calling followed by the downstream analysis and finally visualized by IGV. Results: We observed a large number of SNVs or mutations from an Indian perspective, such as KMT2C, OR4M1, PDPR2P, EPHB1, FAS, OPCML, MGST1, C1QTNF9, HS6ST3, OR4K2, PRPSAP2, KCNJ12, FIGNL1, SFXN1, BAGE2, ARVCF, NAMPT variants of significance, unknown significance and possibly significant are reported. Conclusion: For the first time, KMT2C is observed as a novel potent mutation and pathogenic, showing the variant position at (7:152265091, T>A, SNV 62478356) from the Indian context in CC. Further, we visualized and validated the mutations using the Integrative Genomic Viewer (IGV) browser. Finally, we discuss the inherent challenges through KMT2C mutations.

The molecular mechanisms underlying heat stress tolerance in animals to high temperatures remain unclear. This experiment was performed with Sprague-Dawley rats housed at 22 ℃ (control group; CT), and 42 ℃ for 30 min (H30), 60 min (H60), and 120 min (H120). Transcriptomic analysis using RNA-Sequencing (RNA-Seq) from blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120) was performed. The differentially expressed mRNA isoforms (DEIs) were identified and annotated by the CLC Genomics Workbench. Biological process and metabolic pathway analyses were performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 225, 5,764, and 4,988 DEIs in the blood, liver, and adrenal glands were observed. Furthermore, the number of novel differentially expressed transcript lengths with annotated genes and the novel differentially expressed transcript with non-annotated genes were 136 and 8 in blood, 3,549 and 120 in the liver, as well as 3,078 and 220 in adrenal glands, respectively. A total of 35 genes were involved in the heat stress response, out of which Dnaja1, LOC680121, Chordc1, AABR07011951.1, Hsp90aa1, Hspa1b, Cdkn1a, Hmox1, Bag3 and Dnaja4 were commonly identified in the liver and adrenal glands, suggesting that these genes may regulate heat stress response through interactions between liver and adrenal glands. The results suggest that the identified mRNA isoforms could be considered as potential candidates for selecting mammals to improve thermotolerance.

Here I describe an overlooked form of non-biological paralife (i.e., near-life) that has been evolving on Earth for millions of years, and is currently in the final stages of transitioning into a new form of life. Any consideration of non-biological life or paralife is complicated by the fact that there is no consensus among biologists for the definition of life. This ambiguity has caused disagreement about whether subcellular reproduction systems like viruses are a form of life, despite having genomes, mutations, heritable phenotypes and system-improving evolution. To resolve this problem, I develop a definition of life that is entirely functional and independent of any of the structural idiosyncrasies of biological life on Earth: an order-generating system controlled by internally-encoded information that perpetuates itself by functioning to counteract its entropic decay. Using this definition, subcellular transposons, plasmids, and viruses are paralife because they match the definition of life in all ways except that they induce their order-generating functioning by a living host rather through their own self-sustaining production system. Using this functional definition of life, I show that utility- products (UPs) like fabricated hand tools are part of induced-reproduction systems that have features equivalent to biological genomes, mutations, heritable phenotypes, and a process of system-improving evolution. The perceived benefit of utility-products causes them to induce their reproduction by a biological life-form (humans). For these reasons, human utility products are functionally just as close to being a form of life as subcellular transposons, plasmids, and viruses, i.e., they are Utility-Product paralife (UP-parlife). I also show that some forms of UP-paralife are currently evolving into mechanical life that is capable of both self- sustaining reproduction and system-improving evolution without outside assistance. This transition requires the development of a high level of factory and/or UP automation and artificial intelligence (AI) that is capable of complex reasoning, imagination and creativity. Finally, I consider the influence of UP-life and UP-paralife on the development of the level of structural complexity in the universe, and I briefly speculate about how these non-biological forms of life and paralife will influence the expansion of scientific knowledge about the universe.

When animals evolve sufficient intelligence and dexterity to be able to learn to fabricate utility products (UPs) like tools, the UP's they produce become part of an induced-reproduction system that intrinsically shares many life-like traits with biological organisms, including genome-like fabrication and operation information that is physically-encoded in the animal fabricator’s neural networks. When this set of life- like traits includes a sufficient capacity for system-improving cultural evolution (UP-evolvability), the UPs become ‘para-alive’, i.e., nearly alive, or a form of non-biological UP-paralife that is equivalent to the life- status of biological viruses, plasmids, and transposons. In the companion paper I focus on the evolution of UP-paralife in the context of modern, language-capable humans and its predicted evolution going forward in time (Rice 2022). Here I look backward in time and focus on the origin of UP-paralife and its subsequent coevolution with human intelligence. I begin by determining the pathways leading to the evolution of large brains in the rare lineages of biological life that have sufficient intelligence to learn to fabricate tools –a critical first step in the evolution of UP-paralife. The simplest forms of these learning- based UPs, made by species like chimpanzees and New Caledonian crows, represent only proto-UP- paralife because they lack sufficient UP-evolvability. Expanded UP-evolvability required a combination of three attributes that enabled continuous niche-expansion of the animal fabricator via a new and advanced form of UP-mediated teamwork (TW): i) self-domestication that facilitated TW among low-related individuals, ii) learned volitional words (protolanguage) that represent ephemeral UPs that coordinate TW, and iii) learned fabrication of simple flaked-stone tools with cutting and chopping capabilities (a UP to make other structural UPs) that expanded teammate phenotypes and TW capabilities. This specific triad of attributes is synergistic because each one acts as a TW-enhancer that can gradually erode different components of the three major constraints on TW operation and expansion: too much selfishness, insufficient coordination signals, and insufficient physical traits of teammates. The increase in UP- evolvability was transformative and marked the origin of UP-paralife and the initiation of coevolution between UP-paralife (cultural evolution) and the intelligence of its hominin/human symbiont (genetic evolution) that fostered 2.5 million years of: i) continuous brain size increase and niche-expansion within the genus Homo, and ii) parallel advances in the diversity, complexity and uses of UP-paralife. This coevolution also fostered evolutionary expansion of word-based communication, and eventually language, that acted in a catalyst-like manner to facilitate the evolution of increasingly complex forms of imagination, reasoning, mentalizing, and UP-generating technology. I next focus on the evolution of creativity in the human lineage –in the form of divergent thinking and creative imagination. I conclude that the evolution of this advanced cognitive feature required a preadaptation of sufficient intelligence and is the component of human cognition that was the major causal factor generating the greatly expanded diversity and complexity of UP-paralife currently associated with modern humans. Lastly, I apply my findings to the issue of the prevalence of extraterrestrial intelligent life. I conclude that any exoplanets with detected chemical life will very rarely (e.g., probability ~10-5 for a planet closely matching Earth’s characteristics) have evolved intelligence equalling or exceeding that of humans.

Human pigmentation has been largely associated with different disease prevalence among populations, but most of these studies are observational and inconclusive. Known to be genetically determined, pigmentary traits have been largely studied by GWAS, mostly in Caucasian ancestry cohorts from North Europe, identifying robustly, several loci involved in many of the pigmentary traits. Here, we conduct a detailed analysis of 13 pigmentary-related traits in a South European cohort of Caucasian ancestry (n=20,000). We observed fair phototype strongly associated with non-melanoma skin cancer and other dermatoses and confirmed by PRS-approach the shared genetic basis with skin and eye diseases, such as melanoma (OR=0.95), non-melanoma skin cancer (OR=0.93), basal cell carcinoma (OR=0.97) and darker phototype with vitiligo (OR=1.02), and cataracts (OR=1.04). Detailed genetic analyses revealed 37 risk loci associated with 10 out of 13 analyzed traits, and 16 fine-mapped genes significantly associated with at least two pigmentary traits. Some of them widely reported, such as MC1R, HERC2, OCA2, TYR, TYRP1, SLC45A2, and unveiling three new candidates RP11-1084J3.4, C1QTNF3 and C17orf112, not reported in GWAS Catalog. These results highlight the importance of phototype assessment as a genetic proxy of skin functionality when evaluating disease screening in mixed populations.

Genetic discoveries related to Alzheimer’s disease and other dementias have been performed using either large cohorts of affected subjects or multiple individuals from the same pedigree, therefore disregarding mutations in the context of healthy groups. Moreover, a large portion of studies so far have been performed on individuals of European ancestry, with a remarkable lack of epidemiological and genomic data from underrepresented populations. The present study aims at scanning 70 single-point mutations on the APP gene in a publicly available genetic dataset including 2.504 healthy individuals from 26 populations, and analyzing their distribution. Moreover, after gametic phase reconstruction, a pairwise comparison of the segments surrounding the mutations was performed to reveal patterns of haplotype sharing that could point to specific cross-population and cross-ancestry admixture events. Eight mutations have been detected in the worldwide dataset, with several of them being specific for a single individual, population or macroarea. Patterns of segment sharing reflect recent historical events of migration and admixture possibly linked to colonization campaigns. These observations reveal the population dynamics of the considered APP mutations in worldwide human groups, and support the development of ancestry-informed screening practices for the improvement of precision and personalized approaches to neurodegeneration and dementias.

Abstract: Chloroplast is a new hotspot in the field of plant transformation system of plant genetic engineering. Initially developed in Chlamydomonas and tobacco, it is now feasible in a broad range of species. They exploit the homologous recombination and segregation pathways acting on chloroplast genomes and are based on direct repeats, transient co-integration or co-transformation and segregation of trait and marker genes. Foreign site-specific recombinases and their target sites provide an alternative and effective method for removing marker genes from plastids.Chloroplast genetic engineering has many advantages over nuclear genetic enginering, especially site-specific introduction of foreign genes ,leading to the absence of gene siliency and positon effect,which providing the available to explore the regulation and mechanism of chloroplast genes`expression in vitro .It also can identify the structure and function of the chloroplast genome , the expression of chloroplast affect the nulear genome. In this paper, the basic methodology of chloroplast transformation, the current techniques and applications, and the future possibilities for Chloroplast genetic engineering was reviewed[1-3].

Abstract. The increasing frequency of isolation of methicillin-resistant Staphylococcus aureus (MRSA) limits the chances of effective antibacterial therapy of staphylococcal diseases and results in development of persistent infection such as bacteremia and osteomyelitis. The aim of this study was to identify features of the MRSAST239 0943-1505-2016 (SA943) genome, that contribute to the formation of both acute and chronic musculoskeletal infections. The analysis was performed using comparative genomics data of the dominant epidemic S. aureus lineages namely ST1, ST8, ST30, ST36, ST239. SA943 genome encodes proteins that provide resistance to the host immune system, suppress immunological memory and form biofilms. The molecular mechanisms of adaptation responsible for development of persistent infection were as follows: amino acid substitution in PBP2 and PBP2a, providing resistance to ceftaroline; loss of a large part of prophage DNA and restoration of nucleotide sequence of beta-hemolysin, that greatly facilitates escape of phagocytosed bacteria from phagosome and formation of biofilms; dysfunction of the AgrA system due to the presence of psm-mec and several amino acid substitutions in the AgrC; partial deletion of nucleotide sequence in genomic island vSAβ resulting in the loss of two proteases of Spl - operon; deletion of SD repeats in SdrE amino acid sequence.

Several observational studies have indicated the potential associations between calcium, vitamin D(Vit-D) and irritable bowel syndrome (IBS). However, the causal relationship deduced from these studies is subjected to residual confounding factors and reverse causation. Therefore, we aim to explore the bidirectional causal effects between serum calcium, Vit-D, PTH and IBS at the genetic level by a two-sample Mendelian randomization (MR) analysis. Sensitivity analyses were performed to evaluate the robustness. The estimates were presented as odds ratio (OR) with their 95% confidence intervals (CIs). The results of the inverse-variance-weighted method did not re-veal any causal relationship shared between genetically predisposed calcium (OR = 0.92, 95% CI: 0.80-1.06, P = 0.25) and Vit-D (OR = 0.99, 95% CI: 0.83-1.19, P = 0.94) level and the risk of IBS. The bidirectional analysis demonstrated that genetic predisposition to IBS was associated with a de-creased level of PTH (beta: -0.19, 95%CI: -0.34 to -0.04, P = 0.01). In conclusion, the present study indicates no causal relationship between the serum calcium and Vit-D concentrations and the risk of IBS. The potential mechanisms by which IBS affects serum PTH need to be further investigated.

Delayed cancer detection is one of the common causes of poor prognosis in case of many cancers including the cancers of the oral cavity. Despite improvement and development of new and efficient gene therapy treatments, very little has been done to algorithmically assess the impedance of these carcinomas. In this work, we attempt to annotate viable attributes in oral cancer gene datasets for identification of gingivobuccal cancer (GBC). We further apply supervised and unsupervised machine learning methods to the gene datasets revealing key candidate attributes for GBC prognosis. Our work highlights the importance of automated identification of key genes responsible for GBC that could perhaps be easily replicated to other forms of oral cancer detection.

RNA interference (RNAi) is one of the main mechanisms for disease resistance and small RNA production in plants. The main proteins involved in RNAi include Dicer-like (DCL), RNA-dependent RNA polymerase (RDR), double-stranded RNA-binding (DRB), and Argonaute (AGO). Juglandaceae contains a variety of important woody plants, and walnuts are one of the four major woody plant groups and one of the four major dried fruits in the world. To clarify the evolution and functional differentiation of RNAi-related proteins in the walnut (Juglans regia) genome, this study integrated various web resources from gene family acquisition to structural analysis and transcriptome data to correlate walnuts and their congeners. The walnut genome has 5 DCL, 13 RDR, 15 DRB and 15 AGO genes, similar genes encoding conserved protein structural domains and conserved motifs with similar subcellular localization. Walnut AGO proteins are classified into three classes and seven subclasses. The DCL is divided into four categories, while RDR is mainly divided into four categories, and DRB can be divided into six categories. The exception is that the copy number of walnut RDR1 is 9, in which seven RDR1 are distributed in clusters on chromosome 16. Purifying selection drove the formation of walnut genes, but protein classes were subjected to varying degrees of purifying selection. Additionally, these results showed some similarity in other plants of the walnut family. Moreover, different RNAi-related genes of walnut produced abundant selective expression in response to different tissues and stresses. In this study, DCL, RDR, DRB and AGO gene families were identified and analysed in the genome of the walnut family for the first time and preliminarily examined the evolution, structure and expression characteristics of these families to provide a preliminary basis for the evolution of the walnut RNAi pathway and breeding research.

Sulfate transporters (SULTRs) are responsible for the uptake of the sulfate (SO42−) ions in the rhizosphere by the roots and their distribution in plant organs. In this study, SULTR family members in the genome of the two oilseed crops, Camelina sativa, and Brassica napus, were identified and characterized based on their sequence structure, duplication events, phylogenetic relationships, phosphorylation sites, and expression levels. Herein, 36 and 45 putative SULTR genes were recognized from the genome of C. sativa, and B. napus, respectively. SULTR proteins were predicted as basophilic proteins with low hydrophilicity in both studied species. According to phylogenetic relationships, we divided SULTRs into five groups, in which SULTRs 3 showed highest variation. Besides, several duplication events were observed between SULTRs. The first duplication event was predicted approximately five million years ago between three SULTRs 3.1 in C. sativa. Two subunits were indicated in the 3D structure of SULTRs that the active binding sites differed between C. sativa and B. napus. According to available RNA-seq data, SULTRs showed diverse expression in tissues and response to stimuli. SULTRs 3 showed an expression in all tissues. SULTRs 3.1 were more upregulated in response to abiotic stresses in C. sativa, while SULTRs 3.3, and SULTRs 2.1 showed an upregulation in B. napus. Furthermore, SULTRs 3 and SULTRs 4.1 showed an upregulation in response to biotic stresses in B. napus. Based on the distribution of cis-regulatory elements in the promoter region, we speculated that SULTRs might be controlled by phytohormones such as ABA, and MeJA. Therefore, it seems that SULTR genes in C. sativa have been more influenced by evolutionary processes and have acquired more diversity.

Background: The Ibizan Hound is a canine breed native to the Mediterranean region, where leishmaniasis is an endemic zoonosis. Several studies indicate a low prevalence of this disease in dogs, whereas other canine breeds present a high prevalence. However, the molecular underlaying mechanisms yet remains unknown. Methods: In this study, we analyze the haplotypes of genes encode cytokines related to immune response of Leishmania infantum infection in twenty-four Boxer and twenty-four Ibizan hound apparently healthy using CanineHD DNA Analysis BeadChip including 165,480 mapped positions. Results: The results show that several haplotypes of genes encoding Interleukin 6 (IL6) and Interferon gamma (IFNG) are related to Interferon gamma (IFN-γ), and Interleukins (IL) 2 and 18 serum levels. Our results indicate that the regulation of immune response is different in the two canine breeds analyzed and are related to the haplotype compositions of the genes encoding these cytokines. Conclusions: Future studies are needed to elucidate whether these differences and haplotypes are related to different phenotypes in immune response and expression gene regulation to L. infantum infection in dogs.

Most pathogens are capable of infecting multiple host species and such cross-species transmission (CST) can have dramatic consequences, as highlighted by recent disease emergence events affecting human, animal and plant health. Understanding the ecological and evolutionary factors that constrain or facilitate the ability of disease agents to infect and establish in a novel host is therefore a timely and important research area. Previous work across different pathogens, including rabies virus (RABV), found that increased evolutionary distance between hosts reduces the frequency of cross-species transmission and of permanent host shifts. However, whether this effect of host relatedness still holds for transmission among recently diverged hosts, and the importance of this effect relative to other predictors, is not well understood. We addressed this question by quantifying the CST frequency of RABV between North American bat species within the genus Myotis, using a multi-decade data set containing 128 nucleoprotein (N) sequences from ten host species. For comparison, we also conducted an equivalent analysis of a RABV dataset from North American bat species comprising nine genera. We found that at the within genus scale, host relatedness failed to explain the frequency of CST events. However, CST frequency increased with overlap in species’ host range. Moreover, we found evidence of CST occurring among a higher proportion of species, and CST more frequently resulting in sustained transmission in the novel host in the Myotis dataset compared to the multi-genus dataset. Our results suggest that among recently diverged species, the ability to infect a novel host is no longer restricted by physiological barriers but instead is limited by physical contact. Our results improve predictions of where future CST events for RABV might occur and clarify the relationship between host divergence and pathogen emergence.

Since the discovery of mitochondrial-derived peptides (MDP), their participation in cellular metabolism is no longer considered as the sole function of the mitochondria, but importance was also attached to its role as a source of protective factors of metabolic stress. These peptides are encoded in the mitochondrial genome and translated into the mitochondria or cytoplasm, to signal within the cell or be released and bind to membrane receptors. The objective of this work was explored and compare the frequency of MT-RNR1 and MT-RNR2 variants in sequences obtained from T2D individuals and control population. 213 different mitochondrial polymorphisms previously reported in the literature associated with T2D and cardiovascular diseases were analyzed. We can found three variants in the MT-RNR1 not related with MOTS-c coding sequence: m.1189T>C (rs28358571), m.1420T>C (rs111033356), and m.1438A>G (rs2001030); and secondly, three polymorphisms associated to MT-RNR2 m.2667T>C (rs878870626) related to humanin, m.1811A>G (rs28358576) in SHPL3 and m.3027T>C (rs199838004) in SHPL6 associated with statistical differences between the T2D and control group. All these findings were previously related to cardiovascular complications in literature and, as far as we know, relating for the first time in diabetic patients.

The main objective of this study is to demonstrate proof of principle that computational fluid dynamics (CFD) modeling is a tool for studying the contribution of covert and overt vascular architecture to the risk of cerebrovascular disease in in sickle cell disease (SCD) as well as uncover one or more mechanism of response to therapy such as chronic red blood cell (cRBC) transfusion. We analyzed baseline (screening), pre-randomization and study exit magnetic resonance angiogram (MRA) images from 10 (5 each from the transfusion and observation arms) pediatric sickle SCD participants in the silent cerebral infarct transfusion (SIT) trial, using CFD modeling. We reconstructed the intracranial portion of the internal carotid artery and branches and extracted the geometry using 3D Slicer. We cut specific potions of the large intracranial artery to include segments of the internal carotid, middle, anterior, and posterior cerebral artery such that the vessel segment analyzed extended from the intracranial beginning of the internal carotid artery up to immediately after (~0.25 inches) the middle cerebral artery branching point. Cut models were imported into Ansys 2021R2/2022R1 and laminar and time-dependent flow simulation was performed. Change in time averaged mean velocity, wall shear stress, and vessel tortuosity were compared between the observation and cRBC arm. We did not observe a correlation between time averaged mean velocity (TAMV) and mean transcranial doppler (TCD) velocity at study entry. There was also no difference in change in time average mean velocity, wall shear stress (WSS), and vessel tortuosity between the observation and cRBC transfusion arms. WSS and TAMV were abnormal for 2 (developed TIA) out of the 3 participants (one participant had SCI) that developed neurovascular outcomes. CFD approaches allows for the evaluation of vascular topology and hemodynamics in SCD using MRA images. In this proof of principle study, we show that CFD could be a useful tool and we intend to carry out future studies with a larger sample to enable more robust conclusions.

Systems genetics is key for integrating a large number of variants associated with diseases. Vitamin K (VK) is one of the scarcely studied conditions in lieu of ascertaining either the differentially expressed genes (DEGs) or variants in an individual subpopulation of diseased phenotypes associated with VK, viz. myocardial infarction, renal failure, prostate cancer, thrombosis, thrombocytopenia, coagulation related diseases to name a few. In this work, we have screened characteristic DEGs common to three VK-related diseases, viz. myocardial infarction, renal failure and prostate cancer and asked whether or not any DEGs in addition to pathogenic variants are common to these conditions. We attempt to bridge the gap in finding characteristic biomarkers and discuss the role of long noncoding RNAs (lncRNAs) in the biogenesis of VK deficiencies.

Lifestyle has been associated with breast cancer risk through different pathways, including oxidative stress. Antioxidant enzymes are endogenous defense mechanisms against oxidative stress damage, and this response might be modulated by the genetic variation in these enzyme-codifying genes. This study aimed to analyze the synergistic effect of an antioxidant Healthy Lifestyle Index (HeLiX) composed of principal components of Western dietary pattern, alcohol consumption, smoking and physical activity, and genetic polymorphisms in the first-line antioxidant response family genes SOD, GPX, and CAT on breast cancer risk. We included 176 SNPs, and only CAT rs554576 remained significant after correction for multiple comparisons. Breast cancer odds were reduced at the highest (T3) and medium (T2) tertiles of the HeLiX. When stratified by HeLiX, we observed a reduced risk of breast cancer with at least one T-allele, and the effect increased in a dose-dependent manner. Compared to the reference category (HeLiX T1 and AA genotype), women at the HeLiX T3 with AT and TT genotypes in postmenopausal women showed an OR = 0.15 (95% CI 0.07–0.32). For HeLiX T2 and AT genotype OR = 0.26 (95% CI 0.13–0.54); for TT genotype OR = 0.24 (95% CI 0.12–0.45). For premenopausal women, at the HeLiX T3 and AT genotype OR = 0.29 (95% CI 0.13–0.62); for the TT genotype OR = 0.21 (95% CI 0.08–0.51). We also observed an inverse association for HeLiX T2 and TT genotype (OR 0.39 95% CI 0.17–0.87). Our study shows a significant synergistic gene-environment interaction on an additive scale, contributing to understanding pathways involved in breast cancer etiology and prevention.

Taking a DNA-sequence, a word with letters/bases A, T, G and C, as the relation between the generators of an infinite group $\pi$, one can discriminate two important families: (i) the cardinality structure for conjugacy classes of subgroups of $\pi$ is that of a free group on $1$ to $4$ bases and the DNA word, viewed as a substitution sequence, is aperiodic. (ii) The cardinality structure for conjugacy classes of subgroups of $\pi$ is not that of a free group, the sequence is generally not aperiodic and topological properties of $\pi$ have to be determined differently. The two cases rely on DNA conformations such as A-DNA, B-DNA, Z-DNA, G-quadruplexes, etc. We found a few salient results: Z-DNA, when involved in transcription, replication and regulation in a healthy situation implies (i). The sequence of telomeric repeats comprising $3$ distinct bases, most of the time, satisfies (i). For $2$-base sequences in the free case (i) or non free case (ii), the topology of $\pi$ may be found in terms of the $SL(2,\mathbb{C})$ character variety of $\pi$ and the attached algebraic surfaces. The linking of two unknotted curves --the Hopf link-- may occur in the topology of $\pi$ in cases of biological importance, in telomeres, G-quadruplexes, hairpins and junctions, a feature we already found in the context of models of topological quantum computing. For $3$- and $4$ base sequences, other knotting configurations are noticed and a building block of the topology is the $4$-punctured sphere. Our methods have the potential to discriminate potential diseases associated to the sequences.

Targeting genes involved in sexual determinism for vector or pest control purpose, requires a better understanding of their polymorphism in natural populations in order to ensure a rapid spread of the construct. By using genomic data from An. gambiae s.l., we analyzed the genetic variation and the conservation score of the fru gene in 18 natural populations across Africa. A total of 34339 SNPs were identified including 3.11% non-synonymous segregating sites. Overall, the nucleotide diversity was low and the Tajima's D neutrality test was negative indicating an excess of low frequency SNPs in the fru gene. The allelic frequencies of the non-synonymous SNPs were low (freq < 0.26) except two SNPs identified at high frequencies (freq > 0.8) in the Zinc-finger A and B protein domains. The conservation score was variable throughout the fru gene with maximum values in the exonic compared to the intronic regions. These results showed a low genetic variation in overall the exonic regions especially the male sex-specific exon and the BTB-exon 1 of the fru gene. These findings are crucial for the development of a gene drive construct targeting the fru gene that can rapidly spread without encountering resistance in wild populations.

Free-living amoebae (FLA) are protozoa widely distributed in the environment, found in a great diversity of terrestrial biomes. However, few genera of FLA are linked to human infections. Within these genera, Acanthamoeba spp., classified by genotypes (T1-T23), being T1, T2, T4, T5, T10, T12, and T18 as capable of causing granulomatous amoebic encephalitis (GAE) in immunocompromised patients mostly and Acanthamoeba keratitis related to genotypes T2, T3, T4, T5, T6, T10, T11, T12 and T15 in apparently healthy patients. Meanwhile, Naegleria fowleri is the causative agent of an acute infection called primary amoebic meningoencephalitis (PAM), while Balamuthia mandrillaris, like some Acanthamoeba genotypes, causes GAE, differing from the latter in the description of numerous cases in patients immunocompetent. Finally, other FLA related to the pathologies mentioned above have been reported; Sappinia pedata is responsible for one case of amoebic encephalitis; Vermamoeba vermiformis has been found in cases of ocular damage, and its extraordinary capacity as endo-cytobiont for microorganisms of public health importance such as Legionella pneumophila, Bacillus anthracis, Pseudomonas aeruginosa, among others. In this review, issues related to the epidemiology of each one are addressed, updating their geographic distribution and cases reported in recent years for pathogenic FLA.

Alzheimer´s Disease (AD) has currently no effective treatment; however, preventive measures can significantly delay the progress/onset of the disease. Thus, accurate and early prediction of risk is an important strategy to alleviate the AD burden. Neuroinflammation is a major factor prompting the onset of the disease. Inflammation exerts its toxic effect via multiple mechanisms. Amongst others, it is affecting gene expression via modulation of non-coding RNAs (ncRNAs), such as miRNAs. Recent evidence supports that inflammation can also affect long non-coding RNAs (lncRNAs) expression. While the association between miRNAs and inflammation in AD has been extensively studied, the role of lncRNAs in neurodegenerative diseases has been less explored. In this review, we focus on lncRNAs and inflammation in the context of AD. Furthermore, since plasma-isolated extracellular vesicles (EVs) are increasingly recognized as an effective monitoring strategy of brain pathologies, we have focused on the studies reporting dysregulated lncRNAs in EVs isolated from AD patients and controls. The revised literature shows a positive association between pro-inflammatory lncRNAs and AD. However, the reports evaluating lncRNAs alterations in EVs isolated from plasma of patients and controls, although still limited confirm the value of specific lncRNAs associated with AD as reliable biomarkers. This is an emerging field that will open new avenues to improve risk prediction, patients’ stratification and may lead to the discovery of potential novel therapeutic targets for AD

The article is devoted to algebraic modeling of universal rules of stochastic organization of genomic DNA of higher and lower organisms, previously published by the author. The proposed algebraic apparatus, which uses formalisms of quantum mechanics and quantum informatics and which is based on the so-called tensor-unitary transformations of vectors that generate families of interrelated stochastic-deterministic vectors of increased dimensions. The features of the vectors' interconnections in these families model the stochastic-deterministic properties of the named phenomenological rules. There are new approaches to modeling of developing multi-parameter biosystems, whose number of parameters increases in the course of step-by-step development. In the light of the presented materials, the issues of fractal-like organization in genetically inherited biosystems are considered. The development of the theory of stochastic determinism as an antipode of deterministic chaos is discussed.

Although originally multi-ethnic in its structure, nowadays the Calabria region of southern Italy represents an area with a low genetic heterogeneity and a high level of consanguinity that allows rare mutations to be maintained due to the founder effect. A complex research methodology ranging from clinical activity to genealogical reconstruction of families/populations along the centuries, creation of databases, and molecular/genetic research, has been modelled on the characteristics of the Calabrian population for more than three decades. This methodology allows to the identification of several novel genetic mutations or variants associated with neurodegenerative diseases. In addition, in this population it has been reported a higher prevalence of several hereditary neurodegenerative diseases such as Alzheimer’s disease, Frontotemporal dementia, Parkinson’s disease, Niemann Pick type C disease, Spino-cerebellar ataxia, Creutzfeldt–Jakob disease and Gerstmann Straussler Scheincker disease. Thus, Calabria constitutes a model for the study of neurodegenerative diseases, a sort of "outdoor laboratory" useful for the advancement of knowledge in this field. Here, we summarize and discuss some results of research data supporting the view that Calabria is a genetic isolate and could represent a useful model for the study and characterization of neurodegenerative diseases.

The Peruvian creole cattle (PCC) is a neglected breed, and is an essential livestock resource in the Andean region of Peru. To develop a modern breeding program and conservation strategies for the PCC, a better understanding of the genetics of this breed is needed. We sequenced the whole genome of the PCC using a paired-end 150 strategy on the Illumina HiSeq 2500 platform, obtaining 320 GB of sequencing data. The obtained genome size of the PCC was 2.77 Gb with a contig N50 of 108Mb and 92.59% complete BUSCOs. Also, we identified 40.22% of repetitive DNA of the genome assembly, of which retroelements occupy 32.39% of the total genome. A total of 19,803 protein-coding genes were annotated in the PCC genome. We downloaded proteomes and genomes of the Bovinae subfamily, and conducted a comparative analysis with our draft genome. Phylogenomic analysis showed that PCC is related to Bos indicus. Also, we identified 7,746 family genes shared among the Bovinae subfamily. This first PCC genome is expected to contribute to a better understanding of its genetics to adapt to the tough conditions of the Andean ecosystem, and evolution.

Breast cancer has an important incidence in the worldwide female population. Although alterations in the mitochondrial genome probably play an important role in carcinogenesis, the actual evidence is ambiguous and inconclusive. The purpose of the present work was to explore mitochondrial sequences of clinical cases with breast cancer from different origins and determine the polymorphisms associated. The search for complete and partial mtDNA sequences obtained from breast cancer patients and controls was performed in NCBI Genbank database. We identified 124 mtDNA sequences associated to breast cancer cases of which 86 were complete and 38 partial sequences. Of these 86 complete sequences, 52 belong to patients with a confirmed diagnosis of breast cancer and 34 sequences were obtained from healthy mammary tissue of the same patients used as controls. From mtDNA analysis, two polymorphisms with significative statistical differences were found in D130 in sequences analyzed: m.310del (rs869289246) in 34.6% (27/78) breast cancer cases and 61.7% (21/34) of controls; and m.315dup (rs369786048) in 60.2% (47/78) of breast cancer cases and 38.2% (13/34) of controls. Also, the variant m.16519T>C (rs3937033) was found in 59% of control sequences and 52% of breast cancer sequences with a significant statistical difference. Polymorphic changes are evolutionarily related to haplogroup H of Indo-European and Euro Asiatic origins, however, were found in all non-European sequences with breast cancer.

The ability to directly look into genome sequences has opened great opportunities in plant breeding. Yet, the assembly of full-length chromosomes remains one of the most difficult problems in modern genomics. Genetic maps are commonly used in de novo genome assembly and are constructed on the basis of a statistical analysis of the number of recombination rather than physical distance in base pairs. This may affect the accuracy of the ordering and orientation of scaffolds within the chromosome, especially in the region of recombination suppression. Here we report the use of Tyr-FISH for validation of genetic and pseudochromosome maps. For probe design we developed a pipeline that is based on selection of a unique sequence with minimal potential fluorescent background arising from non-specific in situ hybridization. In total of 24 unique genes were located on physical chromosomes 2 and 6. The order of markers has been corrected by integration genetic and cytogenetic maps. Tyr-FISH mapping showed that the order of 23.1% (chromosome 2) and 27.3% (chromosome 6) of the tested genes differed between physical chromosomes and pseudochromosomes. Also the position of the mlh1 gene, which was not on the genetic map, was defined on physical chromosome 2. The results suggest the possibility of Tyr-FISH mapping of any gene that may not be neither on the genetic map nor in the assembly. Hence, Tyr-FISH provides valuable information for the improvement of the genome assembly.

Genomics involving tens of thousands of genes is a complex system determining phenotype. An interesting and vital issue is that how to integrate highly sparse genetic genomics data with a mass of minor effects into prediction model for improving prediction power. We find that deep learning method can work well to extract features by transforming highly sparse dichotomous data to lower dimensional continuous data in a non-linear way. This idea may provide benefits in risk prediction based on genome-wide data associated e.g. integrating most of the information in the genotype data. Hence, we developed a multi-stage strategy to extract information from highly sparse binary genotype data and applied it for risk prediction. Specifically, we first reduced the number of biomarkers via a univariable regression model to a moderate size. Then a trainable auto-encoder was used to extract compact representations from the reduced data. Next, we performed a LASSO problem process over a grid of tuning parameter values to select the optimal combination of extracted features. Finally, we applied such feature combination to two prognostic models, and evaluated predictive effect of the models. The results of simulation studies and real data applying indicated that these highly compressed transformation features could better improve predictive performance and did not easily lead to over-fitting.

Background: The trends in the role of the germline in epigenetic transgenerational inheritance starts with the environmental factors acting on the first generation of a gestating mother. These factors influence the developing second-generation fetus by altering gonadal development, thereby reprogramming the primordial germ cell DNA methylation and leading to consequences that would be seen along generations. Objective: Despite these epigenetic factors now surfacing, the few available studies are on animal-based experiments and to make follow up on human intergenerational trials might take decades. To this response, this study aimed to determine the influence of parental energy, toxicant exposure, age and nutrient restriction on the early life of offspring growth in Gambia. Method: The study was on population-based observational study on parental energy influence, toxicant exposure, age, and nutrient restriction on offspring growth in Gambia. Results: This study showed that parents who worked in industrial areas were more likely to have offspring with poor psychosocial skills. In addition, mothers who are exposed to oxidative stress and high temperature are more likely to have offspring with poor psychosocial skills. Mothers who consumed a high protein diet were almost three more likely to have infants with good psychosocial skills in their offspring. Furthermore, there was a negative correlation between maternal stress during pregnancy and psychosocial skills of offspring. Conclusion: This study was able ascertain if maternal diet during gestation, toxicant exposure, maternal stress and parental smoking habits have influence on the early life of offspring. While the study is recommending a large sample size study to eliminate selection bias, there should be an increased level of awareness of mothers on their offspring's health and their husbands' lifestyles that might influence the adulthood health of offspring.

Congenital Disorders of Glycosylation (CDG) are a heterogenous group of primarily autosomal recessive mendelian diseases caused by disruptions in the synthesis of lipid linked oligosaccha-rides and their transfer to proteins. CDGs affect multiple organ systems and vary in presentation, even within families. Here we describe a chemically induced mouse mutant, tvrm76, with early onset photoreceptor degeneration. The recessive mutation was mapped to Chromosome 9 and as-sociated with a missense mutation in the Dpagt1 gene encoding UDP-N-acetyl-D-glucosamine:dolichyl-phosphate N-acetyl-D-glucosaminephosphotransferase (EC 2.7.8.15). The mutation is predicted to cause a substitution of aspartic acid with glycine at residue 166 of DPAGT1. Increased expression of Ddit3, and elevated levels of HSPA5 (BiP) sug-gest the presence of early-onset endoplasmic reticulum (ER) stress. These changes were associated with induction of photoreceptor apoptosis in tvrm76 retinas. Mutations in human DPAGT1 cause Myasthenic Syndrome 13 and severe forms of Congenital Disorder of Glycosylation Type Ij. In contrast, Dpagt1tvrm76 homozygous mice present with congenital photoreceptor degeneration without overt muscle or muscular junction involvement. Our results suggest the possibility of DPAGT1 mutations in human patients that present primarily with retinitis pigmentosa with little or no muscle disease. Variants in DPAGT1 should be considered when evaluating cases of non-syndromic retinal degeneration.

Myocilin is an enigmatic glaucoma-associated glycoprotein whose biological role remains incompletely understood. To gain novel insight into its normal function, we used transposon-mediated transgenesis to generate the first zebrafish line stably overexpressing myocilin [Tg(actb1:myoc-2A-mCherry)]. qPCR showed an approximately four-fold increased myocilin expression in transgenic zebrafish embryos (144 hpf). Adult (13 months old) transgenic animals displayed variable and age-dependent ocular anterior segment alterations. Almost 60% of two-years old male, but not female, transgenic zebrafish developed enlarged eyes with severe asymmetrical and variable abnormalities in the anterior segment, characterized by corneal limbus hypertrophy, and thickening of the cornea, iris, annular ligament and lens capsule. The most severe phenotype presented small or absent ocular anterior chamber and pupils, due to iris overgrowth along with dysplastic retinal growth and optic nerve hypertrophy. Immunohistochemistry revealed increased presence of myocilin in most altered ocular tissues of adult transgenic animals, as well as signs of retinal gliosis, and expanded ganglion cells and nerve fibers. The preliminary results indicate that these cells contributed to retinal dysplasia. Visual impairment was demonstrated in all old male transgenic zebrafish. Transcriptomic analysis of the abnormal transgenic eyes identified disrupted expression of genes involved in lens, muscular and extracellular matrix activities, among other processes. In summary, the developed transgenic zebrafish provides a new tool to investigate this puzzling protein and provides evidence for the role of zebrafish myocilin in ocular anterior segment and retinal biology, through the influence of extracellular matrix organization and cellular proliferation.

Pharmacogenetics (PGx) aims to identify the genetic factors that determine inter-individual differences in response to drug treatment maximizing efficacy while decreasing the risk of adverse events. Estimating the prevalence of PGx variants involved in drug response, is a critical preparatory step for large-scale implementation of a personalized medicine program in a target population. Here, we profiled pharmacogenetic variation in fourteen clinically relevant genes in a representative sample set of 1,577 unrelated sequenced Sardinians -- an ancient island population that well accounts for genetic variation in Europe as a whole and, at the same time are enriched in genetic variants that are very rare elsewhere. To this end, we used PGxPOP, a PGx allele caller based on the guidelines created by the Clinical Pharmacogenetics Implementation Consortium (CPIC), to identify the main phenotypes associated with the PGx alleles most represented in Sardinians. We estimate that 99.43% of Sardinian individuals may potentially respond atypically to at least one drug, that on average each individual is expected to have an abnormal response to about 17 drugs, and that for 27 drugs the fraction of the population at risk of atypical responses to therapy is more than 40%. Finally, we highlighted 174 pharmacogenetic variants for which the minor allele frequency is at least 10% higher among Sardinians as compared to other European populations, a fact that may contribute to substantial interpopulation variability in drug response phenotypes. This study provides baseline information for further large-scale pharmacogenomic investigations in the Sardinian population and underlines the importance of the PGx characterization of diverse than European populations as Sardinians.

Lynch syndrome (LS) is an inherited disorder in which affected individuals have a significantly higher-than-average risk of developing colorectal and non-colorectal cancers, often before the age of 50 years. In LS, mutations in DNA repair genes lead to a dysfunctional post-replication repair system. As a result, the unrepaired errors in coding regions of the genome produce novel proteins, called neoantigens. Neoantigens are recognised by the immune system as foreign and trigger an immune response. Due to the invasive nature of cancer screening tests, universal cancer screening guidelines unique for LS (including colonoscopy) are poorly adhered to by LS variant heterozygotes (LSVH). Currently, it is unclear whether immunogenomic components produced as a result of neoantigen formation can be used as novel biomarkers in LS. We hypothesize that: (i) LSVH produce measurable and dynamic immunogenomic components in blood, and (ii) these quantifiable immunogenomic components correlate with cancer onset and stage. Here, we discuss feasibility and propose the potential to: (a) identify personalised novel immunogenomic biomarkers, (b) reduce invasive cancer screening tests and thus increase non-invasive cancer surveillance, (c) improve compliance and adherence to recommended cancer screening guidelines in LSVH.

Cutaneous Melanoma (CM) is the most lethal form of skin cancer if it becomes metastatic, where treatment options and survival chances decrease dramatically. Immunotherapy treatments based on the immunologic checkpoint inhibitors (PD-1 and CTLA4) constituted a main breakthrough in the treatment of metastatic CM, particularly in the long-term benefit. However, several molecular pathways are responsible for the failure of this strategy in about 50-70% of CM patients. Some Long Non-coding RNAs (lncRNAs), and circular RNAs (circRNA) are implicated in triggering pro- and antitumorigenic responses to various cancer treatments. The relationship between lncRNA, circRNA and Immune Checkpoint Blockade (ICB) immunotherapy is not extensively explored in cutaneous metastatic melanoma (CMM). The aim of this study is to evaluate the potential role of both circRNA and lncRNA as a predictive immunotherapy biomarker in CMM. RNA-seq from 12 FFPE samples from the metastatic biopsy of metastatic melanoma patients treated with Nivolumab were analyzed. Our findings indicate that specific lncRNA and circRNA are involved in regulatory networks of the immune response against metastatic melanoma under treatment with nivolumab. Moreover, we have established a risk score that allows the prediction of Overall survival (OS) and Progression-free survival (PFS) of CMM patients with high accuracy. This proof of principle work provides a possible insight on the function of ceRNA, contributing to decipher the complex molecular mechanism of ICB cancer treatment response.

Abstract: Since ancient times aging has also been regarded as a disease, and humankind has always strived to extend the natural lifespan. Analyzing the genes involved in aging and disease allows for finding important indicators and biological markers for pathologies and possible therapeutic targets. An example of the use of omics technologies is the research regarding aging and the rare and fatal premature aging syndrome progeria (Hutchinson-Gilford progeria syndrome, HGPS). In our study, we focused on the in silico analysis of differentially expressed genes (DEGs) in progeria and aging, using a publicly available RNA Seq dataset (GEO dataset GSE113957) and a variety of bioinformatics tools. We identified several genes that appear to be involved both in natural aging and progeria (KRT8, KRT18, ACKR4, CCL2, UCP2, ADAMTS15, ACTN4P1, WNT16, IGFBP2). Further analyzing these genes and the pathways involved confirmed their possible roles in aging, suggesting the need for further in vitro and in vivo research. The graphical abstract illustrates the analysis workflow we used and will introduce in the following as an example to demonstrate the power of omics and bioinformatics.

Recent results have shown that the human malaria-resistant hemoglobin S mutation originates de novo more frequently in the gene and in the population where it is of adaptive significance, namely, in the hemoglobin subunit beta gene compared to the non-resistant but otherwise identical 20A>T mutation in the hemoglobin subunit delta gene, and in sub-Saharan Africans, who have been subject to intense malarial pressure for many generations, compared to Northern Europeans, who have not. This finding raises a fundamental challenge to the traditional notion of accidental mutation. Here we address this finding with the replacement hypothesis, according to which pre-existing genetic interactions can lead directly and mechanistically to mutations that simplify and replace them. Thus, a gradual evolutionary process under selection can hone in on interactions of importance for the currently evolving adaptations, from which large-effect mutations follow that are relevant to these adaptations. We exemplify this hypothesis using multiple types of mutation, including gene fusion mutations, gene duplication mutations, A>G mutations in RNA-edited sites and transcription-associated mutations, and place it in the broader context of a system-level view of mutation origination called Interaction-based Evolution. Potential consequences include that similarity of mutation pressures may contribute to parallel evolution in genetically related species, that the evolution of genome organization may be driven by mutational mechanisms, that transposable element movements may also be explained by replacement, and that long-term directional mutational responses to specific environmental pressures are possible. Such responses need to be further tested by future studies in natural and artificial settings.

The number of children diagnosed with autism spectrum disorder (ASD) has increased substantially over the past two decades with current research unable to fully account for this dramatic increase in prevalence. One explanation proposes that both intrinsic (e.g., genetic) and extrinsic (e.g., environmental) risk factors may be involved in the etiology of ASD. The goal of this review paper is to explore modifiable pathways for intervention in children at risk for ASD, specifically examining how early social experience may be correlated with epigenetic change in genes associated with autism. We present an innovative model which proposes that polygenic risk and social experience (via epigenetic mechanisms) may both contribute to the observed ASD phenotype. Previous research on genetic, environmental, and epigenetic mechanisms implicated in the etiology of ASD will be reviewed, with an emphasis on the oxytocin receptor gene, which is epigenetically altered by early social experience, plays a crucial role in mammalian social and cognitive development, and is associated with both genetic and epigenetic risk for ASD. Identifying intrinsic (e.g., genetic) and extrinsic (e.g., social experience) risk markers for ASD, a combination of which has not previously been examined, would transform our understanding of this condition, facilitate earlier identification of ASD risk, and guide early intervention efforts. This may have a far-reaching impact on individuals with ASD, their families, and society.

Decades of research identified numerous gene biomarkers of cardiac diseases whose restored sequence or/and expression level was hoped to recover the normal cardiac function. However, each human has unique and dynamic pathophysiological characteristics resulting from the unrepeatable combination of favoring factors such are: race, sex, age, medical history, diet, stress, exposure to toxins, habits etc. As such, no treatment fits everybody and finding personalized solutions is a top priority for medicine of 21st century. The Genomic Fabric Paradigm (GFP) provides the most theoretically possible comprehensive characterization of the transcriptome, its alterations in disease and recovery following a treatment. By attaching to each gene the independent average expression level, expression variation and expression coordination with each other gene, GFP delivers thousands times more information than the traditional analysis. This report presents the theoretical bases of the GFP and some applications to our microarray data from mouse models of post ischemic, and constant and intermittent hypoxia-induced heart failure. The GFP analyses revealed novel transcriptomic aspects of the gene expression control and networking under ischemic conditions. Through all-inclusive characterization of the transcriptome and the unrepeatable gene hierarchy in each condition, GFP is an essential avenue towards development of a truly personalized cardiogenomic therapy.

In the past decade, RNA fragments derived from full length small nucleolar RNAs (snoRNAs) have been shown to be specifically excised and functional. These sno-derived RNAs (sdRNAs) have been implicated as gene regulators in a multitude of cancers, controlling a variety of genes post-transcriptionally via association with the RNA-induced silencing complex (RISC). In this review, we have summarized the literature connecting sdRNAs to cancer gene regulation. SdRNAs possess miRNA-like functions, and are able to fill the role of tumor-suppressor or tumor-promoter in a tissue context-dependent manner. Indeed, there are many miRNAs that are actually derived from snoRNA transcripts, meaning that they are truly sdRNAs and as such are included in this review. As sdRNAs are frequently discarded from ncRNA analyses, we emphasize that sdRNAs are functionally relevant gene regulators and likely represent an overlooked subclass of miRNAs. Based on the evidence provided by the papers reviewed here, we propose that sdRNAs deserve more extensive study to better understand their underlying biology and to identify previously overlooked biomarkers and therapeutic targets for a multitude of human cancers.

The Drosophila imaginal disc has been an excellent model for the study of developmental gene regulation. In particular, long non-coding RNAs (lncRNAs) have gained widespread attention in recent years due to their important role in gene regulation. Their specific spatiotemporal expressions further support their role in developmental processes and diseases. In this study, we explored the role of a novel lncRNA in Drosophila leg development by dissecting and dissociating w1118 third-instar larval third leg (L3) discs into single cells and single nuclei, and performing single-cell RNA-sequencing (scRNA-seq) and single-cell assays for transposase-accessible chromatin (scATAC-seq). Single-cell transcriptomics analysis of the L3 discs across three developmental timepoints revealed different cell types and identified lncRNA:CR33938 as a distal specific gene with high expression in late development. This was further validated by fluorescence in-situ hybridization (FISH). The scATAC-seq results reproduced the single-cell transcriptomics landscape and elucidated the distal cell functions at different timepoints. Furthermore, overexpression of lncRNA:CR33938 in the S2 cell line increased the expression of leg development genes, further confirming its important role in development.

We are all similar, but a bit different. These differences are partially due to variations in our genomes and are related to the heterogeneity of symptoms and responses to treatments that patients exhibit. Most animal studies are performed in one single strain with one manipulation. However, due to the lack of variability, therapies are not always reproducible when treatments are translated to humans. Panels of already sequenced organisms are valuable tools for mimicking human phenotypic heterogeneities and gene mapping. This review summarizes the current knowledge of mouse, fly and yeast panels with insightful applications for translational research.

Meier–Gorlin syndrome (MGS) is a rare genetic developmental disorder that causes primordial proportional dwarfism, microtia, absent or hypoplastic patellae and other skeletal anomalies. Overlapping skeletal symptoms make MGS difficult to diagnose clinically. We describe a 3-year-old boy with short stature, recurrent respiratory infections, short-rib dysplasia, tower head and facial dysmorphisms who was admitted to the Tomsk Genetic Clinic to verify a clinical diagnosis of Jeune syndrome. Clinical Exome sequencing revealed two variants (compound heterozygosity) in the ORC6 gene: c.2T&gt;C(p.Met1Thr) and c.449+5G&gt;A. In silico analysis showed the pathogenicity of these two mutations and predicted a decrease in donor splicing site strength for c.449+5G&gt;A. An in vitro minigene assay demonstrated that variant c.449+5G&gt;A causes a complete skipping of exon 4 in ORC6 gene. The parents asked for urgent prenatal testing for MGS for the next pregnancy, but it ended in a miscarriage. Thus, this case report may help to prevent MGS misdiagnosis in the future. We also performed in silico and functional analyses of ORC6 mutations and developed a restriction fragment length polymorphism and haplotype-based short-tandem-repeat assay for prenatal genetic testing for MGS. These findings should elucidate MGS aetiology and improve the quality of genetic counselling for affected families.

A whole-exome capture and next-generation sequencing applied to an 11 y/o patient with a clinical history of congenital hypotonia, generalized motor and cognitive neurodevelopmental delay, severe cognitive deficit, without any identifiable Syndromic pattern, and to her parents, disclosed a de novo heterozygous pathogenic mutation, c.697_699del p.Phe233del (rs786204835)(ACMG classification PS2, PM1, PM2, PP5), harbored in the PURA gene (MIM*600473) (5q31.3), associated to Autosomal Dominant Mental Retardation 31 (MIM # 616158). We used the significant improvement in the accuracy of protein structure prediction recently implemented in AlphaFold that incorporates novel neural network architectures and training procedures based on the evolutionary, physical, and geometric constraints of protein structures. The wild-type (WT) sequence and the mutated one, missing the Phe233, were reconstructed. The predicted local Distance Difference Test (lDDT) for the PURA WT and the PURA-Phe233del showed that the occurrence of the Phe233del affects between 220-320 amino acids. The distortion in the PURA structural confor-mation in the ~5Å surrounding area after the p.Phe233del produces a conspicuous disruption of the repeat III, where the DNA and RNA helix unwinding capability occurs. PURA Protein-DNA Docking corroborated these results in silico Analysis that showed a loss of the contact of the PURA-Phe233del III repeat domain model with the DNA. Together, i) the energetic and stereochemical, ii) the hydropathic indices and polarity surfaces, and iii) the hybrid Quantum Mechanics-Molecular Mechanics (QM-MM) analyses of the PURA molecular models demarcate at the atomic resolution the specific surrounding region affected by these mutations and paves the way for future cell-based functional analysis. To the best of our knowledge, this is the first report of a de novo mutation underpinning a PURA syndrome in a Latin American patient and highlights the importance of predicting the molecular effects in protein structure using artificial intelligence algorithms and molecular and atomic resolution stereochemical analyses.

Background: Immunoglobulin A nephropathy (IgAN) is a complex autoimmune disease, and the exact pathogenesis remains to be elucidated. Methods: We conducted summary data-based Mendelian randomization (SMR) analysis and performed functional mapping and annotation using FUMA to explore genetic loci that are po-tentially involved in the pathogenies of IgAN. Both analyses used summarized data of a recent genome-wide association study (GWAS) on IgANs, which included 477,784 Europeans (15,587 cases and 462,197 controls) and 175,359 East Asians (71 cases and 175,288 controls). We performed separate SMR analysis using CAGE and GTEx eQTL data. Results: Using the CAGE eQTL data, our SMR analysis identified 32 probes tagging 25 unique genes that were pleiotropically/potentially causally associated with IgAN, with the top three probes being ILMN_2150787 (tagging HLA-C, PSMR=2.10×10-18), ILMN_1682717 (tagging IER3, PSMR=1.07×10-16) and ILMN_1661439 (tagging FLOT1, PSMR=1.16×10-14). Using GTEx eQTL data, our SMR analysis identified 24 probes tagging 24 unique genes, with the top three probes being ENSG00000271581.1 (tagging XXbac-BPG248L24.12, PSMR=1.44×10-10), ENSG00000186470.9 (tagging BTN3A2, PSMR=2.28×10-10), and ENSG00000224389.4 (tagging C4B, PSMR=1.23×10-9). FUMA analysis identified 3 independent, significant and lead SNPs, 2 genomic risk loci and 39 genes. Conclusion: We identified many genetic variants/loci that are potentially involved in the patho-genesis of IgAN.

Background: Alcohol consumption is associated with the development of cardiovascular diseases, cancer, and liver disease. The biological mechanisms are still largely unclear. Here, we aimed to use an agnostic approach to identify phenotypes mediating the effect of alcohol on various diseases. Methods: We performed an agnostic association analysis between alcohol consumption (red, and white wine, beer/cider, fortified wine, and spirits) with over 7,800 phenotypes from the UK biobank comprising 223,728 participants. We performed Mendelian randomisation analysis to infer causality. We additionally performed a Phenome-wide association analysis and a mediation analysis between alcohol consumption as exposure, traits in causal relationship with alcohol consumption as mediators, and various diseases as outcome. Results: Of 45 traits in association with alcohol consumption, 20 were in causal relationship with alcohol consumption. Gamma glutamyltransferase (GGT; β=9.44; CI,5.94-12.93; Pfdr=9.04×10-7), mean sphered cell volume (β=0.189; CI,0.11-0.27; Pfdr=1.00×10-4), mean corpuscular volume (β=0.271; CI,0.19-0.35; Pfdr=7.09×10-10) and mean corpuscular haemoglobin (β=0.278; CI,0.19-0.36; Pfdr=1.60×10-6) showed the strongest causal relationships. We also identified GGT and physical activity as mediators causing liver cirrhosis and alcohol dependence. Conclusion: Our study provides evidence of causality between alcohol consumption and 20 traits and a mediation effect for physical activity on health consequences of alcohol consumption.

Sweet cucumber (Solanum muricatum) sect. Basarthrum, is a neglected horticultural crop native of the Andean region. It is naturally distributed very close to potatoes (Solanum sect. Petota) and tomatoes (Solanum sect. Lycopersicon), two groups of high economic importance. To date, molecular tools for this crop are still undetermined. We here obtained the first complete chloroplast (cp) genome of sweet cucumber and compared with seven Solanaceae species. Pair-end clean reads were obtained by PE 150 library and the Illumina HiSeq 2500 platform. The complete cp genome of S. muricatum had a 155,681 bp with typical quadripartite structure, containing a large single copy (LSC) region (86,182 bp) and a small single-copy (SSC) region (18,360 bp), separated by two inverted repeat (IR) regions (25,568 bp). The annotation of chloroplast genome predicted 88 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes, and one pseudogene. A total of 48 perfect microsatellites were identified, divided in mononucleotide repeats (32), followed by tetranucleotide (6) and dinucleotides (5). SSRs with trinucleotides repeats (3), pentanucleotide (1) and hexanucleotide (1) repeats motifs in these genomes were identified in lower quantity. Most of these repeats were distributed in the noncoding regions. Whole chloroplast genome comparison with the other seven Solanaceae species revealed that the small and large single copy regions showed more divergence than inverted regions. Finally, phylogenetic analysis resolved that S. muricatum is a sister species to members of sections Petota + Lycopersicum + Etuberosum. This study reports for the first time the genome organization, gene content, and structural features of the cp genome of S. muricatum. Also, this study may provide the basis for evaluating genetic diversity within Solanum, and will be useful to examine the evolutionary processes in sweet cucumber landraces.

Non-vertebrate species represent about ~95% of known metazoan (animal) diversity. They remain to this day relatively unexplored genetically, but understanding their genome structure and function is pivotal for expanding our current knowledge of evolution, ecology and biodiversity. Following the continuous improvements and decreasing costs of sequencing technologies, many genome assembly tools have been released, leading to a significant amount of genome projects being completed in recent years. In this review, we examine the current state of genome projects of non-vertebrate animal species. We present an overview of available sequencing technologies, assembly approaches, as well as pre and post-processing steps, genome assembly evaluation methods, and their application to non-vertebrate animal genomes.

Standardized gene nomenclature supports unambiguous communication and identification of the scientific literature associated with genes. To support the increasing number of annotated genomes that are now available for comparative studies, gene nomenclature authorities coordinate the assignment of approved gene names that can be readily propagated across species. Theofanopoulou et al. (Theofanopoulou et al. 2021) propose a new nomenclature for the genes encoding oxytocin and arginine vasopressin and their receptors. Rather than changing to a different nomenclature system, we propose minor updates to the current approved nomenclature of these vertebrate genes to better reflect their evolutionary history. We call on authors, journal editors and reviewers to help support communication and indexing of gene-related publications by working with existing gene nomenclature committees and ensuring that standardized gene nomenclature is routinely used.

The classical human satellite DNAs, also referred to as human satellites 1, 2 and 3 (HSat1, HSat2, HSat3, collectively HSat1-3) constitute the largest individual arrays of tandemly repeated DNA sequences in the genome. Even though they were among the first human DNA sequences to be isolated and characterized at the dawn of molecular biology, HSat1-3 have been left behind in the genomics era and remain among the most enigmatic sequences in the human genome. Although HSat1-3 total roughly 3% of the genome on average, they were almost entirely missing from the human genome reference assembly for 20 years. Recently, the Telomere-to-Telomere Consortium produced the first truly complete assembly of a human genome, including the enormous HSat1-3 arrays, opening them up for a new wave of discovery. Towards this end, here, I provide an account of the history and current understanding of HSat1-3 genomics, evolution, and roles in disease.

Mandevilla (Apocynaceae) is a greatly appreciated genus in the world ornamental market. In this study, we attempted to address the poor genetic knowledge and the huge taxonomic gaps existing in this genus by analyzing the germplasm of 55 accessions. After cytometrically determining the triploid genome size (1,512.64 Mb) of a reference sample (variety “Mandevilla 2001”), the plastidial genome (cpDNA, 0,18 Mb) and a draft of the nuclear genome (nuDNA, 207 Mb) were assembled. While cpDNA was effective in reconstructing the phylogenesis of the Apocynaceae family based on a DNA superbarcoding approach, the nuDNA assembly length was found to be only 41% of the haploid genome size (506 Mb, predicted based on the K-mer frequency distribution). Its annotation enabled the prediction of thousands of genes, of which 5,275 resulted in full CDS. Among them, we identified nine genes whose orthologs (in Catharanthus roseus) encode enzymes involved in the biosynthesis of monoterpene indole alkaloids (MIAs), including catharanthine, tabersonine and vincadifformine. The nuclear genome draft was also useful to develop a highly informative (mean PIC=0.62) set of 23 SSR markers that was validated on the Mandevilla collection. These results were integrated with cytometric measurements, nuclear ITS1 haplotyping and chloroplast DNA barcoding analyses to assess the origin, divergence and relationships existing among the 55 accessions object of the study. As expected, based on the scarce information available in the literature, the scenario was extremely intricate. A reasonable hypothesis is that most of the accessions represent interspecific hybrids sharing the same species as maternal parent (i.e., Mandevilla sanderi).

The article continues the author's publications about the matrix-tensor study of universal rules of stochastic (probabilistic) organization of long single-stranded DNA sequences in eukaryotic and prokaryotic genomes. The author reveals that corresponding matrices of probabilities of n-plets in n-textual representations of each genomic DNA are numerically interrelated each with other in such algebraic form, which has analogies with formalisms of the known tensor-matrix theory of digital antenna arrays. These arrays combine many separate antennas into a single coordinated ensemble with unique emergent properties, due to which antenna arrays are widely used in devices of medicine, astrophysics, avionics, etc. The noted analogies allow putting forward the author's hypothesis that the stochastic organization of genomic DNAs is connected with bio-antenna arrays. From the point of view of this hypothesis, many known facts about using principles of antenna arrays in inherited physiological phenomena are collected in a single grouping with genomic DNAs. This new topic about the biological meaning of profitable properties of antenna arrays includes problems of biological evolution, the origin of the genetic code, regenerative medicine, and the development of algebraic biology. These issues are discussed jointly with the author's results of quantum information analysis of stochastic features of genomic DNAs.

Conventional phenotyping breeding approaches for vegetable crops like Solanaceae, Bulb, Root crops, have made a significant contribution by developing many varieties. Despite this, conventional phenotyping approaches are not sufficient due to the longer time taken to develop a variety, low genetic gain, environmental factors and some other externalities that affect the phenotype-based selection. To address the challenges of conventional phenotype, a new recent method of high throughput phenotyping (HTP) is considered a promising tool. The development of high-throughput phenotyping technology began in the preceding decade as advancements in sensor, computer vision, automation, and advanced machine learning technologies. HTP platforms are being utilized to undertake non-destructive assessments of the complete plant system in a range of crops. HTP provides the precise measurements and suggests the collection of high-quality and accurate data which is necessary for standardizing phenotyping for the collection of genetic dissection and genomic assisted breeding such as genome-wide association studies (GWAS), linkage mapping, marker-assisted selection (MAS), genomic selection (GS). The remainder of this chapter discusses how high-throughput phenotyping technologies can be used in genomic-assisted breeding for vegetable crops

The Omicron variant of SARS-CoV-2 is the latest pandemic lineage causing COVID-19. Despite having a vaccination rate ≥ 85% Ecuador recorded a high incidence of Omicron from December 2021 to March 2022. Since Omicron emerged it is evolving into multiple sublineages with distinct prevalence in different regions. In this work, we use all Omicron sequences from Ecuador available at GISAID until March 2022 and the software Nextclade and Pangolin to identify which lineages circulate in this country. We detected 12 different sublineages (BA.1, BA.1.1, BA.1.1.1, BA.1.1.14, BA.1.1.2, BA.1.14, BA.1.15, BA.1.16, BA.1.17, BA.1.6, BA.2, BA.2.3), which has been reported in Africa, America, Europe, and Asia suggesting multiple introduction events. Sublineages BA.1.1 and BA.1 were the most prevalent. Genomic surveillance must continue to evaluate the dynamic of current sublineages, early introduction of new ones and vaccine efficacy against evolving SARS-CoV-2.

Abstract: Due to their rapid evolution and their impact on health care, be-ta-lactamases, protein degrading beta-lactam antibiotics, are used as generic model of protein evolution. Therefore, we investigated the mutation effects in two distant beta-lactamases TEM-1 and CTX-M-15. Interestingly we found a site with a complex pattern of genetic interactions. Mutation G251W in TEM-1 is in-activating the protein’s function just as the reciprocal mutation W251G in CTXM-15. Phylogenetic analysis revealed that mutation G has been entrenched in TEM-1 background: while rarely observed throughout the phylogeny it is es-sential in TEM-1. Using a rescue experiment in TEM-1 G251W mutant, we could identify sites that alleviates the deviation from G to W. While few of these muta-tions could potentially involve local interactions, most of them were found on distant residues in the 3D structure. Many well-known mutations having an impact on protein stability, such as M182T, were recovered. Our results there-fore suggest that entrenchment of an amino acid may rely on diffuse interactions among multiple sites with a major impact on protein stability.

Cattle spread throughout the American continent during the colonization years, originating creole breeds that adapted to a wide range of climate conditions. Population of creole cattle in Peru is decreasing mainly due to the introduction of more productive breeds in recent years. During the last 15 years, there have been a significant progress on cattle genomics. However, little is known about the genetics of the Peruvian creole cattle (PCC) even though its importance to (i) improve productivity in the Andean region, (ii) agricultural labor, and (iii) cultural traditions. In addition, the origin and phylogenetic relationship of the PCC is still unclear. In order to promote the conservation of the PCC, we sequenced for the first time the mitochondrial genome of a creole bull from the highlands of Arequipa, which also possessed exceptional fighting skills and was employed for agricultural tasks. The total mitochondrial genome sequence is 16,339 bp in length with the base composition of 31.43 % for A, 28.64 % for T, 26.81 % for C, and 13.12 % for G. It contains 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes and a control region. Among the 37 genes, 28 were positioned on the H-strand and nine were positioned on the L-strand. The most frequently used codons were CUA (Leucine), AUA (Isoleucine), AUU (Isoleucine), AUC (Isoleucine), y ACA (Threonine). Maximum likelihood reconstruction using complete mitochondrial genome sequences clearly demonstrated that the PCC is strongly related to native African breeds, giving insights into the ancestry of PCC. The annotated mitochondrial genome of PCC would serve as an important genetic data set for further breeding work and conservation strategies.

The cytosolic PNGase (peptide:N-glycanase; Png1 in yeast; NGLY1/Ngly1 in human/mice), also known as peptide-N4-(N-acetyl-beta-glucosaminyl)-asparagine ami-dase, is a well-conserved deglycosylation enzyme (EC 3.5.1.52) which catalyzes the non-lysosomal hydrolysis of an N(4)-(acetyl-β-D-glucosaminyl) asparagine residue into N-acetyl-β-D-glucosaminylamine and a peptide containing an aspartate residue. This enzyme (NGLY1) plays essential roles in clearance of misfolded or unassembled gly-coproteins through a process named ER-associated degradation (ERAD). Accumulating evidence also points out that NGLY1 deficiency can cause an autosomal recessive hu-man genetic disorder associated with abnormal development and congenital disorder of deglycosylation. In addition, the loss of NGLY1 can affect multiple cellular pathways, including but not limited to NFE2L1 pathway, Creb1/Atf1-AQP pathway, BMP path-way, AMPK pathway, and SLC12A2 ion transporter, which might be the underlying reasons for a constellation of clinical phenotypes of NGLY1 deficiency. The current comprehensive review indeed uncovers the detailed NGLY1’s structure and its im-portant roles for participation in ERAD, involvement in CDDG and potential treatment for NGLY1 deficiency.

Chemically modified forms of tetraiodothyroacetic acid (tetrac), an L-thyroxine derivative, have been shown to exert their anticancer activity at plasma membrane integrin αvβ3 of tumor cells. Via a specific hormone receptor on the integrin, tetrac-based therapeutic agents modulate expression of genes relevant to cancer cell proliferation, survival and energy metabolism. P-bi-TAT, a novel bivalent tetrac-containing synthetic compound has anticancer activity in vitro and in vivo against Glioblastoma Multiforme (GBM) and other types of human cancers. In the current study, microarray analysis was carried out on a primary culture of human GBM cells exposed to P-bi-TAT (10-6 tetrac equivalent) for 24 h. P-bi-TAT significantly affected expression of a large panel of genes implicated in cancer cell stemness, growth, survival, and angiogenesis. Recent interest elsewhere in ATP synthase as a target in GBM cells caused us to focus attention on expression of genes involved in energy metabolism. Significantly downregulated transcripts included multiple energy metabolism-related genes: electron transport chain genes ATP5A1 (ATP synthase 1), ATP51, ATP5G2, COX6B1 (cytochrome c oxidase subunit 6B1), NDUFA8 (NADH dehydrogenase [ubiquinone] FA8), NDUFV2I and other NDUF genes. The NDUF and ATP genes are also relevant to control of oxidative phosphorylation and transcription. Qualitatively similar actions of P-bi-TAT on expression of energy metabolism-linked genes were also detected in established human GBM and pancreatic cancer cell lines. In conclusion, acting at αvβ3 integrin, P-bi-TAT caused downregulation in human cancer cells of expression of a large number of genes involved in electron transport and oxidative phosphorylation. These observations suggest that cell surface thyroid hormone receptors on αvβ3 regulate expression of genes relevant to tumor cell stemness and energy metabolism.

Pharmacogenomics (PGx) is the study of how variations of the patient genetic makeup could influence the selection and dosing of medications, including urate-lowering therapies (ULTs). Patients with gout present with several comorbidities, warranting the use of several long-term medications that increase their pill burden and the risk of adverse drug events. PGx testing can identify individuals who are more or less likely to benefit from a given treatment and ultimately improve adherence and reduce polypharmacy. Gout management involves using various ULTs associated with specific genetic variations that could modulate medication safety and efficacy. Individuals carrying the HLA-B*58:01 variant are at a higher risk of serious and life-threatening skin reactions from taking allopurinol. Besides, racial disparities in HLA-B distribution warrant testing for its status in high-risk populations, specifically some Asian subgroups and African Americans. G6PD-deficient individuals are prone to develop hemolytic anemia and methemoglobinemia with pegloticase and probenecid use. Patients with the less active form of the drug-metabolizing CYP2C9 are at higher risk for NSAID-related upper gastrointestinal (GI) bleeding. The drug-gene pair interactions among various gout therapies were also reviewed. Allopurinol-related genes included AOX, ABCG2, and SLC22A12. Febuxostat-related genes included UGT1A1. Benzbromarone- and probenecid- related genes included SLC22A12 and ABCB1. Colchicine-related genes included CYP2D6, ABCB1, rs6916345 G>A, and SEPHS1. Finally, the corticosteroid-related gene included HCG22; and the anakinra-related gene included IL1RA. This review synthesizes the most contemporary level of evidence for using PGx in gout pharmacotherapy, as referenced in the CPIC guidelines, FDA, PharmGKB, and current literature.

The genetic aetiology of familial glioma is largely unknown. To gain further insight into the role of rare disruptive variants we performed whole exome sequencing (WES) of 17 glioma families, identifying two families where loss of function mutations in dynein axonemal heavy chain 11 (DNAH11) co-segregated with glioma (DNAH11 p.Ser1470ArgfsTer6 and p.Thr3900Lys) and single families implicating Telomere Replication Complex Component 1 (CTC1 p.Leu1058ProfsTer32) and filamin alpha (FLNA p.Ser1356Phe). We complemented linkage analysis by WES of 150 additional unrelated familial glioma cases and conducting gene burden tests leveraging WES data on 691 cases of glioma from TCGA. Our analysis does not provide support for a hitherto unidentified major cancer gene for glioma but highlights a series of novel candidate predisposition variants/genes, worthy of further investigation to shed light on glioma risk.

(1) Background: The goal of our study was to identify SNPs, metabolic pathways (KEGG), and gene ontology (GO) terms significantly associated with calving and workability; (2) Methods: Based on the EuroGenomics reference data set, we analyzed direct (DCE) and maternal (MCE) calving ease, direct (DSB), and maternal (MSB) stillbirth, milking speed (MSP), and temperament (TEM). We estimated SNP effects using a multi-SNP mixed-model. Further, SNP positions were mapped to genes, and GO terms/KEGG pathways of the corresponding genes were assigned. The estimation of GO term/KEGG pathway effects was based on a mixed-model using SNP effects as dependent variables; (3) Results: The number of significant SNPs comprised 59 for DCE, 25 for DSB and MSP, 17 for MCE and MSB, and 7 for TEM. Significant KEGG pathways were found for MSB (2), TEM (2), and MSP (1), while 11 GO terms were significant for MSP, 10 for DCE, 8 for DSB and TEM, 5 for MCE, and 3 for MSB.; (4) Conclusions: From the perspective of a better understanding of the genomic background of the phenotypes, traits with low heritabilities suggest that the focus should be moved from single genes to metabolic pathways or gene ontologies significant for the phenotype.

The experiments described in this research article were designed to test the effect of rare variants into genomic prediction in dairy cattle. Common polymorphisms are able to explain only a small proportion of the underlying genetic variation of complex phenotypes. Variants representing functional mutations with large effects on complex phenotypes are expected to be rare due to natural (humans) or artificial (livestock) selection pressure. Therefore, it is important to check whether the use of rare variants could increase the accuracy of ranking of animals by providing the tool for more precise differentiation among the bulls with high additive genetic merit. The goal of our study was to verify whether including rare variants in a genomic selection model allows for a more accurate description of the additive genetic background of traits under selection in dairy cattle. We used the linear mixed model for comparison SNP estimates for Holstein-Friesian cattle of the two data sets – a set containing only single nucleotide polymorphisms defined by minor allele frequency ≥ 0.01, which is routinely used in the Polish genomic evaluation system (46,216 SNPs), and a set containing SNPs selected based only on the call rate (54,378 SNPs). Based on the SNP estimates we also calculated DGV and GEBV and compared them between both data sets. In all the analyses we used production, fertility, conformation and udder health traits. We also assessed the time required for the two most computationally demanding components of genomic selection: preparing genotype data, and estimation of SNP effects between those two data sets. The results of our study indicated that the analysis including rare variants resulted in changes in the individual ranking of the top 100 male and female candidates, but had no effect on the outcome of the quality of EBV prediction as expressed by the Interbull validation test.

Plant geneticists and breeders have used marker technology since the 1980s in quantitative trait locus (QTL) identification. Marker-assisted selection is effective for large-effect QTL but has been challenging to use with quantitative traits controlled by multiple minor effect alleles. Therefore, genomic selection (GS) was proposed to estimate all markers simultaneously, thereby capturing all their effects. However, breeding programs are still struggling to identify the best strategy to implement it into their programs. Traditional breeding programs need to be optimized to implement GS effectively. This review explores the optimization of breeding programs for variety release based on aspects of the breeder’s equation. Optimizations include reorganizing field designs, training populations, increasing the number of lines evaluated, and leveraging the large amount of genomic and phenotypic data collected across different growing seasons and environments to increase heritability estimates, selection intensity, and selection accuracy. Breeding programs can leverage their phenotypic and genotypic data to maximize genetic gain and selection accuracy through GS methods utilizing multi-trait and, multi-environment models, high-throughput phenotyping, and deep learning approaches. Overall, this review describes various methods that plant breeders can utilize to increase genetic gains and effectively implement GS in breeding .

Wheat (Triticum aestivum L.) breeding programs can take over a decade to release a new variety. However, new methods of selection such as genomic selection (GS) must be integrated to decrease the time it takes to release new varieties to meet the demand of a growing population. The implementation of GS into breeding programs is still being explored, with many studies showing its potential to change wheat breeding through achieving higher genetic gain. In this review, we explore the integration of GS for a wheat breeding program by redesigning the traditional breeding pipeline to implement GS. We propose implementing a two-part breeding strategy by differentiating between population improvement and product development. The implementation of GS in the product development pipeline can be integrated into most stages and can predict within and across breeding cycles. Additionally, we explore optimizing the population improvement strategy through GS recurrent selection schemes to reduce crossing cycle time and significantly increase genetic gain. The recurrent selection schemes can be optimized for parental selection, maintenance of genetic variation, and optimal cross-prediction. Overall, we outline the ability to increase the genetic gain of a breeding program by implementing GS and a two-part breeding strategy.

Branched-chain amino acid (BCAA) levels are associated with skeletal muscle cross-sectional area (CSA). Serum BCAA levels are enhanced by whey protein supplementation (WPS), and evidence in clinical populations suggests an association of single nucleotide polymorphisms (SNPs) with BCAA metabolite levels. It is not known whether the same SNPs are associated with the ability to catabolise BCAAs from exogenous sources, such as WPS. The present study investigated whether possessing a higher number of alleles associated with increased BCAA metabolites correlates with muscle fiber CSA of m. vastus lateralis in physically active participants, and whether any relationship is enhanced by WPS. Endurance-trained participants (n=75) were grouped by self-reported habitual WPS consumption and genotyped for 5 SNPs (PPM1K rss1440580, APOA5 rs2072560, CBLN1 rs1420601, DDX19B rs12325419, TRMT61A rs58101275). Body mass, BMI and fat percentage were significantly lower and muscle mass higher in the WPS group (n=22) compared to Non-WPS (n=53). The number of BCAA-increasing alleles was correlated with fiber CSA in the WPS group (r=0.75, p&lt;0.0001) and was stronger for fast-twitch fibers (p=0.001) than slow-twitch fibers (p=0.048). Similar results remained when corrected for multiple covariates (age, physical activity, meat and dairy intake). No correlation was found in the Non-WPS group. This study presents novel evidence of a positive relationship between BCAA-increasing alleles and muscle fiber CSA in athletes habitually consuming WPS. We suggest that a high number of BCAA-increasing alleles improves the efficiency of WPS by stimulation of muscle protein synthesis, and contributes to greater fiber CSA.

Duchenne muscular dystrophy (DMD) is an X-linked recessive neuromuscular disorder with a prevalence of approximately 1 in 3,500-5,000 males. DMD manifests as childhood-onset muscle degeneration, followed by loss of ambulation, cardiomyopathy, and death in early adulthood due to a lack of functional dystrophin protein. Out-of-frame mutations in the dystrophin gene are the most common underlying cause of DMD. Gene editing via the clustered regularly interspaced short palindromic repeats (CRISPR) system is a promising therapeutic for DMD, as it can permanently correct DMD mutations and thus restore the reading frame, allowing for the production of functional dystrophin. The specific mechanism of gene editing can vary based on a variety of factors such as the number of cuts generated by CRISPR, the presence of an exogenous DNA template, or the current cell cycle stage. CRISPR-mediated gene editing for DMD has been tested both in vitro and in vivo, with many of these studies discussed herein. Additionally, novel modifications to the CRISPR system such as base or prime editors allow for more precise gene editing. Despite recent advances, limitations remain including delivery efficiency, off-target mutagenesis, and long-term maintenance of dystrophin. Further studies focusing on safety and accuracy of the CRISPR system are necessary prior to clinical translation.

Transcription factors (TFs) are proteins that recognize specific DNA fragments in order to decode the genome and ensure its optimal functioning. TFs work at the local and global scales by specifying cell type, cell growth and death, cell migration, organization and timely tasks. We investigate the structure of DNA-binding motifs with the theory of finitely generated groups. The DNA ‘word’ in the binding domain -the motif- may be seen as the generator of a finitely generated group Fdna on four letters, the bases A, T, G and C. It is shown that, most of the time, the DNA-binding motifs have subgroup structure close to free groups of rank three or less, a property that we call ‘syntactical freedom’. Such a property is associated to the aperiodicity of the motif when it is seen as a substitution sequence. Examples are provided for the major families of TFs such as leucine zipper factors, zinc finger factors, homeo-domain factors, etc. We also discuss the exceptions to the existence of such a DNA syntactical rule and their functional role. This includes the TATA box in the promoter region of some genes, the single nucleotide markers (SNP) and the motifs of some genes of ubiquitous role in transcription and regulation.

: Using line × tester analysis, the current research analyses parental genotypes and their combinations in normal conditions and identifies the genes influencing yield characteristics. In the present study, 15 diverse genotypes, including 10 lines, 5 testers, and 50 F1s hybrids, were evaluated for 13 morphological and 2 biochemical traits. A suitable location was taken to study the effect of 15 characters. The results exposed that ability mean squares were significant for all studied additive and non-additive components. In this direction, the general combining ability of PBW-343, DBW-39, K-402, K-1317, KRL-210, and K-68 were higher than the remaining parents. For morphological traits like yield, the top five crosses were described based on SCA effects, namely, HD-3086 × HD-3171, K-402 × K-9107, K-1317 × K-9107, HD-2967 × K-0307 and K-402 × K-68 in F1 generation. In addition, the high value of heritability was estimated for plant height (77.32%), spike length (32.26%), biological yield/plant (59.52%), and grain yield/plant (68.76%). However, the moderate values of heritability were estimated for days to maturity (22.78%) and phenol color reaction (18.00%). The higher genetic advance was not found for recorded characters; however, a moderate genetic advance was recorded for grain yield per plant (13.15%) and harvest index (11.72%). High heritability coupled with moderate genetic advance was recorded for two characters grain yield per plant and harvest index in F1 and F2 generations.

B lymphocyte development has two DNA recombination processes: V(D)J recombination of the immunoglobulin (Igh) gene variable region and class switching of the Igh constant regions from IgM to IgG, IgA, or IgE. V(D)J recombination is required for successful maturation of B cells from pro-B to pre-B to immature-B and then to mature B cells in the bone marrow. CSR occurs outside of the bone marrow when mature B cells migrate to peripheral lymphoid organs, such as spleen and lymph nodes. Both V(D)J recombination and CSR depend on an “open chromatin” state that makes DNA accessible to specific enzymes, recombination activating gene (RAG), and activation-induced cytidine deaminase (AID). Acetyltransferases GCN5 and PCAF possess redundant functions acetylating histone H3 lysine 9 (H3K9). Here, we generated a mouse model that lacks both GCN5 and PCAF in B cells. We found that double-deficient mice possess low levels of mature B cells in the bone marrow and peripheral organs, an accumulation of pro-B cells in bone marrow, and reduced CSR levels. We conclude that both GCN5 and PCAF are required for B cell development in vivo.

The domestic pig (Sus scrofa) is a species representative of the Suina, one of the four suborders within Cetartiodactyla. In this paper, we reported our analysis of the pig TRG locus in comparison with the loci of species representative of the Ruminantia, Tylopoda and Cetacea suborders. The pig TRG genomic structure reiterates the peculiarity of the organization of Cetartiodactyla loci in TRGC “cassettes”, each containing the basic V-J-J-C unit. Eighteen genes arranged in four TRGC cassettes, form the pig TRG locus. All the functional TRG genes were expressed, and the TRGV genes preferentially rearrange with the TRGJ genes within their own cassette, which correlates the diversity of the gamma-chain repertoire with the number of cassettes. Among them, the TRGC5, located at the 5’ end of the locus, is the only cassette that retains a marked homology with the corresponding TRGC cassettes of all the analyzed species. The preservation of the TRGC5 cassette for such a long evolutionary time presumes a highly specialized function of its genes, which could be essential for the survival of species. Therefore, the maintenance of this cassette in pigs confirms that it is the most evolutionarily ancient within Cetartiodactyla, and it has undergone a process of duplication to give rise to the other TRGC cassettes in the different artiodactyl species in a lineage-specific manner.

Familial colorectal cancer (CRC) is only partially explained by known germline predisposing genes. We performed whole genome sequencing in 15 Polish families of many affected individuals, without mutations in known CRC predisposing genes. We focused on loss-of-function variants and functionally characterized them. We identified a frameshift variant in the CYBA gene (c.246delC) in one family and a splice site variant in the TRPM4 gene (c.25-1 G>T) in another family. While both variants were absent or extremely rare in gene variant databases, we identified four additional Polish familial CRC cases and two healthy elderly individuals with the CYBA variant (odds ratio 2.46, 95% confidence interval 0.48-12.69). Both variants led to a premature stop codon and to a truncated protein. Functional characterization of the variants showed that knockdown of CYBA or TRPM4 depressed generation of reactive oxygen species (ROS) in LS174T and HT-29 cell lines. Knockdown of TRPM4 resulted in decreased MUC2 protein production. CYBA encodes a component in the NADPH oxidase system which generates ROS and controls, e.g., bacterial colonization in the gut. Germline CYBA variants are associated with early onset inflammatory bowel disease, supported with experimental evidence on loss of intestinal mucus barrier function due to ROS deficiency. TRPM4 encodes a calcium-activated ion channel, which in a human colonic cancer cell line controls calcium-mediated secretion of MUC2, a major component of intestinal mucus barrier. We suggest that the gene defects in CYBA and TRPM4 mechanistically involve intestinal barrier integrity through ROS and mucus biology, which converges in chronic bowel inflammation.

The current and projected climate change that is represented by increasing temperatures, humidity levels and irregular rainfall patterns, promotes the occurrence of preharvest sprouting (PHS) in wheat. PHS results in significant economic losses, globally, which necessitates the need for high-yielding cultivars with increased PHS tolerance, hence this study was conducted. The current study evaluated two doubled-haploid (DH) wheat populations of Tugela-Dn × Elands and Elands × Flamink across six environments in the Free State Province of South Africa to select genotypes with increased PHS tolerance and further map the underlying loci. Significant effects of DH lines (194) and environments (6) were observed for PHS tolerance. The results of this study validate previous findings that PHS is only expressed when environmental conditions are conducive. Quantitative trait loci (QTL) mapping using single nucleotide polymorphism (SNP) and silicoDArT markers revealed three additive QTL with major effects on chromosomes 5B and 7B, and these QTL were detected more than once, when conditions were favourable. These QTL explained a phenotypic variation (PVE) varying between 10.08% and 20.30% (LOD = 2.73 – 3.11). About 16.50% of DH lines performed to the level of Elands (the PHS-tolerant parent) and are recommended for further selection in a pre-breeding or breeding programme. The findings of the study are expected to facilitate the on-going breeding efforts for PHS tolerance in winter wheat.

Adipor1^tm1Dgen and Mfrp^rd6 mutant mice share similar eye disease characteristics. Previously, studies established a functional relationship of ADIPOR1 and MFRP proteins in maintaining retinal lipidome homeostasis and visual function. However, the independent and/or interactive contribution of both the genes to similar disease phenotypes, including fundus spots, decreased axial length and photoreceptor degeneration has yet to be examined. We performed a gene-interaction study where homozygous Adipor1^tm1Dgen and Mfrp^rd6 mice were bred together and the resulting doubly heterozygous F1 offspring were intercrossed to produce 210 F2 progeny. Four-month-old mice from all nine genotypic combinations obtained in the F2 generation were assessed for white spots by fundus photo documentation, for axial length by caliper measurements, and for photoreceptor degeneration by histology. Two-way factorial ANOVA was performed to study individual as well as gene interaction effects on each phenotype. Here, we report the first observation of reduced axial length in Adipor1^tmlDgen homozygotes. We show that while Adipor1 and Mfrp interact to affect spotting and degeneration, they act independently to control axial length, highlighting the complex functional association between these two genes. Further examination of the molecular basis of this interaction may help in uncovering mechanisms by which these genes perturb ocular homeostasis.

Recent climatic disasters (dzudzs) and uncontrolled massive breeding endangered valuable ge-netic resources of Mongolian goats labored by five thousand years of evolution in extremes of Al-tai mountains and Gobi deserts. Meanwhile, Mongolian goats has never been characterized by genomic data. We used Illumina Goat SNP50 to investigate phylogenic relationships and genetic risks in 5 local Mongolian populations: Erchim (N=37) from geographically remote Darkhat Valley in the northern mountains, Ulgii Red (N=35) and Dorgon (N=28) from the western mountain region, Buural (N=34) from the western steppe and mountain areas and Gobi Gurvan Saikhan (N=33) from semi-arid steppe area in the South Gobi region. ROH analyses, estimated populations sizes and Fst values showed South Gobi and Darkhat Valley goats to be of a high risk of inbreed-ing that however appear to be of distinct origin (artificial selection vs geographical isolation). In-dices of genetic differentiations between Mongolian goats were relatively low compared to Euro-pean breeds. Meanwhile, we suggest direct comparison is not fully appropriate given that Mongo-lian breeds were not subjected to intense selection. Darkhat Valley goats were clearly differentiat-ed from other Mongolian breeds according to various types of analyses. Phylogenetic relationships within XXX breeds of the rest of the world plotted Mongolian goats between some Russian (Altai, Orenburg) and Chinese populations (Nanjiang and Qinggeli). Thus, present study (i) highlights demographic history Mongolian goats and (ii)provides unified SNP-data called for support deci-sions in conservational genetics. Finally, our work (iii) raises a question of how exactly these data should be compared to make objective choices.

Avian sex chromosomes evolved after the divergence birds and crocodilians from their common ancestor, so they are much younger than the better-studied chromosomes of mammals. It has long been recognized that there may have been several stages to the evolution of avian sex chromosomes. For example, the CHD1 undergoes recombination in paleognaths but not neognaths. Genome assemblies have suggested there may be variation in the timing of barriers to recombination among Neognathae, but there remains little understanding of the extent of this variability. Here, we look at partial sequences of ATP5F1A, which is on the avian Z and W chromosomes. It is known that recombination of this gene has independently ceased in Galliformes, Anseriformes, and at least five neoavian orders, but whether there are other independent cessations of recombination among Neoaves is not understood. We used a combination of data extracted from published chromosomal-level genomes with data collected using PCR and cloning to identify Z and W copies in 22 orders. Our results suggest there may be at least 19 independent cessations of recombination within Neognathae, and 3 clades that may still be undergoing recombination (or have only recently ceased recombination). Analyses of ATP5F1A protein sequences revealed an increased amino acid substitution rate for W chromosome gametologs, suggesting relaxed purifying selection on the W chromosome. Supporting this hypothesis, we found that the increased substitution rate was particularly pronounced for buried residues, which are expected to be more strongly constrained by purifying selection. This highlights the dynamic nature of avian sex chromosomes, and that this level of variation among clades means they should be a good system to understand sex chromosome evolution.

Heteromorphic sex chromosomes are rarely found in plants. They were observed only in 47 species from phylogenetically distant families, suggesting that the evolution of sex chromosomes was independent in these species. It was shown that DNA repeat sequences are one of the major factors driving sex chromosomes evolution, and an accumulation or elimination of the repetitive DNA elements are closely linked with the formation of differences in the sex chromosomes. The goal of this study was to characterize the transposon composition in male and female plants of Cannabis sativa L., Humulus lupulus L. and Humulus japonicus Siebold &amp; Zucc. For the first time, the male and female genomes of H. japonicus as well as male genomes of H. lupulus and C. sativa have been sequenced (there were no open data about them). The analysis of genome-wide sequencing data with using Repeatexplorer2 and author’s scripts was carried out. It was shown that accumulation of Ty3-gypsy may be associated with speciation in Cannabaceae family which is the opposite of the theory of speciation throw whole-genome duplication. Moreover, the sex-specific DNA repeat clusters in C. sativa and H. japonicus were found. The analysis also revealed that the concentration of Tekay, Retand and Ikeros repeats in the Y chromosome of C. sativa is lower than in the X chromosome and the Angela concentration is higher in the Y chromosome.

Glioblastoma (GB) is the most frequent malignant brain tumor among adults and currently there is no effective treatment. It is a very aggressive tumor that grows fast and spreads through the brain causing the death of patients in 15 months. GB cells mutate frequently and generate a heterogeneous population of tumoral cells genetically distinct. Thus, the contribution of genes and signaling pathways relevant for GB progression is of great relevance. We use a Drosophila model of GB that reproduces the features of human GB, and describe the upregulation of the circadian gene cry in GB patients and in a Drosophila GB model. We study the contribution of cry to the expansion of GB cells, to the neurodegeneration caused by GB, and to premature death and determine that cry is required for GB progression. Moreover, we analyze the mechanisms that regulate cry expression by the PI3K pathway. Finally, we conclude that cry is necessary and sufficient to regulate myc expression in GB. These results contribute to the understanding of the signals that impulse GB malignancy and lethality and open novel opportunities for the treatment of GB patients.

Capirona (Calycophyllum spruceanum Benth.) belongs to subfamily Ixoroideae, one of de major lineages in the Rubiaceae family, and is an important timber tree, with origin in the Amazon Basin and has widespread distribution in Bolivia, Peru, Colombia, and Brazil. In this study, we obtained the first complete chloroplast (cp) genome of capirona from department of Madre de Dios located in the Peruvian Amazon. High-quality genomic DNA was used to construct librar-ies. Pair-end clean reads were obtained by PE 150 library and the Illumina HiSeq 2500 platform. The complete cp genome of C. spruceanum has a 154,480 bp in length with typical quadripartite structure, containing a large single copy (LSC) region (84,813 bp) and a small single-copy (SSC) region (18,101 bp), separated by two inverted repeat (IR) regions (25,783 bp). The annotation of C. spruceanum cp genome predicted 87 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, 37 transfer RNA (tRNA) genes and 01 pseudogene. A total of 41 simple sequence repeats (SSR) of this cp genome were divided into mononucleotides (29), dinucleotides (5), trinucleotides (3), and tetranucleotide (4). Most of these repeats were distributed in the noncoding regions. Whole chloroplast genome comparison with the other six Ixoroideae species revealed that the small single copy and large single copy regions showed more divergence than invert regions. Finally, phylogenetic analysis resolved that C. spruceanum is a sister species to Emmenopterys henryi, and confirms its position within the subfamily Ixoroideae. This study reports for the first time the genome organization, gene content, and structural features of the chloroplast genome of C. spruceanum, providing valuable information for genetic and evolutionary studies in the genus Calycophyllum and beyond.