High-quality breeding group selection method, device and equipment and storage medium

Abstract

The invention relates to the technical field of breeding, and discloses a high-quality breeding group selection method, device and equipment and a storage medium. The method comprises the following steps: uploading genome re-sequencing data to a cloud computing platform, and fully utilizing high-performance computing resources and a distributed storage technology of cloud computing to realize automatic processing and efficient analysis of genome data; a dynamic weighting comparison algorithm based on a graph theory is adopted to perform segmentation comparison and weight optimization on the uploaded genome data to generate a genome sequence map, so that the comparison efficiency and accuracy are greatly improved; in addition, through a genome association analysis technology and a deep learning model, breeding individuals with excellent characters can be accurately identified and screened out; finally, the individuals with the highest comprehensive score are screened out through comprehensive scores, and a high-quality breeding group is formed. According to the method, the accuracy of individual selection is improved, manual intervention is remarkably reduced, and the breeding process is more efficient and reliable.

Description

Technical Field

[0001] The present invention relates to the field of breeding technology, and in particular to a method, device, equipment and storage medium for selecting a high-quality breeding population. Background Art

[0002] The rapid development of genome resequencing technology has provided unprecedented opportunities for agricultural breeding. Although traditional breeding methods, such as phenotypic selection and hybrid breeding, have been successful to a certain extent, they still have many problems such as low efficiency, long cycle, and uncertain results. Through genome resequencing technology, scientists can fully understand the information of crop genomes and identify gene loci associated with excellent traits, thereby providing a scientific basis for breeding decisions. However, with the rapid accumulation of genomic data, how to efficiently process and analyze these massive data has become a huge challenge.

[0003] In traditional breeding methods, the processing and analysis of genomic data usually rely on local computing resources, which not only requires high-performance hardware support, but also requires a lot of time for data processing and analysis. For example, the genome resequencing data of a single breeding individual is usually as high as hundreds of GB. Faced with the data processing tasks of hundreds of breeding individuals, traditional methods often seem to be unable to cope with it. In addition, when processing genomic data, the existing technology often requires manual steps such as quality control, sequence splicing, and genetic variation detection. These steps are cumbersome and error-prone, further limiting the breeding efficiency. Another defect of the existing technology is that it is difficult to accurately identify breeding individuals with excellent traits while processing large-scale data. Traditional genomic association analysis methods mostly rely on linear models or simple statistical methods. These methods are often unable to cope with complex traits and it is difficult to effectively identify gene loci that are significantly associated with the target traits. In addition, traditional methods usually ignore the effects of multi-gene interactions and environmental factors on traits, resulting in low accuracy and reliability of screening results.

[0004] The rise of cloud computing has provided a new way to solve these problems. By uploading genome resequencing data to the cloud platform, we can make full use of the high-performance computing resources and distributed storage technology of cloud computing to quickly complete the processing and analysis of large-scale data. However, there are still many challenges in realizing the automated processing and efficient analysis of genomic data on the cloud computing platform. For example, how to achieve efficient alignment and variation detection of genomic data on the cloud platform, how to use deep learning and artificial intelligence technology for genome association analysis, and how to accurately screen individuals with excellent traits based on the results of association analysis are all hot spots and difficulties in current research.

[0005] Traditional genome alignment algorithms are inefficient when processing large-scale data and are difficult to meet the needs of practical applications. Existing alignment algorithms mostly use fixed weights or simple dynamic programming methods, which make it difficult to fully utilize the parallel computing capabilities of cloud computing. In addition, traditional genetic variation detection methods usually require manual adjustment of parameters and thresholds, which are difficult to adapt to the characteristics of different data sets, affecting the accuracy of the test results. In terms of genome association analysis, existing technologies mainly rely on single statistical methods or simple machine learning models, which perform poorly when dealing with complex traits and multi-gene interactions. Although deep learning technology has made significant progress in fields such as image recognition and natural language processing, its application in genome association analysis is still in the exploratory stage. How to design an effective deep learning model and combine genomic data and phenotypic data for association analysis is still an urgent problem to be solved.

[0006] In summary, existing technologies have many shortcomings in genomic data processing, genetic variation detection and genomic association analysis, and it is difficult to meet the efficiency and precision requirements of modern breeding.

[0007] Therefore, the prior art still needs to be improved and developed. Summary of the invention

[0008] In order to overcome the shortcomings of the prior art, the purpose of the present invention is to provide a method, device, equipment and storage medium for selecting high-quality breeding populations, aiming to solve the many shortcomings of the prior art in genome data processing, genetic variation detection and genome association analysis, and the problem that it is difficult to meet the requirements of high efficiency and accuracy of modern breeding.

[0009] The first aspect of the present invention provides a method for selecting a high-quality breeding population, comprising: performing multiple high-throughput sequencing on each breeding individual in the breeding population to obtain genome sequence data; performing preliminary processing on the genome sequence data to obtain a pre-processed genome sequence; uploading the pre-processed genome sequence to a cloud platform, using a dynamic weighted alignment algorithm based on graph theory to compare and analyze the uploaded genome sequence data to generate an individual genome sequence map; performing genetic variation detection on the genome sequence map to identify genetic variation data, wherein the genetic variation data includes single nucleotide polymorphisms, insertions and deletions, and structural variations; performing genome association analysis based on the genetic variation data to identify key gene sites related to target traits; using a multi-layer convolutional neural network based on deep learning, combined with an attention mechanism, to construct an association model for predicting the association between target traits and key gene sites; and screening out breeding individuals with excellent traits based on the results of the genome association analysis and the output value of the association model.

