Potential DNA Marker Research Articles

Development of DNA Markers for Acute Hepatopancreatic Necrosis Disease Tolerance in Litopenaeus vannamei through a Genome-Wide Association Study.

Shrimp aquaculture is facing a serious disease, acute hepatopancreatic necrosis disease (AHPND), caused by Vibrio paraheamolyticus (VpAPHND). For sustainable shrimp aquaculture, massive losses of shrimp infected with VpAPHND must be prevented. Research and selection of shrimp tolerant to VpAPHND infection is a sustainable approach to reducing the risk of AHPND. This study focused on the identification and development of potential DNA markers associated with AHPND using DArT sequencing (DArTSeq) and a genome-wide association study. Three populations of post-larval Litopenaeus vannamei were immersed in VpAPHND to collect susceptible (D) and tolerant (S) samples. The 45 D and 48 S shrimp had their genotypes analyzed using DArTSeq. A total of 108,983 SNPs and 17,212 InDels were obtained from the DArTseq data, while the biallelic 516 SNPs and 2293 InDels were finally filtered with PIC < 0.1, MAF < 0.05, and a call rate ≥ 80%. The filtered variants were analyzed for their association with AHPND tolerance. Although there were no significantly associated SNPs and InDels above the Bonferroni correction threshold, candidate variants, four SNPs and 17 InDels corresponding to p < 0.01, were provided for further validation of the AHPND tolerance trait. The candidate SNPs are located on an exon of the zinc finger protein 239-like gene, an intron of an uncharacterized gene, and in intergenic regions. Most of the candidate InDels are in the intergenic regions, with fewer in the intronic and exonic regions. This study provides information on SNPs and InDels for white shrimp. These markers will support the variant database of shrimp and be useful in shrimp aquaculture for breeding selection.

Plastome structure and phylogenetic relationships of genus Hydrocotyle (apiales): provide insights into the plastome evolution of Hydrocotyle

BackgroundThe genus Hydrocotyle Tourn. ex L. is a key group for further study on the evolution of Apiales, comprising around 170 species globally. Previous studies mainly focused on separate sections and provided much information about this genus, but its infrageneric relationships are still confusing. In addition, the genetic basis of its adaptive evolution remains poorly understood. To investigate the phylogeny and evolution of the genus, we selected ten representative species covering two of three diversity distribution centers and exhibiting rich morphology diversity. Comparative plastome analysis was conducted to clarify the structural character of Hydrocotyle plastomes. Positive selection analyses were implemented to assess the evolution of the genus. Phylogenetic inferences with protein-coding sequences (CDS) of Hydrocotyle and 17 related species were also performed.ResultsPlastomes within Hydrocotyle were generally conservative in structure, gene order, and size. A total of 14 regions (rps16-trnK, trnQ-rps16, atpI-atpH, trnC-petN-psbM, ycf3-trnS, accD-psaI-ycf4, petA-psbJ, rps12-rpl20, rpl16 intron, rps3-rpl16 intron, rps9-rpl22, ndhF-rpl32, ndhA intron, and ycf1a) were recognized as hotspot regions within the genus, which suggested to be promising DNA barcodes for global phylogenetic analysis of Hydrocotyle. The ycf15 gene was suggested to be a protein-coding gene for Hydrocotyle species, and it could be used as a DNA barcode to identify Hydrocotyle. In phylogenetic analysis, three monophyletic clades (Clade I, II, III) were identified with evidence of rapid radiation speciation within Clade I. The selective pressure analysis detected that six CDS genes (ycf1b, matK, atpF, accD, rps14, and psbB) of Hydrocotyle species were under positive selection. Within the genus, the last four genes were conservative, suggesting a relation to the unique evolution of the genus in Apiales. Seven genes (atpE, matK, psbH, ycf1a, ycf1b, rpoA, and ycf2) were detected to be under some degree of positive selection in different taxa within the genus Hydrocotyle, indicating their role in the adaptive evolution of species.ConclusionsOur study offers new insights into the phylogeny and adaptive evolution of Hydrocotyle. The plastome sequences could significantly enhance phylogenetic resolution and provide genomic resources and potential DNA markers useful for future studies of the genus.

