Molecular Genetics and Genomics

Establishment of adaptability to the northern-limit of rice production Abstract The domestication of cultivated crops from their wild relatives narrowed down their genetic diversity in a bottleneck effect. Subsequently, the cultivation areas of crops have expanded all over the world into various environmental conditions from the original area along with human migration after domestication. Here, we demonstrated the genetic changes in the adaptation of rice to Hokkaido (41°2–45°3N latitude), Japan, from the tropics of their origin in Asian cultivated rice, Oryza sativa L. Although cultivated rice originated from the tropics, Hokkaido is one of the northern-limits of rice cultivation worldwide. Population genomics focusing on the local populations showed the varieties had genetically distinct classes with limited genetic diversity. In addition, some varieties in the class carried unique genotypes for flowering time, exhibiting extremely early flowering time. Certain mutations in unique genotypes can split off the varieties that are able to grow in Hokkaido. Furthermore, the changes in the genotype for flowering time during rice cultivation in Hokkaido demonstrated novel combinations of genes for flowering time owing to the intensive artificial selection on natural variation and rice breeding programs to achieve stable rice production in Hokkaido.

Genotype–phenotype correlation and expansion of orodental anomalies in LTBP3 -related disorders Abstract The latent transforming growth factor-beta-binding protein 3 (LTBP3), encoding extracellular matrix proteins, plays a role in skeletal formation. Mutations in LTBP3 have been associated with various types of skeletal dysplasia. We aimed to characterize clinical and molecular features of more patients with mutations in the gene, which may help suggest genotype–phenotype correlation. The first two East Asian patients with short stature, heart defects, and orodental anomalies having LTBP3 mutations were identified. Whole exome and Sanger sequencing revealed that the one with a novel heterozygous missense (c.2017G>T, p.Gly673Cys) mutation in LTBP3 had clinical features consistent with acromicric dysplasia (ACMICD). The variant was located in the highly conserved EGF-like calcium-binding domain adjacent to the single reported LTBP3 variant associated with ACMICD. This finding supports that LTBP3 is a disease gene for ACMICD. Another patient with a novel homozygous splice site acceptor (c.1721-2A>G) mutation in LTBP3 was affected with dental anomalies and short stature (DASS). Previously undescribed orodental features included multiple unerupted teeth, high-arched palate, and microstomia found in our patient with ACMICD, and extensive dental infection, condensing osteitis, and deviated alveolar bone formation in our patient with DASS. Our results and comprehensive reviews suggest a genotype–phenotype correlation: biallelic loss-of-function mutations cause DASS, monoallelic missense gain-of-function mutations in the EGF-like domain cause ACMICD, and monoallelic missense gain-of-function mutations with more drastic effects on the protein functions cause geleophysic dysplasia (GPHYSD3). In summary, we expand the phenotypic and genotypic spectra of LTBP3-related disorders, support that LTBP3 is a disease gene for ACMICD, and propose the genotype–phenotype correlation of LTBP3 mutations.

Genome-wide association analysis for body weight identifies candidate genes related to development and metabolism in rainbow trout ( Oncorhynchus mykiss ) Abstract Growth is one of the most important traits from both a physiological and economic perspective in aquaculture species. Thus, identifying the genomic regions and genes underpinning genetic variation for this trait is of particular interest in several fish species, including rainbow trout. In this work, we perform a genome-wide association study (GWAS) to identify the genomic regions associated with body weight at tagging (BWT) and at 18 months (BW18M) using a dense SNP panel (57 k) and 4596 genotyped rainbow trout from 105 full-sib families belonging to a Chilean breeding population. Analysis was performed by means of single-step GBLUP approach. Genetic variance explained by 20 adjacent SNP windows across the whole genome is reported. To further explore candidate genes, we focused on windows that explained the highest proportion of genetic variance in the top 10 chromosomes for each trait. The main window from the top 10 chromosomes was explored by BLAST using the first and last SNP position of each window to determine the target nucleotide sequence. As expected, the percentage of genetic variance explained by windows was relatively low, due to the polygenic nature of body weight. The most important genomic region for BWT and BW18M were located on chromosomes 15 and 24 and they explained 2.14% and 3.02% of the genetic variance for each trait, respectively. Candidate genes including several growth factors, genes involved in development of skeletal muscle and bone tissue and nutrient metabolism were identified within the associated regions for both traits BWT and BW18M. These results indicate that body weight is polygenic in nature in rainbow trout, with the most important loci explaining as much as 3% of the genetic variance for the trait. The genes identified here represent good candidates for further functional validation to uncover biological mechanisms underlying variation for growth in rainbow trout.

