Syntenic Regions Research Articles

Assembly and comparative analysis of the first complete mitochondrial genome of Astragalus membranaceus (Fisch.) Bunge: an invaluable traditional Chinese medicine

BackgroundAstragalus membranaceus (Fisch.) Bunge is one of the most well-known tonic herbs in traditional Chinese medicine, renowned for its remarkable medicinal value in various clinical contexts. The corresponding chloroplast (cp) and nuclear genomes have since been accordingly sequenced, providing valuable information for breeding and phylogeny studies. However, the mitochondrial genome (mitogenome) of A. membranaceus remains unexplored, which hinders comprehensively understanding the evolution of its genome.ResultsFor this study, we de novo assembled the mitogenome of A. membranaceus (Fisch.) Bunge var. mongholicus (Bunge) P. K. Hsiao using a strategy integrating Illumina and Nanopore sequencing technology and subsequently performed comparative analysis with its close relatives. The mitogenome has a multi-chromosome structure, consisting of two circular chromosomes with a total length of 398,048 bp and an overall GC content of 45.3%. It encodes 54 annotated functional genes, comprising 33 protein-coding genes (PCGs), 18 tRNA genes, and 3 rRNA genes. An investigation of codon usage in the PCGs revealed an obvious preference for codons ending in A or U (T) bases, given their high frequency. RNA editing identified 500 sites in the coding regions of mt PCGs that exhibit a perfect conversion of the base C to U, a process that tends to lead to the conversion of hydrophilic amino acids into hydrophobic amino acids. From the mitogenome analysis, a total of 399 SSRs, 4 tandem repeats, and 77 dispersed repeats were found, indicating that A. membranaceus possesses fewer repeats compared to its close relatives with similarly sized mitogenomes. Selection pressure analysis indicated that most mt PCGs were purifying selection genes, while only five PCGs (ccmB, ccmFc, ccmFn, nad3, and nad9) were positive selection genes. Notably, positive selection emerged as a critical factor in the evolution of ccmB and nad9 in all the pairwise species comparisons, suggesting the extremely critical role of these genes in the evolution of A. membranaceus. Moreover, we inferred that 22 homologous fragments have been transferred from cp to mitochondria (mt), in which 5 cp-derived tRNA genes remain intact in the mitogenome. Further comparative analysis revealed that the syntenic region and mt gene organization are relatively conserved within the provided legumes. The comparison of gene content indicated that the gene composition of Fabaceae mitogenomes differed. Finally, the phylogenetic tree established from analysis is largely congruent with the taxonomic relationships of Fabaceae species and highlights the close relationship between Astragalus and Oxytropis.ConclusionsWe provide the first report of the assembled and annotated A. membranaceus mitogenome, which enriches the genetic resources available for the Astragalus genus and lays the foundation for comprehensive exploration of this invaluable medicinal plant.

Convergence and divergence of diploid and tetraploid cotton genomes.

Polyploidy is an important driving force in speciation and evolution; however, the genomic basis for parallel selection of a particular trait between polyploids and ancestral diploids remains unexplored. Here we construct graph-based pan-genomes for diploid (A2) and allotetraploid (AD1) cotton species, enabled by an assembly of 50 genomes of genetically diverse accessions. We delineate a mosaic genome map of tetraploid cultivars that illustrates genomic contributions from semi-wild forms into modern cultivars. Pan-genome comparisons identify syntenic and hyper-divergent regions of continued variation between diploid and tetraploid cottons, and suggest an ongoing process of sequence evolution potentially linked to the contrasting genome size change in two subgenomes. We highlight 43% of genetic regulatory relationships for gene expression in diploid encompassing sequence divergence after polyploidy, and specifically characterize six underexplored convergent genetic loci contributing to parallel selection of fiber quality. This study offers a framework for pan-genomic dissection of genetic regulatory components underlying parallel selection of desirable traits in organisms.

Acquisition and evolution of the neurotoxin domoic acid biosynthesis gene cluster in Pseudo-nitzschia species