Characterization of chicken Relaxin3 gene: mRNA expression and response to reproductive hormone treatment in ovarian granulosa cells, and single nucleotide polymorphisms associated with egg laying traits in hens

As a protein structurally similar to insulin, relaxin3 (RLN3) plays a role in promoting arousal, suppressing depressive or anxious behaviors. Two studies revealed the increase of RLN3 expression during chicken follicle selection. In this study, by real-time quantitative PCR and luciferase assay, mRNA expression and single nucleotide polymorphisms (SNPs) of chicken RLN3 were investigated. The mRNA expression of chicken RLN3 was higher in the granulosa cell of hierarchal follicles (Post-GCs) than that of pre-hierarchal follicles (Pre-GCs). In Pre-GCs, the mRNA expression of chicken RLN3 was stimulated by FSH and progesterone; in Post-GCs, it was stimulated by higher concentration of estrogen and FSH, however, was inhibited by progesterone. Four SNPs including g.-655G > C, g-592G > A, g.-372T > A and g.-282G > C were identified in the critical promoter region from −1291 bp to −207 bp of chicken RLN3, among which g.-655G > C, and g-592G > A were associated with age at first laying and clutch size, respectively, in Zaozhuang Sunzhi chickens. At g.-655G > C and g-592G > A, allele C and allele A had higher transcriptional activity, respectively. These data suggest that RLN3 plays an important role in chicken follicle development and SNPs in its promoter region are potential DNA markers for improving egg production traits.

Comparative analysis the chloroplast genomes of Celastrus (Celastraceae) species: Provide insights into molecular evolution, species identification and phylogenetic relationships

BackgroundThe genus Celastrus is an important medicinal plant resource. The similarity of morphology and the lack of complete chloroplast genome analysis have significantly impeded the exploration of species identification, molecular evolution and phylogeny of Celastrus. PurposeIn order to resolve the phylogenic controversy of Celastrus species, the chloroplast genome comparative analysis was performed to provide genetic evidence. MethodsIn this study, we collected and sequenced ten chloroplast genomes of Celastrus species from China and downloaded three chloroplast genomes from the databases. The chloroplast genomes were compared and analyzed to explore their characteristics and evolution. Furthermore, the phylogenetic relationships of Celastrus species were inferred based on the whole chloroplast genomes and protein-coding genes. ResultsAll the 13 Celastrus species chloroplast genomes showed a typical quadripartite structure with genome sizes ranging from 155,113 to 157,366 bp. The intron loss of the rps16 gene occurred in all the 13 Celastrus species. The GC content, gene sequence, repeat types and codon bias pattern were highly conserved. Ten highly variation regions were identified, which can be used as potential DNA markers in molecular identification of Celastrus species. Eight genes, including accD, atp4, ndhB, rpoC1, rbcL, rpl2, rpl20 and ycf1, were detected to experience positive selection. Phylogenetic analysis showed that Celastrus was a monophyletic group and Tripterygium was the closest sister-group. Noteworthy, C. gemmatus Loes. and C. orbiculatus Thunb. can be discriminated using the chloroplast genome as a super barcode. The comparative and phylogenetic analysis results proposed that C. tonkinensis Pitard. was the synonym of C. hindsii Benth. ConclusionThe comparative analysis of the Celastrus chloroplast genomes can provide comprehensive genetic evidence for molecular evolution, species identification and phylogenetic relationships.

Cytochrome b gene as a potential DNA barcoding marker in ecoraces of tropical Tasar silkworm Antheraea mylitta Drury

