Indel Variants Research Articles

Rescue of the endogenous FVIII expression in hemophilia A mice using CRISPR Cas9 mRNA LNPs

Gene editing provides a promising alternative approach that may achieve long-term FVIII expression for hemophilia A (HemA) treatment. In this study, we investigated in vivo correction of a mutant factor VIII (FVIII) gene in HemA mice. We first developed MC3-based LNPs for efficient mRNA delivery into liver sinusoidal endothelial cells (LSECs), the major site of FVIII biosynthesis. To target a five base pair deletion in FVIII exon 1 in a specific HemA mouse strain, we injected LNPs encapsulating Cas9 mRNA and specifically designed sgRNAs intravenously for in vivo gene editing of the mutant FVIII. Indel variants generated at the mutant site contained mostly a single base pair deletion , resulting in frame-shift correction of FVIII gene. Sustained endogenous FVIII activity up to 6% was achieved over 26 weeks in treated HemA mice. Sequencing data indicated an average gene editing rate of 15.3% in LSECs. Our study suggests that optimized MC3 LNP formulations, combined with CRISPR/Cas9 technology, can effectively correct the mutant FVIII gene in LSECs and restore FVIII activity for therapeutic treatment of HemA.

The landscape of sequence variations between resistant and susceptible hot peppers to predict functional candidate genes against bacterial wilt disease

BackgroundBacterial wilt (BW), caused by Ralstonia solanacearum (Ral), results in substantial yield losses in pepper crops. Developing resistant pepper varieties through breeding is the most effective strategy for managing BW. To achieve this, a thorough understanding of the genetic information connected with resistance traits is essential. Despite identifying three major QTLs for bacterial wilt resistance in pepper, Bw1 on chromosome 8, qRRs-10.1 on chromosome 10, and pBWR-1 on chromosome 1, the genetic information of related BW pepper varieties has not been sufficiently studied. Here, we resequenced two pepper inbred lines, C. annuum ‘MC4’ (resistant) and C. annuum ‘Subicho’ (susceptible), and analyzed genomic variations through SNPs and Indels to identify candidate genes for BW resistance.ResultsAn average of 139.5 Gb was generated among the two cultivars, with coverage ranging from 44.94X to 46.13X. A total of 8,815,889 SNPs was obtained between ‘MC4’ and ‘Subicho’. Among them, 31,190 (0.35%) were non-synonymous SNPs (nsSNPs) corresponding to 10,926 genes, and these genes were assigned to 142 Gene Ontology (GO) terms across the two cultivars. We focused on three known BW QTL regions by identifying genes with sequence variants through gene set enrichment analysis and securing those belonging to high significant GO terms. Additionally, we found 310 NLR genes with nsSNP variants between ‘MC4’ (R) and ‘Subicho’ (S) within these regions. Also, we performed an Indel analysis on these genes. By integrating all this data, we identified eight candidate BW resistance genes, including two NLR genes with nsSNPs variations in qRRs-10.1 on chromosome 10.ConclusionWe identified genomic variations in the form of SNPs and Indels by re-sequencing two pepper cultivars with contrasting traits for bacterial wilt. Specifically, the four genes associated with pBWR-1 and Bw1 that exhibit both nsSNP and Indel variations simultaneously in ‘Subicho’, along with the two NLR genes linked to qRRs-10.1, which are known for their direct involvement in immune responses, are identified as most likely BW resistance genes. These variants in leading candidate genes associated with BW resistance can be used as important markers for breeding pepper varieties.

Single Nucleotide Polymorphisms and Insertion/Deletion Variation Analysis of Octoploid and Decaploid Tropical Oil Tea Camellia Populations Based on Whole-Genome Resequencing