Of the hitherto over 60 taxonomically identified species in the genus of Pseudo-nitzschia, 26 have been confirmed to be toxigenic. Nevertheless, the acquisition and evolution of the toxin biosynthesis (dab) genes by this extensive group of Pseudo-nitzschia species remains unclear. Through constructing chromosome-level genomes of three Pseudo-nitzschia species and draft genomes of ten additional Pseudo-nitzschia species, putative genomic integration sites for the dab genes in Pseudo-nitzschia species were explored. A putative breakpoint was observed in syntenic regions in the dab gene cluster-lacking Pseudo-nitzschia species, suggesting potential independent losses of dab genes. The breakpoints between this pair of conserved genes were also identified in some dab genes-possessing Pseudo-nitzschia species, suggesting that the dab gene clusters transposed to other loci after the initial integration. A “single acquisition, multiple independent losses (SAMIL)” model is proposed to explain the acquisition and evolution of the dab gene cluster in Pseudo-nitzschia species.

Regional active transcription associates with homoeologous exchange breakpoints in synthetic Brassica tetraploids.

Polyploidization plays a crucial role in plant evolution and is becoming increasingly important in breeding. Structural variations and epigenomic repatterning have been observed in synthetic polyploidizations. However, the mechanisms underlying the occurrence and their effects on gene expression and phenotype remain unknown. Here, we investigated genome-wide large deletion/duplication regions (DelDups) and genomic methylation dynamics in leaf organs of progeny from the first eight generations of synthetic tetraploids derived from Chinese cabbage (Brassica rapa L. ssp. pekinensis) and cabbage (Brassica oleracea L. var. capitata). One- or two-copy DelDups, with a mean size of 5.70 Mb (400 kb to 65.85 Mb), occurred from the first generation of selfing and thereafter. The duplication of a fragment in one subgenome consistently coincided with the deletion of its syntenic fragment in the other subgenome, and vice versa, indicating that these DelDups were generated by homoeologous exchanges (HEs). Interestingly, the larger the genomic syntenic region, the higher the frequency of DelDups, further suggesting that the pairing of large homoeologous fragments is crucial for HEs. Moreover, we found that the active transcription of continuously distributed genes in local regions is positively associated with the occurrence of HE breakpoints. In addition, the expression of genes within DelDups exhibited a dosage effect, and plants with extra parental genomic fragments generally displayed phenotypes biased toward the corresponding parent. Genome-wide methylation fluctuated remarkably, which did not clearly affect gene expression on a large scale. Our findings provide insights into the early evolution of polyploid genomes, offering valuable knowledge for polyploidization-based breeding.

Molecular mechanism of vivipary as revealed by the genomes of viviparous mangroves and non-viviparous relatives

Vivipary is a prominent feature of mangroves, allowing seeds to complete germination while attached to the mother plant, and equips propagules to endure and flourish in challenging coastal intertidal wetlands. However, vivipary-associated genetic mechanisms remain largely elusive. Genomes of two viviparous mangrove species and a non-viviparous inland relative were sequenced and assembled at the chromosome level. Comparative genomic analyses between viviparous and non-viviparous genomes revealed that DELAY OF GERMINATION 1 (DOG1) family genes (DFGs), the proteins from which are crucial for seed dormancy, germination, and reserve accumulation, are either lost or dysfunctional in the entire lineage of true viviparous mangroves but are present and functional in their inland, non-viviparous relatives. Transcriptome dynamics at key stages of vivipary further highlighted the roles of phytohormonal homeostasis, proteins stored in mature seeds, and proanthocyanidins in vivipary under conditions lacking DFGs. Population genomic analyses elucidate dynamics of syntenic regions surrounding the missing DFGs. Our findings demonstrated the genetic foundation of constitutive vivipary in Rhizophoraceae mangroves.

Dissecting the contribution of human chromosome 21 syntenic regions to recognition memory processes in adult and aged mouse models of Down syndrome.