BackgroundTropical Tasar silkworm Antheraea mylitta Drury plays a very crucial role in the economic development of the country employing around 9.76 million people in rural and semi-urban areas. Since the population of this species has spread over many states with varying environmental factors, the population has been isolated into many pockets, specifically adapted to that niche, and are called as “ecoraces”. Even though phenotypically these ecoraces may look similar to some level but the intra-generic genotypic variability between them is correlated to the economically important yield traits of their cocoons. The purpose of our study was to evaluate the effectiveness of Cytochrome b (Cyt b) gene of mitochondrial genome as a potential DNA barcoding marker for identification of these ecoraces as to are they really ecoraces, and attempts to establish the first country report of DNA barcode reference library of this silkworm species. Methods and resultsDNA was extracted from pupae of pure ecoraces collected from core ecozones. PCR amplicon was bead purified and further subjected to Sanger Sequencing. The parsimony tree was also created using the MEGA 11 software. The genetic distance, gene flow and relation among the ecoraces was evaluated. The conserved regions of Cyt b gene were identified along with the documentation of the mode and tempo of nucleotide variation and evolution. ConclusionThe present study confirms the efficiency of Cyt b gene as a potential DNA marker for identification of ecoraces of Tasar silkworm and is the first country report on barcoding of A. mylitta Drury. Apart from being used as an early-warning system for ecological degradation, this knowledge will help in assessing the evolutionary relationship of these economically important ecoraces and identify their position in the process of evolution aiding in the assessment of conservation priorities.

The Complete Chloroplast Genomes of Bulbophyllum (Orchidaceae) Species: Insight into Genome Structure Divergence and Phylogenetic Analysis.

Bulbophyllum is one of the largest genera and presents some of the most intricate taxonomic problems in the family Orchidaceae, including species of ornamental and medical importance. The lack of knowledge regarding the characterization of Bulbophyllum chloroplast (cp) genomes has imposed current limitations on our study. Here, we report the complete cp genomes of seven Bulbophyllum species, including B. ambrosia, B. crassipes, B. farreri, B. hamatum, B. shanicum, B. triste, and B. violaceolabellum, and compared with related taxa to provide a better understanding of their genomic information on taxonomy and phylogeny. A total of 28 Bulbophyllum cp genomes exhibit typical quadripartite structures with lengths ranging from 145,092 bp to 165,812 bp and a GC content of 36.60% to 38.04%. Each genome contained 125-132 genes, encompassing 74-86 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The genome arrangements, gene contents, and length were similar, with differences observed in ndh gene composition. It is worth noting that there were exogenous fragment insertions in the IR regions of B. crassipes. A total of 18-49 long repeats and 38-80 simple sequence repeats (SSRs) were detected and the single nucleotide (A/T) was dominant in Bulbophyllum cp genomes, with an obvious A/T preference. An analysis of relative synonymous codon usage (RSCU) revealed that leucine (Leu) was the most frequently used codon, while cysteine (Cys) was the least used. Six highly variable regions (rpl32-trnLUAG > trnTUGU-trnLUAA > trnFGAA-ndhJ > rps15-ycf1 > rbcL-accD > psbI-trnSGCU) and five coding sequences (ycf1 > rps12 > matK > psbK > rps15) were identified as potential DNA markers based on nucleotide diversity. Additionally, 31,641 molecular diagnostic characters (MDCs) were identified in complete cp genomes. A phylogenetic analysis based on the complete cp genome sequences and 68 protein-coding genes strongly supported that 28 Bulbophyllum species can be divided into four branches, sects. Brachyantha, Cirrhopetalum, and Leopardinae, defined by morphology, were non-monophyly. Our results enriched the genetic resources of Bulbophyllum, providing valuable information to illustrate the complicated taxonomy, phylogeny, and evolution process of the genus.

Genome-Wide Association Study of Milk Composition in Karachai Goats.

This study is first to perform a genome-wide association study (GWAS) to investigate the milk quality traits in Karachai goats. The objective of the study was to identify candidate genes associated with milk composition traits based on the identification and subsequent analysis of all possible SNPs, both genome-wide (high-confidence) and suggestive (subthreshold significance). To estimate the milk components, 22 traits were determined, including several types of fatty acids. DNA was extracted from ear tissue or blood samples. A total of 167 Karachai goats were genotyped using an Illumina GoatSNP53K BeadChip panel (Illumina Inc., San Diego, CA, USA). Overall, we identified 167 highly significant and subthreshold SNPs associated with the milk components of Karachai goats. A total of 10 SNPs were located within protein-coding genes and 33 SNPs in close proximity to them (±0.2 Mb). The largest number of genome-wide significant SNPs was found on chromosomes 2 and 8 and some of them were associated with several traits. The greatest number of genome-wide significant SNPs was identified for crude protein and lactose (6), and the smallest number-only 1 SNP-for freezing point depression. No SNPs were identified for monounsaturated and polyunsaturated fatty acids. Functional annotation of all 43 SNPs allowed us to identify 66 significant candidate genes on chromosomes 1, 2, 3, 4, 5, 8, 10, 13, 16, 18, 21, 23, 25, 26, and 27. We considered these genes potential DNA markers of the fatty acid composition of Karachai goat milk. Also, we found 12 genes that had a polygenic effect: most of them were simultaneously associated with the dry matter content and fatty acids (METTL, SLC1A 8, PHACTR1, FMO2, ECI1, PGP, ABCA3, AMDHD2). Our results suggest that the genes identified in our study affecting the milk components in Karachai goats differed from those identified in other breeds of dairy goats.