Oil tea camellia (Camellia spp.) is an important woody oil crop with a high nutritional and economic value. Whole-genome resequencing (WGR) technology can provide an in-depth understanding of the genetic background of this plant as well as a reference for breeding research, germplasm resource conservation, and genetic modification. In this study, we analyzed SNP and InDel variations in 49 individual oil tea camellia germplasm samples collected from five populations located in three provinces of China: Hainan, Guangdong and Guangxi. The samples were analyzed through WGR after the ploidy of the samples was determined through flow cytometry. A total of 239,441,603 high-quality single nucleotide polymorphisms (SNPs) and 23,510,374 high-quality insertion/deletion variation sites (InDels) were obtained. The distribution of SNPs and InDels in different functional regions differed significantly, with a high density of variations in non-coding regions, such as intergenic regions and introns, and a relatively low density of variations in coding regions. Transition was the main type of SNP variation. A population genetic diversity analysis revealed that the sampled oil tea camellia populations exhibited a high genetic diversity and extensive genetic variation. The genetic diversity of the oil tea camellia populations in the Hainan region was higher than inland regions. This study also determined the genetic diversity of and variations between octoploid and decaploid oil tea camellia in the tropics and between Hainan-based and inland oil tea camellia. Such findings provide a reference for the conservation of germplasm resources and the genetic modification of oil tea camellia.

Supervised machine learning of outbred mouse genotypes to predict hepatic immunological tolerance of individuals

It is essential to elucidate the molecular mechanisms underlying liver transplant tolerance and rejection. In cases of mouse liver transplantation between inbred strains, immunological rejection of the allograft is reduced with spontaneous apoptosis without immunosuppressive drugs, which differs from the actual clinical result. This may be because inbred strains are genetically homogeneous and less heterogeneous than others. We exploited outbred CD1 mice, which show highly heterogeneous genotypes among individuals, to search for biomarkers related to immune responses and to construct a model for predicting the outcome of liver allografting. Of the 36 mice examined, 18 died within 3 weeks after transplantation, while the others survived for more than 6 weeks. Whole-exome sequencing of the 36 donors revealed more than 9 million variants relative to the C57BL/6 J reference. We selected 6517 single-nucleotide and indel variants and performed machine learning to determine whether or not we could predict the prognosis of each genotype. Models were built by both deep learning with a one-dimensional convolutional neural network and linear classification and evaluated by leave-one-out cross-validation. Given that one short-lived mouse died early in an accident, the models perfectly predicted the outcome of all individuals, suggesting the importance of genotype collection. In addition, linear classification models provided a list of loci potentially responsible for these responses. The present methods as well as results is likely to be applicable to liver transplantation in humans.

High-coverage nanopore sequencing of samples from the 1000 Genomes Project to build a comprehensive catalog of human genetic variation.

Fewer than half of individuals with a suspected Mendelian or monogenic condition receive a precise molecular diagnosis after comprehensive clinical genetic testing. Improvements in data quality and costs have heightened interest in using long-read sequencing (LRS) to streamline clinical genomic testing, but the absence of control datasets for variant filtering and prioritization has made tertiary analysis of LRS data challenging. To address this, the 1000 Genomes Project ONT Sequencing Consortium aims to generate LRS data from at least 800 of the 1000 Genomes Project samples. Our goal is to use LRS to identify a broader spectrum of variation so we may improve our understanding of normal patterns of human variation. Here, we present data from analysis of the first 100 samples, representing all 5 superpopulations and 19 subpopulations. These samples, sequenced to an average depth of coverage of 37x and sequence read N50 of 54 kbp, have high concordance with previous studies for identifying single nucleotide and indel variants outside of homopolymer regions. Using multiple structural variant (SV) callers, we identify an average of 24,543 high-confidence SVs per genome, including shared and private SVs likely to disrupt gene function as well as pathogenic expansions within disease-associated repeats that were not detected using short reads. Evaluation of methylation signatures revealed expected patterns at known imprinted loci, samples with skewed X-inactivation patterns, and novel differentially methylated regions. All raw sequencing data, processed data, and summary statistics are publicly available, providing a valuable resource for the clinical genetics community to discover pathogenic SVs.

Identification of Pathogenic Copy Number Variants in Mexican Patients With Inherited Retinal Dystrophies Applying an Exome Sequencing Data-Based Read-Depth Approach.