[0010] Optionally, in a first implementation of the first aspect of the present invention, the collected genome sequence data Perform preliminary quality control, use quality control software to evaluate the quality of each sequence, remove low-quality sequence data with a quality score lower than Q20, and obtain a high-quality sequence collection , where Q(d ij ) represents the sequence d ij The quality score of i QC It is a set of high-quality sequences retained after quality control, among which, , d in represents the nth DNA sequencing of the i-th individual sample; the software is used to remove the adapter sequence of the high-quality sequence set to generate trimmed sequence data , where m≤n; use the splicing software to assemble the trimmed sequence data Splice to generate the preliminary genome splicing sequence G for each individual i ; For the initially assembled genome sequence G i Remove repetitive sequences, identify and remove redundant and repetitive sequence fragments, and generate pre-processed genome sequences after removing repetitive sequences , where p≤k.

[0011] Optionally, in a second implementation of the first aspect of the present invention, the pre-processed genome sequence is uploaded to a cloud platform, and a dynamic weighted alignment algorithm based on graph theory is used to compare and analyze the uploaded genome sequence data to generate an individual genome sequence map, including the steps of: constructing a sequence fragment graph structure G(V, E) based on graph theory on the cloud platform, where V represents a node set, and each node corresponds to a sequence fragment s ijk , E represents the edge set, each edge corresponds to the comparison relationship between two sequence fragments, and the nodes and edges are constructed as follows: ; ; where N represents the number of samples, m i represents the number of fragments of sample i, v ij represents the fragment j of sample i, w ijk represents the comparison weight between fragments j and k; design an alignment algorithm that combines breeding goals, where each sequence fragment is divided into several sub-segments ijk , each sub-fragment corresponds to a child node v in the graph structure ijk : ,in, Representation fragment The sub-segment between positions a and b, Representation fragment The length of the breeding target is determined by introducing a high-quality trait scoring function specific to the breeding target, and determining the weights w between the child nodes. ijkl, and construct the initial weight matrix W, where w ijkl Represents the child node v ijk and v ijl The quality trait association weights between ; ; ,in, and are the balance coefficients of similarity weight and quality trait weight, is the Kronecker delta function, is the scoring function for high-quality traits, is the weight of gene locus y, Represents sub-fragment s ijk and ijl The specific breeding target similarity at gene locus y, T is the gene locus set associated with high-quality traits; introduce the characteristic information of high-quality breeding individuals and specific breeding targets, dynamically adjust the weight matrix W, and obtain the adjusted weight matrix W', where ; ;in, represents the adjusted weight, is the weight adjustment coefficient, Represents sub-fragment s ijk and ijl The matching function between the gene loci T associated with a specific breeding target; using the optimized dynamic programming algorithm, the optimal alignment path is solved for the adjusted weight matrix W' to generate the individual genome sequence map M i , , where π represents the comparison path, (u, v) represents the edge in the path, Represents the adjusted weight.

[0012] Optionally, in a third implementation of the first aspect of the present invention, genetic variation detection is performed on the genome sequence map to identify genetic variation data, wherein the genetic variation data includes single nucleotide polymorphisms, insertions and deletions, and structural variations, including the steps of: using a cloud platform tool to perform genetic variation detection on the genome sequence map to obtain genetic variation data, wherein the detection method for single nucleotide polymorphisms is: , where p represents the genomic position, b represents the base type, q represents the quality score of the variant, and Q min is the quality score threshold; the detection method for insertion and deletion is: ,in, Indicates the length of the insertion or deletion; the detection method for structural variation is: , where t represents the mutation type, and DUP, DEL, INV, and TRA represent duplication, deletion, inversion, and translocation, respectively.

[0013] Optionally, in a fourth implementation of the first aspect of the present invention, performing genome association analysis based on the genetic variation data to identify key gene loci associated with the target trait comprises the steps of: Annotate detected single nucleotide polymorphisms, insertions and deletions, and structural variations and integrate them into a comprehensive annotation result A i , including the location, type, frequency, functional impact and related annotation information of each variant; comprehensive annotation results A i Perform statistical analysis and generate statistical reports of genetic variation, including the number, distribution, and functional classification of variation: ,in, represents the total number of single nucleotide polymorphisms detected in the genome of individual i, represents the sum of the number of single nucleotide polymorphisms detected in each sample, m i is the number of sequenced fragments in individual i, SNP ij represents the set of SNPs detected in fragment j, represents the size of the set, i.e. the number of SNPs in fragment j, represents the total number of indels detected in the genome of individual i, It means to add up the number of insertions and deletions detected in each sample. represents the set of indels detected in fragment j, represents the size of the set, i.e. the number of indels in fragment j, represents the total number of structural variations detected in the genome of individual i, Indicates the sum of the number of structural variations detected in each sample, SV i represents the set of structural variations detected in fragment j, represents the functional classification of all genetic variants detected in the genome of individual i, It means to classify all detected genetic variations according to functional classification; to associate the statistical report and annotation results of genetic variations with breeding goals, and to identify key gene loci related to target traits: ,in, Indicates whether the variation v is related to the breeding target T.