Analysis for Mitochondrial Encoded Transfer RNA Leucine2 (MT-TL2) Gene in Breast Cancer Patients

Background: Breast cancer (BC) has been linked to a variety of nuclear DNA changes as well as mitochondrial DNA (mtDNA) alterations. The study aims to evaluate/analyze the association of mitochondrial transfer RNA (mt-tRNA) leucine 2 gene with BC patients. Methods: In the current study, 24 samples have been collected from various families in Peshawar. DNA was extracted from blood. Polymerase chain reaction was used to amplify the mt-tRNA MT-TL2 gene, and 20 samples were sequenced. Results: The sequence was compared with accession #NC-012920.1 of the Revised Cambridge Reference Sequence (rCRS). The results (chromatograph, nucleotide sequence, and rCRS alignment) show mutations in mt-tRNA MT-TL2 gene in our participants is not the cause of BC. Conclusion: Yet, a significant number of BC patients must be studied, and their full mtDNA must be analyzed. This will provide an indication of the potential DNA marker that might be used to prevent BC deaths at the earliest stages.

Polymorphism and expression of GLUD1 in relation to reproductive performance in Jining Grey goats.

Understanding the molecular mechanism of mammalian reproduction (puberty and prolificacy) will play a part in improving animal reproductive performance. GLUD1 (glutamate dehydrogenase1) is important for mammalian reproduction, as shown in previous studies; however, its roles in puberty and prolificacy have rarely been reported. In this study, we designed seven pairs of primers (P1 to P7) for cloning and sequencing genomic DNA of Jining Grey goats and Liaoning Cashmere goats. Primer8 (P8) was designed to detect single nucleotide polymorphism (SNP) of the GLUD1 in both sexually precocious and high-fecundity breeds (Jining Grey, Nanjiang Brown and Matou goats) and sexually late-maturing and low-fecundity breeds (Liaoning Cashmere, Inner Mongolia Cashmere and Taihang goats) by PCR-RFLP (restriction fragment length polymorphism). The real-time quantitative polymerase chain reaction (RT-qPCR) technique was used to detect the expression of GLUD1 in a variety of tissues. The results showed that the A197C mutation was only found in the amplification product of P6. For this SNP locus, only two genotypes (AA and AC) were detected in Nanjiang Brown goats, while three genotypes (AA, AC and CC) were detected in the other five breeds. In Jining Grey goats, the frequency of genotypes AA, AC and CC was 0.69, 0.26 and 0.05, respectively. In Jining Grey goats, AA genotype had 0.54 () and 0.3 () more kids than the CC and AC genotype, respectively, and no significant difference () was found in kidding number between the AC and CC genotype. GLUD1 was expressed in five tissues of different developmental stages. The expression level of GLUD1 in the hypothalamus was higher than that in the other four tissues except during puberty of Liaoning Cashmere goats. In puberty in goats, GLUD1 expression was significantly higher in ovaries than that in the juvenile period (). RT-qPCR results showed that the expression of GLUD1 in ovaries may relate to the puberty of goats. The present study preliminarily indicated that there might be an association between the 197 locus of GLUD1 and sexual precocity in goats, and allele A of GLUD1 was a potential DNA marker for improving kidding number in Jining Grey goats.

Microbiome and plant cell transformation trigger insect gall induction in cassava.