Retinal dystrophies (RDs) are the most common cause of inherited blindness worldwide and are caused by genetic defects in about 300 different genes. While targeted next-generation sequencing (NGS) has been demonstrated to be a reliable and efficient method to identify RD disease-causing variants, it doesn't routinely identify pathogenic structural variant as copy number variations (CNVs). Targeted NGS-based CNV detection has become a crucial step for RDs molecular diagnosis, particularly in cases without identified causative single nucleotide or Indels variants. Herein, we report the exome sequencing (ES) data-based read-depth bioinformatic analysis in a group of 30 unrelated Mexican RD patients with a negative or inconclusive genetic result after ES. CNV detection was performed using ExomeDepth software, an R package designed to detect CNVs using exome data. Bioinformatic validation of identified CNVs was conducted through a commercially available CNV caller. All identified candidate pathogenic CNVs were orthogonally verified through quantitative PCR assays. Pathogenic or likely pathogenic CNVs were identified in 6 out of 30 cases (20%), and of them, a definitive molecular diagnosis was reached in 5 cases, for a final diagnostic rate of ~17%. CNV-carrying genes included CLN3 (2 cases), ABCA4 (novel deletion), EYS, and RPGRIP1. Our results indicate that bioinformatic analysis of ES data is a reliable method for pathogenic CNV detection and that it should be incorporated in cases with a negative or inconclusive molecular result after ES.

GL5.2, a Quantitative Trait Locus for Rice Grain Shape, Encodes a RING-Type E3 Ubiquitin Ligase.

Grain weight and grain shape are important traits that determine rice grain yield and quality. Mining more quantitative trait loci (QTLs) that control grain weight and shape will help to further improve the molecular regulatory network of rice grain development and provide gene resources for high-yield and high-quality rice varieties. In the present study, a QTL for grain length (GL) and grain width (GW), qGL5.2, was firstly fine-mapped into a 21.4 kb region using two sets of near-isogenic lines (NILs) derived from the indica rice cross Teqing (TQ) and IRBB52. In the NIL populations, the GL and ratio of grain length to grain width (RLW) of the IRBB52 homozygous lines increased by 0.16-0.20% and 0.27-0.39% compared with the TQ homozygous lines, but GW decreased by 0.19-0.75%. Then, by analyzing the grain weight and grain shape of the knock-out mutant, it was determined that the annotation gene Os05g0551000 encoded a RING-type E3 ubiquitin ligase, which was the cause gene of qGL5.2. The results show that GL and RLW increased by 2.44-5.48% and 4.19-10.70%, but GW decreased by 1.69-4.70% compared with the recipient. Based on the parental sequence analysis and haplotype analysis, one InDel variation located at -1489 in the promoter region was likely to be the functional site of qGL5.2. In addition, we also found that the Hap 5 (IRBB52-type) increased significantly in grain length and grain weight compared with other haplotypes, indicating that the Hap 5 can potentially be used in rice breeding to improve grain yield and quality.

QTL-Seq: Rapid, Cost-Effective, and Reliable Method for QTL Identification

QTL-seq is a powerful method that integrates whole-genome sequencing (WGS) with bulk-segregant analysis to rapidly and reliably identify quantitative trait loci (QTLs) associated with specific traits. This approach significantly advances traditional QTL mapping by eliminating the need for genome wide DNA markers such as SSR, RFLP, and INDELs, which are typically used in linkage-based QTL mapping. Instead, QTL-seq leverages WGS to detect all genetic variations such as SNPs, Indels, and Structural Variants across the entire genome, providing a comprehensive resource for marker development in marker-assisted selection. The QTL-seq process begins with the creation of genetically diverse mapping populations, such as F2 or RILs, followed by detailed phenotypic characterization. DNA from plants exhibiting similar phenotypes is pooled into bulk groups and sequenced, allowing for cost-effective and efficient QTL identification. Identified QTLs can be further validated through fine mapping using recombinant screenings and progeny testing, leading to the identification of candidate genes associated with traits of interest. In this study, we outline a user-friendly QTL-seq pipeline, from sequencing to data visualization, using the methodology and data provided by Takagi et al., 2013, to demonstrate its practical application. While the manuscript primarily focuses on describing the pipeline, we also conducted a case study analysis with real data to showcase its effectiveness. Our work contributes to the broader understanding of QTL-seq applications and offers practical recommendations for optimizing this method in future breeding programs.