[0014] Optionally, in a fifth implementation of the first aspect of the present invention, a multi-layer convolutional neural network based on deep learning is used in combination with an attention mechanism to construct an association model for predicting the association between the target trait and the key gene loci, comprising the steps of: based on the detected comprehensive annotation result A i, perform preliminary screening on each variant and select a set of potential candidate gene loci C associated with the target trait i , , where v represents the variation, T is the breeding target, and P is the characteristic vector of the variation. is a function that evaluates the correlation between the variation v and the breeding target T, is a set threshold; a multi-layer convolutional neural network based on deep learning is used to further analyze the candidate gene loci, and an initial association model for predicting the association between the target trait and the key gene loci is constructed, wherein the multi-layer convolutional neural network structure comprises multiple convolutional layers, pooling layers and fully connected layers, and an attention mechanism is introduced into the convolutional neural network to assign different weights to the characteristics of different gene loci, thereby generating an attention weight matrix A; the predicted value of the association between the key gene loci and the target trait is output through the fully connected layer; the classification loss function is used to train the initial association model, optimize the network parameters, and obtain the association model.

[0015] Optionally, in a sixth implementation of the first aspect of the present invention, breeding individuals with excellent traits are screened out according to the results of the genome association analysis and the output value of the association model, comprising the steps of: obtaining a set of key gene loci for each breeding individual based on the results of the genome association analysis; , where v represents the variation, C i is the set of candidate gene loci, T k is the kth target trait, Represents the variation v and the target trait T k The correlation, p k Represents the target trait T k The weight of The comprehensive score of each breeding individual is calculated based on the weight of each key gene locus and the predicted value output by the association model. , where w v represents the weight of the key gene locus v, Represents the output value of the association model for the gene locus v, Indicates the similarity between breeding individuals i and j; according to the comprehensive score S i Sort the breeding individuals, select the individuals with the highest comprehensive scores from high to low, and build a high-quality breeding population: ,in, is the ratio of the set comprehensive score threshold, and N represents the total number of breeding individuals.

[0016] The second aspect of the present invention provides a high-quality breeding population selection device, comprising: a sequencing module, which is used to perform multiple high-throughput sequencing on each breeding individual in the breeding population to obtain genome sequence data; a preprocessing module, which is used to perform preliminary processing on the genome sequence data to obtain a preprocessed genome sequence; a comparison and analysis module, which is used to upload the preprocessed genome sequence to a cloud platform, and use a dynamic weighted comparison algorithm based on graph theory to compare and analyze the uploaded genome sequence data to generate an individual genome sequence map; a variation detection module, which is used to perform genetic variation detection on the genome sequence map to identify genetic variation data, wherein the genetic variation data includes single nucleotide polymorphisms, insertions and deletions, and structural variations; an association analysis module, which is used to perform genome association analysis based on the genetic variation data to identify key gene sites related to target traits; a model construction module, which is used to use a multi-layer convolutional neural network based on deep learning, combined with an attention mechanism, to construct an association model for predicting the association between target traits and key gene sites; a screening module, which is used to screen out breeding individuals with excellent traits based on the results of the genome association analysis and the output value of the association model.

[0017] The third aspect of the present invention provides a high-quality breeding population selection device, comprising: a memory and at least one processor, wherein the memory stores computer-readable instructions, and the memory and the at least one processor are interconnected through lines; the at least one processor calls the computer-readable instructions in the memory to enable the high-quality breeding population selection device to perform the various steps of the high-quality breeding population selection method as described above.

[0018] A fourth aspect of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores computer-readable instructions, which, when executed on a computer, enable the computer to execute the various steps of the method for selecting a high-quality breeding population as described above.

[0019] Beneficial effects: The present invention uploads genome resequencing data to a cloud computing platform, making full use of the high-performance computing resources and distributed storage technology of cloud computing to achieve automated processing and efficient analysis of genome data. Through distributed storage technology, a large amount of genome data can be efficiently uploaded, stored and backed up to ensure data security and availability. A dynamic weighted alignment algorithm based on graph theory is used to segment and optimize the weight of the uploaded genome data to generate a genome sequence map, which greatly improves the alignment efficiency and accuracy. In addition, through genome association analysis technology and deep learning models, breeding individuals with excellent traits can be accurately identified and screened out. A multi-layer convolutional neural network structure is used to extract features through multiple convolutional layers, pooling layers and fully connected layers, and an attention mechanism is introduced to assign different weights to the features of different gene loci, thereby improving the accuracy of association analysis. Finally, through comprehensive scoring and further verification, the individuals with the highest comprehensive scores are screened out to form a high-quality breeding population. This process not only improves the accuracy of individual selection, but also significantly reduces manual intervention, making the breeding process more efficient and reliable. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] Figure 1 A flow chart of a method for selecting a high-quality breeding population provided in an embodiment of the present invention.

[0021] Figure 2 This is a schematic diagram of the structure of the high-quality breeding population selection device provided by the present invention.

[0022] Figure 3 This is a schematic diagram of the structure of the high-quality breeding population selection equipment provided by the present invention. DETAILED DESCRIPTION