Several specialised insects can manipulate normal plant development to induce a highly organised structure known as a gall, which represents one of the most complex interactions between insects and plants. Thus far, the mechanism for insect-induced plant galls has remained elusive. To study the induction mechanism of insect galls, we selected the gall induced by Iatrophobia brasiliensis (Diptera: Cecidomyiidae) in cassava (Euphorbiaceae: Manihot esculenta Crantz) as our model. PCR-based molecular markers and deep metagenomic sequencing data were employed to analyse the gall microbiome and to test the hypothesis that gall cells are genetically transformed by insect vectored bacteria. A shotgun sequencing discrimination approach was implemented to selectively discriminate between foreign DNA and the reference host plant genome. Several known candidate insertion sequences were identified, the most significant being DNA sequences found in bacterial genes related to the transcription regulatory factor CadR, cadmium-transporting ATPase encoded by the cadA gene, nitrate transport permease protein (nrtB gene), and arsenical pump ATPase (arsA gene). In addition, a DNA fragment associated with ubiquitin-like gene E2 was identified as a potential accessory genetic element involved in gall induction mechanism. Furthermore, our results suggest that the increased quality and rapid development of gall tissue are mostly driven by microbiome enrichment and the acquisition of critical endophytes. An initial gall-like structure was experimentally obtained in M. esculenta cultured tissues through inoculation assays using a Rhodococcus bacterial strain that originated from the inducing insect, which we related to the gall induction process. We provide evidence that the modification of the endophytic microbiome and the genetic transformation of plant cells in M. esculenta are two essential requirements for insect-induced gall formation. Based on these findings and having observed the same potential DNA marker in galls from other plant species (ubiquitin-like gene E2), we speculate that bacterially mediated genetic transformation of plant cells may represent a more widespread gall induction mechanism found in nature.

Detection of SNPs in FABP4 Gene and Its Relationship with Milk Quality Traits in Iraqi Jenoubi Cows

Jenoubi cattle breed classified as a Bos indicus L., is one of Iraqi indigenous cattle that kept for milk production purpose, but the influence of fatty acid-binding protein 4 (FABP4) on milk-related traits remain poorly understood. The current study was conducted to investigate the effect of single nucleotide polymorphisms (SNPs) in the FABP4 gene of Jenoubi cows. The DNA samples were extracted from the blood of 21 cows. Along 565-bp nucleotide (nt) amplicon from the intron 2- exon 3- intron 3 of FABP4 was amplified by polymerase chain reaction (PCR) while forward sequencing was used to detect SNPs and results were analyzed using bioinformatics software. The findings resulted the existence of two SNPs (c.3689G&gt;A and c.3709G&gt;C) in the exon 3 and four SNPs (g.3494A&gt;C, g.3531A&gt;T, g.3743T&gt;C and g.3765T&gt;C) in the intron regions. Therefore, three haplotypes of FABP4 gene were obtained in this study based on three SNPs of g.3494A&gt;C, g.3531A&gt;T and g.3765T&gt;C. Nonetheless, those haplotypes were not significantly influencing to the milk quality traits of Jenoubi cows. However, the SNPs examined in this study might not be used as potential DNA marker to improve milk production traits of the current breed.

Chloroplast genome characteristics and phylogeny of the sinodielsia clade (apiaceae: apioideae)

BackgroundThe Sinodielsia clade of the subfamily Apioideae (Apiacieae) was established in 2008, and it is composed of 37 species from 17 genera. Its circumscription is still poorly delimited and unstable, and interspecific relationships in the clade lack comprehensive analysis. Chloroplast (cp.) genomes provide valuable and informative data sources for evolutionary biology and have been widely used in studies on plant phylogeny. To infer the phylogenetic history of the Sinodielsia clade, we assembled complete cp. genomes of 39 species and then performed phylogenetic analysis based on these cp. genome sequence data combined with 66 published cp. genomes from 16 genera relative to the Sinodielsia clade.ResultsThese 39 newly assembled genomes had a typical quadripartite structure with two inverted repeat regions (IRs: 17,599–31,486 bp) separated by a large single-copy region (LSC: 82,048–94,046 bp) and a small single-copy region (SSC: 16,343–17,917 bp). The phylogenetic analysis showed that 19 species were clustered into the Sinodielsia clade, and they were divided into two subclades. Six mutation hotspot regions were detected from the whole cp. genomes among the Sinodielsia clade, namely, rbcL–accD, ycf4–cemA, petA–psbJ, ycf1–ndhF, ndhF–rpl32 and ycf1, and it was found that ndhF–rpl32 and ycf1 were highly variable in the 105 sampled cp. genomes.ConclusionThe Sinodielsia clade was subdivided into two subclades relevant to geographical distributions, except for cultivated and introduced species. Six mutation hotspot regions, especially ndhF–rpl32 and ycf1, could be used as potential DNA markers in the identification and phylogenetic analyses of the Sinodielsia clade and Apioideae. Our study provided new insights into the phylogeny of the Sinodielsia clade and valuable information on cp. genome evolution in Apioideae.

