SNP and indel frequencies at transcription start sites and at canonical and alternative translation initiation sites in the human genome

Relationship between Indel Variants within the JAK2 Gene and Growth Traits in Goats

Janus kinase 2 (JAK2) plays a critical role in myoblast proliferation and fat deposition in animals. Our previous RNA-Seq analyses identified a close association between the JAK2 gene and muscle development. To date, research delving into the relationship between the JAK2 gene and growth traits has been sparse. In this study, we sought to investigate the relationship between novel mutations within the JAK2 gene and goat growth traits. Herein, two novel InDel (Insertion/Deletion) polymorphisms within the JAK2 gene were detected in 548 goats, and only two genotypes were designated as ID (Insertion/Deletion) and DD (Deletion/Deletion). The results indicate that the two InDels, the del19008 locus in intron 2 and del72416 InDel in intron 6, showed significant associations with growth traits (p < 0.05). Compared to Nubian and Jianzhou Daer goats, the del72416 locus displayed a more pronounced effect in the Fuqing breed group. In the Nubian breed (NB) group, both InDels showed a marked influence on body height (BH). There were strong linkages observed for these two InDels between the Fuqing (FQ) and Jianzhou (JZ) populations. The DD-ID diplotype was associated with inferior growth traits in chest width (ChW) and cannon circumference (CaC) in the FQ goats compared to the other diplotypes. In the NB population, the DD-DD diplotype exhibited a marked negative impact on BH and HuWI (hucklebone width index), in contrast to the other diplotypes. In summary, our findings suggest that the two InDel polymorphisms within the JAK2 gene could serve as valuable molecular markers for enhancing goat growth traits in breeding programs.

Exploiting public databases of genomic variation to quantify evolutionary constraint on the branch point sequence in 30 plant and animal species

The branch point sequence is a degenerate intronic heptamer required for the assembly of the spliceosome during pre-mRNA splicing. Disruption of this motif may promote alternative splicing and eventually cause phenotype variation. Despite its functional relevance, the branch point sequence is not included in most genome annotations. Here, we predict branch point sequences in 30 plant and animal species and attempt to quantify their evolutionary constraints using public variant databases. We find an implausible variant distribution in the databases from 16 of 30 examined species. Comparative analysis of variants from whole-genome sequencing shows that variants submitted from exome sequencing or false positive variants are widespread in public databases and cause these irregularities. We then investigate evolutionary constraint with largely unbiased public variant databases in 14 species and find that the fourth and sixth position of the branch point sequence are more constrained than coding nucleotides. Our findings show that public variant databases should be scrutinized for possible biases before they qualify to analyze evolutionary constraint.

Identification and Functional Analysis of Transcriptome Profiles, Long Non-Coding RNAs, Single-Nucleotide Polymorphisms, and Alternative Splicing from the Oocyte to the Preimplantation Stage of Sheep by Single-Cell RNA Sequencing

Numerous dynamic and complicated processes characterize development from the oocyte to the embryo. However, given the importance of functional transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms, and alternative splicing during embryonic development, the effect that these features have on the blastomeres of 2-, 4-, 8-, 16-cell, and morula stages of development has not been studied. Here, we carried out experiments to identify and functionally analyze the transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms (SNPs), and alternative splicing (AS) of cells from sheep from the oocyte to the blastocyst developmental stages. We found between the oocyte and zygote groups significantly down-regulated genes and the second-largest change in gene expression occurred between the 8- and 16-cell stages. We used various methods to construct a profile to characterize cellular and molecular features and systematically analyze the related GO and KEGG profile of cells of all stages from the oocyte to the blastocyst. This large-scale, single-cell atlas provides key cellular information and will likely assist clinical studies in improving preimplantation genetic diagnosis.

Synonymous alterations of cancer-associated Trp53 CpG mutational hotspots cause fatal developmental jaw malocclusions but no tumors in knock-in mice