Genome-wide association study of inflorescence length of cultivated soybean based on the high-throughout single-nucleotide markers Abstract As an important and complex trait, inflorescence length (IL) of soybean [Glycine max (L.) Merr.] significantly affected seed yields. Therefore, elucidating molecular basis of inflorescence architecture, especially for IL, was important for improving soybean yield potentials. Longer IL meaned to have more pod and seed in soybean. Hence, increasing IL and improving yield are targets for soybean breeding. In this study, a association panel, comprising 283 diverse samples, was used to dissect the genetic basis of IL based on genome-wide association analysis (GWAS) and haplotype analysis. GWAS and haplotype analysis were conducted through high-throughout single-nucleotide polymorphisms (SNP) developed by SLAF-seq methodology. A total of 39, 057 SNPs (minor allele frequency ≥ 0.2 and missing data ≤ 10%) were utilized to evaluate linkage disequilibrium (LD) level in the tested association panel. A total of 30 association signals were identified to be associated with IL via GWAS. Among them, 13 SNPs were novel, and another 17 SNPs were overlapped or located near the linked regions of known quantitative trait nucleotide (QTN) with soybean seed yield or yield component. The functional genes, located in the 200-kb genomic region of each peak SNP, were considered as candidate genes, such as the cell division/ elongation, specific enzymes, and signaling or transport of specific proteins. These genes have been reported to participant in the regulation of IL. Ten typical long-IL lines and ten typical short-IL lines were re-sequencing, and then, six SNPs from five genes were obtained based on candidate gene-based association. In addition, 42 haplotypes were defined based on haplotype analysis. Of them, 11 haplotypes were found to regulate long IL (> 14 mm) in soybean. The identified 30 QTN with beneficial alleles and their candidate genes might be valuable for dissecting the molecular mechanisms of IL and further improving the yield potential of soybean.

Exploring molecular tools for transformation and gene expression in the cultivated edible mushroom Agrocybe aegerita Abstract Agrocybe aegerita is a cultivated edible mushroom in numerous countries, which also serves as a model basidiomycete to study fruiting body formation. Aiming to create an easily expandable customised molecular toolset for transformation and constitutive gene of interest expression, we first created a homologous dominant marker for transformant selection. Progeny monokaryons of the genome-sequenced dikaryon A. aegerita AAE-3 used here were identified as sensitive to the systemic fungicide carboxin. We cloned the wild-type gene encoding the iron-sulphur protein subunit of succinate dehydrogenase AaeSdi1 including its up- and downstream regions, and introduced a single-point mutation (His237 to Leu) to make it confer carboxin resistance. PEG-mediated transformation of protoplasts derived from either oidia or vegetative monokaryotic mycelium with the resulting carboxin resistance marker (CbxR) plasmid pSDI1E3 yielded carboxin-resistant transformants in both cases. Plasmid DNA linearised within the selection marker resulted in transformants with ectopic multiple insertions of plasmid DNA in a head-to-tail repeat-like fashion. When circular plasmid was used, ectopic single integration into the fungal genome was favoured, but also gene conversion at the homologous locus was seen in 1 out of 11 analysed transformants. Employing CbxR as selection marker, two versions of a reporter gene construct were assembled via Golden Gate cloning which allows easy recombination of its modules. These consisted of an eGFP expression cassette controlled by the native promoter PAaeGPDII and the heterologous terminator Tnos, once with and once without an intron in front of the eGFP start codon. After protoplast transformation with either construct as circular plasmid DNA, GFP fluorescence was detected with either transformants, indicating that expression of eGFP is intron-independent in A. aegerita. This paves the way for functional genetics approaches to A. aegerita, e.g., via constitutive expression of fruiting-related genes.