CNN_FunBar: Advanced Learning Technique for Fungi ITS Region Classification.

Fungal species identification from metagenomic data is a highly challenging task. Internal Transcribed Spacer (ITS) region is a potential DNA marker for fungi taxonomy prediction. Computational approaches, especially deep learning algorithms, are highly efficient for better pattern recognition and classification of large datasets compared to in silico techniques such as BLAST and machine learning methods. Here in this study, we present CNN_FunBar, a convolutional neural network-based approach for the classification of fungi ITS sequences from UNITE+INSDC reference datasets. Effects of convolution kernel size, filter numbers, k-mer size, degree of diversity and category-wise frequency of ITS sequences on classification performances of CNN models have been assessed at all taxonomic levels (species, genus, family, order, class and phylum). It is observed that CNN models can produce >93% average accuracy for classifying ITS sequences from balanced datasets with 500 sequences per category and 6-mer frequency features at all levels. The comparative study has revealed that CNN_FunBar can outperform machine learning-based algorithms (SVM, KNN, Naïve-Bayes and Random Forest) as well as existing fungal taxonomy prediction software (funbarRF, Mothur, RDP Classifier and SINTAX). The present study will be helpful for fungal taxonomy classification using large metagenomic datasets.

Microsatellite Analysis Revealed Potential DNA Markers for Gestation Length and Sub-Population Diversity in Kari Sheep

Kari sheep inhabiting the Chitral district of Pakistan show variation in gestation length. In this study, we have analyzed the genetic differences between the three subtypes of Kari sheep (based on variation in gestation length) using microsatellite markers. Kari sheep samples were collected from their breeding tract and were characterized for gestation length and genetic diversity using microsatellite markers. A total of 78 Kari ewes were grouped into three categories based on gestation length (GL), i.e., Kari-S (with a shorter GL), Kari-M (with a medium GL), and Kari-L (with a longer GL). DNA from these samples was used to amplify 31 ovine-specific microsatellite loci through PCR. Of the total 78 Kari specimens, 24 were grouped in Kari-S (GL = 100.7 ± 1.8), 26 were from the Kari-M subtype (GL = 123.1 ± 1.0), and 28 were Kari-L (GL = 143.8 ± 1.5). Microsatellite analysis revealed an association of genotypes at two marker sites (MAF214 and ILSTS5) with variation in GL. A total of 158 alleles were detected across the 22 polymorphic loci with an average of 7.18 alleles per locus. Unique alleles were found in all three subtypes. The highest number of unique alleles was observed in Kari-L (15), followed by Kari-S (10) and Kari-M (8). The results indicated that Kari-S is a genetically distinct subtype (with higher genetic differentiation and distance) from Kari-M and Kari-L. The genetic uniqueness of Kari-S is important for further exploration of the genetic basis for shorter gestation length, and exploitation of their unique values.

Detection of small sequence variations within the goat GHR gene and its effects on growth traits