Intragenic CpG dinucleotides are tightly conserved in evolution yet are also vulnerable to methylation-dependent mutation, raising the question as to why these functionally critical sites have not been deselected by more stable coding sequences. We previously showed in cell lines that altered exonic CpG methylation can modify promoter start sites, and hence protein isoform expression, for the human TP53 tumor suppressor gene. Here we extend this work to the in vivo setting by testing whether synonymous germline modifications of exonic CpG sites affect murine development, fertility, longevity, or cancer incidence. We substituted the DNA-binding exons 5-8 of Trp53, the mouse ortholog of human TP53, with variant-CpG (either CpG-depleted or -enriched) sequences predicted to encode the normal p53 amino acid sequence; a control construct was also created in which all non-CpG sites were synonymously substituted. Homozygous Trp53-null mice were the only genotype to develop tumors. Mice with variant-CpG Trp53 sequences remained tumor-free, but were uniquely prone to dental anomalies causing jaw malocclusion (p < .0001). Since the latter phenotype also characterises murine Rett syndrome due to dysfunction of the trans-repressive MeCP2 methyl-CpG-binding protein, we hypothesise that CpG sites may exert non-coding phenotypic effects via pre-translational cis-interactions of 5-methylcytosine with methyl-binding proteins which regulate mRNA transcript initiation, expression or splicing, although direct effects on mRNA structure or translation are also possible.

Relationships between the Mini-InDel Variants within the Goat CFAP43 Gene and Body Traits

The cilia- and flagella-associated protein 43 (CFAP43) gene encodes a member of the cilia- and flagellum-associated protein family. Cilia on the cell surface influence intercellular signaling and are involved in biological processes such as osteogenesis and energy metabolism in animals. Previous studies have shown that insertion/deletion (InDel) variants in the CFAP43 gene affect litter size in Shaanbei white cashmere (SBWC) goats, and that litter size and body traits are correlated in this breed. Therefore, we hypothesized that there is a significant relationship between InDel variants within the CFAP43 gene and body traits in SBWC goats. Herein, we first investigated the association between three InDel variant loci (L-13, L-16, and L-19 loci) within CFAP43 and body traits in SBWC goats (n = 1827). Analyses revealed that the L-13, L-16, and L-19 loci were significantly associated with chest depth, four body traits, and three body traits, respectively. The results of this study are in good agreement with those previously reported and could provide useful molecular markers for the selection and breeding of goats for body traits.

To Determine the Genotyping of Fc-gamma Receptor FCGR2A Polymorphism as Genetic Susceptibility to Neonatal Sepsis: A Study from a Tertiary Center of North India

Background

Neonatal sepsis term is an infection of newborns <28 days of age. It is a common cause of death in developing countries. The receptor-gamma receptor FCGR2A has been shown to be associated with neonatal sepsis. It is an activating receptor found in many cell types such as monocytes, neutrophils, macrophages, platelets, and others. The receptor has a polymorphism (single-nucleotide polymorphism rs1801274) in its gene (FCGR2A) that encodes either a histidine (H) or arginine (R) at amino acid position 131. There are many studies showing the impact of these FCGR2A polymorphisms on sepsis. Our study aims to determine the prevalence of Fc-gamma receptor FCGR2A (rs1801274) polymorphism in neonatal sepsis and control in Eastern UP populations.

Patients and Methods

We conducted a cross-sectional descriptive study of 590 patients (310 healthy individuals and 280 sepsis patients) to determine polymorphisms in the CD32A coding region in neonates. All individuals were genotyped for a variant at position 131 of the FcγRIIA gene.

Discussion

In our study, the prevalence of FcγRIIa polymorphism is more in neonates with sepsis than in noninfected neonates. It was observed that the heterozygous allele (AG) were significantly increased in septic neonates when compared to the normal.

Conclusion

Our data indicate that FcγRIIA genotyping can be used as a marker of genetic susceptibility to sepsis.

The specific DNA barcodes based on chloroplast genes for species identification of Theaceae plants