Comparative genomics reveals origin of MIR159A – MIR159B paralogy, and complexities of PTGS interaction between miR159 and target GA-MYBs in Brassicaceae Abstract Whole-genome and segmental duplications coupled with sequence and functional diversification are responsible for gene family expansion, and morphological and adaptive diversity. Although broad contours of such processes are understood, detailed investigations on regulatory elements, such as miRNA-transcription factor modules, especially in non-model crop plants with complex genomes, are few. The present study was performed to understand evolutionary history of MIR159 family, and changes in the miRNA-binding site (MBS) of the targets MYB33, MYB65, and MYB101 that may affect post-transcriptional gene silencing. We established orthology and paralogy between members of MIR159 family by reconstructing the phylogeny based on 240 precursor sequences sampled across green plants. An unambiguous paralogous relationship between MIR159Aand MIR159B was observed only in Brassicaceae which prompted us to analyze the origin of this paralogy. Comparative micro-synteny of ca. 100 kb genomic segments surrounding MIR159A, MIR159B, and MIR159C loci across 15 genomes of Brassicaceae revealed segmental duplication that occurred in the common ancestor of Brassicaceae to be responsible for origin of MIR159A–MIR159B paralogy; extensive gene loss and rearrangements were also encountered. The impact of polyploidy was revealed when the three sub-genomes—least fractionated (LF), moderately fractionated (MF1), and most fractionated (MF2) sub-genomes of Brassicaand Camelina sativa—were analyzed. Extensive gene loss was observed among sub-genomes of Brassica, whereas those in Camelina were largely conserved. Analysis of the target MYBs revealed the complete loss of MYB33 homologs in a Brassica lineage-specific manner. Our findings suggest that mature miR159a/b /c are capable of targeting MYB65 across Brassicaceae, MYB33 in all species except Brassica, and MYB101 only in Arabidopsis thaliana. Comparative analysis of the mature miRNA sequence and the miRNA-binding site (MBS) in MYB33, MYB65, and MYB101 showed the complexity of regulatory network that is dependent on strict sequence complementarity potentially leading to regulatory diversity.

Genomic copy number variation of the CHKB gene alters gene expression and affects growth traits of Chinese domestic yak ( Bos grunniens ) breeds Abstract Copy number variation (CNV) influences the mRNA transcription levels and phenotypic traits through gene dosage, position effects, alteration of downstream pathways, and modulation of the structure and position of chromosomes. A previous study using the read depth approach to genome resequencing analysis revealed CNVs of the choline kinase beta (CHKB) gene in the copy number variable regions (CNVRs) of yak breeds may influence muscle development and therefore the phenotypic traits of yak breeds. Further work is required to attain a more complete understanding and validate the importance of the detected CNVR of the CHKB gene found in yak breeds, because there is no association studies of the CHKB gene with yak growth traits that have been reported. The goal of this study was to determine the distribution of CHKB copy numbers in five Chinese domestic yak breeds and evaluate their impact on gene expression and growth traits. The data were analyzed using real-time quantitative PCR. In this study, the normal CNV of the CHKB gene was found to be significantly (p < 0.05) associated with greater chest girth and body weight for three age groups of Datong yaks. Our results indicated that the copy number of the CHKB gene is negatively correlated with the mRNA expression level. From this result, we conclude that CNVs of the CHKB gene could be novel markers for growth traits of Chinese domestic yak breeds and might therefore provide a novel opportunity to utilize data on CNVs in designing molecular markers for the selection of animal breeding programs for larger populations of various yaks.

Galactitol catabolism in Sinorhizobium meliloti is dependent on a chromosomally encoded sorbitol dehydrogenase and a pSymB-encoded operon necessary for tagatose catabolism Abstract The legume endosymbiont Sinorhizobium meliloti can utilize a broad range of carbon compounds to support its growth. The linear, six-carbon polyol galactitol is abundant in vascular plants and is metabolized in S. melilotiby the contribution of two loci SMb21372-SMb21377 and SMc01495-SMc01503 which are found on pSymB and the chromosome, respectively. The data suggest that several transport systems, including the chromosomal ATP-binding cassette (ABC) transporter smoEFGK, contribute to the uptake of galactitol, while the adjacent gene smoS encodes a protein for oxidation of galactitol into tagatose. Subsequently, genes SMb21374 and SMb21373, encode proteins that phosphorylate and epimerize tagatose into fructose-6-phosphate, which is further metabolized by the enzymes of the Entner–Doudoroff pathway. Of note, it was found that SMb21373, which was annotated as a 1,6-bis-phospho-aldolase, is homologous to the E. coli gene gatZ, which is annotated as encoding the non-catalytic subunit of a tagatose-1,6-bisphosphate aldolase heterodimer. When either of these genes was introduced into an Agrobacterium tumefaciens strain that carries a tagatose-6-phosphate epimerase mutation, they are capable of complementing the galactitol growth deficiency associated with this mutation, strongly suggesting that these genes are both epimerases. Phylogenetic analysis of the protein family (IPR012062) to which these enzymes belong, suggests that this misannotation is systemic throughout the family. S. meliloti galactitol catabolic mutants do not exhibit symbiotic deficiencies or the inability to compete for nodule occupancy.