Growth hormone receptor (GHR) gene is considered to be an important candidate gene in growth traits. Therefore, the purpose of this study was to detect whether there were potential indel variations in the GHR gene that were related to the growth traits of the Shaanbei white cashmere goats (SBWC). In this study, genomic DNA from 931 healthy SBWC individuals were used to verify the relationship between the indel of the GHR gene and growth traits. Two indel variants, P49-bp indel in intron 1 and P1410-bp indel in 3’-UTR, were confirmed. Association analyses demonstrated that these two indel polymorphism loci were associated with the chest circumference and chest width of SBWC. Additionally, for the P49-bp and P1410-bp indel loci, the ID and II genotypes were dominant genotypes, respectively. Moreover, the genotypic distributions of these two indel loci in SBWC were significantly different from those in three other Chinese indigenous goat breeds (HNBG, GZDG and IMWC) (p < 0.05). Taken together, two indel loci (P49-bp indel and P1410-bp indel) both significantly affected the growth traits of goats. This illustrated that these two indel loci might be the potential DNA marker for use in improving the selection and breeding of goats.

The complete chloroplast genome of two Firmiana species and comparative analysis with other related species.

Firmiana is a small genus within the subfamily Sterculioideae of the Malvaceae. There are nine Firmiana species distributed in South and South-west China, most of which are endangered. Due to the shortage of plastid genomes data, the phylogenetic relationships and the evolutionary history of this genus remain unclear. Therefore, the complete chloroplast genomes of F. calcarean and F. hainanensis were sequenced using high-throughput sequencing and then compared with the chloroplast genomes of other reported Firmiana species. The genome size of F. calcarean and F. hainanensis is 161,263 and 160,031bp long, respectively, containing a total of 131 genes (including 85 protein coding genes, 37 tRNAs, 8 rRNAs, and one pseudogene). Comparative analysis revealed that the genome structure, GC content, gene content and order, as well as the RNA editing sites within the chloroplast genomes of F. calcarean and F. hainanensis were similar to previously reported Firmiana species. ML phylogenetic analysis revealed that F. danxiaensis, F. hainanensis, F. calcarean, F. simplex, and F. major form a sister group to F. colorata, F. pulcherrima, and F. kwangsiensis. The SSRs, long repeats, and 21 highly divergent regions (Pi > 0.01) identified in this study might provide potential DNA markers for further population genetics and phylogenetic studies of Firmiana. Our findings can help design new species-specific molecular markers and the general framework to further explore the evolutionary history of Firmiana and to address their conservation challenges.

Differentially expressed genes against Colletotrichum lindemuthiamum in a bean genotype carrying the Co-2 gene revealed by RNA-sequencing analysis.

Anthracnose is responsible for large yield losses in common bean crops. RNA-sequencing was used to investigate the differentially expressed genes (DEGs) in response to race 38 of Colletotrichum lindemuthianum in two near-isogenic lines (A25 and A4804) that differ in the presence of a resistance gene located in the cluster Co-2. Their responses were analyzed at different hours after inoculation (0, 24, and 48) and within and between genotypes. In all, 2,850 DEGs were detected, with 2,373 assigned to at least one functional GO term. Enriched GO terms in the resistant genotype were mainly related to functions as a response to stimulus, hormone signaling, cellular component organization, phosphorylation activities, and transcriptional regulation. The region containing the Co-2 cluster was delimited at the end of chromosome Pv11 (46.65–48.65 Mb) through a comparison with the SNP genotypes, obtained using ‘Genotyping by Sequencing,’ among seven resistant lines harboring the Co-2 gene and the susceptible line A25. The delimited region contained 23 DEGs, including 8 typical R genes, that showed higher expression levels in the resistant genotype and non-changes in the susceptible genotype after inoculation. Six R genes encoding protein kinases and an LRR domain formed a cluster in a core region between 46.98 and 47.04 Mb. The alignment of the raw transcriptome reads in the core region revealed structural changes that were used to design four potential breeder-friendly DNA markers, and it revealed some alignments with the intergenic regions, suggesting the presence of genes in addition to those annotated in the reference genome.

Genetic polymorphisms within the ETAA1 gene associated with growth traits in Chinese sheep breeds.