More than 600 species in over 40 genera have been identified in family Theaceae worldwide. The accurate identification of Theaceae plants can ensure the market economic order, and it plays a vital role in achieving the sustainable utilization of germplasm resources. DNA barcoding, one of the most potential species identification technologies at present, has advanced in the rapid, accurate and repetitive discrimination of species. In this study, matK + ndhF + ycf1 was observed as the optimal combined candidate gene sequence of DNA barcodes by analyzing genetic information of four single chloroplast DNA sequences, including matK, rbcL, ndhF and ycf1, as well as six combined gene sequences. Subsequently, the experiments were performed on phylogenetic analysis based on genetic distance to study the phylogenetic relationship of Theaceae plants and evaluate the species identification accuracy of matK + ndhF + ycf1. Lastly, the species-specific DNA barcodes were designed by searching the variable sites (one type of single nucleotide polymorphism sites) for the accurate identification of Camellia amplexicaulis, Franklinia alatamaha, Gordonia brandegeei and Stewartia micrantha. The previous methods of screening and testing candidate gene sequences were optimized, and innovation was made in the above methods. The process of making visual DNA barcodes was standardized. Besides, DNA barcoding technology increased the accuracy of species identification and DNA barcoding was analyzed in accordance with the theories of population genetics (e.g., neutral theory of molecular evolution). The results of the study will lay a basis for the identification and protection of Theaceae species and germplasm resources.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01175-7.

ITGB5 mutation discovered in a Chinese family with blepharophimosis-ptosis-epicanthus inversus syndrome

Blepharophimosis-ptosis-epicanthus inversus syndrome (BPES) is a rare autosomal-dominant genetic disorder, and mutations in the forkhead box L2 (FOXL2) gene are one of the major genetic causes. As this study shows, there are many patients with BPES who do not have FOXL2 mutations, as the screening results in all family members were negative. Using whole-exome sequence analysis, we discovered another possible mutational cause of BPES in integrin subunit beta 5 (ITGB5). The ITGB5 mutation (c.608T>C, p.Ile203Thr) appears in the base sequence of all BPES⁺ patients in this family, and it appears to be a three-generation-inherited mutation. It can cause changes in base sequence and protein function, and there may be cosegregation of disease phenotypes. ITGB5 is located on the long arm of chromosome three (3q21.2) and is close to the known pathogenic gene FOXL2 (3q23). This study is the first to report ITGB5 mutations in BPES, and we speculate that it may be directly involved in the pathogenesis of BPES or indirectly through the regulation of FOXL2.

rs6426881 in the 3'-UTR of PBX1 is involved in breast and gastric cancers via altering the binding potential of miR-522-3p

BACKGROUND

Breast and gastric cancers are the most important diseases that lead to cancer death and social healthcare challenge. Overexpression of PBX1, a proto-oncogene, is correlated with the progression and metastasis of various cancers. For the first time, in this study the researchers evaluated the relationship between rs6426881, affecting miR-522-3p binding to the PBX1, with breast and gastric cancers.

METHODS AND RESULTS

The Microarray analysis was performed for finding the relative expression level of PBX1 and hsa-miR-522-3p, based on high throughput experiments. The GSE54397, GSE112369, GSE10810, GSE241585.ER, GSE24185.PR, GSE68373, and GSE38167 datasets were analyzed. A case-control study was carried out in 123 Iranian suffering from breast cancer and 132 participants as control samples as well as 130 people suffering from gastric cancer and 54 people as control group members. SNP rs6426881 in the 3'-UTR of PBX1 was genotyped by the High-Resolution Melting (HRM) method. Association analysis revealed that rs6426881 is correlated with Estrogen and Progesterone receptors, grade, and stage of breast cancer. Furthermore, a significant relationship was observed between the genotypes and blood groups in gastric cancer, while the distribution of alleles was significantly related to smoking, status of the primary tumor, and metastasis (Chi-Square P < 0.05). Finally, Bioinformatics analyses suggested that rs6426881 contains binding sites for miR-522-3p in the 3'-UTR of PBX1 transcript. The finding suggested that TT genotype is associated with poor prognosis in breast and gastric cancer.

CONCLUSIONS

The rs6426881 T allele at PBX1 3'-UT is significantly related to breast and gastric cancers by altering the regulatory affinity of miR-522-3p to PBX1 3'-UTR and may be suggested as a novel prognostic biomarker for the diseases.