[0023] The embodiments of the present invention provide a method, device, equipment and storage medium for selecting a high-quality breeding population. The terms "first", "second", "third", "fourth", etc. (if any) in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments described here can be implemented in an order other than that illustrated or described here. In addition, the terms "including" or "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0024] See also Figure 1 , Figure 1A flow chart of a high-quality breeding population selection method provided by the present invention, as shown in the figure, comprises the steps of: S10, performing high-throughput sequencing multiple times on each breeding individual in the breeding population to obtain genome sequence data; Specifically, this embodiment collects samples from each breeding individual in the breeding population. Individual samples include plant tissues of leaves, stems and roots. Ensure that each sample is large enough and healthy to avoid collecting damaged or diseased parts. Each sample is uniquely marked and relevant information (such as individual number, collection date, location, etc.) is recorded. The samples are stored in an appropriate storage medium (such as liquid nitrogen or a refrigerator at a specific temperature) to maintain the integrity of their genetic information. A sequencing instrument is used to extract DNA from the collected individual samples; the extracted DNA is sequenced multiple times using sequencing technology, and each individual sample is sequenced at least three times independently to generate genome sequence data. ,In this embodiment, each individual sample is sequenced at least three times independently to reduce the impact of sequencing errors and biases; In the sequencing process of this embodiment, the Illumina platform is used to generate short read length or long read length sequence data, where the read length is defined as the length of a single sequence read generated by the sequencing instrument, which is between 100bp and 300bp. The Illumina platform has the advantages of high throughput, high accuracy and low cost, and is suitable for large-scale genome sequencing. By setting appropriate sequencing parameters, short read length or long read length sequence data that meets research needs can be generated.

[0025] S20, performing preliminary processing on the genome sequence data to obtain a preprocessed genome sequence; In this embodiment, the collected genome data D i Perform preliminary quality control, use quality control software to evaluate the quality of each sequence, remove low-quality sequence data with a quality score lower than Q20, and obtain a high-quality sequence collection , where Q(d ij ) represents the sequence d ij The quality score of i QC It is a set of high-quality sequences retained after quality control, among which, , d in represents the nth DNA sequencing of the i-th individual sample; preliminary quality control can remove low-quality sequence data and improve the accuracy and reliability of subsequent analysis; setting the quality threshold Q20 can ensure that the retained sequence data has a high quality level.

[0026] Then use the software to remove the adapter sequence from the high-quality sequence collection to generate trimmed sequence data , where m≤n, removing the adapter sequence can reduce interference and errors in subsequent analysis; using splicing software to assemble the trimmed sequence data Splice to generate the preliminary genome splicing sequence G for each individual i , the spliced ​​sequence data can splice multiple short sequences into long sequences, providing complete genome sequence information for subsequent genome analysis; finally, the preliminary spliced ​​genome sequence Remove repetitive sequences, identify and remove redundant and repetitive sequence fragments, and generate genome sequences after removing repetitive sequences , where p≤k, removing duplicate sequences can reduce the complexity and computational effort of subsequent analysis.

[0027] S30, uploading the pre-processed genome sequence to a cloud platform, using a dynamic weighted alignment algorithm based on graph theory to compare and analyze the uploaded genome sequence data, and generating an individual genome sequence map; In this embodiment, the pre-processed genome sequence is first uploaded to the cloud platform, and then a sequence fragment graph structure G(V, E) based on graph theory is constructed on the cloud platform, where V represents a node set, and each node corresponds to a sequence fragment s ijk , E represents the edge set, each edge corresponds to the comparison relationship between two sequence fragments, and the nodes and edges are constructed as follows: ; ; where N represents the number of samples, m i represents the number of fragments of sample i, v ij represents the fragment j of sample i, w ijk represents the comparison weight between fragments j and k; Design an alignment algorithm that incorporates breeding objectives, where each sequence fragment is divided into several sub-segments ijk , each sub-fragment corresponds to a child node v in the graph structure ijk : ,in, Representation fragment The sub-segment between positions a and b, Representation fragment Length; By introducing a high-quality trait scoring function specific to the breeding goal, the weights w between the child nodes are determined. ijkl , and construct the initial weight matrix W, where w ijkl Represents the child node v ijk and v ijl The quality trait association weights between ; ; ,in, and are the balance coefficients of similarity weight and quality trait weight, is the Kronecker delta function, is the scoring function for high-quality traits, is the weight of gene locus y, Represents sub-fragment s ijk and ijl The similarity of a specific breeding target at the gene locus y, T is the set of gene loci associated with high-quality traits; Introducing the characteristic information of high-quality breeding individuals and specific breeding goals, dynamically adjusting the weight matrix W, and obtaining the adjusted weight matrix W', where ; ;in, represents the adjusted weight, is the weight adjustment coefficient, Represents sub-fragment s ijk and ijl The matching function between the gene loci T associated with a specific breeding goal; The optimized dynamic programming algorithm is used to solve the optimal alignment path for the adjusted weight matrix W' to generate the individual genome sequence map M i , , where π represents the comparison path, (u, v) represents the edge in the path, Indicates the adjusted weight. i Finally, the map is checked and corrected, and the accuracy and stability of high-quality traits for specific breeding targets are verified.

[0028] This embodiment, by combining a dynamic weighted alignment algorithm for specific breeding targets, not only achieves efficient alignment of genome sequences, but also introduces high-quality trait information of specific breeding targets into the alignment process, ensuring that the generated genome sequence map more accurately reflects the specific high-quality traits of the breeding individuals, and providing more accurate data support for subsequent genetic variation detection and genome association analysis.