Construction of genetic linkage map and genome dissection of domestication-related traits of moth bean ( Vigna aconitifolia ), a legume crop of arid areas Abstract The moth bean (Vigna aconitifolia), possibly the most primitive crop of the genus Vigna, is a highly drought- and heat-resistant legume grown in arid areas. Moth bean domestication involved phenotypic changes, including reduction of seed dormancy and pod shattering, increased organ size, and earlier flowering and maturity. However, the genetics of the domestication process in moth bean is not known. In this study, we constructed a genetic linkage map for moth bean and used the map to identify quantitative trait loci (QTL) for domestication-related traits of an F2 population of 188 individuals produced from a cross of wild moth bean (TN67) and cultivated moth bean (ICPMO056). The genetic linkage map comprised 11 linkage groups (LG) of 172 simple sequence repeat markers and spanned a total length of 1016.8 centiMorgan (cM), with an average marker distance of 7.34 cM. A comparative genome analysis showed high genome synteny between moth bean and mungbean (Vigna radiata), adzuki bean (Vigna angularis), rice bean (Vigna umbellata), and yardlong bean (Vigna unguiculata). In total, 50 QTLs and 3 genes associated with 20 domestication-related traits were identified. Most of the QTLs belonged to five LGs (1, 2, 4, 7, and 10). Key traits related to domestication such as seed dormancy and pod shattering were controlled by large-effect QTLs (PVE > 20%) with one or two minor QTLs, whereas all other traits were controlled by one–seven minor QTLs, apart from seed weight, which was controlled by one major and seven minor QTLs. These results suggest that a small number of mutations with large phenotypic effects have contributed to the domestication of the moth bean. Comparative analysis of QTLs with related Vigna crops revealed that there are several domestication-related large-effect QTLs that had not been used in moth bean domestication. This study provides a basic genetic map and identified genome regions associated with domestication-related traits, which will be useful for the genetic improvement of the moth bean and related Vigna species.

Comparative analysis of the mitochondrial genomes of Orthonectida: insights into the evolution of an invertebrate parasite species Abstract Among invertebrates, only a few groups still have uncertain phylogenetic position, Orthonectida, a small group of rare multi-cellular parasites of marine invertebrates, being one of them. Recent molecular and morphological findings suggest that orthonectids belong to Lophotrochozoa and are close to Annelida. Nevertheless, phylogenetic relationships between orthonectids and annelids are unclear, and the phylogeny within the group itself has never been studied. Sequencing of mitochondrial genomes is used here to clarify this issue. Complete mt genomes of the orthonectids Intoshia variabili and Rhopalura litoralis were characterized and compared with Intoshia linei mt genome. Our results show that Orthonectida mt genomes have undergone reduction and gene loss, and that they have complicated organization revealed in strand asymmetry in nucleotide composition, in some features of intergenic non-coding regions, tRNA duplication and folding. Moreover, all species of Orthonectida have a unique gene order with complicated rearrangement landscape. Significant differences in mitochondrial genomes in the three orthonectid species could be explained by the fact that their host species belong to different taxa (flat worms, nemertines and gastropods). Among the analyzed mt genomes of Orthonectida, I. linei possesses the closest gene order to the ancestral genome. All Orthonectida species are monophyletic, and in the phylogenetic tree are close to Pleistoannelida, and specifically, to Clitellata.

Alexandros Sfakianakis
Anapafseos 5 . Agios Nikolaos
Crete.Greece.72100
2841026182

6948891480

alsfakia@gmail.com