The intronic branch point sequence is under strong evolutionary constraint in the bovine and human genome

The branch point sequence is a cis-acting intronic motif required for mRNA splicing. Despite their functional importance, branch point sequences are not routinely annotated. Here we predict branch point sequences in 179,476 bovine introns and investigate their variability using a catalogue of 29.4 million variants detected in 266 cattle genomes. We localize the bovine branch point within a degenerate heptamer "nnyTrAy". An adenine residue at position 6, that acts as branch point, and a thymine residue at position 4 of the heptamer are more strongly depleted for mutations than coding sequences suggesting extreme purifying selection. We provide evidence that mutations affecting these evolutionarily constrained residues lead to alternative splicing. We confirm evolutionary constraints on branch point sequences using a catalogue of 115 million SNPs established from 3,942 human genomes of the gnomAD database.

Gene Position Index Mutation Detection Algorithm Based on Feedback Fast Learning Neural Network

In the detection of genome variation, the research on the internal correlation of reference genome is deepening; the detection of variation in genome sequence has become the focus of research, and it has also become an effective path to find new genes and new functional proteins. The targeted sequencing sequence is used to sequence the exon region of a specific gene in cancer gene detection, and the sequencing depth is relatively large. Traditional alignment algorithms will lose some sequences, which will lead to inaccurate mutation detection. This paper proposes a mutation detection algorithm based on feedback fast learning neural network position index. By establishing a position index relationship for ACGT in the DNA sequence, the subsequence is decomposed into the position relationship of different subsequences corresponding to the main sequence. The positional relationship of the subsequence in the main sequence is determined by the positional relationship. Analyzing SNP and InDel mutations, even structural mutations, through the position correlation of sequences has the advantages of high precision and easy implementation by personal computers. The feedback fast learning neural network is used to verify whether there is a linear relationship between two or more positions. Experimental results show that the mutation points detected by position index are more than those detected by Bcftools, Freebye, Vanscan2, and Gatk.

Transcription initiation of distant core promoters in a large-sized genome of an insect

Background

Core promoters have a substantial influence on various steps of transcription, including initiation, elongation, termination, polyadenylation, and finally, translation. The characterization of core promoters is crucial for exploring the regulatory code of transcription initiation. However, the current understanding of insect core promoters is focused on those of Diptera (especially Drosophila) species with small genome sizes.

Results

Here, we present an analysis of the transcription start sites (TSSs) in the migratory locust, Locusta migratoria, which has a genome size of 6.5 Gb. The genomic differences, including lower precision of transcription initiation and fewer constraints on the distance from transcription factor binding sites or regulatory elements to TSSs, were revealed in locusts compared with Drosophila insects. Furthermore, we found a distinct bimodal log distribution of the distances from the start codons to the core promoters of locust genes. We found stricter constraints on the exon length of mRNA leaders and widespread expression activity of the distant core promoters in locusts compared with fruit flies. We further compared core promoters in seven arthropod species across a broad range of genome sizes to reinforce our results on the emergence of distant core promoters in large-sized genomes.

Conclusions

In summary, our results provide novel insights into the effects of genome size expansion on distant transcription initiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01004-5.

Epigenomic differences in the human and chimpanzee genomes are associated with structural variation

Structural variation (SV), including insertions and deletions (indels), is a primary mechanism of genome evolution. However, the mechanism by which SV contributes to epigenome evolution is poorly understood. In this study, we characterized the association between lineage-specific indels and epigenome differences between human and chimpanzee to investigate how SVs might have shaped the epigenetic landscape. By intersecting medium-to-large human-chimpanzee indels (20 bp-50 kb) with putative promoters and enhancers in cranial neural crest cells (CNCCs) and repressed regions in induced pluripotent cells (iPSCs), we found that 12% of indels overlap putative regulatory and repressed regions (RRRs), and 15% of these indels are associated with lineage-biased RRRs. Indel-associated putative enhancer and repressive regions are approximately 1.3 times and approximately three times as likely to be lineage-biased, respectively, as those not associated with indels. We found a twofold enrichment of medium-sized indels (20-50 bp) in CpG island (CGI)-containing promoters than expected by chance. Lastly, from human-specific transposable element insertions, we identified putative regulatory elements, including NR2F1-bound putative CNCC enhancers derived from SVAs and putative iPSC promoters derived from LTR5s. Our results show that different types of indels are associated with specific epigenomic diversity between human and chimpanzee.