[0029] S40, performing genetic variation detection on the genome sequence map to identify genetic variation data, wherein the genetic variation data includes single nucleotide polymorphisms, insertions and deletions, and structural variations; In this embodiment, a cloud platform tool is used to perform genetic variation detection on the genome sequence map to obtain genetic variation data, wherein the detection method for single nucleotide polymorphism is: , where p represents the genomic position, b represents the base type, q represents the quality score of the variant, and Q min is the quality score threshold; the detection method for insertion and deletion is: ,in, Indicates the length of the insertion or deletion; The detection method for structural variation is: , where t represents the type of variation, and DUP, DEL, INV, and TRA represent duplication, deletion, inversion, and translocation, respectively. This embodiment uses tools on the cloud platform to detect genetic variation, which can accurately identify the location and type of variation sites, providing a reliable basis for in-depth research. At the same time, the cloud platform also provides visualization tools and rich parameter setting options, which are convenient for users to flexibly adjust and optimize according to research needs.

[0030] S50, performing genome association analysis based on the genetic variation data to identify key gene loci associated with the target trait; In this embodiment, the detected single nucleotide polymorphisms are annotated, and the position, function and potential impact of each single nucleotide polymorphism are annotated with reference to the genome and annotation database, and annotation results are generated; the detected indels are annotated, the position and length of each indel are analyzed, and the potential functional impact of the indel is annotated using the tools of the cloud platform, and annotation results are generated; the detected structural variations are counted and annotated, structural variations including chromosome rearrangement, copy number variation and large indel are identified, and each structural variation is annotated in detail using the tools of the cloud platform, and annotation results are generated; all detected single nucleotide polymorphisms, indels and structural variations are integrated into a comprehensive annotation result A i , including the location, type, frequency, functional impact and related annotation information of each variant; comprehensive annotation results A i Perform statistical analysis and generate statistical reports of genetic variation, including the number, distribution, and functional classification of variation: ,in, represents the total number of single nucleotide polymorphisms detected in the genome of individual i, represents the sum of the number of single nucleotide polymorphisms detected in each sample, m i is the number of sequenced fragments in individual i, SNP ij represents the set of SNPs detected in fragment j, represents the size of the set, i.e. the number of SNPs in fragment j, represents the total number of indels detected in the genome of individual i, It means to add up the number of insertions and deletions detected in each sample. represents the set of indels detected in fragment j, represents the size of the set, i.e. the number of indels in fragment j, represents the total number of structural variations detected in the genome of individual i, Indicates the sum of the number of structural variations detected in each sample, SV i represents the set of structural variations detected in fragment j, represents the functional classification of all genetic variants detected in the genome of individual i, It means to classify all detected genetic variations by functional classification. The statistical report and annotation results of genetic variations are associated with breeding goals to identify key gene loci related to target traits: ,in, Indicates whether the variation v is related to the breeding target T.

[0031] S60. Using a multi-layer convolutional neural network based on deep learning and combined with an attention mechanism, an association model is constructed to predict the association between target traits and key gene loci. In this embodiment, based on the detected comprehensive annotation result A i , perform preliminary screening on each variant and select a set of potential candidate gene loci C associated with the target trait i , , where v represents the variation, T is the breeding target, and P is the characteristic vector of the variation. is a function that evaluates the correlation between the variation v and the breeding target T, is a set threshold; a multi-layer convolutional neural network based on deep learning is used to further analyze the candidate gene loci, and an initial association model for predicting the association between the target trait and the key gene loci is constructed, wherein the multi-layer convolutional neural network structure comprises multiple convolutional layers, pooling layers and fully connected layers, and an attention mechanism is introduced into the convolutional neural network to assign different weights to the characteristics of different gene loci, thereby generating an attention weight matrix A; the predicted value of the association between the key gene loci and the target trait is output through the fully connected layer; the classification loss function is used to train the initial association model, optimize the network parameters, and obtain the association model.

[0032] Specifically, the main purpose of this embodiment is to analyze the specific degree of association between the candidate gene loci and the target trait, and to construct a high-precision association model between the trait and the gene variation. Among them, the convolution layer in deep learning automatically extracts the key features of the candidate gene loci (such as gene variation patterns, functional associations, etc.), and each candidate gene loci is modeled as a feature vector P, which can be a property of genetic variation (such as position, variation frequency, annotation function, etc.). The potential relationship between the candidate gene loci and the target trait is extracted using the convolution operation, especially the association analysis of the candidate loci with the specific attributes of the target trait (such as gene function, impact path). The multi-layer convolutional neural network structure includes multiple convolutional layers, pooling layers, and fully connected layers. The convolutional layers are used to extract higher-level features. In the network training, the attention mechanism is added to give different weights to different gene loci, emphasizing those sites that are more related to the target trait; finally, the prediction value is generated through the fully connected layer to quantify the contribution of each candidate gene locus to the target trait; the network is trained using the classification loss function to optimize the parameters of the neural network to ensure that the analysis results have higher accuracy.

[0033] In this embodiment, the feature vector p of the training data includes: genetic variation information (such as variation type, annotation function) and target traits of the candidate gene loci, wherein the target traits refer to traits that need to be optimized or improved in breeding objectives, such as disease resistance, yield, etc., and these trait values ​​are obtained through experiments or field trials; the association information between the candidate gene loci and the target traits is used as a label.