Effect of VDR and TLR2 gene variants on the clinical course of patients with COVID-19 disease.

The coronavirus disease 2019 (COVID-19) pandemic, which has caused a major global health crisis, primarily targets the upper and lower respiratory tract. But infected individuals may experience different clinical symptoms, ranging from asymptomatic to critical. The vitamin D receptor (VDR) and Toll-like receptor 2 (TLR2) polymorphisms play a role in the immune response. This study aimed to evaluate the effect of VDR Bsml (rs1544410) and TLR2 23bp indel variants on the clinical status of Turkish patients with COVID-19 disease. A total of 312 people, including 106 intensive care unit (ICU) patients, 103 symptomatic hospitalized patients, and 103 healthy controls, were included in the study. The VDR BsmI and TLR2 23bp indel were genotyped using polymerase chain reaction and/or restriction fragment length fraction methods. The VDR BsmI b/b genotype and b allele were higher in symptomatic patients compared to the healthy control group (p = 0.035). The VDR BsmI B/B and B/b genotype distribution did not differ between ICU patients and both symptomatic patients and controls (p > 0.05). We found that B/B:B/b+b/b and B/B+B/b:b/b were significantly different in symptomatic patients compared to controls (p = 0.033 and p = 0.041, respectively). The VDR BsmI b/b genotype distribution was found to be lower in deceased patients than in living patients (p = 0.023). There was no significant difference between the groups in terms of TLR2 23bp indel genotype and allele distribution (p > 0.05). Our study results suggest that the VDR BsmI b allele may have a role in COVID-19 patients with symptomatic findings. These data need to be repeated in different ethnic and larger sample groups.

Whole-Genome Resequencing Reveals Signatures of Adaptive Evolution in Acanthopagrus latus and Rhabdosargus sarba.

Acanthopagrus latus and Rhabdosargus sarba are economically important marine species along the coast of China, with similar external morphological characteristics and living habits, with wide distribution and strong adaptability. To investigate the molecular mechanisms underlying the adaptive evolution of these two species, we conducted whole-genome resequencing of 10 individuals of both species from the coastal waters of Wuyu Island, Fujian, China, using high-throughput sequencing technology. We obtained SNP, InDel, CNV, and SV variation information and annotated these variations, constructing a genomic variation database for both species. By comparing the resequencing data with reference genomes, we identified 9,829,511 SNP loci in the population of A. latus and 34,051,056 SNP loci in the population of R. sarba. Using whole-genome SNP data, we employed Fst and ROD methods to identify candidate genomic regions under selection. Functional annotation and enrichment analysis using GO and KEGG databases revealed potential adaptive evolution in R. sarba associated with immune response, feeding, growth and development, and locomotion, while A. latus showed potential adaptive evolution associated with immune response, nervous system, growth and development, and metabolism.

Genomic insights into prenatal diagnosis of congenital heart defects: value of CNV-seq and WES in clinical practice.

This study aimed to assess the efficiency of CNV-seq and WES in detecting genetic cause of congenital heart disease (CHDs) in prenatal diagnoses and to compare CNV detection rate between isolated and non-isolated CHD cases. We conducted a retrospective study of 118 Chinese fetuses diagnosed with CHD by prenatal ultrasound. Participants underwent CNV-seq and, if necessary, WES to detect chromosomal and single nucleotide variations. The overall detection rate for pathogenic or likely pathogenic chromosomal abnormalities was 16.9%, including 7.6% aneuploidies and 9.3% pathogenic/likely pathogenic copy number variations (CNVs), predominantly 22q11.2 deletion syndrome (54.4%). The sensitivity and specificity of CNV-Seq for detecting P/Lp CNVs were 95% and 100%, respectively. CNV-Seq offered a 6.7% improvement in detecting chromosomal abnormalities over karyotyping. WES further identified significant single nucleotide and small indel variations contributing to CHD in genes such as TMEM67, PLD1, ANKRD11, and PNKP, enhancing diagnostic yield by 14.8% in cases negative for CNVs. Non-isolated CHD cases exhibited higher rates of detectable chromosomal abnormalities compared to isolated cases (32.4% vs. 9.9%, p = 0.005), underlining the genetic complexity of these conditions. The combined use of CNV-seq and WES provides a comprehensive approach to prenatal genetic testing for CHDs, unveiling significant genetic cause that could impact clinical management and parental decision-making. This study supports the integration of these advanced genomic technologies in routine prenatal diagnostics to increase detection diagnostic yields of causal genetic variants associated with CHDs.