Improved genomic resources for the black tiger prawn (Penaeus monodon)

World production of farmed crustaceans was 7.8 million tons in 2016. While only making up approximately 10% of world aquaculture production, crustaceans are generally high-value species and can earn significant export income for producing countries. Viet Nam is a major seafood producing country earning USD 7.3 billion in 2016 in export income with shrimp as a major commodity. However, there is a general lack of genomic resources available for shrimp species, which is challenging to obtain due to the need to deal with large repetitive genomes, which characterize many decapod crustaceans. The first tiger prawn (P. monodon) genome assembly was assembled in 2016 using the standard Illumina PCR-based pair-end reads and a computationally-efficient but relatively suboptimal assembler, SOAPdenovo v2. As a result, the current P. monodon draft genome is highly fragmented (> 2 million scaffolds with N₅₀ length of <1000 bp), exhibiting only moderate genome completeness (< 35% BUSCO complete single-copy genes). We sought to improve upon the recently published P. monodon genome assembly and completeness by generating Illumina PCR-free pair-end sequencing reads to eliminate genomic gaps associated with PCR-bias and performing de novo assembly using the updated MaSurCA de novo assembler. Furthermore, we scaffolded the assembly with low coverage Nanopore long reads and several recently published deep Illumina transcriptome paired-end sequencing data, producing a final genome assembly of 1.6 Gbp (1,211,364 scaffolds; N₅₀ length of 1982 bp) with an Arthropod BUSCO completeness of 96.8%. Compared to the previously published P. monodon genome assembly from China (NCBI Accession Code: NIUS01), this represents an almost 20% increase in the overall BUSCO genome completeness that now consists of more than 90% of Arthropod BUSCO single-copy genes. The revised P. monodon genome assembly (NCBI Accession Code: VIGR01) will be a valuable resource to support ongoing functional genomics and molecular-based breeding studies in Vietnam.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

Genotype Variations of rs13381800 in TCF4 Gene and rs17039988 in NRXN1 Gene among a Sample of Iranian Patients with Schizophrenia

Objective: Schizophrenia is a complicated mental disorder that affects about 1% of the world's population. It is a complex disease and is approximately 80% inherited. One of the candidate genes in schizophrenia is transcription factor 4 (TCF4), which is positioned on chromosome 18 and is a transcription factor that plays a role in the transcription of Neurexin 1(NRXN1) gene, which is one of the candidate genes for developing schizophrenia. This case-control study aimed to investigate the correlation of TCF4 rs13381800 and NRXN1 rs17039988 polymorphisms with the risk of schizophrenia in a sample of Iranian patients with schizophrenia. Method : A total of 200 individuals were included in this study: 100 patients with schizophrenia (65 males and 35 females), with the mean age of 40.80 ± 11.298 years, and 100 as a control group (63 males and 37 females), with the mean age 32.92 ± 7.391 years. Allele specific polymerase chain reaction and restriction fragment length polymorphism (PCR-RFLP) were done, respectively, for genotyping of rs13381800 (T/C) and rs17039988 (A/C) polymorphisms. Results: The results showed that the frequency of C / C genotype in rs13381800 in patients' group was 9%, while it was 13% in the control group. Also, the frequency of C / C genotype in rs17039988 was 9% in patients and 7% in control groups. Statistical analysis of polymorphisms showed no correlation between patients and controls in rs13381800 (OR = 1.51; CI = 95%; P = 0.366) and rs17039988 (OR = 0.76; CI = 95%; P = 0.602). Conclusion: No significant difference was found between rs13381800 and rs17039988 genotypes between patients and control groups in terms of gender, age and education in the patients group. Our study suggests that there was no correlation between desired polymorphisms with schizophrenia in the studied population.