[0034] The final output of the model is the predicted value of the association between the key gene locus and the target trait. The predicted value is not a simple probability or percentage, but a quantitative "association score", which is a scalar (such as a floating point number) that indicates the degree of contribution of a specific gene locus to the target trait.

[0035] S70. Screen out breeding individuals with excellent traits based on the results of genome association analysis and the output values ​​of the association model.

[0036] Based on the results of genome association analysis, obtain the key gene loci set for each breeding individual , where v represents the variation, C i is the set of candidate gene loci, T k is the kth target trait, Represents the variation v and the target trait T k The correlation, p k Represents the target trait T k The weight of is the set threshold; The comprehensive score of each breeding individual is calculated based on the weight of each key gene locus and the predicted value output by the association model , where w v represents the weight of the key gene locus v, Represents the output value of the association model for the gene locus v, represents the similarity between breeding individuals i and j; According to the comprehensive score S i Sort the breeding individuals, select the individuals with the highest comprehensive scores from high to low, and build a high-quality breeding population: ,in, is the ratio of the set comprehensive score threshold, and N represents the total number of breeding individuals.

[0037] In the present invention, the selected high-quality breeding individuals can be further verified, including field trials and gene expression analysis, to confirm their breeding potential and the stability of excellent traits; based on the verification results, the high-quality breeding individuals are finally determined and used for breeding.

[0038] The present invention uploads genome resequencing data to a cloud computing platform, makes full use of the high-performance computing resources and distributed storage technology of cloud computing, realizes the automatic processing and efficient analysis of genome data, and through distributed storage technology, can efficiently upload, store and back up a large amount of genome data to ensure the security and availability of data, adopts a dynamic weighted alignment algorithm based on graph theory, performs segmented alignment and weight optimization on the uploaded genome data to generate a genome sequence map, greatly improving the alignment efficiency and accuracy, in addition, through genome association analysis technology and deep learning model, can accurately identify and screen out breeding individuals with excellent traits, utilize multi-layer convolutional neural network structure, extract features through multiple convolutional layers, pooling layers and fully connected layers, and introduce attention mechanism to give different weights to the features of different gene sites, thereby improving the accuracy of association analysis. Finally, through comprehensive scoring and further verification, the individual with the highest comprehensive score is screened out to form a high-quality breeding population, this process not only improves the accuracy of individual selection, but also significantly reduces manual intervention, making the breeding process more efficient and reliable.

[0039] The above describes the method for selecting a high-quality breeding population in the embodiment of the present invention. The following describes the device for selecting a high-quality breeding population in the embodiment of the present invention. Figure 2 , an embodiment of the high-quality breeding population selection device in the embodiment of the present invention comprises: A sequencing module 10 is used to perform multiple high-throughput sequencing on each breeding individual in the breeding population to obtain genome sequence data; A preprocessing module 20, used to perform preliminary processing on the genome sequence data to obtain a preprocessed genome sequence; The comparison and analysis module 30 is used to upload the pre-processed genome sequence to the cloud platform, and use a dynamic weighted comparison algorithm based on graph theory to compare and analyze the uploaded genome sequence data to generate an individual genome sequence map; A variation detection module 40 is used to perform genetic variation detection on the genome sequence map to identify genetic variation data, wherein the genetic variation data includes single nucleotide polymorphisms, insertions and deletions, and structural variations; An association analysis module 50 is used to perform genome association analysis based on the genetic variation data to identify key gene loci associated with the target trait; A model building module 60 is used to construct an association model for predicting the association between the target trait and the key gene loci using a multi-layer convolutional neural network based on deep learning and in combination with an attention mechanism; The screening module 70 is used to screen out breeding individuals with excellent traits according to the results of the genome association analysis and the output value of the association model.

[0040] Based on the same idea as the method in the above-mentioned embodiment, the device provided in the present application can implement the method in the above-mentioned embodiment. For the convenience of explanation, the structural diagram of the device embodiment only shows the parts related to the embodiment of the present application. Those skilled in the art can understand that the illustrated structure does not constitute a limitation on the device, and may include more or fewer modules than shown in the figure, or a combination of certain modules, or different module arrangements.

[0041] Figure 2 The high-quality breeding population selection device in the embodiment of the present invention is described in detail from the perspective of modular functional entities, and the high-quality breeding population selection device in the embodiment of the present invention is described in detail from the perspective of hardware processing.

[0042] Figure 3 1 is a schematic diagram of the structure of a high-quality breeding population selection device provided by an embodiment of the present invention. The high-quality breeding population selection device 100 may have relatively large differences due to different configurations or performances, and may include one or more processors (central processing units, CPU) 11 (for example, one or more processors) and a memory 12, and one or more storage media 13 (for example, one or more mass storage devices) storing application programs 133 or data 132. Among them, the memory 12 and the storage medium 13 can be short-term storage or permanent storage. The program stored in the storage medium 13 may include one or more modules (not shown in the figure), and each module may include a series of instruction operations in the high-quality breeding population selection device 100. Furthermore, the processor 11 can be configured to communicate with the storage medium 13 to execute a series of instruction operations in the storage medium 13 on the high-quality breeding population selection device 100.