Trisomy of human chromosome 21 (Hsa21) results in a constellation of features known as Down syndrome (DS), the most common genetic form of intellectual disability. Hsa21 is orthologous to three regions in the mouse genome on mouse chromosome 16 (Mmu16), Mmu17 and Mmu10. We investigated genotype-phenotype relationships by assessing the contribution of these three regions to memory function and age-dependent cognitive decline, using three mouse models of DS, Dp1Tyb, Dp(17)3Yey, Dp(10)2Yey, that carry an extra copy of the Hsa21-orthologues on Mmu16, Mmu17 and Mmu10, respectively. Prior research on cognitive function in DS mouse models has largely focused on models with an extra copy of the Mmu16 region and relatively little is known about the effects of increased copy number on Mmu17 and Mmu10 on cognition and how this interacts with the effects of aging. As aging is is a critical contributor to cognitive and psychiatric changes in DS, we hypothesised that ageing would differentially impact memory function in Dp1Tyb, Dp(17)3Yey, and Dp(10)2Yey, models of DS. Young (12-13 months and old (18-20 months mice Dp1Tyb, Dp(17)3Yey and Dp(10)2Yey mice were tested on a battery of object recognition memory test that assessed object novelty detection, novel location detection and associative object-in place memory. Following behavioral testing, hippocampal and frontal cortical tissue was analysed for expression of glutamatergic receptor proteins using standard immunoblot techniques. Young (12-13 months and old (18-20 months mice Dp1Tyb, Dp(17)3Yey and Dp(10)2Yey mice were tested on a battery of object recognition memory test that assessed object novelty detection, novel location detection and associative object-in place memory. Following behavioral testing, hippocampal and frontal cortical tissue was analysed for expression of glutamatergic receptor proteins using standard immunoblot techniques. Our results show that distinct Hsa21-orthologous regions contribute differentially to cognitive dysfunction in DS mouse models and that aging interacts with triplication of Hsa21-orthologous genes on Mmu10.

The recombination landscape of introgression in yeast.

Meiotic recombination is an evolutionary force that acts by breaking up genomic linkage, increasing the efficacy of selection. Recombination is initiated with a double-strand break which is resolved via a crossover, which involves the reciprocal exchange of genetic material between homologous chromosomes, or a non-crossover, which results in small tracts of non-reciprocal exchange of genetic material. Crossover and non-crossover rates vary between species, populations, individuals, and across the genome. In recent years, recombination rate has been associated with the distribution of ancestry derived from past interspecific hybridization (introgression) in a variety of species. We explore this interaction of recombination and introgression by sequencing spores and detecting crossovers and non-crossovers from two crosses of the yeast Saccharomyces uvarum. One cross is between strains which each contain introgression from their sister species, S. eubayanus, while the other cross has no introgression present. We find that the recombination landscape is significantly different between S. uvarum crosses, and that some of these differences can be explained by the presence of introgression in one cross. Crossovers are reduced and non-crossovers are increased in heterozygous introgression compared to syntenic regions in the cross without introgression. This translates to reduced allele shuffling within introgressed regions, and an overall reduction of shuffling on most chromosomes with introgression compared to the syntenic regions and chromosomes without introgression. Our results suggest that hybridization can significantly influence the recombination landscape, and that the reduction in allele shuffling contributes to the initial purging of introgression in the generations following a hybridization event.

The TIR1/AFB Family in Solanum melongena: Genome-Wide Identification and Expression Profiling under Stresses and Picloram Treatment

TIR1/AFB proteins are a class of auxin receptors with key roles in plant development and biotic and abiotic stress responses; several have been identified as targets of the auxin-mimicking herbicide picloram. In this study, we identified five putative TIR1/AFB gene family members in the important vegetable crop Solanum melongena (eggplant) and characterized them using bioinformatics tools and gene expression analyses. Phylogenetic analysis of the TIR1/AFBs classified them into three subgroups based on their Arabidopsis and Solanum lycopersicum homologs. AFB6 homologs were present only in S. melongena and S. lycopersicum, whereas AFB2/3 homologs were found only in Arabidopsis. One pair of S. melongena TIR1 homologs were located in syntenic regions in the genome and appeared to have arisen by segmental duplication. Promoter analysis revealed 898 cis-elements in the TIR1/AFB promoters, 125 of which were related to hormones, stress, light, or growth responses, but only SmAFB5 had a cis-acting regulatory element involved in auxin responsiveness (AuxRR-core). RNA sequencing and expression profiling showed that the TIR1/AFB genes were differentially expressed at different growth stages and in response to light, temperature, and drought. Only SmTIR1A expression was significantly induced by picloram treatment and different growth stages. TIR1/AFB expression is regulated by microRNAs (miRNAs) in other plant species, and we identified 6 or 29 miRNAs that potentially targeted the five TIR1/AFB genes on the basis of comparisons with S. lycopersicum and S. tuberosum miRNAs, respectively. Three-dimensional protein structure predictions revealed that all the TIR1/AFB proteins were very similar in structure, differing only in the numbers of alpha helices and in one angle linking an α helix and a β sheet. For measuring the function of TIR1/AFB genes in response to drought, SmAFB5 was selected, and knockdown by virus-induced gene silence (VIGS) 35S::SmAFB5 lines showed resistance to drought compared to controls. These analyses provide insight into the potential functions of TIR1/AFBs during growth and in response to stress; they highlight differences among the SmTIR1/AFBs that may be useful for eggplant breeding.

Chromosome-level genome assembly provides insights into the genetic diversity, evolution, and flower development of Prunus conradinae

Prunus conradinae, a valuable flowering cherry belonging to the Rosaceae family subgenus Cerasus and endemic to China, has high economic and ornamental value. However, a high-quality P. conradinae genome is unavailable, which hinders our understanding of its genetic relationships and phylogenesis, and ultimately, the possibility of mining of key genes for important traits. Herein, we have successfully assembled a chromosome-scale P. conradinae genome, identifying 31,134 protein-coding genes, with 98.22% of them functionally annotated. Furthermore, we determined that repetitive sequences constitute 46.23% of the genome. Structural variation detection revealed some syntenic regions, inversions, translocations, and duplications, highlighting the genetic diversity and complexity of Cerasus. Phylogenetic analysis demonstrated that P. conradinae is most closely related to P. campanulata, from which it diverged ~ 19.1 million years ago (Mya). P. avium diverged earlier than P. cerasus and P. conradinae. Similar to the other Prunus species, P. conradinae underwent a common whole-genome duplication event at ~ 138.60 Mya. Furthermore, 79 MADS-box members were identified in P. conradinae, accompanied by the expansion of the SHORT VEGETATIVE PHASE subfamily. Our findings shed light on the complex genetic relationships, and genome evolution of P. conradinae and will facilitate research on the molecular breeding and functions of key genes related to important horticultural and economic characteristics of subgenus Cerasus.

The Chromatin Organization Close to SNP rs12913832, Involved in Eye Color Variation, Is Evolutionary Conserved in Vertebrates.

The most significant genetic influence on eye color pigmentation is attributed to the intronic SNP rs12913832 in the HERC2 gene, which interacts with the promoter region of the contiguous OCA2 gene. This interaction, through the formation of a chromatin loop, modulates the transcriptional activity of OCA2, directly affecting eye color pigmentation. Recent advancements in technology have elucidated the precise spatial organization of the genome within the cell nucleus, with chromatin architecture playing a pivotal role in regulating various genome functions. In this study, we investigated the organization of the chromatin close to the HERC2/OCA2 locus in human lymphocyte nuclei using fluorescence in situ hybridization (FISH) and high-throughput chromosome conformation capture (Hi-C) data. The 3 Mb of genomic DNA that belonged to the chromosomal region 15q12-q13.1 revealed the presence of three contiguous chromatin loops, which exhibited a different level of compaction depending on the presence of the A or G allele in the SNP rs12913832. Moreover, the analysis of the genomic organization of the genes has demonstrated that this chromosomal region is evolutionarily highly conserved, as evidenced by the analysis of syntenic regions in species from other Vertebrate classes. Thus, the role of rs12913832 variant is relevant not only in determining the transcriptional activation of the OCA2 gene but also in the chromatin compaction of a larger region, underscoring the critical role of chromatin organization in the proper regulation of the involved genes. It is crucial to consider the broader implications of this finding, especially regarding the potential regulatory role of similar polymorphisms located within intronic regions, which do not influence the same gene by modulating the splicing process, but they regulate the expression of adjacent genes. Therefore, caution should be exercised when utilizing whole-exome sequencing for diagnostic purposes, as intron sequences may provide valuable gene regulation information on the region where they reside. Thus, future research efforts should also be directed towards gaining a deeper understanding of the precise mechanisms underlying the role and mode of action of intronic SNPs in chromatin loop organization and transcriptional regulation.

Apoptotic chromatin condensation inducer in the nucleus: Genome-wide analysis in plants and expression profile during Cowpea Severe Mosaic Virus infection in Vigna unguiculata [L.] Walp

Apoptotic Chromatin Condensation Inducer in the Nucleus (ACIN1) is a scaffold protein that was first described as a complex component responsible for triggering apoptosis in human cells. In plants, ACIN1 participates in silencing of Flowering Locus C (FLC), involved in vernalization in Arabidopsis thaliana. Contrary to what has been observed for humans, there are no reports on ACIN1 linked to programmed cell death (PCD) in plants. Actually, the function of ACIN1 in plants is still poorly understood. In the present study, a genome-wide analysis of the ACIN1 gene family in plants identified 27 ACIN1 orthologs from 19 species belonging to 12 plant families. The phylogenetic relationships, physicochemical properties, gene structure, conserved motifs, promoter cis-elements, chromosomal localization, syntenic regions, and protein network were investigated. Altogether, these analyzes revealed highly conserved domains in the structure of the ACIN1 proteins, as well as putative metacaspase cleavage sites, which suggest that they play a conserved function probably associated with the programmed cell death in plants. For instance, differential expression pattern and modulation of ACIN1 were noticed after inoculation of cowpea with Cowpea severe mosaic virus (CPSMV). Therefore, this study was conducted to provide, for the first time, information on the evolutionary, structural, and functional characteristics of the ACIN1 gene family as an initial effort towards understanding the role of these proteins in studied plant development and stress responses.

Increased hippocampal epigenetic age in the Ts65Dn mouse model of Down Syndrome.

Down syndrome (DS) is a segmental progeroid genetic disorder associated with multi-systemic precocious aging phenotypes, which are particularly evident in the immune and nervous systems. Accordingly, people with DS show an increased biological age as measured by epigenetic clocks. The Ts65Dn trisomic mouse, which harbors extra-numerary copies of chromosome 21 (Hsa21)-syntenic regions, was shown to recapitulate several progeroid features of DS, but no biomarkers of age have been applied to it so far. In this pilot study, we used a mouse-specific epigenetic clock to measure the epigenetic age of hippocampi from Ts65Dn and euploid mice at 20 weeks. Ts65Dn mice showed an increased epigenetic age in comparison with controls, and the observed changes in DNA methylation partially recapitulated those observed in hippocampi from people with DS. Collectively, our results support the use of the Ts65Dn model to decipher the molecular mechanisms underlying the progeroid DS phenotypes.

A de novo genome assembly of Solanum bulbocastanum Dun., a Mexican diploid species reproductively isolated from the A-genome species, including cultivated potatoes.

Potato and its wild relatives are distributed mainly in the Mexican highlands and central Andes of South America. The South American A-genome species, including cultivated potatoes, are reproductively isolated from Mexican diploid species. Whole-genome sequencing has disclosed genome structure and similarity, mostly in cultivated potatoes and their closely related species. In this study, we generated a chromosome-scale assembly of the genome of a Mexican diploid species, Solanum bulbocastanum Dun., using PacBio long-read sequencing, optical mapping, and Hi-C scaffolding technologies. The final sequence assembly consisted of 737.9 Mb, among which 647.0 Mb were anchored to the 12 chromosomes. Compared with chromosome-scale assemblies of S. lycopersicum (tomato), S. etuberosum (non-tuber-bearing species with E-genome), S. verrucosum, S. chacoense, S. multidissectum, and S. phureja (all four are A-genome species), the S. bulbocastnum genome was the shortest. It contained fewer transposable elements (56.2%) than A-genome species. A cluster analysis was performed based on pairwise ratios of syntenic regions among the seven chromosome-scale assemblies, showing that the A-genome species were first clustered as a distinct group. Then, this group was clustered with S. bulbocastanum. Sequence similarity in 1,624 single-copy orthologous gene groups among 36 Solanum species and clones separated S. bulbocastanum as a specific group, including other Mexican diploid species, from the A-genome species. Therefore, the S. bulbocastanum genome differs in genome structure and gene sequences from the A-genome species. These findings provide important insights into understanding and utilizing the genetic diversity of S. bulbocastanum and the other Mexican diploid species in potato breeding.

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence.

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Epigenetic Regulation of Cardiomyocyte Maturation by Arginine Methyltransferase CARM1.

During the neonatal stage, the cardiomyocyte undergoes a constellation of molecular, cytoarchitectural, and functional changes known collectively as cardiomyocyte maturation to increase myocardial contractility and cardiac output. Despite the importance of cardiomyocyte maturation, the molecular mechanisms governing this critical process remain largely unexplored. We leveraged an in vivo mosaic knockout system to characterize the role of Carm1, the founding member of protein arginine methyltransferase, in cardiomyocyte maturation. Using a battery of assays, including immunohistochemistry, immuno-electron microscopy imaging, and action potential recording, we assessed the effect of loss of Carm1 function on cardiomyocyte cell growth, myofibril expansion, T-tubule formation, and electrophysiological maturation. Genome-wide transcriptome profiling, H3R17me2a chromatin immunoprecipitation followed by sequencing, and assay for transposase-accessible chromatin with high-throughput sequencing were used to investigate the mechanisms by which CARM1 (coactivator-associated arginine methyltransferase 1) regulates cardiomyocyte maturation. Finally, we interrogated the human syntenic region to the H3R17me2a chromatin immunoprecipitation followed by sequencing peaks for single-nucleotide polymorphisms associated with human heart diseases. We report that mosaic ablation of Carm1 disrupts multiple aspects of cardiomyocyte maturation cell autonomously, leading to reduced cardiomyocyte size and sarcomere thickness, severe loss and disorganization of T tubules, and compromised electrophysiological maturation. Genomics study demonstrates that CARM1 directly activates genes that underlie cardiomyocyte cytoarchitectural and electrophysiological maturation. Moreover, our study reveals significant enrichment of human heart disease-associated single-nucleotide polymorphisms in the human genomic region syntenic to the H3R17me2a chromatin immunoprecipitation followed by sequencing peaks. This study establishes a critical and multifaceted role for CARM1 in regulating cardiomyocyte maturation and demonstrates that deregulation of CARM1-dependent cardiomyocyte maturation gene expression may contribute to human heart diseases.

MYH13, a superfast myosin expressed in extraocular, laryngeal and syringeal muscles.

MYH13 is a unique type of sarcomeric myosin heavy chain (MYH) first detected in mammalian extraocular (EO) muscles and later also in vocal muscles, including laryngeal muscles of some mammals and syringeal muscles of songbirds. All these muscles are specialized in generating very fast contractions while producing relatively low force, a design appropriate for muscles acting against a much lower load than most skeletal muscles inserting into the skeleton. The definition of the physiological properties of muscle fibres containing MYH13 has been complicated by the mixed fibre type composition of EO muscles and the coexistence of different MYH types within the same fibre. A major advance in this area came from studies on isolated recombinant myosin motors and the demonstration that the affinity of actin-bound human MYH13 for ADP is much weaker than those of fast-type MYH1 (type 2X) and MYH2 (type 2A). This property is consistent with a very fast detachment of myosin from actin, a major determinant of shortening velocity. The MYH13 gene arose early during vertebrate evolution but was characterized only in mammals and birds and appears to have been lost in some teleost fish. The MYH13 gene is located at the 3' end of the mammalian fast/developmental gene cluster and in a similar position to the orthologous cluster in syntenic regions of the songbird genome. MYH13 gene regulation is controlled by a super-enhancer in the mammalian locus and deletion of the neighbouring fast MYH1 and MYH4 genes leads to abnormal MYH13 expression in mouse leg muscles.

TRNA-Cys gene clusters exhibit high variability in Arabidopsis thaliana

Although most of the genes encoding tRNAs in plants are dispersed throughout the genome, a fraction of them form tRNA gene clusters. In Arabidopsis thaliana, the smallest of tRNA clusters on chromosome 5 consists of four tRNA-Cys-GCA genes placed within repeating units of 0.4 kbp. A systematic analysis of the genomic sequences of syntenic regions from various ecotypes of A. thaliana showed that the general structure of the cluster, consisting of a tRNA-Cys pseudogene followed by repeating units containing tRNA-Cys genes, is well conserved. However, there is significant heterogeneity in the number of repeating units between different ecotypes. A unique feature of this cluster is the presence of putative transposable elements (Helitron). In addition, two further tRNA-Cys gene mini-clusters (gene pairs) in A. thaliana were identified. RNA-seq-based evaluation of expression of tRNA-Cys-GCA genes showed a positive signal for 11 out of 13 unique transcripts. An analysis of the conservation of the tRNA-Cys clusters from A. thaliana with the corresponding regions from four other Arabidopsis species suggests a sequence of events that led to the divergence of these regions.

Identification of QTLs for resistance to 10 pathotypes of Plasmodiophora brassicae in Brassica oleracea cultivar ECD11 through genotyping-by-sequencing

Key messageTwo major quantitative trait loci (QTLs) and five minor QTLs for 10 pathotypes were identified on chromosomes C01, C03, C04 and C08 through genotyping-by-sequencing from Brassica oleracea.Clubroot caused by Plasmodiophora brassicae is an important disease in brassica crops. Managing clubroot disease of canola on the Canadian prairie is challenging due to the continuous emergence of new pathotypes. Brassica oleracea is considered a major source of quantitative resistance to clubroot. Genotyping-by-sequencing (GBS) was performed in the parental lines; T010000DH3 (susceptible), ECD11 (resistant) and 124 BC1 plants. A total of 4769 high-quality polymorphic SNP loci were obtained and distributed on 9 chromosomes of B. oleracea. Evaluation of 124 BC1S1 lines for resistance to 10 pathotypes: 3A, 2B, 5C, 3D, 5G, 3H, 8J, 5K, 5L and 3O of P. brassicae, was carried out. Seven QTLs, 5 originating from ECD11 and 2 from T010000DH3, were detected. One major QTL designated as Rcr_C03-1 on C03 contributed 16.0–65.6% of phenotypic variation explained (PVE) for 8 pathotypes: 2B, 5C, 5G, 3H, 8J, 5K, 5L and 3O. Another major QTL designated as Rcr_C08-1 on C08 contributed 8.3 and 23.5% PVE for resistance to 8J and 5K, respectively. Five minor QTLs designated as Rcr_C01-1, Rcr_C03-2, Rcr_C03-3, Rcr_C04-1 and Rcr_C08-2 were detected on chromosomes C01, C03, C04 and C08 that contributed 8.3–23.5% PVE for 5 pathotypes each of 3A, 2B, 3D, 8J and 5K. There were 1, 10 and 4 genes encoding TIR-NBS-LRR/CC-NBS-LRR class disease resistance proteins in the Rcr_C01-1, Rcr_C03-1 and Rcr_C08-1 flanking regions. The syntenic regions of the two major QTLs Rcr_C03-1 and Rcr_C08-1 in the B. rapa genome ‘Chiifu’ were searched.

Integrated Genomic Analysis Identifies a Downstream Regulatory Element Oncogenically Activating Hematopoietically Expressed Homeobox ( HHEX) in ETP-ALL

Hematopoietically-expressed Homeobox ( HHEX) is required for the maintenance of hematopoietic stem cells (HSPCs) and common lymphoid progenitor cells. HHEX is also the second most frequent integration site in retroviral insertional mutagenesis screens of leukemias and lymphomas arising in AKXD recombinant inbred mouse strains, implying that it is an important oncogene. Enforced expression of HHEX induces T-cell acute lymphoblastic leukemias (T-ALL) in murine bone marrow transplantation models. In human studies, HHEXmRNA expression is upregulated in human T-ALL studies, especially in Early T-cell Precursor (ETP-) ALL subtypes where it is concordantly expressed with LMO2. The HHEX locus is not rearranged in ETP-ALL where it is upregulated so the mechanism of HHEX activation is not clear. We analyzed T-ALL induced by retroviral mutagenesis and found an intergenic site of frequent integration, 65kb 3‘ of the Hhex coding exons. Integrations clustered within a 3 kb genomic area that is also an open chromatin region (OCR) by ATAC-seq analysis and that induced Hhex upregulation. The syntenic human region was highly enriched for H3K27 acetylation and for occupancy by multiple transcription factors, notably GATA1 and TAL1, per ENCODE. Our prior data had shown that HHEX was a downstream target of the LMO2 complex comprised of GATA1/2, TAL1 or LYL1, and nucleated by the scaffolding protein, LIM domain binding protein 1 (LDB1). We had previously shown by ChIP-PCR that the LMO2 complex was bound to the intron 1 enhancer. Thus, we performed ChIP-exonuclease analysis of LDB1 in the human ETP-ALL model cell line, LOUCY. We confirmed its occupancy at the intron 1 enhancer but we also observed LDB1 occupancy 68kb 3‘ of the HHEX coding exons, in the exact region syntenic to the murine common insertion sites. ChIP-exo analysis allowed us to pinpoint a core element with composite GATA/E box sites that was highly conserved across multiple mammalian species. LDB1 occupancy at intron 1 and at +68kb of the HHEX locus was reminiscent of LDB1's occupancy at the beta globin locus where LDB1's occupancy and homodimerization mediates chromatin looping between composite E box/GATA sites in the locus control region (LCR) and the beta globin proximal promoters. To test whether a similar interaction was occurring at the HHEX locus, we analyzed virtual 4C (chromatin conformation capture) data. We confirmed looping between the +68kb element and the proximal promoter of HHEX in human ETP-ALL primary samples and in human HSPCs. The experiments showed that looping occurs in certain developmental contexts in HSPCs and is recapitulated in ETP-ALL and mediated by the LMO2/LDB1 protein complex. HHEX is a downstream oncogene of the LMO2/LDB1 complex and is activated through this distal +68kb regulatory element. Our studies raise the possibility that dissection of this regulatory element and interference with chromatin looping are potential therapeutic mechanisms that could disrupt oncogene expression.

Identification of Genomic Regions Implicated in Susceptibility to Schistosoma mansoni Infection in a Murine Backcross Genetic Model.

Only a small number of infected people are highly susceptible to schistosomiasis, showing high levels of infection or severe liver fibrosis. The susceptibility to schistosome infection is influenced by genetic background. To assess the genetic basis of susceptibility and identify the chromosomal regions involved, a backcross strategy was employed to generate high variation in schistosomiasis susceptibility. This strategy involved crossing the resistant C57BL/6J mouse strain with the susceptible CBA/2J strain. The resulting F1 females (C57BL/6J × CBA/2J) were then backcrossed with CBA/2J males to generate the backcross (BX) cohort. The BX mice exhibited a range of phenotypes, with disease severity varying from mild to severe disease, lacking a fully resistant group. We observed four levels of infection intensity using cluster and principal component analyses and K-means based on parasitological, pathological, and immunological trait measurements. The mice were genotyped with 961 informative SNPs, leading to the identification of 19 new quantitative trait loci (QTL) associated with parasite burden, liver lesions, white blood cell populations, and antibody responses. Two QTLs located on chromosomes 15 and 18 were linked to the number of granulomas, liver lesions, and IgM levels. The corresponding syntenic human regions are located in chromosomes 8 and 18. None of the significant QTLs had been reported previously.