Diversity of ribosomes at the level of rRNA variation associated with human health and disease

With hundreds of copies of rDNA, it is unknown whether they possess sequence variations that form different types of ribosomes. Here, we developed an algorithm for long-read variant calling, termed RGA, which revealed that variations in human rDNA loci are predominantly insertion-deletion (indel) variants. We developed full-length rRNA sequencing (RIBO-RT) and in situ sequencing (SWITCH-seq), which showed that translating ribosomes possess variation in rRNA. Over 1,000 variants are lowly expressed. However, tens of variants are abundant and form distinct rRNA subtypes with different structures near indels as revealed by long-read rRNA structure probing coupled to dimethyl sulfate sequencing. rRNA subtypes show differential expression in endoderm/ectoderm-derived tissues, and in cancer, low-abundance rRNA variants can become highly expressed. Together, this study identifies the diversity of ribosomes at the level of rRNA variants, their chromosomal location, and unique structure as well as the association of ribosome variation with tissue-specific biology and cancer.

Identification and validation of novel breed-specific biomarker for the purpose of village chicken authentication using genomics approaches

Local village chicken, or "Ayam kampung" as it's known in Malaysia, is considered a premium chicken breed with a higher price than other chicken breeds. As a result of their comparable appearances and sizes, colored broiler chickens are often sold as village chickens, which is a form of food fraud that can result in a 3- to 4-fold rise in profit. Therefore, developing a breed-specific authentication method is crucial for preventing food fraud in the poultry industry. This study aims to investigate the genetic diversity of village chickens from other commercial chicken breed populations available in the market (broiler [Cobb], colored broiler [Hubbard], and layer [DeKalb]) to identify breed-specific DNA fragments as biomarkers for village chicken authentication. The Whole-genome sequencing and mutation calling of 12 chickens (3 chickens/breed) led to the identification of a total of 73,454,654 single nucleotide polymorphisms (SNP) and 8,762,338 insertion and deletions (InDel) variants, with more variants detected in the village chicken population (6,346,704 SNPs; 752,408 InDels) compared to commercial breeds. Therefore, this study revealed that village chickens were more genetically variable compared to other breeds in Malaysia. Furthermore, the breed-specific genomic region located on chromosome 1 (1:84,405,652) harboring SNP (C-T) with high discrimination power was discovered and validated which can be considered as a novel breed-specific biomarker to develop a method for accurate authentication of village chickens in Malaysia. This authentication method offers potentialw applications in the chicken industry and food safety.

Identification of truncated variants in GLI family zinc finger 3 (GLI3) associated with polydactyly

BackgroundPolydactyly is a prevalent congenital anomaly with an incidence of 2.14 per 1000 live births in China. GLI family zinc finger 3 (GLI3) is a classical causative gene of polydactyly, and serves as a pivotal transcription factor in the hedgehog signaling pathway, regulating the development of the anterior-posterior axis in limbs.MethodsThree pedigrees of polydactyly patients were enrolled from Hunan Province, China. Pathogenic variants were identified by whole-exome sequencing (WES) and Sanger sequencing.ResultsThree variants in GLI3 were identified in three unrelated families, including a novel deletion variant (c.1372del, p.Thr458GlnfsTer44), a novel insertion-deletion (indel) variant (c.1967_1968delinsAA, p.Ser656Ter), and a nonsense variant (c.2374 C > T, p.Arg792Ter). These variants were present exclusively in patients but not in healthy individuals.ConclusionsWe identified three pathogenic GLI3 variants in polydactyly patients, broadening the genetic spectrum of GLI3 and contributing significantly to genetic counseling and diagnosis for polydactyly.