[0043] The high-quality breeding population selection device 100 may also include one or more power supplies 14, one or more wired or wireless network interfaces 15, one or more input and output interfaces 16, and / or one or more operating systems 131, such as Windows Serve, Mac OS X, Unix, Linux, FreeBSD, etc. Those skilled in the art will appreciate that Figure 3 The device structure shown does not constitute a limitation on the high-quality breeding population selection device 100, and may include more or less components than shown in the figure, or combine certain components, or arrange the components differently.

[0044] The present invention also provides a computer-readable storage medium, which may be a non-volatile computer-readable storage medium or a volatile computer-readable storage medium. Instructions are stored in the computer-readable storage medium. When the instructions are executed on a computer, the computer executes the steps of the method for selecting a high-quality breeding population.

[0045] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the above-described system, device, or unit can refer to the corresponding process in the aforementioned method embodiment and will not be repeated here.

[0046] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the whole or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM), random access memory (RAM), disk or optical disk and other media that can store program code.

[0047] As described above, the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features thereof may be replaced by equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for selecting a high-quality breeding population, characterized in that: Includes steps: Performing multiple high-throughput sequencing on each breeding individual in the breeding population to obtain genome sequence data; After preliminary processing of the genome sequence data, a preprocessed genome sequence is obtained; The pre-processed genome sequence is uploaded to a cloud platform, and a dynamic weighted alignment algorithm based on graph theory is used to compare and analyze the uploaded genome sequence data to generate an individual genome sequence map; Performing genetic variation detection on the genome sequence map to identify genetic variation data, wherein the genetic variation data includes single nucleotide polymorphisms, insertions and deletions, and structural variations; Performing genome association analysis based on the genetic variation data to identify key gene loci associated with the target trait; Using a multi-layer convolutional neural network based on deep learning and combined with an attention mechanism, an association model is constructed to predict the association between target traits and key gene loci; Breeding individuals with excellent traits are screened out based on the results of genome association analysis and the output values ​​of the association model.

2. The method for selecting a high-quality breeding population according to claim 1, characterized in that: The genome sequence data is preliminarily processed, comprising the steps of: The collected genome sequence data D i Perform preliminary quality control, use quality control software to evaluate the quality of each sequence, remove low-quality sequence data with a quality score lower than Q20, and obtain a high-quality sequence collection , where Q(d ij ) represents the sequence d ij The quality score of i QC It is a set of high-quality sequences retained after quality control, among which, , d in represents the nth DNA sequencing of the i-th individual sample; The software is used to remove the adapter sequence from the high-quality sequence set to generate trimmed sequence data. , where m≤n; Use stitching software to assemble the trimmed sequence data Splice to generate the preliminary genome splicing sequence G for each individual i ; The initial assembled genome sequence G i Remove repetitive sequences, identify and remove redundant and repetitive sequence fragments, and generate pre-processed genome sequences after removing repetitive sequences , where p≤k.

3. The method for selecting a high-quality breeding population according to claim 2, characterized in that: The pre-processed genome sequence is uploaded to a cloud platform, and a dynamic weighted alignment algorithm based on graph theory is used to compare and analyze the uploaded genome sequence data to generate an individual genome sequence map, including the following steps: A graph structure G(V, E) based on graph theory is constructed on the cloud platform, where V represents a set of nodes and each node corresponds to a sequence fragment s. ijk , E represents the edge set, each edge corresponds to the comparison relationship between two sequence fragments, and the nodes and edges are constructed as follows: ; ; where N represents the number of samples, m i represents the number of fragments of sample i, v ij represents the fragment j of sample i, w ijk represents the comparison weight between fragments j and k; Design an alignment algorithm that incorporates breeding objectives, where each sequence fragment is divided into several sub-segments ijk , each sub-fragment corresponds to a child node v in the graph structure ijk : ,in, Representation fragment The sub-segment between positions a and b, Representation fragment Length; By introducing a high-quality trait scoring function specific to the breeding goal, the weights w between the child nodes are determined. ijkl , and construct the initial weight matrix W, where w ijkl Represents the child node v ijk and v ijl The quality trait association weights between ; ; ,in, and are the balance coefficients of similarity weight and quality trait weight, is the Kronecker delta function, is the scoring function for high-quality traits, is the weight of gene locus y, Represents sub-fragment s ijk and ijl The similarity of a specific breeding target at the gene locus y, T is the set of gene loci associated with high-quality traits; Introducing the characteristic information of high-quality breeding individuals and specific breeding goals, dynamically adjusting the weight matrix W, and obtaining the adjusted weight matrix W', where ; ;in, represents the adjusted weight, is the weight adjustment coefficient, Represents sub-fragment s ijk and ijl The matching function between the gene loci T associated with a specific breeding goal; The optimized dynamic programming algorithm is used to solve the optimal alignment path for the adjusted weight matrix W' to generate the individual genome sequence map M i , , where π represents the comparison path, (u, v) represents the edge in the path, Represents the adjusted weight.

4. The method for selecting a high-quality breeding population according to claim 1, characterized in that: Performing genetic variation detection on the genome sequence map to identify genetic variation data, wherein the genetic variation data includes single nucleotide polymorphisms, insertions and deletions, and structural variations, comprises the steps of: The cloud platform tool is used to perform genetic variation detection on the genome sequence map to obtain genetic variation data, wherein the detection method for single nucleotide polymorphism is: , where p represents the genomic position, b represents the base type, q represents the quality score of the variant, and Q min is the quality score threshold; The detection method for insertion and deletion is: ,in, Indicates the length of the insertion or deletion; The detection method for structural variation is: , where t represents the mutation type, and DUP, DEL, INV, and TRA represent duplication, deletion, inversion, and translocation, respectively.