Ewing tumor-associated antigen 1 (ETAA1) is an important candidate gene for fat deposition in sheep. This study aimed to investigate the indel variations in the ETAA1gene and their associations with growth traits in the 1619sheep. The polymorphism information content of this indel ranged from 0.308 to 0.375, with medium genetic diversity (0.25 ≤ polymorphism information content ≤ 0.5). The correlation analysis showed that an 8-bp insertion within the ETAA1gene was significantly associated with growth traits in Luxi Blackhead sheep (LXBH), Lanzhou fat-tailed sheep, Hu sheep, Tong sheep, and Tan sheep (p<0.05). Furthermore, the ETAA1gene was significantly associated with several growth traits (p<0.01), such as chest width and paunch girth of LXBH, Tong, and Tan rams. It was significantly related to the body morphometric traits of LXBH, Lanzhou fat-tailed sheep, Hu, Tong, and Tan ewes. In conclusion, the detected 8-bp indel within the ETAA1gene was confirmed in sheep, significantly affecting the growth traits, and might be a potential DNA marker for the selection of high-quality individuals in marker-assisted selection for sheep breeding.

Two strongly linked blocks within the KIF16B gene significantly influence wool length and greasy yield in fine wool sheep (Ovis aries)

A previous genome-wide association study (GWAS) identified the kinesin family member 16B (KIF16B) as a candidate gene related to sheep wool production. In this work, DNA pool sequencing and SNPscanTM high-throughput genotyping methods were used to detect single-nucleotide polymorphisms (SNPs) in the sheep KIF16B gene. The correlations between the SNPs and wool length and greasy wool yield were systematically assessed. Forty-five SNPs were identified and 37 of them were genotyped, including 10 exon mutations, 26 intron mutations, and 1 promoter region mutation. Most of the SNPs were of medium genetic diversity and at Hardy-Weinberg equilibrium (HWE). Among them, 10 SNPs were associated with greasy wool yield and 28 SNPs impact the wool length. Five specific SNPs were found to exert significant effects on the wool length in all body parts analyzed in this study. Furthermore, linkage disequilibrium (LD) analysis was conducted among SNP loci and they were found to be significantly associated with economically important traits. Two strongly linked SNP blocks were identified within these SNPs and they might exert significant impacts on the greasy wool yield and wool length. The identified SNPs exert significant effects on wool production and could be considered as potential DNA markers for selecting the individuals with superior phenotypes. How to cite: Zhao H, Hu R, Li F, et al. Two strongly linked blocks within the KIF16B gene significantly influence wool length and greasy yield in fine wool sheep (Ovis aries). Electron J Biotechnol 2021:53. https://doi.org/10.1016/j.ejbt.2021.05.003

Association between single nucleotide polymorphisms of nLvALF1 and PEN2-1 genes and resistance to Vibrio parahaemolyticus in the Pacific white shrimp Litopenaeusvannamei

Vibrio parahaemolyticus carrying a plasmid encoding the PirAvp and PirBvp toxins is the causative agent of acute hepatopancreatic necrosis disease (AHPND) in penaeid shrimps. In the Pacific white shrimp farming industry, one possible strategy to reduce economic loss due to AHPND is the development of a shrimp line resistant to the disease. In this study, we identified single nucleotide polymorphisms (SNPs) in the Litopenaeus vannamei anti-lipopolysaccharide factor 1 (nLvALF1) and penaeidin 2-1 (PEN2-1) genes, and we analyzed the associations between these SNPs and resistance/susceptibility to V. parahaemolyticus infection in the Pacific white shrimp. Postlarvae (PL20) shrimp from a local hatchery in Prachuap Khiri Khan Province were challenged with an isolate of VPAHPND and mortality was observed for 14 days. DNA was extracted from susceptible (died within 6 days) and resistant (survived the challenge) shrimp (45 individuals/group) and used for PCR amplifications of nLvALF1 (397 bp) and PEN2-1 (637 bp) gene fragments. PCR products were sequenced by direct sequencing and SNPs were identified from sequencing chromatograms. Nine and seven SNPs were identified in nLvALF1 and PEN2-1 gene fragments, respectively. Analyses of allele frequencies in susceptible and resistant samples using Chi-square tests revealed that four and six SNPs in nLvALF1 and PEN2-1, respectively, were associated with resistance/susceptibility to V. parahaemolyticus infection (P < 0.05). The SNPs in the candidate genes identified here are potential DNA markers for breeding V. parahaemolyticus-resistant Pacific white shrimp in the study population; however, further validation will be required if these SNPs are to be used across populations.