Association analysis of indel variants and gene expression identifies MDM4 as a novel locus for skeletal muscle hypertrophy and power athlete status.

Insertions and deletions (indels) are the second most common type of variation in the human genome. However, limited data on their associations with exercise-related phenotypes have been documented. The aim of the present study was to examine the association between 18,370 indel variants and power athlete status, followed by additional studies in 357,246 individuals. In the discovery phase, the D allele of the MDM4 gene rs35493922 I/D polymorphism was over-represented in power athletes compared with control subjects (P=7.8×10-9) and endurance athletes (P=0.0012). These findings were replicated in independent cohorts, showing a higher D allele frequency in power athletes compared with control subjects (P=0.016) and endurance athletes (P=0.031). Furthermore, the D allele was positively associated (P=0.0013) with greater fat-free mass in the UK Biobank. MDM4 encodes a protein that inhibits the activity of p53, which induces muscle fibre atrophy. Accordingly, we found that MDM4 expression was significantly higher in the vastus lateralis of power athletes compared with endurance athletes (P=0.0009) and was positively correlated with the percentage of fast-twitch muscle fibres (P=0.0062) and the relative area occupied by fast-twitch muscle fibres (P=0.0086). The association between MDM4 gene expression and an increased proportion of fast-twitch muscle fibres was confirmed in two additional cohorts. Finally, we found that the MDM4 DD genotype was associated with increased MDM4 gene expression in vastus lateralis and greater cross-sectional area of fast-twitch muscle fibres. In conclusion, MDM4 is suggested to be a potential regulator of muscle fibre specification and size, with its indel variant being associated with power athlete status. HIGHLIGHTS: What is the central question of this study? Which indel variants are functional and associated with sport- and exercise-related traits? What is the main finding and its importance? Out of 18,370 tested indels, the MDM4 gene rs35493922 I/D polymorphism was found to be the functional variant (affecting gene expression) and the most significant, with the deletion allele showing associations with power athlete status, fat-free mass and cross-sectional area of fast-twitch muscle fibres. Furthermore, the expression of MDM4 was positively correlated with the percentage of fast-twitch muscle fibres and the relative area occupied by fast-twitch muscle fibres.

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.

A novel three-layer module BoMYB1R1-BoMYB4b/BoMIEL1-BoDFR1 regulates anthocyanin accumulation in kale.

Anthocyanin is an important pigment responsible for plant coloration and beneficial to human health. Kale (Brassica oleracea var. acephala), a primary cool-season flowers and vegetables, is an ideal material to study anthocyanin biosynthesis and regulation mechanisms due to its anthocyanin-rich leaves. However, the underlying molecular mechanism of anthocyanin accumulation in kale remains poorly understood. Previously, we demonstrated that BoDFR1 is a key gene controlling anthocyanin biosynthesis in kale. Here, we discovered a 369-bp InDel variation in the BoDFR1 promoter between the two kale inbred lines with different pink coloration, which resulted in reduced transcriptional activity of the BoDFR1 gene in the light-pink line. With the 369-bp insertion as a bait, an R2R3-MYB repressor BoMYB4b was identified using the yeast one-hybrid screening. Knockdown of the BoMYB4b gene led to increased BoDFR1 expression and anthocyanin accumulation. An E3 ubiquitin ligase, BoMIEL1, was found to mediate the degradation of BoMYB4b, thereby promoting anthocyanin biosynthesis. Furthermore, the expression level of BoMYB4b was significantly reduced by light signals, which was attributed to the direct repression of the light-signaling factor BoMYB1R1 on the BoMYB4b promoter. Our study revealed that a novel regulatory module comprising BoMYB1R1, BoMIEL1, BoMYB4b, and BoDFR1 finely regulates anthocyanin accumulation in kale. The findings aim to establish a scientific foundation for genetic improvement of leaf color traits in kale, meanwhile, providing a reference for plant coloration studies.

Indel variation in the mitochondrial ND5 region supports monophyly of the tribe Hippoglossoidini (sensu Vinnikov et al. 2018) within the family Pleuronectidae