5. The method for selecting a high-quality breeding population according to claim 1, characterized in that: Performing genome association analysis based on the genetic variation data to identify key gene loci associated with the target trait includes the following steps: Annotate detected single nucleotide polymorphisms, insertions and deletions, and structural variations and integrate them into a comprehensive annotation result A i , including the location, type, frequency, functional impact and related annotation information of each variant; Comprehensive annotation results A i Perform statistical analysis and generate statistical reports of genetic variation, including the number, distribution, and functional classification of variation: ,in, represents the total number of single nucleotide polymorphisms detected in the genome of individual i, represents the sum of the number of single nucleotide polymorphisms detected in each sample, m i is the number of sequenced fragments in individual i, SNP ij represents the set of SNPs detected in fragment j, represents the size of the set, i.e. the number of SNPs in fragment j, represents the total number of indels detected in the genome of individual i, It means to add up the number of insertions and deletions detected in each sample. represents the set of indels detected in fragment j, represents the size of the set, i.e. the number of indels in fragment j, represents the total number of structural variations detected in the genome of individual i, Indicates the sum of the number of structural variations detected in each sample, SV i represents the set of structural variations detected in fragment j, represents the functional classification of all genetic variants detected in the genome of individual i, It means that all detected genetic variants are classified into functional categories; Statistical reports and annotations of genetic variation are analyzed in association with breeding goals to identify key gene loci associated with target traits: ,in, Indicates whether the variation v is related to the breeding target T.

6. The method for selecting a high-quality breeding population according to claim 1, characterized in that: Using a multi-layer convolutional neural network based on deep learning and combined with an attention mechanism, an association model is constructed to predict the association between the target trait and key gene loci, including the following steps: Based on the comprehensive annotation results detected i , perform preliminary screening on each variant and select a set of potential candidate gene loci C associated with the target trait i , , where v represents the variation, T is the breeding target, and P is the characteristic vector of the variation. is a function that evaluates the correlation between the variation v and the breeding target T, is the set threshold; A multi-layer convolutional neural network based on deep learning is used to further analyze the candidate gene loci and construct an initial association model for predicting the association between the target trait and the key gene loci. The multi-layer convolutional neural network structure includes multiple convolutional layers, pooling layers, and fully connected layers. An attention mechanism is introduced into the convolutional neural network to assign different weights to the characteristics of different gene loci and generate an attention weight matrix A. Output the predicted value of the association between the key gene loci and the target trait through the fully connected layer; Use the classification loss function to train the initial association model, optimize the network parameters, and obtain the association model.

7. The method for selecting a high-quality breeding population according to claim 1, characterized in that: The method includes the following steps: selecting breeding individuals with excellent traits according to the results of genome association analysis and the output value of the association model. Based on the results of genome association analysis, obtain the key gene loci set for each breeding individual , where v represents the variation, C i is the set of candidate gene loci, T k is the kth target trait, Represents the variation v and the target trait T k The correlation, p k Represents the target trait T k The weight of is the set threshold; The comprehensive score of each breeding individual is calculated based on the weight of each key gene locus and the predicted value output by the association model , where w v represents the weight of the key gene locus v, Represents the output value of the association model for the gene locus v, represents the similarity between breeding individuals i and j; According to the comprehensive score S i Sort the breeding individuals, select the individuals with the highest comprehensive scores from high to low, and build a high-quality breeding population: ,in, is the ratio of the set comprehensive score threshold, and N represents the total number of breeding individuals.

8. A high-quality breeding population selection device, characterized in that: include: A sequencing module is used to perform multiple high-throughput sequencing on each breeding individual in the breeding population to obtain genome sequence data; A preprocessing module, used to perform preliminary processing on the genome sequence data to obtain a preprocessed genome sequence; A comparison and analysis module is used to upload the pre-processed genome sequence to the cloud platform, and use a dynamic weighted comparison algorithm based on graph theory to compare and analyze the uploaded genome sequence data to generate an individual genome sequence map; A variation detection module, used to perform genetic variation detection on the genome sequence map and identify genetic variation data, wherein the genetic variation data includes single nucleotide polymorphisms, insertions and deletions, and structural variations; An association analysis module, used to perform genome association analysis based on the genetic variation data to identify key gene loci associated with the target trait; A model building module is used to build an association model that predicts the association between target traits and key gene loci using a multi-layer convolutional neural network based on deep learning and combined with an attention mechanism; The screening module is used to screen out breeding individuals with excellent traits based on the results of genome association analysis and the output values ​​of the association model.

9. A high-quality breeding population selection device, characterized in that: comprising a memory and at least one processor, wherein the memory has computer-readable instructions stored therein; The at least one processor calls the computer-readable instructions in the memory to execute the various steps of the method for selecting a high-quality breeding population as described in any one of claims 1-7.

10. A computer-readable storage medium having computer-readable instructions stored thereon, characterized in that: When the computer-readable instructions are executed by a processor, the various steps of the method for selecting a high-quality breeding population as described in any one of claims 1 to 7 are implemented.

Images