AbstractIndel (insertion–deletion) events observed in the genome represent irreversible mutational processes, making indel regions crucial characteristics for discussing phylogenetic relationships. The tribe Hippoglossoidini is a recently proposed taxonomic group based on the molecular phylogenetic analyses of both mitochondrial (mt) and nuclear DNA sequences. However, no synapomorphic characteristics have been identified within this tribe, either morphologically or molecularly. In the present study, we sequenced the ND5 region of mtDNA in the righteye flounder species and conducted interspecific comparisons. We found a 12 bp indel immediately upstream of the stop codon in the ND5 region. A comparative analysis of this region with outgroup species from the Paralichthyidae revealed that the indel was a unique insertion shared by the common ancestor of the Hippoglossoidini species, providing irreversible evidence to support the monophyly of this taxonomic group (synapomorphic characters).

Causative role of a novel intronic indel variant in FBN1 and maternal germinal mosaicism in Marfan syndrome

BackgroundMarfan syndrome (MFS) is an autosomal dominant connective tissue disease with wide clinical heterogeneity, and mainly caused by pathogenic variants in fibrillin-1 (FBN1).MethodsA Chinese 4-generation MFS pedigree with 16 family members was recruited and exome sequencing (ES) was performed in the proband. Transcript analysis (patient RNA and minigene assays) and in silico structural analysis were used to determine the pathogenicity of the variant. In addition, germline mosaicism in family member (Ι:1) was assessed using quantitative fluorescent polymerase chain reaction (QF-PCR) and short tandem repeat PCR (STR) analyses.ResultsTwo cis-compound benign intronic variants of FBN1 (c.3464–4 A > G and c.3464-5G > A) were identified in the proband by ES. As a compound variant, c.3464-5_3464-4delGAinsAG was found to be pathogenic and co-segregated with MFS. RNA studies indicated that aberrant transcripts were found only in patients and mutant-type clones. The variant c.3464-5_3464-4delGAinsAG caused erroneous integration of a 3 bp sequence into intron 28 and resulted in the insertion of one amino acid in the protein sequence (p.Ile1154_Asp1155insAla). Structural analyses suggested that p.Ile1154_Asp1155insAla affected the protein’s secondary structure by interfering with one disulfide bond between Cys1140 and Cys1153 and causing the extension of an anti-parallel β sheet in the calcium-binding epidermal growth factor-like (cbEGF)13 domain. In addition, the asymptomatic family member Ι:1 was deduced to be a gonadal mosaic as assessed by inconsistent results of sequencing and STR analysis.ConclusionsTo our knowledge, FBN1 c.3464-5_3464-4delGAinsAG is the first identified pathogenic intronic indel variant affecting non-canonical splice sites in this gene. Our study reinforces the importance of assessing the pathogenic role of intronic variants at the mRNA level, with structural analysis, and the occurrence of mosaicism.

The quantitative genetics of gene expression in Mimulus guttatus.

Gene expression can be influenced by genetic variants that are closely linked to the expressed gene (cis eQTLs) and variants in other parts of the genome (trans eQTLs). We created a multiparental mapping population by sampling genotypes from a single natural population of Mimulus guttatus and scored gene expression in the leaves of 1,588 plants. We find that nearly every measured gene exhibits cis regulatory variation (91% have FDR < 0.05). cis eQTLs are usually allelic series with three or more functionally distinct alleles. The cis locus explains about two thirds of the standing genetic variance (on average) but varies among genes and tends to be greatest when there is high indel variation in the upstream regulatory region and high nucleotide diversity in the coding sequence. Despite mapping over 10,000 trans eQTL / affected gene pairs, most of the genetic variance generated by trans acting loci remains unexplained. This implies a large reservoir of trans acting genes with subtle or diffuse effects. Mapped trans eQTLs show lower allelic diversity but much higher genetic dominance than cis eQTLs. Several analyses also indicate that trans eQTLs make a substantial contribution to the genetic correlations in expression among different genes. They may thus be essential determinants of "gene expression modules," which has important implications for the evolution of gene expression and how it is studied by geneticists.